Cofactor-Mediated Conformational Dynamics Promote Product Release From Escherichia coli Dihydrofolate Reductase via an Allosteric Pathway

What is allosteric regulation? Exploring the exceptions that prove the rule!

"Allosteric" was first introduced to mean the other site (i.e., a site distinct from the active or orthosteric site), an adjective for "regulation" to imply a regulatory outcome resulting from ligand binding at another site. That original idea outlines a system with two ligand-binding events at two distinct locations on a macromolecule (originally a protein system), which defines a four-state energy cycle. An allosteric energy cycle provides a quantifiable allosteric coupling constant and focuses our attention on the unique properties of the four equilibrated protein complexes that constitute the energy cycle. Because many observed phenomena have been referenced as "allosteric regulation" in the literature, the goal of this work is to use literature examples to explore which systems are and are not consistent with the two-ligand thermodynamic energy cycle-based definition of allosteric regulation. We emphasize the need for consistent language so comparisons can be made among the ever-increasing number of allosteric systems. Building on the mutually exclusive natures of an energy cycle definition of allosteric regulation versus classic two-state models, we conclude our discussion by outlining how the often-proposed Rube-Goldberg-like mechanisms are likely inconsistent with an energy cycle definition of allosteric regulation.

Perturbative diffraction methods resolve a conformational switch that facilitates a two-step enzymatic mechanism

Enzymes catalyze biochemical reactions through precise positioning of substrates, cofactors, and amino acids to modulate the transition-state free energy. However, the role of conformational dynamics remains poorly understood due to poor experimental access. This shortcoming is evident with Escherichia coli dihydrofolate reductase (DHFR), a model system for the role of protein dynamics in catalysis, for which it is unknown how the enzyme regulates the different active site environments required to facilitate proton and hydride transfer. Here, we describe ligand-, temperature-, and electric-field-based perturbations during X-ray diffraction experiments to map the conformational dynamics of the Michaelis complex of DHFR. We resolve coupled global and local motions and find that these motions are engaged by the protonated substrate to promote efficient catalysis. This result suggests a fundamental design principle for multistep enzymes in which pre-existing dynamics enable intermediates to drive rapid electrostatic reorganization to facilitate subsequent chemical steps.

Modifying the amino acids in conformational motion pathway of the α-amylase of Geobacillus stearothermophilus improved its activity and stability

Amino acids along the conformational motion pathway of the enzyme molecule correlated to its flexibility and rigidity. To enhance the enzyme activity and thermal stability, the motion pathway of Geobacillus stearothermophilus α-amylase has been identified and molecularly modified by using the neural relational inference model and deep learning tool. The significant differences in substrate specificity, enzymatic kinetics, optimal temperature, and thermal stability were observed among the mutants with modified amino acids along the pathway. Mutants especially the P44E demonstrated enhanced hydrolytic activity and catalytic efficiency (kcat/KM) than the wild-type enzyme to 95.0% and 93.8% respectively, with the optimum temperature increased to 90°C. This mutation from proline to glutamic acid has increased the number and the radius of the bottleneck of the channels, which might facilitate transporting large starch substrates into the enzyme. The mutation could also optimize the hydrogen bonding network of the catalytic center, and diminish the spatial hindering to the substrate entry and exit from the catalytic center.

New Insights into the Cooperativity and Dynamics of Dimeric Enzymes

A survey of protein databases indicates that the majority of enzymes exist in oligomeric forms, with about half of those found in the UniProt database being homodimeric. Understanding why many enzymes are in their dimeric form is imperative. Recent developments in experimental and computational techniques have allowed for a deeper comprehension of the cooperative interactions between the subunits of dimeric enzymes. This review aims to succinctly summarize these recent advancements by providing an overview of experimental and theoretical methods, as well as an understanding of cooperativity in substrate binding and the molecular mechanisms of cooperative catalysis within homodimeric enzymes. Focus is set upon the beneficial effects of dimerization and cooperative catalysis. These advancements not only provide essential case studies and theoretical support for comprehending dimeric enzyme catalysis but also serve as a foundation for designing highly efficient catalysts, such as dimeric organic catalysts. Moreover, these developments have significant implications for drug design, as exemplified by Paxlovid, which was designed for the homodimeric main protease of SARS-CoV-2.

Resolving conformational changes that mediate a two-step catalytic mechanism in a model enzyme

Enzymes catalyze biochemical reactions through precise positioning of substrates, cofactors, and amino acids to modulate the transition-state free energy. However, the role of conformational dynamics remains poorly understood due to lack of experimental access. This shortcoming is evident with E. coli dihydrofolate reductase (DHFR), a model system for the role of protein dynamics in catalysis, for which it is unknown how the enzyme regulates the different active site environments required to facilitate proton and hydride transfer. Here, we present ligand-, temperature-, and electric-field-based perturbations during X-ray diffraction experiments that enable identification of coupled conformational changes in DHFR. We identify a global hinge motion and local networks of structural rearrangements that are engaged by substrate protonation to regulate solvent access and promote efficient catalysis. The resulting mechanism shows that DHFR's two-step catalytic mechanism is guided by a dynamic free energy landscape responsive to the state of the substrate.

