Dynamics of the flexible loop of triosephosphate isomerase: the loop motion is not ligand gated

Reactivity Tracking of an Enzyme Progress Coordinate

The reactivity of individual solvent-coupled protein configurations is used to track and resolve the progress coordinate for the core reaction sequence of substrate radical rearrangement and hydrogen atom transfer in the ethanolamine ammonia-lyase (EAL) enzyme from Salmonella enterica. The first-order decay of the substrate radical intermediate is the monitored reaction. Heterogeneous confinement from sucrose hydrates in the mesophase solvent surrounding the cryotrapped protein introduces distributed kinetics in the non-native decay of the substrate radical pair capture substate, which arise from an ensemble of configurational microstates. Reaction rates increase by >103-fold across the distribution to approach that for the native enabled substate for radical rearrangement, which reacts with monotonic kinetics. The native progress coordinate thus involves a collapse of the configuration space to generate optimized reactivity. Reactivity tracking reveals fundamental features of solvent-protein-reaction configurational coupling and leads to a model that refines the ensemble paradigm of enzyme catalysis for strongly adiabatic chemical steps.

Loop dynamics and the evolution of enzyme activity

In the early 2000s, Tawfik presented his 'New View' on enzyme evolution, highlighting the role of conformational plasticity in expanding the functional diversity of limited repertoires of sequences. This view is gaining increasing traction with increasing evidence of the importance of conformational dynamics in both natural and laboratory evolution of enzymes. The past years have seen several elegant examples of harnessing conformational (particularly loop) dynamics to successfully manipulate protein function. This Review revisits flexible loops as critical participants in regulating enzyme activity. We showcase several systems of particular interest: triosephosphate isomerase barrel proteins, protein tyrosine phosphatases and β-lactamases, while briefly discussing other systems in which loop dynamics are important for selectivity and turnover. We then discuss the implications for engineering, presenting examples of successful loop manipulation in either improving catalytic efficiency, or changing selectivity completely. Overall, it is becoming clearer that mimicking nature by manipulating the conformational dynamics of key protein loops is a powerful method of tailoring enzyme activity, without needing to target active-site residues.

Modulation of Flexible Loops in Catalytic Cavities Reveals the Thermal Activation Mechanism of a Glycogen-Debranching Enzyme

Some thermophilic enzymes not only exhibit high thermostability at high temperatures but also have an activation effect by thermal incubation. However, the correlations between temperature-induced structural modulation and thermal activation are still unclear. In this study, we selected a thermophilic glycogen-debranching enzyme from Saccharolobus solfataricus STB09 (SsGDE), which was a promising starch-debranching enzyme with a thermal activation property at temperatures ranging from 50 to 70 °C, to explore the thermal activation mechanism. Molecular dynamics simulations were performed for SsGDE at 30, 50, or 70 °C to reveal the temperature dependence of structure modulation and catalytic function. The results revealed that four loops (loop1 313-337, loop2 399-418, loop3 481-513, and loop4 540-574) in SsGDE were reshaped, which made the catalytic cavity more open. The internal residues, including the catalytic triad Asp3631, Glu399, and Asp471, could be exposed, due to the structural modulation, to exert catalytic functions. We proposed that the thermal activation effect of SsGDE was closely associated with the temperature-induced modulation of the catalytic cavity, which paved the way for further engineering enzymes to achieve higher catalytic performance and stability.

Ligand-Based Virtual Screening and Molecular Docking of Benzimidazoles as Potential Inhibitors of Triosephosphate Isomerase Identified New Trypanocidal Agents

Trypanosoma cruzi (T. cruzi) is a parasite that affects humans and other mammals. T. cruzi depends on glycolysis as a source of adenosine triphosphate (ATP) supply, and triosephosphate isomerase (TIM) plays a key role in this metabolic pathway. This enzyme is an attractive target for the design of new trypanocidal drugs. In this study, a ligand-based virtual screening (LBVS) from the ZINC15 database using benzimidazole as a scaffold was accomplished. Later, a molecular docking on the interface of T. cruzi TIM (TcTIM) was performed and the compounds were grouped by interaction profiles. Subsequently, a selection of compounds was made based on cost and availability for in vitro evaluation against blood trypomastigotes. Finally, the compounds were analyzed by molecular dynamics simulation, and physicochemical and pharmacokinetic properties were determined using SwissADME software. A total of 1604 molecules were obtained as potential TcTIM inhibitors. BP2 and BP5 showed trypanocidal activity with half-maximal lytic concentration (LC₅₀) values of 155.86 and 226.30 µM, respectively. Molecular docking and molecular dynamics simulation analyzes showed a favorable docking score of BP5 compound on TcTIM. Additionally, BP5 showed a low docking score (−5.9 Kcal/mol) on human TIM compared to the control ligand (−7.2 Kcal/mol). Both compounds BP2 and BP5 showed good physicochemical and pharmacokinetic properties as new anti-T. cruzi agents.

Sampling of Protein Conformational Space Using Hybrid Simulations: A Critical Assessment of Recent Methods

Recent years have seen several hybrid simulation methods for exploring the conformational space of proteins and their complexes or assemblies. These methods often combine fast analytical approaches with computationally expensive full atomic molecular dynamics (MD) simulations with the goal of rapidly sampling large and cooperative conformational changes at full atomic resolution. We present here a systematic comparison of the utility and limits of four such hybrid methods that have been introduced in recent years: MD with excited normal modes (MDeNM), collective modes-driven MD (CoMD), and elastic network model (ENM)-based generation, clustering, and relaxation of conformations (ClustENM) as well as its updated version integrated with MD simulations (ClustENMD). We analyzed the predicted conformational spaces using each of these four hybrid methods, applied to four well-studied proteins, triosephosphate isomerase (TIM), 3-phosphoglycerate kinase (PGK), HIV-1 protease (PR) and HIV-1 reverse transcriptase (RT), which provide extensive ensembles of experimental structures for benchmarking and comparing the methods. We show that a rigorous multi-faceted comparison and multiple metrics are necessary to properly assess the differences between conformational ensembles and provide an optimal protocol for achieving good agreement with experimental data. While all four hybrid methods perform well in general, being especially useful as computationally efficient methods that retain atomic resolution, the systematic analysis of the same systems by these four hybrid methods highlights the strengths and limitations of the methods and provides guidance for parameters and protocols to be adopted in future studies.

Flexible Loop in Carbohydrate-Binding Module 48 Allosterically Modulates Substrate Binding of the 1,4-α-Glucan Branching Enzyme

The 1,4-α-glucan branching enzyme (GBE, EC 2.4.1.18) catalyzes the formation of α-1,6 branching points in starch and plays a key role in synthesis. To obtain mechanistic insights into the catalytic action of the enzyme, we first determined the crystal structure of GBE from Rhodothermus obamensis STB05 (RoGBE) to a resolution of 2.39 Å (PDB ID: 6JOY). The structure consists of three domains: domain A, domain C, and the carbohydrate-binding module 48 (CBM48). An engineered truncated mutant lacking the CBM48 domain (ΔCBM48) showed significantly reduced ligand binding affinity and enzyme activity. Comparison of the structures of RoGBE with other GBEs showed that CBM48 of RoGBE had a longer flexible loop. Truncation of the flexible loops resulted in reduced binding affinity and activity, thereby substantiating the importance of the optimum loop structure for catalysis. In essence, our study shows that CBM48, especially the flexible loop, plays an important role in substrate binding and enzymatic activity of RoGBE. Further, based on the structural analysis, kinetics, and activity assays on wild type and mutants, as well as homology modeling, we proposed a mechanistic model (called the "lid model") to illustrate how the flexible loop triggers substrate binding, ultimately leading to catalysis.

