Induced-fit Conformational Change Research Articles

CHARMM-GUI-Based Induced Fit Docking Workflow to Generate Reliable Protein-Ligand Binding Modes.

Molecular docking is a preferred method to predict ligand binding modes and their binding energy to target protein receptors, which is critical in early phase structure-based drug discovery. However, there is a persistent challenge in docking that can be attributed to the induced fit effect, as receptor binding sites undergo induced fit conformational changes upon ligand binding to achieve better binding modes. In this work, based on CHARMM-GUI LBS Finder& Refiner and High-Throughput Simulator, we present a straightforward CHARMM-GUI induced fit docking (CGUI-IFD) workflow to generate reliable protein-ligand binding modes. The CGUI-IFD workflow generates an ensemble of receptor binding site conformations through ligand-binding site (LBS) refinement, runs rigid receptor docking, and performs high-throughput molecular dynamics (MD) simulations of protein-ligand complex structures in explicit solvents. The results are evaluated based on the ligand root-mean-square deviation (RMSD)-based binding stability and the molecular mechanics generalized Born surface area binding energy. For a benchmark test, we used 258 cross-docking protein-ligand pairs across 41 target proteins from the Schrodinger IFD-MD data set. The application of CGUI-IFD on this data set shows 80% success rate (within 2.5 Å RMSD from the experimental structures). We expect that the CGUI-IFD workflow can be useful to generate reliable ligand binding modes for cross-docking cases.

Induced fit docking workflow based on charmm-gui modules: LBS Finder & Refiner and High-Throughput Simulator

Molecular docking is a preferred method to predict ligand binding modes and their binding energy to target protein receptors, which is a critical first step in the early phase structure-based drug discovery. However, there is a persistent problem in docking that can be attributed to the induced fit effect, where the receptor binding site undergoes induced fit conformational changes upon ligand binding to achieve the correct binding mode. To solve this problem, we present a reliable stepwise docking workflow based on CHARMM-GUI LBS Finder & Refiner and High-Throughput Simulator.

Photooxidase Mimicking with Adaptive Coordination Molecular Capsules.