Engineering the Active Site Lid Dynamics to Improve the Catalytic Efficiency of Yeast Cytosine Deaminase

Conformational dynamics is important for enzyme catalysis. However, engineering dynamics to achieve a higher catalytic efficiency is still challenging. In this work, we develop a new strategy to improve the activity of yeast cytosine deaminase (yCD) by engineering its conformational dynamics. Specifically, we increase the dynamics of the yCD C-terminal helix, an active site lid that controls the product release. The C-terminal is extended by a dynamical single α-helix (SAH), which improves the product release rate by up to ~8-fold, and the overall catalytic rate kcat by up to ~2-fold. It is also shown that the kcat increase is due to the favorable activation entropy change. The NMR H/D exchange data indicate that the conformational dynamics of the transition state analog complex increases as the helix is extended, elucidating the origin of the enhanced catalytic entropy. This study highlights a novel dynamics engineering strategy that can accelerate the overall catalysis through the entropy-driven mechanism.

Single-Molecule Sampling of Dihydrofolate Reductase Shows Kinetic Pauses and an Endosteric Effect Linked to Catalysis

The ability to sample multiple reactions on the same single enzyme is important to link rare intermediates with catalysis and to unravel the role of conformational changes. Despite decades of efforts, however, the single-molecule characterization of nonfluorogenic enzymes during multiple catalytic turnovers has been elusive. Here, we show that nanopore currents allow sampling the dynamic exchange between five structural intermediates during E. coli dihydrofolate reductase (DHFR) catalysis. We found that an endosteric effect promotes the binding of the substrate to the enzyme with a specific hierarchy. The chemical step then switched the enzyme from the closed to the occluded conformation, which in turn promotes the release of the reduced cofactor NADP⁺. Unexpectedly, only a few reactive complexes lead to catalysis. Furthermore, second-long catalytic pauses were observed, possibly reflecting an off-path conformation generated during the reaction. Finally, the free energy from multiple cofactor binding events were required to release the product and switch DHFR back to the reactive conformer. This catalytic fueled concerted mechanism is likely to have evolved to improve the catalytic efficiency of DHFR under the high concentrations of NADP⁺ in E. coli and might be a general feature for complex enzymatic reactions where the binding and release of the products must be tightly controlled.

The role of conformational dynamics on the activity of polymer-conjugated CalB in organic solvents

A perennial interest in enzyme catalysis has been expanding its applicability from aqueous phase where enzymes are naturally evolved to organic solvents in which the majority of industrial chemical synthesis...

Role of Active Site Loop Dynamics in Mediating Ligand Release from E. coli Dihydrofolate Reductase

Conformational fluctuations from ground-state to sparsely populated but functionally important excited states play a key role in enzyme catalysis. For Escherichia coli dihydrofolate reductase (DHFR), the release of the product tetrahydrofolate (THF) and oxidized cofactor NADP⁺ occurs through exchange between closed and occluded conformations of the Met20 loop. A "dynamic knockout" mutant of E. coli DHFR, where the E. coli sequence in the Met20 loop is replaced by the human sequence (N23PP/S148A), models human DHFR and is incapable of accessing the occluded conformation. ¹H and ¹⁵N CPMG relaxation dispersion analysis for the ternary product complex of the mutant enzyme with NADP⁺ and the product analogue 5,10-dideazatetrahydrofolate (ddTHF) (E:ddTHF:NADP⁺) reveals the mechanism by which NADP⁺ is released when the Met20 loop cannot undergo the closed-to-occluded conformational transition. Two excited states were observed: one related to a faster, relatively high-amplitude conformational fluctuation in areas near the active site, associated with the shuttling of the nicotinamide ring of the cofactor out of the active site, and the other to a slower process where ddTHF undergoes small-amplitude motions within the binding site that are consistent with disorder observed in a room-temperature X-ray crystal structure of the N23PP/S148A mutant protein. These motions likely arise due to steric conflict of the pterin ring of ddTHF with the ribose-nicotinamide moiety of NADP⁺ in the closed active site. These studies demonstrate that site-specific kinetic information from relaxation dispersion experiments can provide intimate details of the changes in catalytic mechanism that result from small changes in local amino acid sequence.

Influence of solution ionic strength on the stabilities of M20 loop conformations in apo E. coli dihydrofolate reductase