Loop Dynamics and Enzyme Catalysis in Protein Tyrosine Phosphatases

Protein tyrosine phosphatases (PTPs) play an important role in cellular signaling and have been implicated in human cancers, diabetes, and obesity. Despite shared catalytic mechanisms and transition states for the chemical steps of catalysis, catalytic rates within the PTP family vary over several orders of magnitude. These rate differences have been implied to arise from differing conformational dynamics of the closure of a protein loop, the WPD-loop, which carries a catalytically critical residue. The present work reports computational studies of the human protein tyrosine phosphatase 1B (PTP1B) and YopH from Yersinia pestis, for which NMR has demonstrated a link between their respective rates of WPD-loop motion and catalysis rates, which differ by an order of magnitude. We have performed detailed structural analysis, both conventional and enhanced sampling simulations of their loop dynamics, as well as empirical valence bond simulations of the chemical step of catalysis. These analyses revealed the key residues and structural features responsible for these differences, as well as the residues and pathways that facilitate allosteric communication in these enzymes. Curiously, our wild-type YopH simulations also identify a catalytically incompetent hyper-open conformation of its WPD-loop, sampled as a rare event, previously only experimentally observed in YopH-based chimeras. The effect of differences within the WPD-loop and its neighboring loops on the modulation of loop dynamics, as revealed in this work, may provide a facile means for the family of PTP enzymes to respond to environmental changes and regulate their catalytic activities.

Backbone assignments and conformational dynamics in the S. typhimurium tryptophan synthase α-subunit from solution-state NMR

Backbone assignments for the isolated α-subunit of Salmonella typhimurium tryptophan synthase (TS) are reported based on triple resonance solution-state NMR experiments on a uniformly 2H,13C,15N-labeled sample. From the backbone chemical shifts, secondary structure and random coil index order parameters (RCI-S2) are predicted. Titration with the 3-indole-D-glycerol 3'-phosphate analog, N-(4'-trifluoromethoxybenzenesulfonyl)-2-aminoethyl phosphate (F9), leads to chemical shift perturbations indicative of conformational changes from which an estimate of the dissociation constant is obtained. Comparisons of the backbone chemical-shifts, RCI-S2 values, and site-specific relaxation times with and without F9 reveal allosteric changes including modulation in secondary structures and loop rigidity induced upon ligand binding. A comparison is made to the X-ray crystal structure of the α-subunit in the full TS αββα bi-enzyme complex and to two new X-ray crystal structures of the isolated TS α-subunit reported in this work.

Harnessing Conformational Plasticity to Generate Designer Enzymes

Recent years have witnessed an explosion of interest in understanding the role of conformational dynamics both in the evolution of new enzymatic activities from existing enzymes and in facilitating the emergence of enzymatic activity de novo on scaffolds that were previously non-catalytic. There are also an increasing number of examples in the literature of targeted engineering of conformational dynamics being successfully used to alter enzyme selectivity and activity. Despite the obvious importance of conformational dynamics to both enzyme function and evolvability, many (although not all) computational design approaches still focus either on pure sequence-based approaches or on using structures with limited flexibility to guide the design. However, there exist a wide variety of computational approaches that can be (re)purposed to introduce conformational dynamics as a key consideration in the design process. Coupled with laboratory evolution and more conventional existing sequence- and structure-based approaches, these techniques provide powerful tools for greatly expanding the protein engineering toolkit. This Perspective provides an overview of evolutionary studies that have dissected the role of conformational dynamics in facilitating the emergence of novel enzymes, as well as advances in computational approaches that allow one to target conformational dynamics as part of enzyme design. Harnessing conformational dynamics in engineering studies is a powerful paradigm with which to engineer the next generation of designer biocatalysts.

The role of ligand-gated conformational changes in enzyme catalysis

Structural and biochemical studies on diverse enzymes have highlighted the importance of ligand-gated conformational changes in enzyme catalysis, where the intrinsic binding energy of the common phosphoryl group of their substrates is used to drive energetically unfavorable conformational changes in catalytic loops, from inactive open to catalytically competent closed conformations. However, computational studies have historically been unable to capture the activating role of these conformational changes. Here, we discuss recent experimental and computational studies, which can remarkably pinpoint the role of ligand-gated conformational changes in enzyme catalysis, even when not modeling the loop dynamics explicitly. Finally, through our joint analyses of these data, we demonstrate how the synergy between theory and experiment is crucial for furthering our understanding of enzyme catalysis.

Dynamic Nuclear Polarization Magic-Angle Spinning Nuclear Magnetic Resonance Combined with Molecular Dynamics Simulations Permits Detection of Order and Disorder in Viral Assemblies

We report dynamic nuclear polarization (DNP)-enhanced magic-angle spinning (MAS) NMR spectroscopy in viral capsids from HIV-1 and bacteriophage AP205. Viruses regulate their life cycles and infectivity through modulation of their structures and dynamics. While static structures of capsids from several viruses are now accessible with near-atomic-level resolution, atomic-level understanding of functionally important motions in assembled capsids is lacking. We observed up to 64-fold signal enhancements by DNP, which permitted in-depth analysis of these assemblies. For the HIV-1 CA assemblies, a remarkably high spectral resolution in the 3D and 2D heteronuclear data sets permitted the assignment of a significant fraction of backbone and side-chain resonances. Using an integrated DNP MAS NMR and molecular dynamics (MD) simulation approach, the conformational space sampled by the assembled capsid at cryogenic temperatures was mapped. Qualitatively, a remarkable agreement was observed for the experimental ¹³C/¹⁵N chemical shift distributions and those calculated from substructures along the MD trajectory. Residues that are mobile at physiological temperatures are frozen out in multiple conformers at cryogenic conditions, resulting in broad experimental and calculated chemical shift distributions. Overall, our results suggest that DNP MAS NMR measurements in combination with MD simulations facilitate a thorough understanding of the dynamic signatures of viral capsids.

Active-Site Glu165 Activation in Triosephosphate Isomerase and Its Deprotonation Kinetics

Triosephosphate isomerase (TIM) catalyzes the interconversion between dihydroxyacetone phosphate (DHAP) and d-glyceraldehyde 3-phosphate (GAP) via an enediol(ate) intermediate. The active-site residue Glu165 serves as the catalytic base during catalysis. It abstracts a proton from C1 carbon of DHAP to form the reaction intermediate and donates a proton to C2 carbon of the intermediate to form product GAP. Our difference Fourier transform infrared spectroscopy studies on the yeast TIM (YeTIM)/phosphate complex revealed a C═O stretch band at 1706 cm^-1 from the protonated Glu165 carboxyl group at pH 7.5, indicating that the p K_a of the catalytic base is increased by >3.0 pH units upon phosphate binding, and that the Glu165 carboxyl environment in the complex is still hydrophilic in spite of the increased p K_a. Hence, the results show that the binding of the phosphodianion group is part of the activation mechanism which involves the p K_a elevation of the catalytic base Glu165. The deprotonation kinetics of Glu165 in the μs to ms time range were determined via infrared (IR) T-jump studies on the YeTIM/phosphate and ("heavy enzyme") [U-¹³C,-¹⁵N]YeTIM/phosphate complexes. The slower deprotonation kinetics in the ms time scale is due to phosphate dissociation modulated by the loop motion, which slows down by enzyme mass increase to show a normal heavy enzyme kinetic isotope effect (KIE) ∼1.2 (i.e., slower rate in the heavy enzyme). The faster deprotonation kinetics in the tens of μs time scale is assigned to temperature-induced p K_a decrease, while phosphate is still bound, and it shows an inverse heavy enzyme KIE ∼0.89 (faster rate in the heavy enzyme). The IR static and T-jump spectroscopy provides atomic-level resolution of the catalytic mechanism because of its ability to directly observe the bond breaking/forming process.

