No abstract

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.

The vibrational Stark effect provides insight into the roles of hydrogen bonding, electrostatics, and conformational motions in enzyme catalysis. In a recent application of this approach to the enzyme ketosteroid isomerase (KSI), thiocyanate probes were introduced in site-specific positions throughout the active site. This paper implements a quantum mechanical/molecular mechanical (QM/MM) approach for calculating the vibrational shifts of nitrile (CN) probes in proteins. This methodology is shown to reproduce the experimentally measured vibrational shifts upon binding of the intermediate analog equilinen to KSI for two different nitrile probe positions. Analysis of the molecular dynamics simulations provides atomistic insight into the roles that key residues play in determining the electrostatic environment and hydrogen-bonding interactions experienced by the nitrile probe. For the M116C-CN probe, equilinen binding reorients an active site water molecule that is directly hydrogen bonded to the nitrile probe, resulting in a more linear CNH angle and increasing the CN frequency upon binding. For the F86C-CN probe, equilinen binding orients the Asp103 residue, decreasing the hydrogen-bonding distance between the Asp103 backbone and the nitrile probe and slightly increasing the CN frequency. This QM/MM methodology is applicable to a wide range of biological systems and has the potential to assist in the elucidation of the fundamental principles underlying enzyme catalysis.

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 ...

A quantitative connection between molecular dynamics simulations and vibrational spectroscopy of probe-labeled systems would enable direct translation of experimental data into structural and dynamical information. To constitute this connection, all-atom molecular dynamics (MD) simulations were performed for two SCN probe sites (solvent-exposed and buried) in a calmodulin-target peptide complex. Two frequency calculation approaches with substantial nonelectrostatic components, a quantum mechanics/molecular mechanics (QM/MM)-based technique and a solvatochromic fragment potential (SolEFP) approach, were used to simulate the infrared probe line shapes. While QM/MM results disagreed with experiment, SolEFP results matched experimental frequencies and line shapes and revealed the physical and dynamic bases for the observed spectroscopic behavior. The main determinant of the CN probe frequency is the exchange repulsion between the probe and its local structural neighbors, and there is a clear dynamic explanation for the relatively broad probe line shape observed at the "buried" probe site. This methodology should be widely applicable to vibrational probes in many environments.