Escherichia coli dihydrofolate reductase catalyzed proton and hydride transfers: temporal order and the roles of Asp27 and Tyr100

Statistical Coupling Analysis Predicts Correlated Motions in Dihydrofolate Reductase

Dihydrofolate reductase (DHFR), due to its universality and the depth with which it has been studied, is a model system in the study of protein dynamics. Myriad previous works have identified networks of residues in positions near to and remote from the active site that are involved in the dynamics. For example, specific mutations on the Met20 loop in Escherichia coli DHFR (N23PP/S148A) are known to disrupt millisecond-time scale motions as well as reduce catalytic activity. However, how and if networks of dynamically coupled residues influence the evolution of DHFR is still an unanswered question. In this study, we first identify, by statistical coupling analysis and molecular dynamic simulations, a network of coevolving residues that possesses increased correlated motions. We then go on to show that allosteric communication in this network is knocked down in N23PP/S148A mutant E. coli DHFR. We also identify two sites in the human DHFR sector which may accommodate the Met20 loop double proline motif. Finally, we demonstrate a concerted evolutionary change in the human DHFR allosteric networks, which maintains dynamic communication. These findings strongly implicate protein dynamics as a driving force for evolution.

Statistical Coupling Analysis Predicts Correlated Motions in Dihydrofolate Reductase

The role of dynamics in enzymatic function is a highly debated topic. Dihydrofolate reductase (DHFR), due to its universality and the depth with which it has been studied, is a model system in this debate. Myriad previous works have identified networks of residues in positions near to and remote from the active site that are involved in dynamics and others that are important for catalysis. For example, specific mutations on the Met20 loop in E. coli DHFR (N23PP/S148A) are known to disrupt millisecond-timescale motions and reduce catalytic activity. However, how and if networks of dynamically coupled residues influence the evolution of DHFR is still an unanswered question. In this study, we first identify, by statistical coupling analysis and molecular dynamic simulations, a network of coevolving residues, which possess increased correlated motions. We then go on to show that allosteric communication in this network is selectively knocked down in N23PP/S148A mutant E. coli DHFR. Finally, we identify two sites in the human DHFR sector which may accommodate the Met20 loop double proline mutation while preserving dynamics. These findings strongly implicate protein dynamics as a driving force for evolution.

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.

Mutexa: A Computational Ecosystem for Intelligent Protein Engineering

Protein engineering holds immense promise in shaping the future of biomedicine and biotechnology. This Review focuses on our ongoing development of Mutexa, a computational ecosystem designed to enable "intelligent protein engineering". In this vision, researchers will seamlessly acquire sequences of protein variants with desired functions as biocatalysts, therapeutic peptides, and diagnostic proteins through a finely-tuned computational machine, akin to Amazon Alexa's role as a versatile virtual assistant. The technical foundation of Mutexa has been established through the development of a database that combines and relates enzyme structures and their respective functions (e.g., IntEnzyDB), workflow software packages that enable high-throughput protein modeling (e.g., EnzyHTP and LassoHTP), and scoring functions that map the sequence-structure-function relationship of proteins (e.g., EnzyKR and DeepLasso). We will showcase the applications of these tools in benchmarking the convergence conditions of enzyme functional descriptors across mutants, investigating protein electrostatics and cavity distributions in SAM-dependent methyltransferases, and understanding the role of nonelectrostatic dynamic effects in enzyme catalysis. Finally, we will conclude by addressing the future steps and fundamental challenges in our endeavor to develop new Mutexa applications that assist the identification of beneficial mutants in protein engineering.

X-ray-driven chemistry and conformational heterogeneity in atomic resolution crystal structures of bacterial dihydrofolate reductases

Dihydrofolate reductase (DHFR) catalyzes the NADPH-dependent reduction of dihydrofolate to tetrahydrofolate. Bacterial DHFRs are targets of several important antibiotics as well as model enzymes for the role of protein conformational dynamics in enzyme catalysis. We collected 0.93 Å resolution X-ray diffraction data from both Bacillus subtilis (Bs) and E. coli (Ec) DHFRs bound to folate and NADP+. These oxidized ternary complexes should not be able to perform chemistry, however electron density maps suggest hydride transfer is occurring in both enzymes. Comparison of low- and high-dose EcDHFR datasets show that X-rays drive partial production of tetrahydrofolate. Hydride transfer causes the nicotinamide moiety of NADP+ to move towards the folate as well as correlated shifts in nearby residues. Higher radiation dose also changes the conformational heterogeneity of Met20 in EcDHFR, supporting a solvent gating role during catalysis. BsDHFR has a different pattern of conformational heterogeneity and an unexpected disulfide bond, illustrating important differences between bacterial DHFRs. This work demonstrates that X-rays can drive hydride transfer similar to the native DHFR reaction and that X-ray photoreduction can be used to interrogate catalytically relevant enzyme dynamics in favorable cases.

From a binding module to essential catalytic activity: how nature stumbled on a good thing

Enzymes are complex macromolecules capable of catalyzing a wide variety of chemical reactions with high efficiency. Nonetheless, biological catalysis can be rudimentary. Here, we describe an enzyme that is built from a simple protein fold. This short protein sequence - almost a peptide - belongs to the ancient SH3 family of binding modules. Surprisingly, this binding module catalyzes the specific reduction of dihydrofolate using NADPH as a reducing cofactor, making this a dihydrofolate reductase. Too small to provide all the required binding and catalytic machinery on its own, it homotetramerizes, thus creating a large, central active site environment. Remarkably, none of the active site residues is essential to the catalytic function. Instead, backbone interactions juxtapose the reducing cofactor proximal to the target imine of the folate substrate, and a specific motion of the substrate promotes formation of the transition state. In this feature article, we describe the features that make this small protein a functional enzyme capable of catalyzing a metabolically essential reaction, highlighting the characteristics that make it a model for the evolution of primitive enzymes from binding modules.

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.

Is Posttranslational Folding More Efficient Than Refolding from a Denatured State: A Computational Study

The folding of proteins into their native conformation is a complex process that has been extensively studied over the past half-century. The ribosome, the molecular machine responsible for protein synthesis, is known to interact with nascent proteins, adding further complexity to the protein folding landscape. Consequently, it is unclear whether the folding pathways of proteins are conserved on and off the ribosome. The main question remains: to what extent does the ribosome help proteins fold? To address this question, we used coarse-grained molecular dynamics simulations to compare the mechanisms by which the proteins dihydrofolate reductase, type III chloramphenicol acetyltransferase, and d-alanine-d-alanine ligase B fold during and after vectorial synthesis on the ribosome to folding from the full-length unfolded state in bulk solution. Our results reveal that the influence of the ribosome on protein folding mechanisms varies depending on the size and complexity of the protein. Specifically, for a small protein with a simple fold, the ribosome facilitates efficient folding by helping the nascent protein avoid misfolded conformations. However, for larger and more complex proteins, the ribosome does not promote folding and may contribute to the formation of intermediate misfolded states cotranslationally. These misfolded states persist posttranslationally and do not convert to the native state during the 6 μs runtime of our coarse-grain simulations. Overall, our study highlights the complex interplay between the ribosome and protein folding and provides insight into the mechanisms of protein folding on and off the ribosome.