Loop Motion in Triosephosphate Isomerase Is Not a Simple Open and Shut Case

Conformational changes are crucial for the catalytic action of many enzymes. A prototypical and well-studied example is loop opening and closure in triosephosphate isomerase (TIM), which is thought to determine the rate of catalytic turnover in many circumstances. Specifically, TIM loop 6 "grips" the phosphodianion of the substrate and, together with a change in loop 7, sets up the TIM active site for efficient catalysis. Crystal structures of TIM typically show an open or a closed conformation of loop 6, with the tip of the loop moving ∼7 Å between conformations. Many studies have interpreted this motion as a two-state, rigid-body transition. Here, we use extensive molecular dynamics simulations, with both conventional and enhanced sampling techniques, to analyze loop motion in apo and substrate-bound TIM in detail, using five crystal structures of the dimeric TIM from Saccharomyces cerevisiae. We find that loop 6 is highly flexible and samples multiple conformational states. Empirical valence bond simulations of the first reaction step show that slight displacements away from the fully closed-loop conformation can be sufficient to abolish most of the catalytic activity; full closure is required for efficient reaction. The conformational change of the loops in TIM is thus not a simple "open and shut" case and is crucial for its catalytic action. Our detailed analysis of loop motion in a highly efficient enzyme highlights the complexity of loop conformational changes and their role in biological catalysis.

Conformational Sampling of the Intrinsically Disordered C-Terminal Tail of DERA Is Important for Enzyme Catalysis

2-Deoxyribose-5-phosphate aldolase (DERA) catalyzes the reversible conversion of acetaldehyde and glyceraldehyde-3-phosphate into deoxyribose-5-phosphate. DERA is used as a biocatalyst for the synthesis of drugs such as statins and is a promising pharmaceutical target due to its involvement in nucleotide catabolism. Despite previous biochemical studies suggesting the catalytic importance of the C-terminal tyrosine residue found in several bacterial DERAs, the structural and functional basis of its participation in catalysis remains elusive because the electron density for the last eight to nine residues (i.e., the C-terminal tail) is absent in all available crystal structures. Using a combination of NMR spectroscopy and molecular dynamics simulations, we conclusively show that the rarely studied C-terminal tail of E. coli DERA (ecDERA) is intrinsically disordered and exists in equilibrium between open and catalytically relevant closed states, where the C-terminal tyrosine (Y259) enters the active site. Nuclear Overhauser effect distance restraints, obtained due to the presence of a substantial closed state population, were used to derive the solution-state structure of the ecDERA closed state. Real-time NMR hydrogen/deuterium exchange experiments reveal that Y259 is required for efficiency of the proton abstraction step of the catalytic reaction. Phosphate titration experiments show that, in addition to the phosphate-binding residues located near the active site, as observed in the available crystal structures, ecDERA contains previously unknown auxiliary phosphate-binding residues on the C-terminal tail which could facilitate in orienting Y259 in an optimal position for catalysis. Thus, we present significant insights into the structural and mechanistic importance of the ecDERA C-terminal tail and illustrate the role of conformational sampling in enzyme catalysis.

Asymmetric Anchoring Is Required for Efficient Ω-Loop Opening and Closing in Cytosolic Phosphoenolpyruvate Carboxykinase

Mobile Ω-loops play essential roles in the function of many enzymes. Here we investigated the importance of a residue lying outside of the mobile Ω-loop element in the catalytic function of an H477R variant of cytosolic phosphoenolpyruvate carboxykinase using crystallographic, kinetic, and computational analysis. The crystallographic data suggest that the efficient transition of the Ω-loop to the closed conformation requires stabilization of the N-terminus of the loop through contacts between R461 and E588. In contrast, the C-terminal end of the Ω-loop undergoes changing interactions with the enzyme body through contacts between H477 at the C-terminus of the loop and E591 located on the enzyme body. Potential of mean force calculations demonstrated that altering the anchoring of the C-terminus of the Ω-loop via the H477R substitution results in the destabilization of the closed state of the Ω-loop by 3.4 kcal mol^-1. The kinetic parameters for the enzyme were altered in an asymmetric fashion with the predominant effect being observed in the direction of oxaloacetate synthesis. This is exemplified by a reduction in k_cat for the H477R mutant by an order of magnitude in the direction of OAA synthesis, while in the direction of PEP synthesis, it decreased by a factor of only 2. The data are consistent with a mechanism for loop conformational exchange between open and closed states in which a balance between fixed anchoring of the N-terminus of the Ω-loop and a flexible, unattached C-terminus drives the transition between a disordered (open) state and an ordered (closed) state.

Comparative Study of the Collective Dynamics of Proteins and Inorganic Nanoparticles

Molecular dynamics simulations of ubiquitin in water/glycerol solutions are used to test the suggestion by Karplus and coworkers that proteins in their biologically active state should exhibit a dynamics similar to 'surface-melted' inorganic nanoparticles (NPs). Motivated by recent studies indicating that surface-melted inorganic NPs are in a 'glassy' state that is an intermediate dynamical state between a solid and liquid, we probe the validity and significance of this proposed analogy. In particular, atomistic simulations of ubiquitin in solution based on CHARMM36 force field and pre-melted Ni NPs (Voter-Chen Embedded Atom Method potential) indicate a common dynamic heterogeneity, along with other features of glass-forming (GF) liquids such as collective atomic motion in the form of string-like atomic displacements, potential energy fluctuations and particle displacements with long range correlations ('colored' or 'pink' noise), and particle displacement events having a power law scaling in magnitude, as found in earthquakes. On the other hand, we find the dynamics of ubiquitin to be even more like a polycrystalline material in which the α-helix and β-sheet regions of the protein are similar to crystal grains so that the string-like collective atomic motion is concentrated in regions between the α-helix and β-sheet domains.

Protein dynamics and function from solution state NMR spectroscopy

It is well-established that dynamics are central to protein function; their importance is implicitly acknowledged in the principles of the Monod, Wyman and Changeux model of binding cooperativity, which was originally proposed in 1965. Nowadays the concept of protein dynamics is formulated in terms of the energy landscape theory, which can be used to understand protein folding and conformational changes in proteins. Because protein dynamics are so important, a key to understanding protein function at the molecular level is to design experiments that allow their quantitative analysis. Nuclear magnetic resonance (NMR) spectroscopy is uniquely suited for this purpose because major advances in theory, hardware, and experimental methods have made it possible to characterize protein dynamics at an unprecedented level of detail. Unique features of NMR include the ability to quantify dynamics (i) under equilibrium conditions without external perturbations, (ii) using many probes simultaneously, and (iii) over large time intervals. Here we review NMR techniques for quantifying protein dynamics on fast (ps-ns), slow (μs-ms), and very slow (s-min) time scales. These techniques are discussed with reference to some major discoveries in protein science that have been made possible by NMR spectroscopy.

ClustENM: ENM-Based Sampling of Essential Conformational Space at Full Atomic Resolution

