Electrostatic interactions play an important role in enzyme catalysis by guiding ligand binding and facilitating chemical reactions. These electrostatic interactions are modulated by conformational changes occurring over the catalytic cycle. Herein, the changes in active site electrostatic microenvironments are examined for all enzyme complexes along the catalytic cycle of Escherichia coli dihydrofolate reductase (ecDHFR) by incorporation of thiocyanate probes at two site-specific locations in the active site. The electrostatics and degree of hydration of the microenvironments surrounding the probes are investigated with spectroscopic techniques and mixed quantum mechanical/molecular mechanical (QM/MM) calculations. Changes in the electrostatic microenvironments along the catalytic environment lead to different nitrile (CN) vibrational stretching frequencies and 13C NMR chemical shifts. These environmental changes arise from protein conformational rearrangements during catalysis. The QM/MM calculations reproduce the experimentally measured vibrational frequency shifts of the thiocyanate probes across the catalyzed hydride transfer step, which spans the closed and occluded conformations of the enzyme. Analysis of the molecular dynamics trajectories provides insight into the conformational changes occurring between these two states and the resulting changes in classical electrostatics and specific hydrogen-bonding interactions. The electric fields along the CN axes of the probes are decomposed into contributions from specific residues, ligands, and solvent molecules that make up the microenvironments around the probes. Moreover, calculation of the electric field along the hydride donor–acceptor axis, along with decomposition of this field into specific contributions, indicates that the cofactor and substrate, as well as the enzyme, impose a substantial electric field that facilitates hydride transfer. Overall, experimental and theoretical data provide evidence for significant electrostatic changes in the active site microenvironments due to conformational motion occurring over the catalytic cycle of ecDHFR.

With the rapidly growing wealth of genomic data, experimental inquiries on the functional significance of important divergence sites in protein evolution are becoming more accessible. Here we trace the evolution of dihydrofolate reductase (DHFR) and identify multiple key divergence sites among 233 species between humans and bacteria. We connect these sites, experimentally and computationally, to changes in the enzyme's binding properties and catalytic efficiency. One of the identified evolutionarily important sites is the N23PP modification (∼mid-Devonian, 415-385 Mya), which alters the conformational states of the active site loop in Escherichia coli dihydrofolate reductase and negatively impacts catalysis. This enzyme activity was restored with the inclusion of an evolutionarily significant lid domain (G51PEKN in E. coli enzyme; ∼2.4 Gya). Guided by this evolutionary genomic analysis, we generated a human-like E. coli dihydrofolate reductase variant through three simple mutations despite only 26% sequence identity between native human and E. coli DHFRs. Molecular dynamics simulations indicate that the overall conformational motions of the protein within a common scaffold are retained throughout evolution, although subtle changes to the equilibrium conformational sampling altered the free energy barrier of the enzymatic reaction in some cases. The data presented here provide a glimpse into the evolutionary trajectory of functional DHFR through its protein sequence space that lead to the diverged binding and catalytic properties of the E. coli and human enzymes.phylogenetic | EVB