Modeling Ensembles of Enzyme Reaction Pathways with Hi-MSM Reveals the Importance of Accounting for Pathway Diversity

Conventional quantum mechanical-molecular mechanics (QM/MM) simulation approaches for modeling enzyme reactions often assume that there is one dominant reaction pathway and that this pathway can be sampled starting from an X-ray structure of the enzyme. These assumptions reduce computational cost; however, their validity has not been extensively tested. This is due in part to the lack of a rigorous formalism for integrating disparate pathway information from dynamical QM/MM calculations. Here, we present a way to model ensembles of reaction pathways efficiently using a divide-and-conquer strategy through Hierarchical Markov State Modeling (Hi-MSM). This approach allows information on multiple, distinct pathways to be incorporated into a chemical kinetic model, and it allows us to test these two assumptions. Applying Hi-MSM to the reaction carried out by dihydrofolate reductase (DHFR) we find (i) there are multiple, distinct pathways significantly contributing to the overall flux of the reaction that the conventional approach does not identify and (ii) that the conventional approach does not identify the dominant reaction pathway. Thus, both assumptions underpinning the conventional approach are violated. Since DHFR is a relatively small enzyme, and configuration space scales exponentially with protein size, accounting for multiple reaction pathways is likely to be necessary for most enzymes.

Convergence in determining enzyme functional descriptors across Kemp eliminase variants

Molecular simulations have been extensively employed to accelerate biocatalytic discoveries. Enzyme functional descriptors derived from molecular simulations have been leveraged to guide the search for beneficial enzyme mutants. However, the ideal active-site region size for computing the descriptors over multiple enzyme variants remains untested. Here, we conducted convergence tests for dynamics-derived and electrostatic descriptors on 18 Kemp eliminase variants across six active-site regions with various boundary distances to the substrate. The tested descriptors include the root-mean-square deviation of the active-site region, the solvent accessible surface area ratio between the substrate and active site, and the projection of the electric field (EF) on the breaking C-H bond. All descriptors were evaluated using molecular mechanics methods. To understand the effects of electronic structure, the EF was also evaluated using quantum mechanics/molecular mechanics methods. The descriptor values were computed for 18 Kemp eliminase variants. Spearman correlation matrices were used to determine the region size condition under which further expansion of the region boundary does not substantially change the ranking of descriptor values. We observed that protein dynamics-derived descriptors, including RMSD_{active_site} and SASA_ratio, converge at a distance cutoff of 5 Å from the substrate. The electrostatic descriptor, EF_C-H, converges at 6 Å using molecular mechanics methods with truncated enzyme models and 4 Å using quantum mechanics/molecular mechanics methods with whole enzyme model. This study serves as a future reference to determine descriptors for predictive modeling of enzyme engineering.

Structure-guided functional studies of plasmid-encoded dihydrofolate reductases reveal a common mechanism of trimethoprim resistance in Gram-negative pathogens

Two plasmid-encoded dihydrofolate reductase (DHFR) isoforms, DfrA1 and DfrA5, that give rise to high levels of resistance in Gram-negative bacteria were structurally and biochemically characterized to reveal the mechanism of TMP resistance and to support phylogenic groupings for drug development against antibiotic resistant pathogens. Preliminary screening of novel antifolates revealed related chemotypes that showed high levels of inhibitory potency against Escherichia coli chromosomal DHFR (EcDHFR), DfrA1, and DfrA5. Kinetics and biophysical analysis, coupled with crystal structures of trimethoprim bound to EcDHFR, DfrA1 and DfrA5, and two propargyl-linked antifolates (PLA) complexed with EcDHFR, DfrA1 and DfrA5, were determined to define structural features of the substrate binding pocket and guide synthesis of pan-DHFR inhibitors.

Critical residue variations in two of the most clinically prevalent DHFR isoforms are identified as a common structural element in trimethoprim-resistant DHFR which impose changes on enzyme catalysis and ligand-cofactor cooperativity.

Molecular Dynamics-Derived Descriptor Informs the Impact of Mutation on the Catalytic Turnover Number in Lactonase Across Substrates

Molecular dynamics simulations have been extensively employed to reveal the roles of protein dynamics in mediating enzyme catalysis. However, simulation-derived predictive descriptors that inform the impacts of mutations on catalytic turnover numbers remain largely unexplored. In this work, we report the identification of molecular modeling-derived descriptors to predict mutation effect on the turnover number of lactonase SsoPox with both native and non-native substrates. The study consists of 10 enzyme-substrate complexes resulting from a combination of five enzyme variants with two substrates. For each complex, we derived 15 descriptors from molecular dynamics simulations and applied principal component analysis to rank the predictive capability of the descriptors. A top-ranked descriptor was identified, which is the solvent-accessible surface area (SASA) ratio of the substrate to the active site pocket. A uniform volcano-shaped plot was observed in the distribution of experimental activation free energy against the SASA ratio. To achieve efficient lactonase hydrolysis, a non-native substrate-bound enzyme variant needs to involve a similar range of the SASA ratio to the native substrate-bound wild-type enzyme. The descriptor reflects how well the enzyme active site pocket accommodates a substrate for reaction, which has the potential of guiding optimization of enzyme reaction turnover for non-native chemical transformations.

Evolution of Optimized Hydride Transfer Reaction and Overall Enzyme Turnover in Human Dihydrofolate Reductase

Evolution of dihydrofolate reductase (DHFR) has been studied using the enzyme from Escherichia coli DHFR (ecDHFR) as a model, but less studies have used the enzyme from Homo sapiens DHFR (hsDHFR). Each enzyme maintains a short and narrow distribution of hydride donor-acceptor distances (DAD) at the tunneling ready state (TRS). Evolution of the enzyme was previously studied in ecDHFR where three key sites were identified as important to the catalyzed reaction. The corresponding sites in hsDHFR are F28, 62-PEKN, and 26-PPLR. Each of these sites was studied here through the creation of mutant variants of the enzyme and measurements of the temperature dependence of the intrinsic kinetic isotope effects (KIEs) on the reaction. F28 is mutated first to M (F28M) and then to the L of the bacterial enzyme (F28L). The KIEs of the F28M variant are larger and more temperature-dependent than wild-type (WT), suggesting a broader and longer average DAD at the TRS. To more fully mimic ecDHFR, we also study a triple mutant of the human enzyme (F32L-PP26N-PEKN62G). Remarkably, the intrinsic KIEs, while larger in magnitude, are temperature-independent like the WT enzymes. We also construct deletion mutations of hsDHFR removing both the 62-PEKN and 26-PPLR sequences. The results mirror those described previously for insertion mutants of ecDHFR. Taken together, these results suggest a balancing act during DHFR evolution between achieving an optimal TRS for hydride transfer and preventing product inhibition arising from the different intercellular pools of NADPH and NADP⁺ in prokaryotic and eukaryotic cells.