Accurate sampling of conformational space and, in particular, the transitions between functional substates has been a challenge in molecular dynamic (MD) simulations of large biomolecular systems. We developed an Elastic Network Model (ENM)-based computational method, ClustENM, for sampling large conformational changes of biomolecules with various sizes and oligomerization states. ClustENM is an iterative method that combines ENM with energy minimization and clustering steps. It is an unbiased technique, which requires only an initial structure as input, and no information about the target conformation. To test the performance of ClustENM, we applied it to six biomolecular systems: adenylate kinase (AK), calmodulin, p38 MAP kinase, HIV-1 reverse transcriptase (RT), triosephosphate isomerase (TIM), and the 70S ribosomal complex. The generated ensembles of conformers determined at atomic resolution show good agreement with experimental data (979 structures resolved by X-ray and/or NMR) and encompass the subspaces covered in independent MD simulations for TIM, p38, and RT. ClustENM emerges as a computationally efficient tool for characterizing the conformational space of large systems at atomic detail, in addition to generating a representative ensemble of conformers that can be advantageously used in simulating substrate/ligand-binding events.

Connecting Active-Site Loop Conformations and Catalysis in Triosephosphate Isomerase: Insights from a Rare Variation at Residue 96 in the Plasmodial Enzyme

Despite extensive research into triosephosphate isomerases (TIMs), there exists a gap in understanding of the remarkable conjunction between catalytic loop-6 (residues 166-176) movement and the conformational flip of Glu165 (catalytic base) upon substrate binding that primes the active site for efficient catalysis. The overwhelming occurrence of serine at position 96 (98% of the 6277 unique TIM sequences), spatially proximal to E165 and the loop-6 residues, raises questions about its role in catalysis. Notably, Plasmodium falciparum TIM has an extremely rare residue--phenylalanine--at this position whereas, curiously, the mutant F96S was catalytically defective. We have obtained insights into the influence of residue 96 on the loop-6 conformational flip and E165 positioning by combining kinetic and structural studies on the PfTIM F96 mutants F96Y, F96A, F96S/S73A, and F96S/L167V with sequence conservation analysis and comparative analysis of the available apo and holo structures of the enzyme from diverse organisms.

Using NMR spectroscopy to elucidate the role of molecular motions in enzyme function

Conformational motions play an essential role in enzyme function, often facilitating the formation of enzyme-substrate complexes and/or product release. Although considerable debate remains regarding the role of molecular motions in the conversion of enzymatic substrates to products, numerous examples have found motions to be crucial for optimization of enzyme scaffolds, effective substrate binding, and product dissociation. Conformational fluctuations are often rate-limiting to enzyme catalysis, primarily through product release, with the chemical reaction occurring much more quickly. As a result, the direct involvement of motions at various stages along the enzyme reaction coordinate remains largely unknown and untested. In the following review, we describe the use of solution NMR techniques designed to probe various timescales of molecular motions and detail examples in which motions play a role in propagating catalytic effects from the active site and directly participate in essential aspects of enzyme function.

Utilization of Substrate Intrinsic Binding Energy for Conformational Change and Catalytic Function in Phosphoenolpyruvate Carboxykinase

Phosphoenolpyruvate carboxykinase (PEPCK) is an essential metabolic enzyme operating in the gluconeogenesis and glyceroneogenesis pathways. Previous work has demonstrated that the enzyme cycles between a catalytically inactive open state and a catalytically active closed state. The transition of the enzyme between these states requires the transition of several active site loops to shift from mobile, disordered structural elements to stable ordered states. The mechanism by which these disorder-order transitions are coupled to the ligation state of the active site however is not fully understood. To further investigate the mechanisms by which the mobility of the active site loops is coupled to enzymatic function and the transitioning of the enzyme between the two conformational states, we have conducted structural and functional studies of point mutants of E89. E89 is a proposed key member of the interaction network of mobile elements as it resides in the R-loop region of the enzyme active site. These new data demonstrate the importance of the R-loop in coordinating interactions between substrates at the OAA/PEP binding site and the mobile R- and Ω-loop domains. In turn, the studies more generally demonstrate the mechanisms by which the intrinsic ligand binding energy can be utilized in catalysis to drive unfavorable conformational changes, changes that are subsequently required for both optimal catalytic activity and fidelity.

Role of Loop-Clamping Side Chains in Catalysis by Triosephosphate Isomerase

The side chains of Y208 and S211 from loop 7 of triosephosphate isomerase (TIM) form hydrogen bonds to backbone amides and carbonyls from loop 6 to stabilize the caged enzyme–substrate complex. The effect of seven mutations [Y208T, Y208S, Y208A, Y208F, S211G, S211A, Y208T/S211G] on the kinetic parameters for TIM catalyzed reactions of the whole substrates dihydroxyacetone phosphate and d-glyceraldehyde 3-phosphate [(k_cat/K_m)_GAP and (k_cat/K_m)_DHAP] and of the substrate pieces glycolaldehyde and phosphite dianion (k_cat/K_HPiK_GA) are reported. The linear logarithmic correlation between these kinetic parameters, with slope of 1.04 ± 0.03, shows that most mutations of TIM result in an identical change in the activation barriers for the catalyzed reactions of whole substrate and substrate pieces, so that the transition states for these reactions are stabilized by similar interactions with the protein catalyst. The second linear logarithmic correlation [slope = 0.53 ± 0.16] between k_cat for isomerization of GAP and K_d^⧧ for phosphite dianion binding to the transition state for wildtype and many mutant TIM-catalyzed reactions of substrate pieces shows that ca. 50% of the wildtype TIM dianion binding energy, eliminated by these mutations, is expressed at the wildtype Michaelis complex, and ca. 50% is only expressed at the wildtype transition state. Negative deviations from this correlation are observed when the mutation results in a decrease in enzyme reactivity at the catalytic site. The main effect of Y208T, Y208S, and Y208A mutations is to cause a reduction in the total intrinsic dianion binding energy, but the effect of Y208F extends to the catalytic site.

Dynamic allostery governs cyclophilin A-HIV capsid interplay

Host factor protein Cyclophilin A (CypA) regulates HIV-1 viral infectivity through direct interactions with the viral capsid, by an unknown mechanism. CypA can either promote or inhibit viral infection, depending on host cell type and HIV-1 capsid (CA) protein sequence. We have examined the role of conformational dynamics on the nanosecond to millisecond timescale in HIV-1 CA assemblies in the escape from CypA dependence, by magic-angle spinning (MAS) NMR and molecular dynamics (MD). Through the analysis of backbone (1)H-(15)N and (1)H-(13)C dipolar tensors and peak intensities from 3D MAS NMR spectra of wild-type and the A92E and G94D CypA escape mutants, we demonstrate that assembled CA is dynamic, particularly in loop regions. The CypA loop in assembled wild-type CA from two strains exhibits unprecedented mobility on the nanosecond to microsecond timescales, and the experimental NMR dipolar order parameters are in quantitative agreement with those calculated from MD trajectories. Remarkably, the CypA loop dynamics of wild-type CA HXB2 assembly is significantly attenuated upon CypA binding, and the dynamics profiles of the A92E and G94D CypA escape mutants closely resemble that of wild-type CA assembly in complex with CypA. These results suggest that CypA loop dynamics is a determining factor in HIV-1's escape from CypA dependence.

Unraveling the complexity of protein backbone dynamics with combined (13)C and (15)N solid-state NMR relaxation measurements