ConspectusEnzymes are essential for all living organisms, and their effectiveness as chemical catalysts has driven more than a half century of research seeking to understand the enormous rate enhancements they provide. Nevertheless, a complete understanding of the factors that govern the rate enhancements and selectivities of enzymes remains elusive, due to the extraordinary complexity and cooperativity that are the hallmarks of these biomolecules. We have used a combination of site-directed mutagenesis, pre-steady-state kinetics, X-ray crystallography, nuclear magnetic resonance (NMR), vibrational and fluorescence spectroscopies, resonance energy transfer, and computer simulations to study the implications of conformational motions and electrostatic interactions on enzyme catalysis in the enzyme dihydrofolate reductase (DHFR).We have demonstrated that modest equilibrium conformational changes are functionally related to the hydride transfer reaction. Results obtained for mutant DHFRs illustrated that reductions in hydride transfer rates are correlated with altered conformational motions, and analysis of the evolutionary history of DHFR indicated that mutations appear to have occurred to preserve both the hydride transfer rate and the associated conformational changes. More recent results suggested that differences in local electrostatic environments contribute to finely tuning the substrate pKa in the initial protonation step. Using a combination of primary and solvent kinetic isotope effects, we demonstrated that the reaction mechanism is consistent across a broad pH range, and computer simulations suggested that deprotonation of the active site Tyr100 may play a crucial role in substrate protonation at high pH.Site-specific incorporation of vibrational thiocyanate probes into the ecDHFR active site provided an experimental tool for interrogating these microenvironments and for investigating changes in electrostatics along the DHFR catalytic cycle. Complementary molecular dynamics simulations in conjunction with mixed quantum mechanical/molecular mechanical calculations accurately reproduced the vibrational frequency shifts in these probes and provided atomic-level insight into the residues influencing these changes. Our findings indicate that conformational and electrostatic changes are intimately related and functionally essential. This approach can be readily extended to the study of other enzyme systems to identify more general trends in the relationship between conformational fluctuations and electrostatic interactions. These results are relevant to researchers seeking to design novel enzymes as well as those seeking to develop therapeutic agents that function as enzyme inhibitors.

The reaction catalyzed by Escherichia coli dihydrofolate reductase (ecDHFR) has become a model for understanding enzyme catalysis, and yet several details of its mechanism are still unresolved. Specifically, the mechanism of the chemical step, the hydride transfer reaction, is not fully resolved. We found, unexpectedly, the presence of two reactive ternary complexes [enzyme:NADPH:7,8-dihydrofolate (E:NADPH:DHF)] separated by one ionization event. Furthermore, multiple kinetic isotope effect (KIE) studies revealed a stepwise mechanism in which protonation of the DHF precedes the hydride transfer from the nicotinamide cofactor (NADPH) for both reactive ternary complexes of the WT enzyme. This mechanism was supported by the pH-and temperature-independent intrinsic KIEs for the C-H→C hydride transfer between NADPH and the preprotonated DHF. Moreover, we showed that active site residues D27 and Y100 play a synergistic role in facilitating both the proton transfer and subsequent hydride transfer steps. Although D27 appears to have a greater effect on the overall rate of conversion of DHF to tetrahydrofolate, Y100 plays an important electrostatic role in modulating the pK a of the N5 of DHF to enable the preprotonation of DHF by an active site water molecule.kinetic isotope effect | mechanism | dihydrofolate reductase | synergism E scherichia coli dihydrofolate reductase (ecDHFR) catalyzes the reduction of 7,8-dihydrofolate (DHF) to 5,6,7,8-tetrahydrofolate (THF) through the transfer of a proR hydride from the C4 atom of NADPH to the C6 position of the dihydropterin ring of DHF (1). This enzyme is critical in maintaining the intercellular pool of DHF, which is subsequently used in the biosynthesis of purine nucleotides and some amino acids. Given its biological importance, DHFR is an important drug target (2, 3), and its function has been extensively studied (4-12). However, several fundamental mechanistic details remained incomplete. One unresolved issue is the chronological order of the DHF protonation and hydride transfer steps. Quantum mechanical/molecular mechanical (QM/MM) calculations favor a stepwise protonation-hydride transfer reaction mechanism (10, 13, 14), whereas recent experimental data were interpreted as a change in the reaction mechanism between a concerted process and a stepwise process as a function of pH (15).The nature of DHF N5 protonation is also not fully understood. D27 has been implicated in the protonation of the N5 position of NADPH (4,12,14,(16)(17)(18)) through a linked water molecule. However, it is not clear how the water molecule-promoted DHF protonation at the N5 of the pterin can change in response to protein conformation (19) and to different electrostatic environments. Recently, the Y100 residue, which is located only ca. 3.5 Å from both the amide of the nicotinamide and N8 of the pterin ring (Fig. 1), has been shown to play an important role in electrostatically facilitating the ecDHFR-catalyzed reaction (20, 21). The location of Y100 suggests that it may form hydrogen bonds with ...

No abstract