Interactions among ions and their specific interactions with macromolecular solutes are known to play a central role in biomolecular stability. However, similar effects in the conformational stability of protein loops that play functional roles, such as binding ligands, proteins, and DNA/RNA molecules, remain relatively unexplored. A well-characterized enzyme that has such a functional loop is Escherichia coli dihydrofolate reductase (ecDHFR), whose so-called M20 loop has been observed in three ordered conformations in crystal structures. To explore how solution ionic strengths may affect the M20 loop conformation, we proposed a reaction coordinate that could quantitatively describe the loop conformation and used it to classify the loop conformations in representative ecDHFR x-ray structures crystallized in varying ionic strengths. The Protein Data Bank survey indicates that at ionic strengths (I) below the intracellular ion concentration-derived ionic strength in E. coli (I ≤ 0.237M), the ecDHFR M20 loop tends to adopt open/closed conformations, and rarely an occluded loop state, but when I is >0.237M, the loop tends to adopt closed/occluded conformations. Distance-dependent electrostatic potentials around the most mobile M20 loop region from molecular dynamics simulations of ecDHFR in equilibrated CaCl₂ solutions of varying ionic strengths show that high ionic strengths (I = 0.75/1.5M) can preferentially stabilize the loop in closed/occluded conformations. These results nicely correlate with conformations derived from ecDHFR structures crystallized in varying ionic strengths. Altogether, our results suggest caution in linking M20 loop conformations derived from crystal structures solved at ionic strengths beyond that tolerated by E. coli to the ecDHFR function.

Distal Regions Regulate Dihydrofolate Reductase-Ligand Interactions

Protein motions play a fundamental role in enzyme catalysis and ligand binding. The relationship between protein motion and function has been extensively investigated in the model enzyme dihydrofolate reductase (DHFR). DHFR is an essential enzyme that catalyzes the reduction of dihydrofolate to tetrahydrofolate. Numerous experimental and computational methods have been used to probe the motions of DHFR through the catalytic cycle and to investigate the effect of distal mutations on DHFR motions and ligand binding. These experimental investigations have pushed forward the study of protein motions and their role in protein-ligand interactions. The introduction of mutations distal to the active site has been shown to have profound effects on ligand binding, hydride transfer rates and catalytic efficacy and these changes are captured by enzyme kinetics measurements. Distal mutations have been shown to exert their effects through a network of correlated amino acids and these effects have been investigated by NMR, protein dynamics, and analysis of coupled amino acids. The experimental methods and the findings that are reviewed here have broad implications for our understanding of enzyme mechanisms, ligand binding and for the future design and discovery of enzyme inhibitors.

Dynamic Preference for NADP/H Cofactor Binding/Release in E. coli YqhD Oxidoreductase

YqhD, an E. coli alcohol/aldehyde oxidoreductase, is an enzyme able to produce valuable bio-renewable fuels and fine chemicals from a broad range of starting materials. Herein, we report the first computational solution-phase structure-dynamics analysis of YqhD, shedding light on the effect of oxidized and reduced NADP/H cofactor binding on the conformational dynamics of the biocatalyst using molecular dynamics (MD) simulations. The cofactor oxidation states mainly influence the interdomain cleft region conformations of the YqhD monomers, involved in intricate cofactor binding and release. The ensemble of NADPH-bound monomers has a narrower average interdomain space resulting in more hydrogen bonds and rigid cofactor binding. NADP-bound YqhD fluctuates between open and closed conformations, while it was observed that NADPH-bound YqhD had slower opening/closing dynamics of the cofactor-binding cleft. In the light of enzyme kinetics and structural data, simulation findings have led us to postulate that the frequently sampled open conformation of the cofactor binding cleft with NADP leads to the more facile release of NADP while increased closed conformation sampling during NADPH binding enhances cofactor binding affinity and the aldehyde reductase activity of the enzyme.

Site-Specific Tryptophan Labels Reveal Local Microsecond-Millisecond Motions of Dihydrofolate Reductase

Many enzymes are known to change conformations during their catalytic cycle, but the role of these protein motions is not well understood. Escherichia coli dihydrofolate reductase (DHFR) is a small, flexible enzyme that is often used as a model system for understanding enzyme dynamics. Recently, native tryptophan fluorescence was used as a probe to study micro- to millisecond dynamics of DHFR. Yet, because DHFR has five native tryptophans, the origin of the observed conformational changes could not be assigned to a specific region within the enzyme. Here, we use DHFR mutants, each with a single tryptophan as a probe for temperature jump fluorescence spectroscopy, to further inform our understanding of DHFR dynamics. The equilibrium tryptophan fluorescence of the mutants shows that each tryptophan is in a different environment and that wild-type DHFR fluorescence is not a simple summation of all the individual tryptophan fluorescence signatures due to tryptophan–tryptophan interactions. Additionally, each mutant exhibits a two-phase relaxation profile corresponding to ligand association/dissociation convolved with associated conformational changes and a slow conformational change that is independent of ligand association and dissociation, similar to the wild-type enzyme. However, the relaxation rate of the slow phase depends on the location of the tryptophan within the enzyme, supporting the conclusion that the individual tryptophan fluorescence dynamics do not originate from a single collective motion, but instead report on local motions throughout the enzyme.

Mechanism of allosteric activation of human mRNA cap methyltransferase (RNMT) by RAM: insights from accelerated molecular dynamics simulations