Pressure Adaptations in Deep-Sea Moritella Dihydrofolate Reductases: Compressibility versus Stability

Simple Summary

Deep-sea organisms must have proteins that function under high hydrostatic pressure to survive. Adaptations used in proteins from “pressure-loving” piezophiles may include greater compressibility or greater stability against pressure-induced destabilization. However, while greater compressibility can be accomplished by greater void volume, larger cavities in a protein have been associated with greater destabilization and even unfolding as pressure is increased. Here, computer simulations of dihydrofolate reductase from a moderate piezophile and a hyperpiezophile were performed to understand the balance between adaptations for greater compressibility and those against pressure destabilization and unfolding. The results indicate that while compressibility appears to be important for deep-sea microbes, adaptation for the greatest depths may be to prevent water penetration into the interior.

Abstract

Proteins from “pressure-loving” piezophiles appear to adapt by greater compressibility via larger total cavity volume. However, larger cavities in proteins have been associated with lower unfolding pressures. Here, dihydrofolate reductase (DHFR) from a moderate piezophile Moritella profunda (Mp) isolated at ~2.9 km in depth and from a hyperpiezophile Moritella yayanosii (My) isolated at ~11 km in depth were compared using molecular dynamics simulations. Although previous simulations indicate that MpDHFR is more compressible than a mesophile DHFR, here the average properties and a quasiharmonic analysis indicate that MpDHFR and MyDHFR have similar compressibilities. A cavity analysis also indicates that the three unique mutations in MyDHFR are near cavities, although the cavities are generally similar in size in both. However, while a cleft overlaps an internal cavity, thus forming a pathway from the surface to the interior in MpDHFR, the unique residue Tyr103 found in MyDHFR forms a hydrogen bond with Leu78, and the sidechain separates the cleft from the cavity. Thus, while Moritella DHFR may generally be well suited to high-pressure environments because of their greater compressibility, adaptation for greater depths may be to prevent water entry into the interior cavities.

Effects of Distal Mutations on Ligand-Binding Affinity in E. coli Dihydrofolate Reductase

Mutations far from the center of chemical activity in dihydrofolate reductase (DHFR) can affect several steps in the catalytic cycle. Mutations at highly conserved positions and the distal distance of the catalytic center (Met-42, Thr-113, and Gly-121) were designed, including single-point and double-point mutations. Upon ligand binding, the fluorescence of the intrinsic optical probe, tryptophan, decreases due to either fluorescence quenching or energy transfer. We demonstrated an optical approach in measuring the equilibrium dissociation constant for enzyme-cofactor, enzyme-substrate, and enzyme-product complexes in wildtype ecDHFR and each mutant. We propose that the effects of these distal mutations on ligand-binding affinity stem from the spatial steric hindrance, the disturbance on the hydrogen network, or the modification of the protein flexibility. The modified N-terminus tag in DHFR acts as a cap on the entrance of the substrate-binding cavity, squeezes the adenosine binding subdomain, and influences the binding of NADPH in some mutants. If the mutation positions are away from the N-terminus tag and the adenosine binding subdomain, the additive effects due to the N-terminus tag were not observed. In the double-mutant-cycle analysis, double mutations show nonadditive properties upon either cofactor or substrate binding. Also, in general, the first point mutation strongly affects the ligand binding compared to the second one.

Capturing the Catalytic Proton of Dihydrofolate Reductase: Implications for General Acid-Base Catalysis

Acid–base catalysis, which involves one or more proton transfer reactions, is a chemical mechanism commonly employed by many enzymes. The molecular basis for catalysis is often derived from structures determined at the optimal pH for enzyme activity. However, direct observation of protons from experimental structures is quite difficult; thus, a complete mechanistic description for most enzymes remains lacking. Dihydrofolate reductase (DHFR) exemplifies general acid–base catalysis, requiring hydride transfer and protonation of its substrate, DHF, to form the product, tetrahydrofolate (THF). Previous X-ray and neutron crystal structures coupled with theoretical calculations have proposed that solvent mediates the protonation step. However, visualization of a proton transfer has been elusive. Based on a 2.1 Å resolution neutron structure of a pseudo-Michaelis complex of E. coli DHFR determined at acidic pH, we report the direct observation of the catalytic proton and its parent solvent molecule. Comparison of X-ray and neutron structures elucidated at acidic and neutral pH reveals dampened dynamics at acidic pH, even for the regulatory Met20 loop. Guided by the structures and calculations, we propose a mechanism where dynamics are crucial for solvent entry and protonation of substrate. This mechanism invokes the release of a sole proton from a hydronium (H₃O⁺) ion, its pathway through a narrow channel that sterically hinders the passage of water, and the ultimate protonation of DHF at the N5 atom.

Temperature-Dependent Kinetic Isotope Effects in R67 Dihydrofolate Reductase from Path-Integral Simulations

Calculation of temperature-dependent kinetic isotope effects (KIE) in enzymes presents a significant theoretical challenge. Additionally, it is not trivial to identify enzymes with available experimental accurate intrinsic KIEs in a range of temperatures. In the current work, we present a theoretical study of KIEs in the primitive R67 dihydrofolate reductase (DHFR) enzyme and compare with experimental work. The advantage of R67 DHFR is its significantly lower kinetic complexity compared to more evolved DHFR isoforms. We employ mass-perturbation-based path-integral simulations in conjunction with umbrella sampling and a hybrid quantum mechanics-molecular mechanics Hamiltonian. We obtain temperature-dependent KIEs in good agreement with experiments and ascribe the temperature-dependent KIEs primarily to zero-point energy effects. The active site in the primitive enzyme is found to be poorly preorganized, which allows excessive water access to the active site and results in loosely bound reacting ligands.

Modulating Enzyme Function via Dynamic Allostery within Biliverdin Reductase B