Typically, protein dynamics involve a complex hierarchy of motions occurring on different time scales between conformations separated by a range of different energy barriers. NMR relaxation can in principle provide a site-specific picture of both the time scales and amplitudes of these motions, but independent relaxation rates sensitive to fluctuations in different time scale ranges are required to obtain a faithful representation of the underlying dynamic complexity. This is especially pertinent for relaxation measurements in the solid state, which report on dynamics in a broader window of time scales by more than 3 orders of magnitudes compared to solution NMR relaxation. To aid in unraveling the intricacies of biomolecular dynamics we introduce (13)C spin-lattice relaxation in the rotating frame (R1ρ) as a probe of backbone nanosecond-microsecond motions in proteins in the solid state. We present measurements of (13)C'R1ρ rates in fully protonated crystalline protein GB1 at 600 and 850 MHz (1)H Larmor frequencies and compare them to (13)C'R1, (15)N R1 and R1ρ measured under the same conditions. The addition of carbon relaxation data to the model free analysis of nitrogen relaxation data leads to greatly improved characterization of time scales of protein backbone motions, minimizing the occurrence of fitting artifacts that may be present when (15)N data is used alone. We also discuss how internal motions characterized by different time scales contribute to (15)N and (13)C relaxation rates in the solid state and solution state, leading to fundamental differences between them, as well as phenomena such as underestimation of picosecond-range motions in the solid state and nanosecond-range motions in solution.

How an Inhibitor Bound to Subunit Interface Alters Triosephosphate Isomerase Dynamics

The tunnel region at triosephosphate isomerase (TIM)'s dimer interface, distant from its catalytic site, is a target site for certain benzothiazole derivatives that inhibit TIM's catalytic activity in Trypanosoma cruzi, the parasite that causes Chagas disease. We performed multiple 100-ns molecular-dynamics (MD) simulations and elastic network modeling (ENM) on both apo and complex structures to shed light on the still unclear inhibitory mechanism of one such inhibitor, named bt10. Within the time frame of our MD simulations, we observed stabilization of aromatic clusters at the dimer interface and enhancement of intersubunit hydrogen bonds in the presence of bt10, which point to an allosteric effect rather than destabilization of the dimeric structure. The collective dynamics dictated by the topology of TIM is known to facilitate the closure of its catalytic loop over the active site that is critical for substrate entrance and product release. We incorporated the ligand's effect on vibrational dynamics by applying mixed coarse-grained ENM to each one of 54,000 MD snapshots. Using this computationally efficient technique, we observed altered collective modes and positive shifts in eigenvalues due to the constraining effect of bt10 binding. Accordingly, we observed allosteric changes in the catalytic loop's dynamics, flexibility, and correlations, as well as the solvent exposure of catalytic residues. A newly (to our knowledge) introduced technique that performs residue-based ENM scanning of TIM revealed the tunnel region as a key binding site that can alter global dynamics of the enzyme.

Triosephosphate isomerase I170V alters catalytic site, enhances stability and induces pathology in a Drosophila model of TPI deficiency

Triosephosphate isomerase (TPI) is a glycolytic enzyme which homodimerizes for full catalytic activity. Mutations of the TPI gene elicit a disease known as TPI Deficiency, a glycolytic enzymopathy noted for its unique severity of neurological symptoms. Evidence suggests that TPI Deficiency pathogenesis may be due to conformational changes of the protein, likely affecting dimerization and protein stability. In this report, we genetically and physically characterize a human disease-associated TPI mutation caused by an I170V substitution. Human TPI(I170V) elicits behavioral abnormalities in Drosophila. An examination of hTPI(I170V) enzyme kinetics revealed this substitution reduced catalytic turnover, while assessments of thermal stability demonstrated an increase in enzyme stability. The crystal structure of the homodimeric I170V mutant reveals changes in the geometry of critical residues within the catalytic pocket. Collectively these data reveal new observations of the structural and kinetic determinants of TPI Deficiency pathology, providing new insights into disease pathogenesis.

Protein dynamics control the progression and efficiency of the catalytic reaction cycle of the Escherichia coli DNA-repair enzyme AlkB

A central goal of enzymology is to understand the physicochemical mechanisms that enable proteins to catalyze complex chemical reactions with high efficiency. Recent methodological advances enable the contribution of protein dynamics to enzyme efficiency to be explored more deeply. Here, we utilize enzymological and biophysical studies, including NMR measurements of conformational dynamics, to develop a quantitative mechanistic scheme for the DNA repair enzyme AlkB. Like other iron/2-oxoglutarate-dependent dioxygenases, AlkB employs a two-step mechanism in which oxidation of 2-oxoglutarate generates a highly reactive enzyme-bound oxyferryl intermediate that, in the case of AlkB, slowly hydroxylates an alkylated nucleobase. Our results demonstrate that a microsecond-to-millisecond time scale conformational transition facilitates the proper sequential order of substrate binding to AlkB. Mutations altering the dynamics of this transition allow generation of the oxyferryl intermediate but promote its premature quenching by solvent, which uncouples 2-oxoglutarate turnover from nucleobase oxidation. Therefore, efficient catalysis by AlkB depends upon the dynamics of a specific conformational transition, establishing another paradigm for the control of enzyme function by protein dynamics.

Dynamics in the solid-state: perspectives for the investigation of amyloid aggregates, membrane proteins and soluble protein complexes

Aggregates formed by amyloidogenic peptides and proteins and reconstituted membrane protein preparations differ significantly in terms of the spectral quality that they display in solid-state NMR experiments. Structural heterogeneity and dynamics can both in principle account for that observation. This perspectives article aims to point out challenges and limitations, but also potential opportunities in the investigation of these systems.

Mechanistic Imperatives for Deprotonation of Carbon Catalyzed by Triosephosphate Isomerase: Enzyme-Activation by Phosphite Dianion

The mechanistic imperatives for catalysis of deprotonation of α-carbonyl carbon by triosephosphate isomerase (TIM) are discussed. There is a strong imperative to reduce the large thermodynamic barrier for deprotonation of carbon to form an enediolate reaction intermediate; and, a strong imperative for specificity in the expression of the intrinsic phosphodianion binding energy at the transition state for the enzyme-catalyzed reaction. Binding energies of 2 and 6 kcal/mol, respectively, have been determined for formation of phosphite dianion complexes to TIM and to the transition state for TIM-catalyzed deprotonation of the truncated substrate glycolaldehyde [T. L. Amyes, J. P. Richard, Biochemistry2007, 46, 5841]. We propose that the phosphite dianion binding energy, which is specifically expressed at the transition state complex, is utilized to stabilize a rare catalytically active loop-closed form of TIM. The results of experiments to probe the role of the side chains of Ile172 and Leu232 in activating the loop-closed form of TIM for catalysis of substrate deprotonation are discussed. Evidence is presented that the hydrophobic side chain of Ile172 assists in activating TIM for catalysis of substrate deprotonation through an enhancement of the basicity of the carboxylate side-chain of Glu167. Our experiments link the two imperatives for TIM-catalyzed deprotonation of carbon by providing evidence that the phosphodianion binding energy is utilized to drive an enzyme conformational change, which results in a reduction in the thermodynamic barrier to deprotonation of the carbon acid substrate at TIM compared with the barrier for deprotonation in water. The effects of a P168A mutation on the kinetic parameters for the reactions of whole and truncated substrates are discussed.

High-resolution crystal structure and redox properties of chloroplastic triosephosphate isomerase from Chlamydomonas reinhardtii