The RNA guanine-N7 methyltransferase (RNMT) in complex with RNMT-activating miniprotein (RAM) catalyses the formation of a N7-methylated guanosine cap structure on the 5' end of nascent RNA polymerase II transcripts. The mRNA cap protects the primary transcript from exonucleases and recruits cap-binding complexes that mediate RNA processing, export and translation. By using microsecond standard and accelerated molecular dynamics simulations, we provide for the first time a detailed molecular mechanism of allosteric regulation of RNMT by RAM. We show that RAM selects the RNMT active site conformations that are optimal for binding of substrates (AdoMet and the cap), thus enhancing their affinity. Furthermore, our results strongly suggest the likely scenario in which the cap binding promotes the subsequent AdoMet binding, consistent with the previously suggested cooperative binding model. By employing the network community analyses, we revealed the underlying long-range allosteric networks and paths that are crucial for allosteric regulation by RAM. Our findings complement and explain previous experimental data on RNMT activity. Moreover, this study provides the most complete description of the cap and AdoMet binding poses and interactions within the enzyme's active site. This information is critical for the drug discovery efforts that consider RNMT as a promising anti-cancer target.

Tuning the "violin" of protein kinases: The role of dynamics-based allostery

The intricacies of allosteric regulation of protein kinases continue to engage the research community. Allostery, or control from a distance, is seen as a fundamental biomolecular mechanism for proteins. From the traditional methods of conformational selection and induced fit, the field has grown to include the role of protein motions in defining a dynamics-based allosteric approach. Harnessing of these continuous motions in the protein to exert allosteric effects can be defined by a "violin" model that focuses on distributions of protein vibrations as opposed to concerted pathways. According to this model, binding of an allosteric modifier causes global redistribution of dynamics in the protein kinase domain that leads to changes in its catalytic properties. This model is consistent with the "entropy-driven allostery" mechanism proposed by Cooper and Dryden in 1984 and does not require, but does not exclude, any major structural changes. We provide an overview of practical implementation of the violin model and how it stands amidst the other known models of protein allostery. Protein kinases have been described as the biomolecules of interest. © 2019 IUBMB Life, 71(6):685-696, 2019.

Engineering of lysine cyclodeaminase conformational dynamics for relieving substrate and product inhibitions in the biosynthesis of l-pipecolic acid

Efficient biocatalytic process construction by relieving substrate and product inhibitions via identification and engineering of enzyme conformational dynamics.

The crystal structure of a tetrahydrofolate-bound dihydrofolate reductase reveals the origin of slow product release

Dihydrofolate reductase (DHFR) catalyzes the stereospecific reduction of 7,8-dihydrofolate (FH2) to (6s)-5,6,7,8-tetrahydrofolate (FH4) via hydride transfer from NADPH. The consensus Escherichia coli DHFR mechanism involves conformational changes between closed and occluded states occurring during the rate-limiting product release step. Although the Protein Data Bank (PDB) contains over 250 DHFR structures, the FH4 complex structure responsible for rate-limiting product release is unknown. We report to our knowledge the first crystal structure of an E. coli. DHFR:FH4 complex at 1.03 Å resolution showing distinct stabilizing interactions absent in FH2 or related (6R)-5,10-dideaza-FH4 complexes. We discover the time course of decay of the co-purified endogenous FH4 during crystal growth, with conversion from FH4 to FH2 occurring in 2–3 days. We also determine another occluded complex structure of E. coli DHFR with a slow-onset nanomolar inhibitor that contrasts with the methotrexate complex, suggesting a plausible strategy for designing DHFR antibiotics by targeting FH4 product conformations.

Hongnan Cao et al. present the X-ray crystal structure of E. coli dihydrofolate reductase (DHFR) in complex with its reduced substrate, (6s)-5,6,7,8-tetrahydrofolate (FH4). This structure provides the first glimpse of the rate-limiting product release step of the DHFR mechanism and suggests a strategy for designing DHFR-targeting antibiotics.

High-pressure protein crystal structure analysis of Escherichia coli dihydrofolate reductase complexed with folate and NADP. Acta Crystallogr D Struct Biol 2018;74:895-905. [PMID: 30198899 PMCID: PMC6130465 DOI: 10.1107/s2059798318009397]

A high-pressure crystallographic study was conducted on Escherichia coli dihydrofolate reductase (ecDHFR) complexed with folate and NADP+ in crystal forms containing both the open and closed conformations of the M20 loop under high-pressure conditions of up to 800 MPa. At pressures between 270 and 500 MPa the crystal form containing the open conformation exhibited a phase transition from P21 to C2. Several structural changes in ecDHFR were observed at high pressure that were also accompanied by structural changes in the NADP+ cofactor and the hydration structure. In the crystal form with the closed conformation the M20 loop moved as the pressure changed, with accompanying conformational changes around the active site, including NADP+ and folate. These movements were consistent with the suggested hypothesis that movement of the M20 loop was necessary for ecDHFR to catalyze the reaction. In the crystal form with the open conformation the nicotinamide ring of the NADP+ cofactor undergoes a large flip as an intermediate step in the reaction, despite being in a crystalline state. Furthermore, observation of the water molecules between Arg57 and folate elucidated an early step in the substrate-binding pathway. These results demonstrate the possibility of using high-pressure protein crystallography as a method to capture high-energy substates or transient structures related to the protein reaction cycle.