The biliverdin reductase B (BLVRB) class of enzymes catalyze the NADPH-dependent reduction of multiple flavin substrates and are emerging as critical players in cellular redox regulation. However, the role of dynamics and allostery have not been addressed, prompting studies here that have revealed a position 15 Å away from the active site within human BLVRB (T164) that is inherently dynamic and can be mutated to control global micro-millisecond motions and function. By comparing the inherent dynamics through nuclear magnetic resonance (NMR) relaxation approaches of evolutionarily distinct BLVRB homologues and by applying our previously developed Relaxation And Single Site Multiple Mutations (RASSMM) approach that monitors both the functional and dynamic effects of multiple mutations to the single T164 site, we have discovered that the most dramatic mutagenic effects coincide with evolutionary changes and these modulate coenzyme binding. Thus, evolutionarily changing sites distal to the active site serve as dynamic "dials" to globally modulate motions and function. Despite the distal dynamic and functional coupling modulated by this site, micro-millisecond motions span an order of magnitude in their apparent kinetic rates of motions. Thus, global dynamics within BLVRB are a collection of partially coupled motions tied to catalytic function.

Structure-based mechanisms: On the way to apply alcohol dehydrogenases/reductases to organic-aqueous systems

Alcohol dehydrogenases/reductases catalyze enantioselective syntheses of versatile chiral compounds relying on direct hydride transfer from cofactor to substrates, or to an intermediate and then to substrates. Since most of the substrates catalyzed by alcohol dehydrogenases/reductases are insoluble in aqueous solutions, increasing interest has been turning to organic-aqueous systems. However, alcohol dehydrogenases/reductases are normally instable in organic solvents, leading to the unsatisfied enantioselective synthesis efficiency. The behaviors of these enzymes in organic solvents at an atomic level are unclear, thus it is of great importance to understand its structure-based mechanisms in organic-aqueous systems to improve their relative stability. Here, we summarized the accessible structures of alcohol dehydrogenases/reductases in Protein Data Bank crystallized in organic-aqueous systems, and compared the structures of alcohol dehydrogenases/reductases which have different tolerance towards organic solvents. By understanding the catalytic behaviors and mechanisms of these enzymes in organic-aqueous systems, the efficient enantioselective syntheses mediated by alcohol dehydrogenases/reductases and further challenges are also discussed through solvent engineering and enzyme-immobilization in the last decade.

Ligand Conjugation Directs the Formation of a 1,3-Dihydropyridinate Regioisomer

The selective formation of the 1,4-dihydropyridine isomer of NAD(P)H is mirrored by the selective formation of 1,4-dihydropyridinate ligand-metal complexes in synthetic systems. Here we demonstrate that ligand conjugation can be used to promote selective 1,3-dihydropyridinate formation. This represents an advance toward controlling and tuning the selectivity in dihydropyridinate formation chemistry. The reaction of (I₂P^2-)Al(THF)Cl [1; I₂P = bis(imino)pyridine; THF = tetrahydrofuran] with the one-electron oxidant (2,2,6,6-tetramethylpiperidin-1-yl)oxyl (TEMPO) afforded (I₂P^-)Al(TEMPO)Cl (2), which can be reduced with sodium to the twice-reduced ligand complex (I₂P^2-)Al(TEMPO) (3). Compounds 2 and 3 serve as precursors for high-yielding and selective routes to an aluminum-supported 1,3-dihydropyridinate complex via the reaction of 2 with 3 equiv of potassium metal or the reaction of 3 with KH.

Mutational analysis confirms the presence of distal inhibitor-selectivity determining residues in B. stearothermophilus dihydrofolate reductase

Many antibacterial and antiparasitic drugs work by competitively inhibiting dihydrofolate reductase (DHFR), a vital enzyme in folate metabolism. The interactions between inhibitors and DHFR active site residues are known in many homologs but the contributions from distal residues are less understood. Identifying distal residues that aid in inhibitor binding can improve targeted drug development programs by accounting for distant influences that may be less conserved and subject to frequent resistance causing mutations. Previously, a novel, homology-based, computational approach that mines ligand inhibition data was used to predict residues involved in inhibitor selectivity in the DHFR family. Expectedly, some inhibitor selectivity determining residue positions were predicted to lie in the active site and coincide with experimentally known inhibitor selectivity determining positions. However, other residues that group spatially in clusters distal to the active site have not been previously investigated. In this study, the effect of introducing amino acid substitutions at one of these predicted clusters (His38-Ala39-Ile40) on the inhibitor selectivity profile in Bacillus stearothermophilus dihydrofolate reductase (Bs DHFR) was investigated. Mutations were introduced into these cluster positions to change sidechain chemistry and size. We determined kcat and KM values and measured KD values at equilibrium for two competitive DHFR inhibitors, trimethoprim (TMP) and pyrimethamine (PYR). Mutations in the His38-Ala39-Ile40 cluster significantly impacted inhibitor binding and TMP/PYR selectivity - seven out of nine mutations resulted in tighter binding to PYR when compared to TMP. These data suggest that the His38-Ala39-Ile40 cluster is a distal inhibitor selectivity determining region that favors PYR binding in Bs DHFR and, possibly, throughout the DHFR family.

Enzymatic and Mutational Analysis of the PruA Pteridine Reductase Required for Pterin-Dependent Control of Biofilm Formation in Agrobacterium tumefaciens

Pterins are ubiquitous biomolecules with diverse functions including roles as cofactors, pigments, and redox mediators. Recently, a novel pterin-dependent signaling pathway that controls biofilm formation was identified in the plant pathogen, Agrobacterium tumefaciens A key player in this pathway is a pteridine reductase termed PruA, where its enzymatic activity has been shown to control surface attachment and limit biofilm formation. Here, we biochemically characterize PruA to investigate the catalytic properties and substrate specificity of this pteridine reductase. PruA demonstrates maximal catalytic efficiency with dihydrobiopterin and comparable activities with the stereoisomers dihydromonapterin and dihydroneopterin. Since A. tumefaciens does not synthesize or utilize biopterins, the likely physiological substrate is dihydromonapterin or dihydroneopterin, or both. Notably, PruA does not exhibit pteridine reductase activity with dihydrofolate or fully oxidized pterins. Site-directed mutagenesis studies of a conserved tyrosine residue, the key component of a putative catalytic triad, indicate that this tyrosine is not directly involved in PruA catalysis but may be important for substrate or cofactor binding. Additionally, mutagenesis of the arginine residue in the N-terminal TGX3RXG motif significantly reduces the catalytic efficiency of PruA, supporting its proposed role in pterin binding and catalysis. Finally, we report the enzymatic characterization of PruA homologs from Pseudomonas aeruginosa and Brucella abortus, thus expanding the roles and potential significance of pteridine reductases in diverse bacteria.Importance Biofilms are complex multicellular communities that are formed by diverse bacteria. In the plant pathogen, Agrobacterium tumefaciens, the transition from a free-living motile state to a non-motile biofilm state is governed by a novel signaling pathway involving small molecules called pterins. The involvement of pterins in biofilm formation is unexpected and prompts many questions about the molecular details of this pathway. This work biochemically characterizes the PruA pteridine reductase involved in the signaling pathway to reveal its enzymatic properties and substrate preference, thus providing important insight into pterin biosynthesis and its role in A. tumefaciens biofilm control. Additionally, the enzymatic characteristics of related pteridine reductases from mammalian pathogens are examined to uncover potential roles of these enzymes in other bacteria.