Triosephosphate isomerase (TPI) catalyzes the interconversion of glyceraldehyde-3-phosphate to dihydroxyacetone phosphate. Photosynthetic organisms generally contain two isoforms of TPI located in both cytoplasm and chloroplasts. While the cytoplasmic TPI is involved in the glycolysis, the chloroplastic isoform participates in the Calvin-Benson cycle, a key photosynthetic process responsible for carbon fixation. Compared with its cytoplasmic counterpart, the functional features of chloroplastic TPI have been poorly investigated and its three-dimensional structure has not been solved. Recently, several studies proposed TPI as a potential target of different redox modifications including dithiol/disulfide interchanges, glutathionylation, and nitrosylation. However, neither the effects on protein activity nor the molecular mechanisms underlying these redox modifications have been investigated. Here, we have produced recombinantly and purified TPI from the unicellular green alga Chlamydomonas reinhardtii (Cr). The biochemical properties of the enzyme were delineated and its crystallographic structure was determined at a resolution of 1.1 Å. CrTPI is a homodimer with subunits containing the typical (β/α)8-barrel fold. Although no evidence for TRX regulation was obtained, CrTPI was found to undergo glutathionylation by oxidized glutathione and trans-nitrosylation by nitrosoglutathione, confirming its sensitivity to multiple redox modifications.

What's in your buffer? Solute altered millisecond motions detected by solution NMR

To date, little work has been conducted on the relationship between solute and buffer molecules and conformational exchange motion in enzymes. This study uses solution NMR to examine the effects of phosphate, sulfate, and acetate in comparison to MES- and HEPES-buffered references on the chemical shift perturbation and millisecond, chemical, or conformational exchange motions in the enzyme ribonuclease A (RNase A), triosephosphate isomerase (TIM) and HisF. The results indicate that addition of these solutes has a small effect on (1)H and (15)N chemical shifts for RNase A and TIM but a significant effect for HisF. For RNase A and TIM, Carr-Purcell-Meiboom-Gill relaxation dispersion experiments, however, show significant solute-dependent changes in conformational exchange motions. Some residues show loss of millisecond motions relative to the reference sample upon addition of solute, whereas others experience an enhancement. Comparison of exchange parameters obtained from fits of dispersion data indicates changes in either or both equilibrium populations and chemical shifts between conformations. Furthermore, the exchange kinetics are altered in many cases. The results demonstrate that common solute molecules can alter observed enzyme millisecond motions and play a more active role than what is routinely believed.

Preparation of uniformly isotope labeled KcsA for solid state NMR: expression, purification, reconstitution into liposomes and functional assay

We report the expression, purification, liposome reconstitution and functional validation of uniformly (13)C and (15)N isotope labeled KcsA, a bacterial potassium channel that has high homology with mammalian channels, for solid-state NMR studies. The expression and purification is optimized for an average yield of ∼35-40mg/L of M9 media in a time-efficient way. The protein purity is confirmed by gel electrophoresis and the protein concentration is quantified by UV-vis absorption spectroscopy. Protocols to efficiently reconstitute KcsA into liposomes are also presented. The presence of liposomes is confirmed by cryo-electron microscopy images and the effect of magic angle spinning on liposome packing is shown. High-resolution solid-state NMR spectra of uniformly isotope labeled KcsA in these liposomes reveal that our protocol yields to a very homogenous KcsA sample with high signal to noise and several well-resolved residues in NMR spectra. Electrophysiology of our samples before and after solid-state NMR show that channel function and selectivity remain intact after the solid-state NMR.

Temporally overlapped but uncoupled motions in dihydrofolate reductase catalysis

Temporal correlations between protein motions and enzymatic reactions are often interpreted as evidence for catalytically important motions. Using dihydrofolate reductase as a model system, we show that there are many protein motions that temporally overlapped with the chemical reaction, and yet they do not exhibit the same kinetic behaviors (KIE and pH dependence) as the catalyzed chemical reaction. Thus, despite the temporal correlation, these motions are not directly coupled to the chemical transformation, and they might represent a different part of the catalytic cycle or simply be the product of the intrinsic flexibility of the protein.

Structural and functional perturbation of Giardia lamblia triosephosphate isomerase by modification of a non-catalytic, non-conserved region

BACKGROUND

We have previously proposed triosephosphate isomerase of Giardia lamblia (GlTIM) as a target for rational drug design against giardiasis, one of the most common parasitic infections in humans. Since the enzyme exists in the parasite and the host, selective inhibition is a major challenge because essential regions that could be considered molecular targets are highly conserved. Previous biochemical evidence showed that chemical modification of the non-conserved non-catalytic cysteine 222 (C222) inactivates specifically GlTIM. The inactivation correlates with the physicochemical properties of the modifying agent: addition of a non-polar, small chemical group at C222 reduces the enzyme activity by one half, whereas negatively charged, large chemical groups cause full inactivation.

RESULTS

In this work we used mutagenesis to extend our understanding of the functional and structural effects triggered by modification of C222. To this end, six GlTIM C222 mutants with side chains having diverse physicochemical characteristics were characterized. We found that the polarity, charge and volume of the side chain in the mutant amino acid differentially alter the activity, the affinity, the stability and the structure of the enzyme. The data show that mutagenesis of C222 mimics the effects of chemical modification. The crystallographic structure of C222D GlTIM shows the disruptive effects of introducing a negative charge at position 222: the mutation perturbs loop 7, a region of the enzyme whose interactions with the catalytic loop 6 are essential for TIM stability, ligand binding and catalysis. The amino acid sequence of TIM in phylogenetic diverse groups indicates that C222 and its surrounding residues are poorly conserved, supporting the proposal that this region is a good target for specific drug design.

CONCLUSIONS

The results demonstrate that it is possible to inhibit species-specifically a ubiquitous, structurally highly conserved enzyme by modification of a non-conserved, non-catalytic residue through long-range perturbation of essential regions.

Ultrafast NMR T1 relaxation measurements: probing molecular properties in real time

The longitudinal relaxation properties of NMR active nuclei carry useful information about the site-specific chemical environments and about the mobility of molecular fragments. Molecular mobility is in turn a key parameter reporting both on stable properties, such as size, as well as on dynamic ones, such as transient interactions and irreversible aggregation. In order to fully investigate the latter, a fast sampling of the relaxation parameters of transiently formed molecular species may be needed. Nevertheless, the acquisition of longitudinal relaxation data is typically slow, being limited by the requirement that the time for which the nucleus relaxes be varied incrementally until a complete build-up curve is generated. Recently, a number of single-shot-inversion-recovery methods have been developed capable of alleviating this need; still, these may be challenged by either spectral resolution restrictions or when coping with very fast relaxing nuclei. Here, we present a new experiment to measure the T1s of multiple nuclear spins that experience fast longitudinal relaxation, while retaining full high-resolution chemical shift information. Good agreement is observed between T1s measured with conventional means and T1s measured using the new technique. The method is applied to the real-time investigation of the reaction between D-xylose and sodium borate, which is in turn elucidated with the aid of ancillary ultrafast and conventional 2D TOCSY measurements.

Ultrafast real-time visualization of active site flexibility of flavoenzyme thymidylate synthase ThyX