Folate biosynthesis pathway: mechanisms and insights into drug design for infectious diseases

Folate pathway is a key target for the development of new drugs against infectious diseases since the discovery of sulfa drugs and trimethoprim. The knowledge about this pathway has increased in the last years and the catalytic mechanism and structures of all enzymes of the pathway are fairly understood. In addition, differences among enzymes from prokaryotes and eukaryotes could be used for the design of specific inhibitors. In this review, we show a panorama of progress that has been achieved within the folate pathway obtained in the last years. We explored the structure and mechanism of enzymes, several genetic features, strategies, and approaches used in the design of new inhibitors that have been used as targets in pathogen chemotherapy.

Expanding the Paradigm: Intrinsically Disordered Proteins and Allosteric Regulation

Allosteric regulatory processes are implicated at all levels of biological function. Recent advances in our understanding of the diverse and functionally significant class of intrinsically disordered proteins have identified a multitude of ways in which disordered proteins function within the confines of the allosteric paradigm. Allostery within or mediated by intrinsically disordered proteins ensures robust and efficient signal integration through mechanisms that would be extremely unfavorable or even impossible for globular protein interaction partners. Here, we highlight recent examples that indicate the breadth of biological outcomes that can be achieved through allosteric regulation by intrinsically disordered proteins. Ongoing and future work in this rapidly evolving area of research will expand our appreciation of the central role of intrinsically disordered proteins in ensuring the fidelity and efficiency of cellular regulation.

Measuring the signs of the methyl 1H chemical shift differences between major and 'invisible' minor protein conformational states using methyl 1H multi-quantum spectroscopy

Carr-Purcell-Meiboom-Gill (CPMG) type relaxation dispersion experiments are now routinely used to characterise protein conformational dynamics that occurs on the μs to millisecond (ms) timescale between a visible major state and 'invisible' minor states. The exchange rate(s) ([Formula: see text]), population(s) of the minor state(s) and the absolute value of the chemical shift difference [Formula: see text] (ppm) between different exchanging states can be extracted from the CPMG data. However the sign of [Formula: see text] that is required to reconstruct the spectrum of the 'invisible' minor state(s) cannot be obtained from CPMG data alone. Building upon the recently developed triple quantum (TQ) methyl [Formula: see text] CPMG experiment (Yuwen in Angew Chem 55:11490-11494, 2016) we have developed pulse sequences that use carbon detection to generate and evolve single quantum (SQ), double quantum (DQ) and TQ coherences from methyl protons in the indirect dimension to measure the chemical exchange-induced shifts of the SQ, DQ and TQ coherences from which the sign of [Formula: see text] is readily obtained for two state exchange. Further a combined analysis of the CPMG data and the difference in exchange induced shifts between the SQ and DQ resonances and between the SQ and TQ resonances improves the estimates of exchange parameters like the population of the minor state. We demonstrate the use of these experiments on two proteins undergoing exchange: (1) the ~ 18 kDa cavity mutant of T4 Lysozyme ([Formula: see text]) and (2) the [Formula: see text] kDa Peripheral Sub-unit Binding Domain (PSBD) from the acetyl transferase of Bacillus stearothermophilus ([Formula: see text]).

Slow Dynamics of Tryptophan-Water Networks in Proteins

Water has a profound effect on the dynamics of biomolecules and governs many biological processes, leading to the concept that function is slaved to solvent dynamics within and surrounding the biomolecule. Protein conformational changes on μs-ms time scales are frequently associated with protein function, but little is known about the behavior of protein-bound water on these time scales. Here we have used NMR relaxation dispersion measurements to probe the tryptophan indoles in the enzyme dihydrofolate reductase (DHFR). We find that during structural changes on the μs-ms time scale, large chemical shift changes are often observed for the NH proton on the indole ring, while relatively smaller chemical shift changes are observed for the ring nitrogen atom. Comparison with experimental chemical shifts and density functional theory-based chemical shift predictions show that during the structural change the tryptophan indole NHs remain bound to water, but the geometry of the protein-bound water networks changes. These results establish that relaxation dispersion measurements can indirectly probe water dynamics and indicate that water can influence, or be influenced by, protein conformational changes on the μs-ms time scale. Our data show that structurally conserved bound water molecules can play a critical role in transmitting information between functionally important regions of the protein and provide evidence that internal protein motions can be coupled through the mediation of hydrogen-bonded water bound in the protein structure.

Measurement of protein backbone 13CO and 15N relaxation dispersion at high resolution