What are the signatures of tunnelling in enzyme-catalysed reactions?

Computed tunnelling contributions and correlations between apparent activation enthalpy and entropy are explored for the interpretation of enzyme-catalysed H-transfer reactions.

Dual and selective inhibitors of pteridine reductase 1 (PTR1) and dihydrofolate reductase-thymidylate synthase (DHFR-TS) from Leishmania chagasi

Leishmaniasis is a tropical disease found in more than 90 countries. The drugs available to treat this disease have nonspecific action and high toxicity. In order to develop novel therapeutic alternatives to fight this ailment, pteridine reductase 1 (PTR1) and dihydrofolate reductase-thymidylate synthase (DHF-TS) have been targeted, once Leishmania is auxotrophic for folates. Although PTR1 and DHFR-TS from other protozoan parasites have been studied, their homologs in Leishmania chagasi have been poorly characterized. Hence, this work describes the optimal conditions to express the recombinant LcPTR1 and LcDHFR-TS enzymes, as well as balanced assay conditions for screening. Last but not the least, we show that 2,4 diaminopyrimidine derivatives are low-micromolar competitive inhibitors of both enzymes (LcPTR1 Ki = 1.50-2.30 µM and LcDHFR Ki = 0.28-3.00 µM) with poor selectivity index. On the other hand, compound 5 (2,4-diaminoquinazoline derivative) is a selective LcPTR1 inhibitor (Ki = 0.47 µM, selectivity index = 20).

Loss of Hyperconjugative Effects Drives Hydride Transfer during Dihydrofolate Reductase Catalysis

Hydride transfer is widespread in nature and has an essential role in applied research. However, the mechanisms of how this transformation occurs in living organisms remain a matter of vigorous debate. Here, we examined dihydrofolate reductase (DHFR), an enzyme that catalyzes hydride from C4′ of NADPH to C6 of 7,8-dihydrofolate (H₂F). Despite many investigations of the mechanism of this reaction, the contribution of polarization of the π-bond of H₂F in driving hydride transfer remains unclear. H₂F was stereospecifically labeled with deuterium β to the reacting center, and β-deuterium kinetic isotope effects were measured. Our experimental results combined with analysis derived from QM/MM simulations reveal that hydride transfer is triggered by polarization at the C6 of H₂F. The σ C_β–H bonds contribute to the buildup of the cationic character during the chemical transformation, and hyperconjugation influences the formation of the transition state. Our findings provide key insights into the hydride transfer mechanism of the DHFR-catalyzed reaction, which is a target for antiproliferative drugs and a paradigmatic model in mechanistic enzymology.

Exploring PfDHFR reaction surface: A combined molecular dynamics and QM/MM analysis

The substrate to the enzyme PfDHFR (Plasmodium falciparum Dihydrofolate Reductase) is a small molecule dihydrofolate (DHF), it gets converted to tetrahydrofolate (THF) in the active site of the enzyme. The PfDHFR reaction surface involves the protonation of DHF to DHFP as an initial step before the catalytic conversion. The binding affinities of all these species (DHF, DHFP and THF) contribute to the mechanism of DHFR catalytic action. Molecular dynamics (MD) simulations and Quantum Mechanics/Molecular Mechanics (QM/MM) analysis were performed to evaluate the binding affinity and molecular recognition interactions of the substrate DHF/DHFP and the product THF, in the active site of wild-type PfDHFR (wtPfDHFR). The binding affinities of the cofactor NADPH/NADP⁺ were also estimated in all the three complexes. The molecular dynamics (MD) simulations of the substrate, product and cofactor in the cavities of wtPfDHFR revealed the variation of the atomic level interactions during the course of the catalytic conversion. It was found that the DHFP binds very strongly to the PfDHFR active site and pulls the cofactor NADPH closer to itself. The QM/MM analysis revealed that the binding energy of DHFP (-59.82 kcal/mol) and NADPH (-100.24 kcal/mol) in DHFP-wtPfDHFR complex, is higher in comparison to the binding energy of DHF (-38.67 kcal/mol) and NADPH (-77.53 kcal/mol) in DHF-wtPfDHFR complex and the binding energy of THF (-30.72 kcal/mol) and NADP⁺ (-73.72 kcal/mol) in THF-wtPfDHFR complex. The hydride ion donor-acceptor distance (DAD) analysis was also carried out. This combined MD and QM/MM analysis revealed that the protonation of DHF increases the proximity between the substrate and the cofactor, thus facilitates the reaction profile of PfDHFR.

Effect of Asp122 Mutation on the Hydride Transfer in E. coli DHFR Demonstrates the Goldilocks of Enzyme Flexibility

Dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate (THF) in the presence of NADPH. The key hydride transfer step in the reaction is facilitated by a combination of enzyme active site preorganization and correlated protein motions in the Michaelis-Menten (E:NADPH:DHF) complex. The present theoretical study employs mutagenesis to examine the relation between structural and functional properties of the enzyme. We mutate Asp122 in Escherichia coli DHFR, which is a conserved amino acid in the DHFR family. The consequent effect of the mutation on enzyme catalysis is examined from an energetic, structural and short-time dynamic perspective. Our investigations suggest that the structural and short-time dynamic perturbations caused by Asp122X mutations (X = Asn, Ser, Ala) are along the reaction coordinate and lower the rate of hydride transfer. Importantly, analysis of the correlated and principle component motions in the enzyme suggest that the mutation alters the coupled motions that are present in the wild-type enzyme. In the case of D122N and D122S, the mutations inhibit coupled motion, whereas in the case of D122A, the mutation enhances coupled motion, although all mutations result in similar rate reduction. These results emphasize a Goldilocks principle of enzyme flexibility, that is, enzymes should neither be too rigid nor too flexible.

Evolutionarily Related Dihydrofolate Reductases Perform Coequal Functions Yet Show Divergence in Their Trajectories

The enzyme dihydrofolate reductase (DHFR) catalyzes NADPH dependent reduction of dihydrofolate to tetrahydrofolate. It plays a crucial role in the DNA synthesis. The investigation of evolution of DHFR generates immense curiosity. It aids in predicting how the enzyme has adapted to the surroundings of various cell types. In spite of great similarity in the structure of E. coli DHFR and human DHFR, their primary sequences are divergent to a great extent, which is evident in variations in the kinetics mechanism of their catalysis. In presence of physiological levels of ligands, they possess distinct kinetics and different rate limiting steps. We have reviewed the process of their unfolding and refolding, their behaviour in denaturing conditions and in presence of various chaperones. Although there is structural similarity between these two homologous enzymes yet they have established distinct mechanisms to accomplish the coequal functions.