In many bacteria the flavoenzyme thymidylate synthase ThyX produces the DNA nucleotide deoxythymidine monophosphate from dUMP, using methylenetetrahydrofolate as carbon donor and NADPH as hydride donor. Because all three substrates bind in close proximity to the catalytic flavin adenine dinucleotide group, substantial flexibility of the ThyX active site has been hypothesized. Using femtosecond time-resolved fluorescence spectroscopy, we have studied the conformational heterogeneity and the conformational interconversion dynamics in real time in ThyX from the hyperthermophilic bacterium Thermotoga maritima. The dynamics of electron transfer to excited flavin adenine dinucleotide from a neighboring tyrosine residue are used as a sensitive probe of the functional dynamics of the active site. The fluorescence decay spanned a full three orders of magnitude, demonstrating a very wide range of conformations. In particular, at physiological temperatures, multiple angstrom cofactor-residue displacements occur on the picoseconds timescale. These experimental findings are supported by molecular dynamics simulations. Binding of the dUMP substrate abolishes this flexibility and stabilizes the active site in a configuration where dUMP closely interacts with the flavin cofactor and very efficiently quenches fluorescence itself. Our results indicate a dynamic selected-fit mechanism where binding of the first substrate dUMP at high temperature stabilizes the enzyme in a configuration favorable for interaction with the second substrate NADPH, and more generally have important implications for the role of active site flexibility in enzymes interacting with multiple poly-atom substrates and products. Moreover, our data provide the basis for exploring the effect of inhibitor molecules on the active site dynamics of ThyX and other multisubstrate flavoenzymes.

Specificity in transition state binding: the Pauling model revisited

Linus Pauling proposed that the large rate accelerations for enzymes are caused by the high specificity of the protein catalyst for binding the reaction transition state. The observation that stable analogues of the transition states for enzymatic reactions often act as tight-binding inhibitors provided early support for this simple and elegant proposal. We review experimental results that support the proposal that Pauling's model provides a satisfactory explanation for the rate accelerations for many heterolytic enzymatic reactions through high-energy reaction intermediates, such as proton transfer and decarboxylation. Specificity in transition state binding is obtained when the total intrinsic binding energy of the substrate is significantly larger than the binding energy observed at the Michaelis complex. The results of recent studies that aimed to characterize the specificity in binding of the enolate oxygen at the transition state for the 1,3-isomerization reaction catalyzed by ketosteroid isomerase are reviewed. Interactions between pig heart succinyl-coenzyme A:3-oxoacid coenzyme A transferase (SCOT) and the nonreacting portions of coenzyme A (CoA) are responsible for a rate increase of 3 × 10(12)-fold, which is close to the estimated total 5 × 10(13)-fold enzymatic rate acceleration. Studies that partition the interactions between SCOT and CoA into their contributing parts are reviewed. Interactions of the protein with the substrate phosphodianion group provide an ~12 kcal/mol stabilization of the transition state for the reactions catalyzed by triosephosphate isomerase, orotidine 5'-monophosphate decarboxylase, and α-glycerol phosphate dehydrogenase. The interactions of these enzymes with the substrate piece phosphite dianion provide a 6-8 kcal/mol stabilization of the transition state for reaction of the appropriate truncated substrate. Enzyme activation by phosphite dianion reflects the higher dianion affinity for binding to the enzyme-transition state complex compared with that of the free enzyme. Evidence is presented that supports a model in which the binding energy of the phosphite dianion piece, or the phosphodianion group of the whole substrate, is utilized to drive an enzyme conformational change from an inactive open form E(O) to an active closed form E(C), by closure of a phosphodianion gripper loop. Members of the enolase and haloalkanoic acid dehalogenase superfamilies use variable capping domains to interact with nonreacting portions of the substrate and sequester the substrate from interaction with bulk solvent. Interactions of this capping domain with the phenyl group of mandelate have been shown to activate mandelate racemase for catalysis of deprotonation of α-carbonyl carbon. We propose that an important function of these capping domains is to utilize the binding interactions with nonreacting portions of the substrate to activate the enzyme for catalysis.

Hydrogen tunneling links protein dynamics to enzyme catalysis

The relationship between protein dynamics and function is a subject of considerable contemporary interest. Although protein motions are frequently observed during ligand binding and release steps, the contribution of protein motions to the catalysis of bond making/breaking processes is more difficult to probe and verify. Here, we show how the quantum mechanical hydrogen tunneling associated with enzymatic C-H bond cleavage provides a unique window into the necessity of protein dynamics for achieving optimal catalysis. Experimental findings support a hierarchy of thermodynamically equilibrated motions that control the H-donor and -acceptor distance and active-site electrostatics, creating an ensemble of conformations suitable for H-tunneling. A possible extension of this view to methyl transfer and other catalyzed reactions is also presented. The impact of understanding these dynamics on the conceptual framework for enzyme activity, inhibitor/drug design, and biomimetic catalyst design is likely to be substantial.

Kinetic mechanism of indole-3-glycerol phosphate synthase

The (βα)(8)-barrel enzyme indole-3-glycerol phosphate synthase (IGPS) catalyzes the multistep transformation of 1-(o-carboxyphenylamino)-1-deoxyribulose 5-phosphate (CdRP) into indole-3-glycerol phosphate (IGP) in tryptophan biosynthesis. Mutagenesis data and crystal structure analysis of IGPS from Sulfolobus solfataricus (sIGPS) allowed for the formulation of a plausible chemical mechanism of the reaction, and molecular dynamics simulations suggested that flexibility of active site loops might be important for catalysis. Here we developed a method that uses extrinsic fluorophores attached to active site loops to connect the kinetic mechanism of sIGPS to structure and conformational motions. Specifically, we elucidated the kinetic mechanism of sIGPS and correlated individual steps in the mechanism to conformational motions of flexible loops. Pre-steady-state kinetic measurements of CdRP to IGP conversion monitoring changes in intrinsic tryptophan and IGP fluorescence provided a minimal three-step kinetic model in which fast substrate binding and chemical transformation are followed by slow product release. The role of sIGPS loop conformational motion during substrate binding and catalysis was examined via variants that were covalently labeled with fluorescent dyes at the N-terminal extension of the enzyme and mobile active site loop β1α1. Analysis of kinetic data monitoring dye fluorescence revealed a conformational change that follows substrate binding, suggesting an induced-fit-type binding mechanism for the substrate CdRP. Global fitting of all kinetic results obtained with wild-type sIGPS and the labeled variants was best accommodated by a four-step kinetic model. In this model, both the binding of CdRP and its on-enzyme conversion to IGP are accompanied by conformational transitions. The liberation of the product from the active site is the rate-limiting step of the overall reaction. Our results confirm the importance of flexible active loops for substrate binding and catalysis by sIGPS.

The Ω-loop lid domain of phosphoenolpyruvate carboxykinase is essential for catalytic function

Phosphoenolpyruvate carboxykinase (PEPCK) is an essential metabolic enzyme operating in the gluconeogenesis and glyceroneogenesis pathways. Recent studies have demonstrated that the enzyme contains a mobile active site lid domain that undergoes a transition between an open, disorded conformation and a closed, ordered conformation as the enzyme progresses through the catalytic cycle. The understanding of how this mobile domain functions in catalysis is incomplete. Previous studies showed that the closure of the lid domain stabilizes the reaction intermediate and protects the reactive intermediate from spurious protonation and thus contributes to the fidelity of the enzyme. To more fully investigate the roles of the lid domain in PEPCK function, we introduced three mutations that replaced the 11-residue lid domain with one, two, and three glycine residues. Kinetic analysis of the mutant enzymes demonstrates that none of the enzyme constructs exhibit any measurable kinetic activity, resulting in a decrease in the catalytic parameters of at least 10(6). Structural characterization of the mutants in complexes representing the catalytic cycle suggests that the inactivity is due to a role for the lid domain in the formation of the fully closed state of the enzyme that is required for catalytic function. In the absence of the lid domain, the enzyme is unable to achieve the fully closed state and is rendered inactive despite possessing all of the residues and substrates required for catalytic function. This work demonstrates how enzyme catalytic function can be abolished through the alteration of conformational equilibria despite all the elements required for chemical conversion of substrates to products remaining intact.