Peak overlap in crowded regions of two-dimensional spectra prevents characterization of dynamics for many sites of interest in globular and intrinsically disordered proteins. We present new three-dimensional pulse sequences for measurement of Carr-Purcell-Meiboom-Gill relaxation dispersions at backbone nitrogen and carbonyl positions. To alleviate increase in the measurement time associated with the additional spectral dimension, we use non-uniform sampling in combination with two distinct methods of spectrum reconstruction: compressed sensing and co-processing with multi-dimensional decomposition. The new methodology was validated using disordered protein CD79A from B-cell receptor and an SH3 domain from Abp1p in exchange between its free form and bound to a peptide from the protein Ark1p. We show that, while providing much better resolution, the 3D NUS experiments give the similar accuracy and precision of the dynamic parameters to ones obtained using traditional 2D experiments. Furthermore, we show that jackknife resampling of the spectra yields robust estimates of peak intensities errors, eliminating the need for recording duplicate data points.

Defining the Structural Basis for Allosteric Product Release from E. coli Dihydrofolate Reductase Using NMR Relaxation Dispersion

The rate-determining step in the catalytic cycle of E. coli dihydrofolate reductase is tetrahydrofolate (THF) product release, which can occur via an allosteric or an intrinsic pathway. The allosteric pathway, which becomes accessible when the reduced cofactor NADPH is bound, involves transient sampling of a higher energy conformational state, greatly increasing the product dissociation rate as compared to the intrinsic pathway that obtains when NADPH is absent. Although the kinetics of this process are known, the enzyme structure and the THF product conformation in the transiently formed excited state remain elusive. Here, we use side-chain proton NMR relaxation dispersion measurements, X-ray crystallography, and structure-based chemical shift predictions to explore the structural basis of allosteric product release. In the excited state of the E:THF:NADPH product release complex, the reduced nicotinamide ring of the cofactor transiently enters the active site where it displaces the pterin ring of the THF product. The p-aminobenzoyl-l-glutamate tail of THF remains weakly bound in a widened binding cleft. Thus, through transient entry of the nicotinamide ring into the active site, the NADPH cofactor remodels the enzyme structure and the conformation of the THF to form a weakly populated excited state that is poised for rapid product release.

In Silico Studies of Small Molecule Interactions with Enzymes Reveal Aspects of Catalytic Function

Small molecules, such as solvent, substrate, and cofactor molecules, are key players in enzyme catalysis. Computational methods are powerful tools for exploring the dynamics and thermodynamics of these small molecules as they participate in or contribute to enzymatic processes. In-depth knowledge of how small molecule interactions and dynamics influence protein conformational dynamics and function is critical for progress in the field of enzyme catalysis. Although numerous computational studies have focused on enzyme-substrate complexes to gain insight into catalytic mechanisms, transition states and reaction rates, the dynamics of solvents, substrates, and cofactors are generally less well studied. Also, solvent dynamics within the biomolecular solvation layer play an important part in enzyme catalysis, but a full understanding of its role is hampered by its complexity. Moreover, passive substrate transport has been identified in certain enzymes, and the underlying principles of molecular recognition are an area of active investigation. Enzymes are highly dynamic entities that undergo different conformational changes, which range from side chain rearrangement of a residue to larger-scale conformational dynamics involving domains. These events may happen nearby or far away from the catalytic site, and may occur on different time scales, yet many are related to biological and catalytic function. Computational studies, primarily molecular dynamics (MD) simulations, provide atomistic-level insight and site-specific information on small molecule interactions, and their role in conformational pre-reorganization and dynamics in enzyme catalysis. The review is focused on MD simulation studies of small molecule interactions and dynamics to characterize and comprehend protein dynamics and function in catalyzed reactions. Experimental and theoretical methods available to complement and expand insight from MD simulations are discussed briefly.

Mispacking and the Fitness Landscape of the Green Fluorescent Protein Chromophore Milieu

The autocatalytic maturation of the chromophore in green fluorescent protein (GFP) was thought to require the precise positioning of the side chains surrounding it in the core of the protein, many of which are strongly conserved among homologous fluorescent proteins. In this study, we screened for green fluorescence in an exhaustive set of point mutations of seven residues that make up the chromophore microenvironment, excluding R96 and E222 because mutations at these positions have been previously characterized. Contrary to expectations, nearly all amino acids were tolerated at all seven positions. Only four point mutations knocked out fluorescence entirely. However, chromophore maturation was found to be slower and/or fluorescence reduced in several cases. Selected combinations of mutations showed nonadditive effects, including cooperativity and rescue. The results provide guidelines for the computational engineering of GFPs.