Modulating Enzyme Activity by Altering Protein Dynamics with Solvent

Optimal enzyme activity depends on a number of factors, including structure and dynamics. The role of enzyme structure is well recognized; however, the linkage between protein dynamics and enzyme activity has given rise to a contentious debate. We have developed an approach that uses an aqueous mixture of organic solvent to control the functionally relevant enzyme dynamics (without changing the structure), which in turn modulates the enzyme activity. Using this approach, we predicted that the hydride transfer reaction catalyzed by the enzyme dihydrofolate reductase (DHFR) from Escherichia coli in aqueous mixtures of isopropanol (IPA) with water will decrease by ∼3 fold at 20% (v/v) IPA concentration. Stopped-flow kinetic measurements find that the pH-independent k_hydride rate decreases by 2.2 fold. X-ray crystallographic enzyme structures show no noticeable differences, while computational studies indicate that the transition state and electrostatic effects were identical for water and mixed solvent conditions; quasi-elastic neutron scattering studies show that the dynamical enzyme motions are suppressed. Our approach provides a unique avenue to modulating enzyme activity through changes in enzyme dynamics. Further it provides vital insights that show the altered motions of DHFR cause significant changes in the enzyme's ability to access its functionally relevant conformational substates, explaining the decreased k_hydride rate. This approach has important implications for obtaining fundamental insights into the role of rate-limiting dynamics in catalysis and as well as for enzyme engineering.

Increased substrate affinity in the Escherichia coli L28R dihydrofolate reductase mutant causes trimethoprim resistance

Dihydrofolate reductase (DHFR) is a ubiquitous enzyme with an essential role in cell metabolism. DHFR catalyzes the reduction of dihydrofolate to tetrahydrofolate, which is a precursor for purine and thymidylate synthesis. Several DHFR targeting antifolate drugs including trimethoprim, a competitive antibacterial inhibitor, have therefore been developed and are clinically used. Evolution of resistance against antifolates is a common public health problem rendering these drugs ineffective. To combat the resistance problem, it is important to understand resistance-conferring changes in the DHFR structure and accordingly develop alternative strategies. Here, we structurally and dynamically characterize Escherichia coli DHFR in its wild type (WT) and trimethoprim resistant L28R mutant forms in the presence of the substrate and its inhibitor trimethoprim. We use molecular dynamics simulations to determine the conformational space, loop dynamics and hydrogen bond distributions at the active site of DHFR for the WT and the L28R mutant. We also report their experimental kcat, Km, and Ki values, accompanied by isothermal titration calorimetry measurements of DHFR that distinguish enthalpic and entropic contributions to trimethoprim binding. Although mutations that confer resistance to competitive inhibitors typically make enzymes more promiscuous and decrease affinity to both the substrate and the inhibitor, strikingly, we find that the L28R mutant has a unique resistance mechanism. While the binding affinity differences between the WT and the mutant for the inhibitor and the substrate are small, the newly formed extra hydrogen bonds with the aminobenzoyl glutamate tail of DHF in the L28R mutant leads to increased barriers for the dissociation of the substrate and the product. Therefore, the L28R mutant indirectly gains resistance by enjoying prolonged binding times in the enzyme-substrate complex. While this also leads to slower product release and decreases the catalytic rate of the L28R mutant, the overall effect is the maintenance of a sufficient product formation rate. Finally, the experimental and computational analyses together reveal the changes that occur in the energetic landscape of DHFR upon the resistance-conferring L28R mutation. We show that the negative entropy associated with the binding of trimethoprim in WT DHFR is due to water organization at the binding interface. Our study lays the framework to study structural changes in other trimethoprim resistant DHFR mutants.

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.

Review of NAD(P)H-dependent oxidoreductases: Properties, engineering and application

NAD(P)H-dependent oxidoreductases catalyze the reduction or oxidation of a substrate coupled to the oxidation or reduction, respectively, of a nicotinamide adenine dinucleotide cofactor NAD(P)H or NAD(P)⁺. NAD(P)H-dependent oxidoreductases catalyze a large variety of reactions and play a pivotal role in many central metabolic pathways. Due to the high activity, regiospecificity and stereospecificity with which they catalyze redox reactions, they have been used as key components in a wide range of applications, including substrate utilization, the synthesis of chemicals, biodegradation and detoxification. There is great interest in tailoring NAD(P)H-dependent oxidoreductases to make them more suitable for particular applications. Here, we review the main properties and classes of NAD(P)H-dependent oxidoreductases, the types of reactions they catalyze, some of the main protein engineering techniques used to modify their properties and some interesting examples of their modification and application.

Chemoenzymatic Assembly of Isotopically Labeled Folates

Pterin-containing natural products have diverse functions in life, but an efficient and easy scheme for their in vitro synthesis is not available. Here we report a chemoenzymatic 14-step, one-pot synthesis that can be used to generate ¹³C- and ¹⁵N-labeled dihydrofolates (H₂F) from glucose, guanine, and p-aminobenzoyl-l-glutamic acid. This synthesis stands out from previous approaches to produce H₂F in that the average yield of each step is >91% and it requires only a single purification step. The use of a one-pot reaction allowed us to overcome potential problems with individual steps during the synthesis. The availability of labeled dihydrofolates allowed the measurement of heavy-atom isotope effects for the reaction catalyzed by the drug target dihydrofolate reductase and established that protonation at N5 of H₂F and hydride transfer to C6 occur in a stepwise mechanism. This chemoenzymatic pterin synthesis can be applied to the efficient production of other folates and a range of other natural compounds with applications in nutritional, medical, and cell-biological research.

Measurement of Enzyme Isotope Effects

Enzyme isotope effects, or the kinetic effects of "heavy" enzymes, refer to the effect of isotopically labeled protein residues on the enzyme's activity or physical properties. These effects are increasingly employed in the examination of the possible contributions of protein dynamics to enzyme catalysis. One hypothesis assumed that isotopic substitution of all ¹²C, ¹⁴N, and nonexchangeable ¹H by ¹³C, ¹⁵N, and ²H, would slow down protein picosecond to femtosecond dynamics without any effect on the system's electrostatics following the Born-Oppenheimer approximation. It was suggested that reduced reaction rates reported for several "heavy" enzymes accords with that hypothesis. However, numerous deviations from the predictions of that hypothesis were also reported. Current studies also attempt to test the role of individual residues by site-specific labeling or by labeling a pattern of residues on activity. It appears that in several systems the protein's fast dynamics are indeed reduced in "heavy" enzymes in a way that reduces the probability of barrier crossing of its chemical step. Other observations, however, indicated that slower protein dynamics are electrostatically altered in isotopically labeled enzymes. Interestingly, these effects appear to be system dependent, thus it might be premature to suggest a general role of "heavy" enzymes' effect on catalysis.