One important feature of enzyme catalysis is the induced-fit conformational change after binding substrates. Herein, we report a biomimetic water-soluble molecular capsule featuring adaptive structural change toward substrate binding, which offers an ideal platform for efficient photocatalysis. The molecular capsule was coordination-assembled from three anthracene-bridged bis-TPT [TPT = 2,4,6-tris(4-pyridyl)-1,3,5-triazine] ligands and six (bpy)Pd(NO3)2 (bpy = 2,2'-bipyridine). Once substrates bind to its hydrophobic cavity, this capsule would undergo quantitative capsule-to-bowl transformation. Visible-light absorption brought about by both the anthracene units and the charge-transfer absorption on the late-formed quintuple π-π stacked host-guest complex efficiently facilitates aerobic photooxidation for the sulfide guests by visible-light irradiation under mild conditions. Desired turnover numbers and product selectivity (sulfoxide over sulfone) have been achieved by the transformable nature of the catalyst and the hydrophilicity of the sulfoxide product. Such a photocatalytic process enabled by an adaptive coordination capsule and substrates as the allosteric effector paves the way for constructing artificial systems to mimic enzyme catalysis.

Bdellovibrio bacteriovorus phosphoglucose isomerase structures reveal novel rigidity in the active site of a selected subset of enzymes upon substrate binding.

Glycolysis and gluconeogenesis are central pathways of metabolism across all domains of life. A prominent enzyme in these pathways is phosphoglucose isomerase (PGI), which mediates the interconversion of glucose-6-phosphate and fructose-6-phosphate. The predatory bacterium Bdellovibrio bacteriovorus leads a complex life cycle, switching between intraperiplasmic replicative and extracellular 'hunter' attack-phase stages. Passage through this complex life cycle involves different metabolic states. Here we present the unliganded and substrate-bound structures of the B. bacteriovorus PGI, solved to 1.74 Å and 1.67 Å, respectively. These structures reveal that an induced-fit conformational change within the active site is not a prerequisite for the binding of substrates in some PGIs. Crucially, we suggest a phenylalanine residue, conserved across most PGI enzymes but substituted for glycine in B. bacteriovorus and other select organisms, is central to the induced-fit mode of substrate recognition for PGIs. This enzyme also represents the smallest conventional PGI characterized to date and probably represents the minimal requirements for a functional PGI.

Hairpin RNA-induced conformational change of a eukaryotic-specific lysyl-tRNA synthetase extension and role of adjacent anticodon-binding domain

Human lysyl-tRNA synthetase (hLysRS) is essential for aminoacylation of tRNALys Higher eukaryotic LysRSs possess an N-terminal extension (Nterm) previously shown to facilitate high-affinity tRNA binding and aminoacylation. This eukaryote-specific appended domain also plays a critical role in hLysRS nuclear localization, thus facilitating noncanonical functions of hLysRS. The structure is intrinsically disordered and therefore remains poorly characterized. Findings of previous studies are consistent with the Nterm domain undergoing a conformational transition to an ordered structure upon nucleic acid binding. In this study, we used NMR to investigate how the type of RNA, as well as the presence of the adjacent anticodon-binding domain (ACB), influences the Nterm conformation. To explore the latter, we used sortase A ligation to produce a segmentally labeled tandem-domain protein, Nterm-ACB. In the absence of RNA, Nterm remained disordered regardless of ACB attachment. Both alone and when attached to ACB, Nterm structure remained unaffected by titration with single-stranded RNAs. The central region of the Nterm domain adopted α-helical structure upon titration of Nterm and Nterm-ACB with RNA hairpins containing double-stranded regions. Nterm binding to the RNA hairpins resulted in CD spectral shifts consistent with an induced helical structure. NMR and fluorescence anisotropy revealed that Nterm binding to hairpin RNAs is weak but that the binding affinity increases significantly upon covalent attachment to ACB. We conclude that the ACB domain facilitates induced-fit conformational changes and confers high-affinity RNA hairpin binding, which may be advantageous for functional interactions of LysRS with a variety of different binding partners.

Probing remote residues important for catalysis in Escherichia coli ornithine transcarbamoylase.

Understanding how enzymes achieve their tremendous catalytic power is a major question in biochemistry. Greater understanding is also needed for enzyme engineering applications. In many cases, enzyme efficiency and specificity depend on residues not in direct contact with the substrate, termed remote residues. This work focuses on Escherichia coli ornithine transcarbamoylase (OTC), which plays a central role in amino acid metabolism. OTC has been reported to undergo an induced-fit conformational change upon binding its first substrate, carbamoyl phosphate (CP), and several residues important for activity have been identified. Using computational methods based on the computed chemical properties from theoretical titration curves, sequence-based scores derived from evolutionary history, and protein surface topology, residues important for catalytic activity were predicted. The roles of these residues in OTC activity were tested by constructing mutations at predicted positions, followed by steady-state kinetics assays and substrate binding studies with the variants. First-layer mutations R57A and D231A, second-layer mutation H272L, and third-layer mutation E299Q, result in 57- to 450-fold reductions in kcat/KM with respect to CP and 44- to 580-fold reductions with respect to ornithine. Second-layer mutations D140N and Y160S also reduce activity with respect to ornithine. Most variants had decreased stability relative to wild-type OTC, with variants H272L, H272N, and E299Q having the greatest decreases. Variants H272L, E299Q, and R57A also show compromised CP binding. In addition to direct effects on catalytic activity, effects on overall protein stability and substrate binding were observed that reveal the intricacies of how these residues contribute to catalysis.

The non-glycosylated N-terminal domain of human thrombopoietin is a molten globule under native conditions.

Human thrombopoietin (hTPO) is a primary hematopoietic growth factor that regulates megakaryocytopoiesis and platelet production. The non-glycosylated form of 1-163 residues of hTPO (hTPO163 ) including the N-terminal active site domain (1-153 residues) is a candidate for treating thrombocytopenia. However, the autoantigenicity level of hTPO163 is higher than that of the full-length glycosylated hTPO (ghTPO332 ). In order to clarify the structural and physicochemical properties of hTPO163 , circular dichroism (CD) and differential scanning calorimetry (DSC) analyses were performed. CD analysis indicated that hTPO163 undergoes an induced-fit conformational change (+19.0% for helix and -16.7% for β-strand) upon binding to the neutralizing antibody TN1 in a manner similar to the coupled folding and binding mechanism. Moreover, DSC analysis showed that the thermal transition process of hTPO163 is a multistate transition; hTPO163 is thermally stabilized upon receptor (c-Mpl) binding, as indicated with raising the midpoint (Tm ) temperature of the transition by at least +9.5 K. The conformational variability and stability of hTPO163 indicate that hTPO163 exists as a molten globule under native conditions, which may enable the induced-fit conformational change according to the type of ligands (antibodies and receptor). Additionally, CD and computational analyses indicated that the C-terminal domain (154-332 residues) and glycosylation assists the folding of the N-terminal domain. These observations suggest that the antibody affinity and autoantigenicity of hTPO163 might be reduced, if the conformational variability of hTPO163 is restricted by mutation and/or by the addition of C-terminal domain with glycosylation to keep its conformation suitable for the c-Mpl recognition.

Structure of the heterophilic interaction between the nectin-like 4 and nectin-like 1 molecules

Nectin-like (Necl) molecules are Ca2+-independent Ig-like transmembrane cell adhesion molecules that participate in junctions between different cell types. The specific cell-cell adhesions mediated by Necl proteins are important in neural development and have been implicated in neurodegenerative diseases. Here, we present the crystal structure of the mouse Necl-4 full ectodomain and the structure of the heterophilic Necl ectodomain complex formed by the mNecl-4 and mNecl-1 ectodomains. We demonstrate that, while the ectodomain of mNecl-4 is monomeric, it forms a stable heterodimer with Ig1 of mNecl-1, with an affinity significantly higher than that observed for self-dimerization of the mNecl-1 ectodomain. We validated our structural characterizations by performing a surface plasmon resonance assay and an Fc fusion protein binding assay in mouse primary dorsal root ganglia neurites and Schwann cells and identified a selection of residues important for heterophilic interactions. Finally, we proposed a model of Necl binding specificity that involves an induced-fit conformational change at the dimerization interface.

Transferase Versus Hydrolase: The Role of Conformational Flexibility in Reaction Specificity

Active in the aqueous cellular environment where a massive excess of water is perpetually present, enzymes that catalyze the transfer of an electrophile to a non-water nucleophile (transferases) require specific strategies to inhibit mechanistically related hydrolysis reactions. To identify principles that confer transferase versus hydrolase reaction specificity, we exploited two enzymes that use highly similar catalytic apparatuses to catalyze the transglycosylation (a transferase reaction) or hydrolysis of α-1,3-glucan linkages in the cyclic tetrasaccharide cycloalternan (CA). We show that substrate binding to non-catalytic domains and a conformationally stable active site promote CA transglycosylation, whereas a distinct pattern of active site conformational change is associated with CA hydrolysis. These findings defy the classic view of induced-fit conformational change and illustrate a mechanism by which a stable hydrophobic binding site can favor transferase activity and disfavor hydrolysis. Application of these principles could facilitate the rational reengineering of transferases with desired catalytic properties.

A Network of Conformational Transitions Revealed by Molecular Dynamics Simulations of the Binary Complex of Escherichia coli 6-Hydroxymethyl-7,8-dihydropterin Pyrophosphokinase with MgATP.

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the first reaction in the folate biosynthetic pathway. Comparison of its X-ray and nuclear magnetic resonance structures suggests that the enzyme undergoes significant conformational change upon binding to its substrates, especially in three catalytic loops. Experimental research has shown that, in its binary form, even bound by analogues of MgATP, loops 2 and 3 remain rather flexible; this raises questions about the putative large-scale induced-fit conformational change of the HPPK-MgATP binary complex. In this work, long-time all-atomic molecular dynamics simulations were conducted to investigate the loop dynamics in this complex. Our simulations show that, with loop 3 closed, multiple conformations of loop 2, including the open, semiopen, and closed forms, are all accessible to the binary complex. These results provide valuable structural insights into the details of conformational changes upon 6-hydroxymethyl-7,8-dihydropterin (HP) binding and biological activities of HPPK. Conformational network analysis and principal component analysis related to the loops are also discussed.

Structural and Functional Influence of the Glycine-Rich Loop G302GGGY on the Catalytic Tyrosine of Histone Deacetylase 8.

Histone deacetylase 8 (HDAC8) catalyzes the hydrolysis of acetyl-l-lysine to yield products l-lysine and acetate through a mechanism in which a nucleophilic water molecule is activated by a histidine general base and a catalytic metal ion (Zn2+ or Fe2+). Acetyl-l-lysine also requires activation by metal coordination and a hydrogen bond with catalytic tyrosine Y306, which also functions in transition state stabilization. Interestingly, Y306 is located in the conserved glycine-rich loop G302GGGY. The potential flexibility afforded by the tetraglycine segment may facilitate induced-fit conformational changes in Y306 between "in" and "out" positions, as observed in related deacetylases. To probe the catalytic importance of the glycine-rich loop in HDAC8, we rigidified this loop by preparing the G302A, G303A, G304A, and G305A mutants and measured their steady state kinetics and determined their X-ray crystal structures. Substantial losses of catalytic efficiency are observed (10-500-fold based on kcat/KM), particularly for G304A HDAC8 and G305A HDAC8. These mutants also exhibit the greatest structural changes for catalytic tyrosine Y306 (1.3-1.7 Å shifts of the phenolic hydroxyl group). Molecular dynamics simulations further indicate that G304 and G305 undergo pronounced structural changes as residue 306 undergoes a transition between "in" and "out" conformations. Thus, the G304A and G305A substitutions likely compromise the position and conformational changes of Y306 required for substrate activation and transition state stabilization. The G302A and G303A substitutions have less severe catalytic consequences, and these substitutions may influence an internal channel through which product acetate is believed to exit.

An insight into the thermodynamic characteristics of human thrombopoietin complexation with TN1 antibody.

Human thrombopoietin (hTPO) primarily stimulates megakaryocytopoiesis and platelet production and is neutralized by the mouse TN1 antibody. The thermodynamic characteristics of TN1 antibody-hTPO complexation were analyzed by isothermal titration calorimetry (ITC) using an antigen-binding fragment (Fab) derived from the TN1 antibody (TN1-Fab). To clarify the mechanism by which hTPO is recognized by TN1-Fab the conformation of free TN1-Fab was determined to a resolution of 2.0 Å using X-ray crystallography and compared with the hTPO-bound form of TN1-Fab determined by a previous study. This structural comparison revealed that the conformation of TN1-Fab does not substantially change after hTPO binding and a set of 15 water molecules is released from the antigen-binding site (paratope) of TN1-Fab upon hTPO complexation. Interestingly, the heat capacity change (ΔCp) measured by ITC (-1.52 ± 0.05 kJmol(-1) K(-1) ) differed significantly from calculations based upon the X-ray structure data of the hTPO-bound and unbound forms of TN1-Fab (-1.02 ∼ 0.25 kJmol(-1) K(-1) ) suggesting that hTPO undergoes an induced-fit conformational change combined with significant desolvation upon TN1-Fab binding. The results shed light on the structural biology associated with neutralizing antibody recognition.

Structure of neprilysin in complex with the active metabolite of sacubitril.

Sacubitril is an ethyl ester prodrug of LBQ657, the active neprilysin (NEP) inhibitor, and a component of LCZ696 (sacubitril/valsartan). We report herein the three-dimensional structure of LBQ657 in complex with human NEP at 2 Å resolution. The crystal structure unravels the binding mode of the compound occupying the S1, S1’ and S2’ sub-pockets of the active site, consistent with a competitive inhibition mode. An induced fit conformational change upon binding of the P1’-biphenyl moiety of the inhibitor suggests an explanation for its selectivity against structurally homologous zinc metallopeptidases.

Molecular dynamics simulations of the Escherichia coli HPPK apo-enzyme reveal a network of conformational transitions.

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the first reaction in the folate biosynthetic pathway. Comparison of its X-ray and nuclear magnetic resonance structures suggests that the enzyme undergoes significant conformational change upon binding to its substrates, especially in three catalytic loops. Experimental research has shown that even when confined by crystal contacts, loops 2 and 3 remain rather flexible when the enzyme is in its apo form, raising questions about the putative large-scale induced-fit conformational change of HPPK. To investigate the loop dynamics in a crystal-free environment, we performed conventional molecular dynamics simulations of the apo-enzyme at two different temperatures (300 and 350 K). Our simulations show that the crystallographic B-factors considerably underestimate the loop dynamics; multiple conformations of loops 2 and 3, including the open, semi-open, and closed conformations that an enzyme must adopt throughout its catalytic cycle, are all accessible to the apo-enzyme. These results revise our previous view of the functional mechanism of conformational change upon MgATP binding and offer valuable structural insights into the workings of HPPK. In this paper, conformational network analysis and principal component analysis related to the loops are discussed to support the presented conclusions.

Rational Modulation of the Induced-Fit Conformational Change for Slow-Onset Inhibition in Mycobacterium tuberculosis InhA.

Slow-onset enzyme inhibitors are the subject of considerable interest as an approach to increasing the potency of pharmaceutical compounds by extending the residence time of the inhibitor on the target (the lifetime of the drug-receptor complex). However, rational modulation of residence time presents significant challenges because it requires additional mechanistic insight, such as the nature of the transition state for postbinding isomerization. Our previous work, based on X-ray crystallography, enzyme kinetics, and molecular dynamics simulation, suggested that the slow step in inhibition of the Mycobacterium tuberculosis enoyl-ACP reductase InhA involves a change in the conformation of the substrate binding loop from an open state in the initial enzyme-inhibitor complex to a closed state in the final enzyme-inhibitor complex. Here, we use multidimensional free energy landscapes for loop isomerization to obtain a computational model for the transition state. The results suggest that slow-onset inhibitors crowd key side chains on helices that slide past each other during isomerization, resulting in a steric clash. The landscapes become significantly flatter when residues involved in the steric clash are replaced with alanine. Importantly, this lower barrier can be increased by rational inhibitor redesign to restore the steric clash. Crystallographic studies and enzyme kinetics confirm the predicted effects on loop structure and flexibility, as well as inhibitor residence time. These loss and regain of function studies validate our mechanistic hypothesis for interactions controlling substrate binding loop isomerization, providing a platform for the future design of inhibitors with longer residence times and better in vivo potency. Similar opportunities for slow-onset inhibition via the same mechanism are identified in other pathogens.

Remote Residues Affect Stability of Ornithine Transcarbamylase

Ornithine transcarbamylase (OTC) is known to undergo induced-fit conformational changes upon the binding of carbamoyl phosphate (CP), which allows the sequential binding of L-ornithine (ORN) to produce citrulline (CIT) and inorganic phosphate. Computational tools THEMATICS and POOL predict that catalytic residues of OTC extend into the second and third spheres around the active site of the protein. To test these predictions, mutations of POOL-predicted residues were tested with a kinetics assay to determine relative catalytic efficiencies. The stabilities of these mutations are being characterized in this study by a ThermoFluor assay that measures the melting temperature of OTC variants. OTC variants C273A, D140N, Y229F, H272L, H272N, E299D, E299Q and R57A showed 5-fold decrease or greater in catalytic efficiency with ORN and/or CP. Mutation of these residues also has been found to decrease the melting temperature of OTC up to 18 ᵒC. Thus these variants show less stability compared to wild-type OTC. It is also observed that when OTC is bound to CP, the stability is increased compared to the apoenzyme, or to OTC with ORN- or CIT-bound. Supported by NSF-MCB-1158176.

Ornithine Transcarbamylase Has a Spatially Extended Active Site as Predicted

Partial Order Optimum Likelihood (POOL), developed at Northeastern University, is a machine learning technique to predict catalytically important residues based on the tertiary structure of a protein and on computed electrostatic properties. POOL has been shown to predict accurately the catalytic residues and discern between a compact and extended active site. Classic studies of enzymes typically identify catalytic residues that are in direct contact with the substrate. POOL predicts a spatially extended active site for ornithine transcarbamylase (OTC), for which dynamic processes are believed to play a role in its catalytic mechanism. OTC is reported to undergo induced-fit conformational changes upon binding carbamoyl phosphate, followed by binding of ornithine. POOL predicts OTC to have an extended triple-layer active site. Kinetics assay of Asp140, His272, Glu299 and Arg57 variants show these POOL-predicted remote residues, located in the second and third layers, are important for catalysis. Specifically, variants D140N, H272L, E299Q and R57A (positive control), respectively have catalytic efficiency relative to wild type of 10%, 3%, 6% and 2% for ornithine and 21%, 45%, 66% and 54% for carbamoyl phosphate. Negative controls not predicted by POOL to be catalytically important were S61A and Q104L with relative catalytic efficiencies of 405% and 205% for ornithine and 88% and 64% for carbamoyl phosphate respectively. Our observations indicate the importance of remote residues for OTC activity and POOL's power to accurately predict catalytically important residues. Work supported by NSF-MCB-1158176.

Structure of inorganic pyrophosphatase from Staphylococcus aureus reveals conformational flexibility of the active site.

Cytoplasmic inorganic pyrophosphatase (PPiase) is an enzyme essential for survival of organisms, from bacteria to human. PPiases are divided into two structurally distinct families: family I PPiases are Mg(2+)-dependent and present in most archaea, eukaryotes and prokaryotes, whereas the relatively less understood family II PPiases are Mn(2+)-dependent and present only in some archaea, bacteria and primitive eukaryotes. Staphylococcus aureus (SA), a dangerous pathogen and a frequent cause of hospital infections, contains a family II PPiase (PpaC), which is an attractive potential target for development of novel antibacterial agents. We determined a crystal structure of SA PpaC in complex with catalytic Mn(2+) at 2.1Å resolution. The active site contains two catalytic Mn(2+) binding sites, each half-occupied, reconciling the previously observed 1:1 Mn(2+):enzyme stoichiometry with the presence of two divalent metal ion sites in the apo-enzyme. Unexpectedly, despite the absence of the substrate or products in the active site, the two domains of SA PpaC form a closed active site, a conformation observed in structures of other family II PPiases only in complex with substrate or product mimics. A region spanning residues 295-298, which contains a conserved substrate binding RKK motif, is flipped out of the active site, an unprecedented conformation for a PPiase. Because the mutant of Arg295 to an alanine is devoid of activity, this loop likely undergoes an induced-fit conformational change upon substrate binding and product dissociation. This closed conformation of SA PPiase may serve as an attractive target for rational design of inhibitors of this enzyme.

Energetic dissection of Gleevec's selectivity toward human tyrosine kinases.

Protein kinases are obvious drug targets against cancer due to their central role in cellular regulation. Since the discovery of Gleevec, a potent and specific inhibitor of Abl kinase, as a highly successful cancer therapeutic, the ability of this drug to distinguish between Abl and other tyrosine kinases like Src has been intensely investigated, but without much success. Using NMR and fast kinetics, we establish a novel model that solves this longstanding question of two tyrosine kinases adopting almost identical structures when bound to Gleevec, yet having vastly different affinities. In contrast to all other proposed models we show that the origin of Abl’s high affinity lies predominantly in a conformational change after binding. An energy landscape that provides tight affinity via an induced-fit and binding plasticity via conformational selection mechanism is likely to be general for many inhibitors.

Crystal Structures of Human Transthyretin Complexed with Glabridin

Transthyretin (TTR) is a plasma protein implicated in human amyloid diseases. Several small molecules that bind to the thyroxine-binding site of TTR have been shown to stabilize the TTR tetramer and to inhibit amyloid fibril formation of TTR. Herein, we demonstrated that glabridin (Glab), a prenylated isoflavan isolated from Glycyrrhiza glabra L., inhibited aggregation of TTR in a thioflavin assay. The TTR-Glab complex structure revealed a novel binding mode including a CH-π interaction with A108 and a hydrogen bond with K15. A structural comparison with the wild type-apo structure revealed that the CH-π interaction with A108 was strengthened by the induced-fit conformational change upon Glab binding. Furthermore, the binding of Glab induced a rotation of the T119 side chain, and the inclusion of a water molecule, leading to stabilization of the dimer-dimer interface. These results demonstrate that Glab is a novel inhibitor of TTR fibrillization and suggest the molecular mechanism by which Glab binding stabilizes the tetramer.