Multi-probe relaxation dispersion measurements increase sensitivity to protein dynamics

Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion measurements are a valuable tool for the characterization of structural transitions on the micro-millisecond timescale. While the measurement of (15)N relaxation dispersion is now routine, the measurements with alternative nuclei remain limited. Here we report (15)N as well as (1)H R2 relaxation dispersion measurements of the N23PP/S148A "dynamic knockout" mutant of dihydrofolate reductase. The (1)H dispersion measurements are complementary to (15)N data as many additional residues are observed to have dispersive behavior for the (1)H nucleus. Simultaneous fitting of the dispersion profiles for the two nuclei increases the accuracy of exchange parameters determined for individual residues and clustered groups of residues. The different sensitivity of the two nuclei to changes in backbone torsional angles, ring currents, and hydrogen bonding effects provides important insights into the nature of the structural changes that take place during the exchange process. We observe clear evidence of direct and indirect hydrogen bond effects for the (15)N and (1)H chemical shift changes in the active-site, modulation of ring current shielding in the CD-loop and backbone torsional changes in a cluster of residues associated with the C-terminus. This work demonstrates the power of combined (1)H and (15)N probes for the study of backbone dynamics on the micro-millisecond timescale though the analysis of chemical shift changes.

Probing conformational dynamics of an enzymatic active site by an in situ single fluorogenic probe under piconewton force manipulation

Unraveling the conformational details of an enzyme during the essential steps of a catalytic reaction (i.e., enzyme-substrate interaction, enzyme-substrate active complex formation, nascent product formation, and product release) is challenging due to the transient nature of intermediate conformational states, conformational fluctuations, and the associated complex dynamics. Here we report our study on the conformational dynamics of horseradish peroxidase using single-molecule multiparameter photon time-stamping spectroscopy with mechanical force manipulation, a newly developed single-molecule fluorescence imaging magnetic tweezers nanoscopic approach. A nascent-formed fluorogenic product molecule serves as a probe, perfectly fitting in the enzymatic reaction active site for probing the enzymatic conformational dynamics. Interestingly, the product releasing dynamics shows the complex conformational behavior with multiple product releasing pathways. However, under magnetic force manipulation, the complex nature of the multiple product releasing pathways disappears and more simplistic conformations of the active site are populated.

Role of Proteome Physical Chemistry in Cell Behavior

We review how major cell behaviors, such as bacterial growth laws, are derived from the physical chemistry of the cell’s proteins. On one hand, cell actions depend on the individual biological functionalities of their many genes and proteins. On the other hand, the common physics among proteins can be as important as the unique biology that distinguishes them. For example, bacterial growth rates depend strongly on temperature. This dependence can be explained by the folding stabilities across a cell’s proteome. Such modeling explains how thermophilic and mesophilic organisms differ, and how oxidative damage of highly charged proteins can lead to unfolding and aggregation in aging cells. Cells have characteristic time scales. For example, E. coli can duplicate as fast as 2–3 times per hour. These time scales can be explained by protein dynamics (the rates of synthesis and degradation, folding, and diffusional transport). It rationalizes how bacterial growth is slowed down by added salt. In the same way that the behaviors of inanimate materials can be expressed in terms of the statistical distributions of atoms and molecules, some cell behaviors can be expressed in terms of distributions of protein properties, giving insights into the microscopic basis of growth laws in simple cells.

Accurate determination of rates from non-uniformly sampled relaxation data

The application of non-uniform sampling (NUS) to relaxation experiments traditionally used to characterize the fast internal motion of proteins is quantitatively examined. Experimentally acquired Poisson-gap sampled data reconstructed with iterative soft thresholding are compared to regular sequentially sampled (RSS) data. Using ubiquitin as a model system, it is shown that 25 % sampling is sufficient for the determination of quantitatively accurate relaxation rates. When the sampling density is fixed at 25 %, the accuracy of rates is shown to increase sharply with the total number of sampled points until eventually converging near the inherent reproducibility of the experiment. Perhaps contrary to some expectations, it is found that accurate peak height reconstruction is not required for the determination of accurate rates. Instead, inaccuracies in rates arise from inconsistencies in reconstruction across the relaxation series that primarily manifest as a non-linearity in the recovered peak height. This indicates that the performance of an NUS relaxation experiment cannot be predicted from comparison of peak heights using a single RSS reference spectrum. The generality of these findings was assessed using three alternative reconstruction algorithms, eight different relaxation measurements, and three additional proteins that exhibit varying degrees of spectral complexity. From these data, it is revealed that non-linearity in peak height reconstruction across the relaxation series is strongly correlated with errors in NUS-derived relaxation rates. Importantly, it is shown that this correlation can be exploited to reliably predict the performance of an NUS-relaxation experiment by using three or more RSS reference planes from the relaxation series. The RSS reference time points can also serve to provide estimates of the uncertainty of the sampled intensity, which for a typical relaxation times series incurs no penalty in total acquisition time.

Role of Conformational Motions in Enzyme Function: Selected Methodologies and Case Studies

It is now common knowledge that enzymes are mobile entities relying on complex atomic-scale dynamics and coordinated conformational events for proper ligand recognition and catalysis. However, the exact role of protein dynamics in enzyme function remains either poorly understood or difficult to interpret. This mini-review intends to reconcile biophysical observations and biological significance by first describing a number of common experimental and computational methodologies employed to characterize atomic-scale residue motions on various timescales in enzymes, and second by illustrating how the knowledge of these motions can be used to describe the functional behavior of enzymes and even act upon it. Two biologically relevant examples will be highlighted, namely the HIV-1 protease and DNA polymerase β enzyme systems.