Preparation of RNA samples with narrow line widths for solid state NMR investigations

Solid state NMR can provide detailed structural and dynamic information on biological systems that cannot be studied under solution conditions, and can investigate motions which occur with rates that cannot be fully studied by solution NMR. This approach has successfully been used to study proteins, but the application of multidimensional solid state NMR to RNA has been limited because reported line widths have been too broad to execute most multidimensional experiments successfully. A reliable method to generate spectra with narrow line widths is necessary to apply the full range of solid state NMR spectroscopic approaches to RNA. Using the HIV-1 transactivation response (TAR) RNA as a model, we present an approach based on precipitation with polyethylene glycol that improves the line width of (13)C signals in TAR from >6 ppm to about 1 ppm, making solid state 2D NMR studies of selectively enriched RNAs feasible at ambient temperature.

Coupling between catalytic loop motions and enzyme global dynamics

Catalytic loop motions facilitate substrate recognition and binding in many enzymes. While these motions appear to be highly flexible, their functional significance suggests that structure-encoded preferences may play a role in selecting particular mechanisms of motions. We performed an extensive study on a set of enzymes to assess whether the collective/global dynamics, as predicted by elastic network models (ENMs), facilitates or even defines the local motions undergone by functional loops. Our dataset includes a total of 117 crystal structures for ten enzymes of different sizes and oligomerization states. Each enzyme contains a specific functional/catalytic loop (10–21 residues long) that closes over the active site during catalysis. Principal component analysis (PCA) of the available crystal structures (including apo and ligand-bound forms) for each enzyme revealed the dominant conformational changes taking place in these loops upon substrate binding. These experimentally observed loop reconfigurations are shown to be predominantly driven by energetically favored modes of motion intrinsically accessible to the enzyme in the absence of its substrate. The analysis suggests that robust global modes cooperatively defined by the overall enzyme architecture also entail local components that assist in suitable opening/closure of the catalytic loop over the active site.

Protein loops have critical roles in ligand binding and catalysis. An unresolved issue in this context is the extent to which the intrinsic dynamics of proteins predispose loops to perform their molecular function. In this work, we (i) critically examine the structural changes undergone by functional/catalytic loops based on a set of enzyme crystal structures in the presence/absence of a ligand, and (ii) examine to what extent those motions are correlated with, or driven by, the global modes that are predictable using simplified, physics-based models. Using a dataset of 117 structures for ten enzymes of different sizes and oligomerization states, we show that the collective modes defined by the protein topology favor loop rearrangements in reasonable agreement with those experimentally observed upon activation. These results suggest that simple but robust motions encoded by the entire architecture, not the local binding site only, assist in binding of the ligand, positioning of the catalytic loop, and/or sequestration of the catalytic site, which in turn, enable efficient catalysis.

Loop residues and catalysis in OMP synthase

Residue-to-alanine mutations and a two-amino acid deletion have been made in the highly conserved catalytic loop (residues 100-109) of Salmonella typhimurium OMP synthase (orotate phosphoribosyltransferase, EC 2.4.2.10). As described previously, the K103A mutant enzyme exhibited a 10(4)-fold decrease in k(cat)/K(M) for PRPP; the K100A enzyme suffered a 50-fold decrease. Alanine mutations at His105 and Glu107 produced 40- and 7-fold decreases in k(cat)/K(M), respectively, and E101A, D104A, and G106A were slightly faster than the wild-type (WT) in terms of k(cat), with minor effects on k(cat)/K(M). Equilibrium binding of OMP or PRPP in binary complexes was affected little by loop mutation, suggesting that the energetics of ground-state binding have little contribution from the catalytic loop, or that a favorable binding energy is offset by costs of loop reorganization. Pre-steady-state kinetics for mutants showed that K103A and E107A had lost the burst of product formation in each direction that indicated rapid on-enzyme chemistry for WT, but that the burst was retained by H105A. Δ102Δ106, a loop-shortened enzyme with Ala102 and Gly106 deleted, showed a 10(4)-fold reduction of k(cat) but almost unaltered K(D) values for all four substrate molecules. The 20% (i.e., 1.20) intrinsic [1'-(3)H]OMP kinetic isotope effect (KIE) for WT is masked because of high forward and reverse commitment factors. K103A failed to express intrinsic KIEs fully (1.095 ± 0.013). In contrast, H105A, which has a smaller catalytic lesion, gave a [1'-(3)H]OMP KIE of 1.21 ± 0.0005, and E107A (1.179 ± 0.0049) also gave high values. These results are interpreted in the context of the X-ray structure of the complete substrate complex for the enzyme [Grubmeyer, C., Hansen, M. R., Fedorov, A. A., and Almo, S. C. (2012) Biochemistry 51 (preceding paper in this issue, DOI 10.1021/bi300083p )]. The full expression of KIEs by H105A and E107A may result from a less secure closure of the catalytic loop. The lower level of expression of the KIE by K103A suggests that in these mutant proteins the major barrier to catalysis is successful closure of the catalytic loop, which when closed, produces rapid and reversible catalysis.

Cooperativity in monomeric enzymes with single ligand-binding sites

Cooperativity is widespread in biology. It empowers a variety of regulatory mechanisms and impacts both the kinetic and thermodynamic properties of macromolecular systems. Traditionally, cooperativity is viewed as requiring the participation of multiple, spatially distinct binding sites that communicate via ligand-induced structural rearrangements; however, cooperativity requires neither multiple ligand binding events nor multimeric assemblies. An underappreciated manifestation of cooperativity has been observed in the non-Michaelis-Menten kinetic response of certain monomeric enzymes that possess only a single ligand-binding site. In this review, we present an overview of kinetic cooperativity in monomeric enzymes. We discuss the primary mechanisms postulated to give rise to monomeric cooperativity and highlight modern experimental methods that could offer new insights into the nature of this phenomenon. We conclude with an updated list of single subunit enzymes that are suspected of displaying cooperativity, and a discussion of the biological significance of this unique kinetic response.

A normal mode-based geometric simulation approach for exploring biologically relevant conformational transitions in proteins

A three-step approach for multiscale modeling of protein conformational changes is presented that incorporates information about preferred directions of protein motions into a geometric simulation algorithm. The first two steps are based on a rigid cluster normal-mode analysis (RCNMA). Low-frequency normal modes are used in the third step (NMSim) to extend the recently introduced idea of constrained geometric simulations of diffusive motions in proteins by biasing backbone motions of the protein, whereas side-chain motions are biased toward favorable rotamer states. The generated structures are iteratively corrected regarding steric clashes and stereochemical constraint violations. The approach allows performing three simulation types: unbiased exploration of conformational space; pathway generation by a targeted simulation; and radius of gyration-guided simulation. When applied to a data set of proteins with experimentally observed conformational changes, conformational variabilities are reproduced very well for 4 out of 5 proteins that show domain motions, with correlation coefficients r > 0.70 and as high as r = 0.92 in the case of adenylate kinase. In 7 out of 8 cases, NMSim simulations starting from unbound structures are able to sample conformations that are similar (root-mean-square deviation = 1.0-3.1 Å) to ligand bound conformations. An NMSim generated pathway of conformational change of adenylate kinase correctly describes the sequence of domain closing. The NMSim approach is a computationally efficient alternative to molecular dynamics simulations for conformational sampling of proteins. The generated conformations and pathways of conformational transitions can serve as input to docking approaches or as starting points for more sophisticated sampling techniques.