The role of the Met20 loop in the hydride transfer in Escherichia coli dihydrofolate reductase

A key question concerning the catalytic cycle of Escherichia coli dihydrofolate reductase (ecDHFR) is whether the Met²⁰ loop is dynamically coupled to the chemical step during catalysis. A more basic, yet unanswered question is whether the Met²⁰ loop adopts a closed conformation during the chemical hydride transfer step. To examine the most likely conformation of the Met²⁰ loop during the chemical step, we studied the hydride transfer in wild type (WT) ecDHFR using hybrid quantum mechanics-molecular mechanics free energy simulations with the Met²⁰ loop in a closed and disordered conformation. Additionally, we investigated three mutant forms (I14X; X = Val, Ala, Gly) of the enzyme that have increased active site flexibility and donor-acceptor distance dynamics in closed and disordered Met²⁰ loop states. We found that the conformation of the Met²⁰ loop has a dramatic effect on the ordering of active site hydration, although the Met²⁰ loop conformation only has a moderate effect on the hydride transfer rate and donor-acceptor distance dynamics. Finally, we evaluated the pK_a of the substrate N5 position in closed and disordered Met²⁰ loop states and found a strong correlation between N5 basicity and the conformation of the Met²⁰ loop.

Improved Sugar Puckering Profiles for Nicotinamide Ribonucleoside for Hybrid QM/MM Simulations

The coenzyme nicotinamide adenine dinucleotide (NAD⁺) and its reduced form (NADH) play ubiquitous roles as oxidizing and reducing agents in nature. The binding, and possibly the chemical redox step, of NAD⁺/NADH may be influenced by the cofactor conformational distribution and, in particular, by the ribose puckering of its nicotinamide-ribonucleoside (NR) moiety. In many hybrid quantum mechanics-molecular mechanics (QM/MM) studies of NAD⁺/NADH dependent enzymes, the QM region is treated by semiempirical (SE) methods. Recent work suggests that SE methods do not adequately describe the ring puckering in sugar molecules. In the present work we adopt an efficient and practical strategy to correct for this deficiency for NAD⁺/NADH. We have implemented a cost-effective correction to a SE Hamiltonian by adding a correction potential, which is defined as the difference between an accurate benchmark density functional theory (DFT) potential energy surface (PES) and the SE PES. In practice, this is implemented via a B-spline interpolation scheme for the grid-based potential energy difference surface. We find that the puckering population distributions obtained from free energy QM(SE)/MM simulations are in good agreement with DFT and in fair accord with experimental results. The corrected PES should facilitate a more accurate description of the ribose puckering in the NAD⁺/NADH cofactor in simulations of biological systems.

Examinations of the Chemical Step in Enzyme Catalysis

Advances in computational and experimental methods in enzymology have aided comprehension of enzyme-catalyzed chemical reactions. The main difficulty in comparing computational findings to rate measurements is that the first examines a single energy barrier, while the second frequently reflects a combination of many microscopic barriers. We present here intrinsic kinetic isotope effects and their temperature dependence as a useful experimental probe of a single chemical step in a complex kinetic cascade. Computational predictions are tested by this method for two model enzymes: dihydrofolate reductase and thymidylate synthase. The description highlights the significance of collaboration between experimentalists and theoreticians to develop a better understanding of enzyme-catalyzed chemical conversions.

The Effect of Protein Mass Modulation on Human Dihydrofolate Reductase

Dihydrofolate reductase (DHFR) from Escherichia coli has long served as a model enzyme with which to elucidate possible links between protein dynamics and the catalyzed reaction. Such physical properties of its human counterpart have not been rigorously studied so far, but recent computer-based simulations suggest that these two DHFRs differ significantly in how closely coupled the protein dynamics and the catalyzed C-H → C hydride transfer step are. To test this prediction, two contemporary probes for studying the effect of protein dynamics on catalysis were combined here: temperature dependence of intrinsic kinetic isotope effects (KIEs), which are sensitive to the physical nature of the chemical step, and protein mass modulation, which slows down fast dynamics (femto- to picosecond time scale) throughout the protein. The intrinsic H/T KIEs of human DHFR, like those of E. coli DHFR, are shown to be temperature-independent in the range from 5 to 45 °C, indicating fast sampling of donor and acceptor distances (DADs) at the reaction's transition state (or tunneling ready state, TRS). Mass modulation of these enzymes through isotopic labeling with (13)C, (15)N, and (2)H at nonexchangeable hydrogens yields an 11% heavier enzyme. The additional mass has no effect on the intrinsic KIEs of the human enzyme. This finding indicates that the mass modulation of the human DHFR affects neither DAD distribution nor the DAD's conformational sampling dynamics. Furthermore, reduction in the enzymatic turnover number and the dissociation rate constant for the product indicate that the isotopic substitution affects kinetic steps that are not the catalyzed C-H → C hydride transfer. The findings are discussed in terms of fast dynamics and their role in catalysis, the comparison of calculations and experiments, and the interpretation of isotopically modulated heavy enzymes in general.

Hydride Transfer in DHFR by Transition Path Sampling, Kinetic Isotope Effects, and Heavy Enzyme Studies

Escherichia coli dihydrofolate reductase (ecDHFR) is used to study fundamental principles of enzyme catalysis. It remains controversial whether fast protein motions are coupled to the hydride transfer catalyzed by ecDHFR. Previous studies with heavy ecDHFR proteins labeled with (13)C, (15)N, and nonexchangeable (2)H reported enzyme mass-dependent hydride transfer kinetics for ecDHFR. Here, we report refined experimental and computational studies to establish that hydride transfer is independent of protein mass. Instead, we found the rate constant for substrate dissociation to be faster for heavy DHFR. Previously reported kinetic differences between light and heavy DHFRs likely arise from kinetic steps other than the chemical step. This study confirms that fast (femtosecond to picosecond) protein motions in ecDHFR are not coupled to hydride transfer and provides an integrative computational and experimental approach to resolve fast dynamics coupled to chemical steps in enzyme catalysis.

Dihydrofolate reductase as a model for studies of enzyme dynamics and catalysis

Dihydrofolate reductase from Escherichia coli (ecDHFR) serves as a model system for investigating the role of protein dynamics in enzyme catalysis. We discuss calculations predicting a network of dynamic motions that is coupled to the chemical step catalyzed by this enzyme. Kinetic studies testing these predictions are presented, and their potential use in better understanding the role of these dynamics in enzyme catalysis is considered. The cumulative results implicate motions across the entire protein in catalysis.