Advances in NMR Methods To Map Allosteric Sites: From Models to Translation

The last five years have witnessed major developments in the understanding of the allosteric phenomenon, broadly defined as coupling between remote molecular sites. Such advances have been driven not only by new theoretical models and pharmacological applications of allostery, but also by progress in the experimental approaches designed to map allosteric sites and transitions. Among these techniques, NMR spectroscopy has played a major role given its unique near-atomic resolution and sensitivity to the dynamics that underlie allosteric couplings. Here, we highlight recent progress in the NMR methods tailored to investigate allostery with the goal of offering an overview of which NMR approaches are best suited for which allosterically relevant questions. The picture of the allosteric "NMR toolbox" is provided starting from one of the simplest models of allostery (i.e., the four-state thermodynamic cycle) and continuing to more complex multistate mechanisms. We also review how such an "NMR toolbox" has assisted the elucidation of the allosteric molecular basis for disease-related mutations and the discovery of novel leads for allosteric drugs. From this overview, it is clear that NMR plays a central role not only in experimentally validating transformative theories of allostery, but also in tapping the full translational potential of allosteric systems.

NMR Methods to Study Dynamic Allostery

Nuclear magnetic resonance (NMR) spectroscopy provides a unique toolbox of experimental probes for studying dynamic processes on a wide range of timescales, ranging from picoseconds to milliseconds and beyond. Along with NMR hardware developments, recent methodological advancements have enabled the characterization of allosteric proteins at unprecedented detail, revealing intriguing aspects of allosteric mechanisms and increasing the proportion of the conformational ensemble that can be observed by experiment. Here, we present an overview of NMR spectroscopic methods for characterizing equilibrium fluctuations in free and bound states of allosteric proteins that have been most influential in the field. By combining NMR experimental approaches with molecular simulations, atomistic-level descriptions of the mechanisms by which allosteric phenomena take place are now within reach.

Ligand-Dependent Conformational Dynamics of Dihydrofolate Reductase

Enzymes are known to change among several conformational states during turnover. The role of such dynamic structural changes in catalysis is not fully understood. The influence of dynamics in catalysis can be inferred, but not proven, by comparison of equilibrium structures of protein variants and protein-ligand complexes. A more direct way to establish connections between protein dynamics and the catalytic cycle is to probe the kinetics of specific protein motions in comparison to progress along the reaction coordinate. We have examined the enzyme model system dihydrofolate reductase (DHFR) from Escherichia coli with tryptophan fluorescence-probed temperature-jump spectroscopy. We aimed to observe the kinetics of the ligand binding and ligand-induced conformational changes of three DHFR complexes to establish the relationship among these catalytic steps. Surprisingly, in all three complexes, the observed kinetics do not match a simple sequential two-step process. Through analysis of the relationship between ligand concentration and observed rate, we conclude that the observed kinetics correspond to the ligand binding step of the reaction and a noncoupled enzyme conformational change. The kinetics of the conformational change vary with the ligand's identity and presence but do not appear to be directly related to progress along the reaction coordinate. These results emphasize the need for kinetic studies of DHFR with highly specific spectroscopic probes to determine which dynamic events are coupled to the catalytic cycle and which are not.

Simple activation by acid of latent Ru-NHC-based metathesis initiators bearing 8-quinolinolate co-ligands

A straightforward synthesis utilizing the ring-opening metathesis polymerization (ROMP) reaction is described for acid-triggered N,O-chelating ruthenium-based pre-catalysts bearing one or two 8-quinolinolate ligands. The innovative pre-catalysts were tested regarding their behavior in ROMP and especially for their use in the synthesis of poly(dicyclopentadiene) (pDCPD). Bearing either the common phosphine leaving ligand in the first and second Grubbs olefin metathesis catalysts, or the Ru-O bond cleavage for the next Hoveyda-type catalysts, this work is a step forward towards the control of polymer functionalization and living or switchable polymerizations.

Protein Allostery and Conformational Dynamics

The functions of many proteins are regulated through allostery, whereby effector binding at a distal site changes the functional activity (e.g., substrate binding affinity or catalytic efficiency) at the active site. Most allosteric studies have focused on thermodynamic properties, in particular, substrate binding affinity. Changes in substrate binding affinity by allosteric effectors have generally been thought to be mediated by conformational transitions of the proteins or, alternatively, by changes in the broadness of the free energy basin of the protein conformational state without shifting the basin minimum position. When effector binding changes the free energy landscape of a protein in conformational space, the change affects not only thermodynamic properties but also dynamic properties, including the amplitudes of motions on different time scales and rates of conformational transitions. Here we assess the roles of conformational dynamics in allosteric regulation. Two cases are highlighted where NMR spectroscopy and molecular dynamics simulation have been used as complementary approaches to identify residues possibly involved in allosteric communication. Perspectives on contentious issues, for example, the relationship between picosecond-nanosecond local and microsecond-millisecond conformational exchange dynamics, are presented.