Vacuum-Ultraviolet Circular Dichroism Spectra of Escherichia coli Dihydrofolate Reductase and Its Mutants: Contributions of Phenylalanine and Tyrosine Side Chains and Exciton Coupling of Two Tryptophan Side Chains

Vacuum-ultraviolet (VUV) circular dichroism (CD) spectroscopy has recently been used for secondary structure analysis of proteins; however, the contribution of aromatic side chains to protein VUV CD spectra is unresolved. In this report, VUV CD spectra of 10 Escherichia coli dihydrofolate reductase (DHFR) mutants, in which each phenylalanine or tyrosine residue was mutated to leucine, were measured down to 175 nm at 25 °C and pH 8.0 to elucidate the contributions of these aromatic side chains to the high-energy transitions of peptide bonds. The VUV CD spectra of these mutants were different from the spectrum of the wild-type protein, indicating that the contribution of the phenylalanine and tyrosine side chains of DHFR extends to the VUV region. Furthermore, the VUV CD spectrum and the folate- or NADP(+)-induced spectral change of F103L mutant DHFR indicated a modification and regeneration of exciton coupling between the Trp47 and Trp74 side chains, respectively, suggesting that exciton coupling may also contribute to the CD spectrum of DHFR in the VUV region. These results should be useful for theoretically characterizing the contribution of aromatic side chains to protein CD spectra, leading to the improvement of protein secondary-structure analysis by VUV CD spectroscopy.

Structural Characterization of Dihydrofolate Reductase Complexes by Top-Down Ultraviolet Photodissociation Mass Spectrometry

The stepwise reduction of dihydrofolate to tetrahydrofolate entails significant conformational changes of dihydrofolate reductase (DHFR). Binary and ternary complexes of DHFR containing cofactor NADPH, inhibitor methotrexate (MTX), or both NADPH and MTX were characterized by 193 nm ultraviolet photodissociation (UVPD) mass spectrometry. UVPD yielded over 80% sequence coverage of DHFR and resulted in production of fragment ions that revealed the interactions between DHFR and each ligand. UVPD of the binary DHFR·NADPH and DHFR·MTX complexes led to an unprecedented number of fragment ions containing either an N- or C-terminal protein fragment still bound to the ligand via retention of noncovalent interactions. In addition, holo-fragments retaining both ligands were observed upon UVPD of the ternary DHFR·NADPH·MTX complex. The combination of extensive holo and apo fragment ions allowed the locations of the NADPH and MTX ligands to be mapped, with NADPH associated with the adenosine binding domain of DHFR and MTX interacting with the loop domain. These findings are consistent with previous crystallographic evidence. Comparison of the backbone cleavage propensities for apo DHFR and its holo counterparts revealed significant variations in UVPD fragmentation in the regions expected to experience conformational changes upon binding NADPH, MTX, or both ligands. In particular, the subdomain rotation and loop movements, which are believed to occur upon formation of the transition state of the ternary complex, are reflected in the UVPD mass spectra. The UVPD spectra indicate enhanced backbone cleavages in regions that become more flexible or show suppressed backbone cleavages for those regions either shielded by the ligand or involved in new intramolecular interactions. This study corroborates the versatility of 193 nm UVPD mass spectrometry as a sensitive technique to track enzymatic cycles that involve conformational rearrangements.

Network of remote and local protein dynamics in dihydrofolate reductase catalysis

Molecular dynamics calculations and bionformatic studies of dihydrofolate reductase (DHFR) have suggested a network of coupled motions across the whole protein that is correlated to the reaction coordinate. Experimental studies demonstrated that distal residues G121, M42 and F125 in E. coli DHFR participate in that network. The missing link in our understanding of DHFR catalysis is the lack of a mechanism by which such remote residues can affect the catalyzed chemistry at the active site. Here, we present a study of the temperature dependence of intrinsic kinetic isotope effects (KIEs) that indicates synergism between a remote residue in that dynamic network, G121, and the active site's residue I14. The intrinsic KIEs for the I14A-G121V double mutant showed steeper temperature dependence (ΔE_a(T-H)) than expected from comparison of the wild type and two single mutants. That effect was non-additive, i.e., ΔE_a(T-H)G121V +ΔE_{a(T-H) I14A} < ΔE_{a(T-H) double mutant}, which indicates a synergism between the two residues. This finding links the remote residues in the network under investigation to the enzyme's active site, providing a mechanism by which these residues can be coupled to the catalyzed chemistry. This experimental evidence validates calculations proposing that both remote and active site residues constitute a network of coupled promoting motions correlated to the bond activation step (C-H→C hydride transfer in this case). Additionally, the effect of I14A and G121V mutations on single turnover rates was additive rather than synergistic. Although single turnover rate measurements are more readily available and thus more popular than assessing intrinsic kinetic isotope effects, the current finding demonstrates that for these rates, which in DHFR reflect several microscopic rate constants, can fall short of revealing the nature of the C-H bond activation per se.

Insights into the slow-onset tight-binding inhibition of Escherichia coli dihydrofolate reductase: detailed mechanistic characterization of pyrrolo [3,2-f] quinazoline-1,3-diamine and its derivatives as novel tight-binding inhibitors

Dihydrofolate reductase (DHFR) is a pivotal enzyme involved in the de novo pathway of purine synthesis, and hence, represents an attractive target to disrupt systems that require rapid DNA turnover. The enzyme acquires resistance to available drugs by various molecular mechanisms, which necessitates the continuous discovery of novel antifolates. Previously, we identified a set of novel molecules that showed binding to E. coli DHFR by means of a thermal shift without establishing whether they inhibited the enzyme. Here, we show that a fraction of those molecules represent potent and novel inhibitors of DHFR activity. 7-[(4-aminophenyl)methyl]-7H-pyrrolo [3,2-f] quinazoline-1,3-diamine, a molecule with no reported inhibition of DHFR, potently inhibits the enzyme with a Ki value of 7.42 ± 0.92 nm by competitive displacement of the substrate dihydrofolic acid. It shows uncompetitive inhibition vis-à-vis NADPH, indicating that the inhibitor has markedly increased affinity for the NADPH-bound form of the enzyme. Further, we demonstrate that the mode of binding of the inhibitor to the enzyme-NADPH binary complex conforms to the slow-onset, tight-binding model. By contrast, mechanistic characterization of the parent molecule 7H-pyrrolo [3,2-f] quinazoline-1,3-diamine shows that lack of (4-aminophenyl)-methyl group at the seventh position abolishes the slow onset of inhibition. This finding provides novel insights into the role of substitutions on inhibitors of E. coli DHFR and represents the first detailed kinetic investigation of a novel diaminopyrroloquinazoline derivative on a prokaryotic DHFR. Furthermore, marked differences in the potency of inhibition for E. coli and human DHFR makes this molecule a promising candidate for development as an antibiotic.