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.

Caught in Action: X-ray Structure of Thymidylate Synthase with Noncovalent Intermediate Analog.

Methylation of 2-deoxyuridine-5'-monophosphate (dUMP) at the C5 position by the obligate dimeric thymidylate synthase (TSase) in the sole de novo biosynthetic pathway to thymidine 5'-monophosphate (dTMP) proceeds by forming a covalent ternary complex with dUMP and cosubstrate 5,10-methylenetetrahydrofolate. The crystal structure of an analog of this intermediate gives important mechanistic insights but does not explain the half-of-the-sites activity of the enzyme. Recent experiments showed that the C5 proton and the catalytic Cys are eliminated in a concerted manner from the covalent ternary complex to produce a noncovalent bisubstrate intermediate. Here, we report the crystal structure of TSase with a close synthetic analog of this intermediate in which it has partially reacted with the enzyme but in only one protomer, consistent with the half-of-the-sites activity of this enzyme. Quantum mechanics/molecular mechanics simulations confirmed that the analog could undergo catalysis. The crystal structure shows a new water 2.9 Å from the critical C5 of the dUMP moiety, which in conjunction with other residues in the network, may be the elusive general base that abstracts the C5 proton of dUMP during the reaction.

Rigid Double-Stranded DNA Linkers for Single Molecule Enzyme-Drug Interaction Measurements Using Molecular Recognition Force Spectroscopy.

Single-molecule studies can reveal the distribution of states and interactions between ligand-enzyme complexes not accessible for most studies that measure a large ensemble average response of many molecules. Furthermore, in some biological applications, the information regarding the outliers, not the average of measured properties, can be more important. The high spatial and force resolution provided by atomic force microscopy (AFM) under physiological conditions has been utilized in this study to quantify the force-distance relations of enzyme-drug interactions. Different immobilization techniques of the protein to a surface and the drug to AFM tip were quantitatively compared to improve the accuracy and precision of the measurement. Protein that is directly bound to the surface, forming a monolayer, was compared to enzyme molecules bound to the surface with rigid double-stranded (ds) DNA spacers. These surfaces immobilization techniques were studied with the drug bound directly to the AFM tip and drug bound via flexible poly(ethylene glycol) and rigid dsDNA linkers. The activity of the enzyme was found to be not significantly altered by immobilization methods relative to its activity in solution. The findings indicate that the approach for studying drug-enzyme interaction based on rigid dsDNA linker on the surface and either flexible or rigid linker on the tip affords straightforward, highly specific, reproducible, and accurate force measurements with a potential for single-molecule level studies. The method could facilitate in-depth examination of a broad spectrum of biological targets and potential drugs.

Isotopic Labeling of Formate Dehydrogenase Perturbs the Protein Dynamics.

Isotope substitution of enzymes has become a means of addressing the participation of protein motions in enzyme-catalyzed reactions. The idea is that only the enzyme mass will be altered and not the electrostatics, so that the protein dynamics are essentially the same but at lower frequencies because of the mass change. In this study, we variably label all carbon atoms in formate dehydrogenase (FDH) with 13C, all nitrogen atoms with 15N, and all nonexchangeable hydrogen atoms with deuterium and investigate the impact that isotopic substitution has on the dynamics at the active site by two-dimensional infrared spectroscopy and compare with the measurements of the temperature dependence of the intrinsic kinetic isotope effects (KIEs). We show that 15N labeling of FDH has the largest effect and makes the active site more heterogeneous, whereas the addition of nonexchangeable deuterium appears to have the opposite effect of 15N on active-site dynamics, resulting in a behavior similar to that of native FDH. Nevertheless, the temperature dependence of the KIEs shows a monotonic trend with protein mass that does not correspond with the changes in dynamics. These results suggest that isotope labeling has more than just a mass effect on enzyme dynamics and may influence electrostatics in ways that complicate the interpretation of the protein isotope effect.

Evolution Conserves the Network of Coupled Residues in Dihydrofolate Reductase.

Understanding protein motions and their role in enzymatic reactions is an important and timely topic in enzymology. Protein motions that are involved in the chemical step of catalysis are particularly intriguing but difficult to identify. A global network of coupled residues in Escherichia coli dihydrofolate reductase (E. coli DHFR), which assists in catalyzing the chemical step, has previously been demonstrated through quantum mechanical/molecular mechanical and molecular dynamics simulations as well as bioinformatic analyses. A few specific residues (M42, G121, F125, and I14) were shown to function synergistically with measurements of single-turnover rates and the temperature dependence of intrinsic kinetic isotope effects (KIEsint) of site-directed mutants. This study hypothesizes that the global network of residues involved in the chemical step is evolutionarily conserved and probes homologous residues of the potential global network in human DHFR through measurements of the temperature dependence of KIEsint and computer simulations based on the empirical valence bond method. We study mutants M53W and S145V. Both of these remote residues are homologous to network residues in E. coli DHFR. Non-additive isotope effects on activation energy are observed between M53 and S145, indicating their synergistic effect on the chemical step in human DHFR, which suggests that both of these residues are part of a network affecting the chemical step in enzyme catalysis. This finding supports the hypothesis that human and E. coli DHFR share similar networks, consistent with evolutionary preservation of such networks.

Oscillatory Active-site Motions Correlate with Kinetic Isotope Effects in Formate Dehydrogenase.

Thermal motions of enzymes have been invoked to explain the temperature dependence of kinetic isotope effects (KIE) in enzyme-catalyzed hydride transfers. Formate dehydrogenase (FDH) from Candida boidinii exhibits a temperature independent KIE that becomes temperature dependent upon mutation of hydrophobic residues in the active site. Ternary complexes of FDH that mimic the transition state structure allow investigation of how these mutations influence active-site dynamics. A combination of X-ray crystallography, two-dimensional infrared (2D IR) spectroscopy, and molecular dynamic simulations characterize the structure and dynamics of the active site. FDH exhibits oscillatory frequency fluctuations on the picosecond timescale, and the amplitude of these fluctuations correlates with the temperature dependence of the KIE. Both the kinetic and dynamic phenomena can be reproduced computationally. These results provide experimental evidence for a connection between the temperature dependence of KIEs and motions of the active site in an enzyme-catalyzed reaction consistent with activated tunneling models of the hydride transfer reaction.

Evolutionary Effects on Bound Substrate p Ka in Dihydrofolate Reductase.

In the present study, we address the effect of active site structure and dynamics of different dihydrofolate reductase (DHFR) isoforms on the p Ka of the bound substrate 7,8-dihydrofolate, in an attempt to understand possible evolutionary trends. We apply a hybrid QM/MM free energy perturbation method to estimate the p Ka of the N5 position of the bound substrate. We observe a gradual increase in N5 basicity as we move from primitive to more evolved DHFR isoforms. Structural analysis of these isoforms reveals a gradual sequestering of water molecules from the active site in the more evolved enzymes, thereby modulating the local dielectric environment near the substrate. Furthermore, the present study reveals a clear correlation between active site hydration and the N5 p Ka of the substrate. We emphasize the role of the M20 loop in controlling the active site hydration level, via a preorganized active site with a more hydrophobic environment and reduced loop flexibility as evolution progresses from bacterial to the human enzyme.

Parallel reaction pathways and noncovalent intermediates in thymidylate synthase revealed by experimental and computational tools.

Thymidylate synthase was one of the most studied enzymes due to its critical role in molecular pathogenesis of cancer. Nevertheless, many atomistic details of its chemical mechanism remain unknown or debated, thereby imposing limits on design of novel mechanism-based anticancer therapeutics. Here, we report unprecedented isolation and characterization of a previously proposed intact noncovalent bisubstrate intermediate formed in the reaction catalyzed by thymidylate synthase. Free-energy surfaces of the bisubstrate intermediates interconversions computed with quantum mechanics/molecular mechanics (QM/MM) methods and experimental assessment of the corresponding kinetics indicate that the species is the most abundant productive intermediate along the reaction coordinate, whereas accumulation of the covalent bisubstrate species largely occurs in a parallel nonproductive pathway. Our findings not only substantiate relevance of the previously proposed noncovalent intermediate but also support potential implications of the overstabilized covalent intermediate in drug design targeting DNA biosynthesis.

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.

Bacterial versus human thymidylate synthase: Kinetics and functionality.

Thymidylate Synthase (TSase) is a highly conserved enzyme that catalyzes the production of the DNA building block thymidylate. Structurally, functionally and mechanistically, bacterial and mammalian TSases share remarkable similarities. Because of this closeness, bacterial enzymes have long been used as model systems for human TSase. Furthermore, while TSase inhibitors have long served as chemotherapeutic drugs, no TSase inhibitor serves as an antibiotic. Despite their high resemblance, the mammalian TSases are distinct in a few known aspects, such as having a N-terminal tail and two insertions in the primary sequence and active/inactive conformations. Here, we aim to comprehensively characterize human (hs) TSase and delineate its contrasts and the similarities to the well-studied Escherichia coli (ec) TSase. We found that, in contrast to ecTSase, Mg2+ does not enhance reaction rates for hsTSase. The temperature dependence of intrinsic kinetic isotope effects (KIEs), on the other hand, suggests that Mg2+ has little or no impact on the transition state of hydride transfer in either enzyme, and that the transition state for the hydride transfer in hsTSase is looser than in ecTSase. Additionally, the substrates' binding order is strictly ordered for ecTSase but slightly less ordered for hsTSase. The observed kinetic and functional differences between bacterial and human enzymes may aid in the development of antibiotic drugs with reduced toxicity.

Chemo-enzymatic synthesis of the exocyclic olefin isomer of thymidine monophosphate.

Exocyclic olefin variants of thymidylate (dTMP) recently have been proposed as reaction intermediates for the thymidyl biosynthesis enzymes found in many pathogenic organisms, yet synthetic reports on these materials are lacking. Here we report two strategies to prepare the exocyclic olefin isomer of dTMP, which is a putative reaction intermediate in pathogenic thymidylate biosynthesis and a novel nucleotide analog. Our most effective strategy involves preserving the existing glyosidic bond of thymidine and manipulating the base to generate the exocyclic methylene moiety. We also report a successful enzymatic deoxyribosylation of a non-aromatic nucleobase isomer of thymine, which provides an additional strategy to access nucleotide analogs with disrupted ring conjugation or with reduced heterocyclic bases. The strategies reported here are straightforward and extendable towards the synthesis of various pyrimidine nucleotide analogs, which could lead to compounds of value in studies of enzyme reaction mechanisms or serve as templates for rational drug design.

Experimental and Computational Studies Delineate the Role of Asparagine 177 in Hydride Transfer for E. coli Thymidylate Synthase.

Thymidylate synthase (TSase), an enzyme responsible for the de novo biosynthesis of 2'-deoxythymidine 5'-monophosphate (thymidylate, dTMP) necessary for DNA synthesis, has been a drug target for decades. TSase is a highly conserved enzyme across species ranging from very primitive organisms to mammals. Among the many conserved active site residues, an asparagine (N177, using Escherichia coli residues numbering) appears to make direct hydrogen bonds with both the C4=O4 carbonyl of the 2'-deoxyuridine 5'-monophosphate (uridylate, dUMP) substrate and its pyrimidine ring's N3. Recent studies have reassessed the TSase catalytic mechanism, focusing on the degree of negative charge accumulation at the O4 carbonyl of the substrate during two critical H-transfers - a proton abstraction and a hydride transfer. To obtain insights into the role of this conserved N177 on the hydride transfer, we examined its aspartic acid (D) and serine (S) mutants - each of which is expected to alter hydrogen bonding and charge stabilization around the C4=O4 carbonyl of the 2'-deoxyuridine 5'-monophosphate (uridylate, dUMP) substrate. Steady-state kinetics, substrate binding order studies and temperature-dependency analysis of intrinsic KIEs for the hydride transfer step of the TSase catalytic cycle suggest the active site of N177D is not precisely organized for that step. A smaller disruption was observed for N177S, which could be rationalized by partial compensation by water molecules and rearrangement of other residues toward preparation of the system for the hydride transfer under study. These experimental findings are qualitatively mirrored by QM/MM computational simulations, thereby shedding light on the sequence and synchronicity of steps in the TSase-catalyzed reaction. This information could potentially inform the design of mechanism-based drugs targeting this enzyme.

Protein Mass Effects on Formate Dehydrogenase.

Isotopically labeled enzymes (denoted as "heavy" or "Born-Oppenheimer" enzymes) have been used to test the role of protein dynamics in catalysis. The original idea was that the protein's higher mass would reduce the frequency of its normal-modes without altering its electrostatics. Heavy enzymes have been used to test if the vibrations in the native enzyme are coupled to the chemistry it catalyzes, and different studies have resulted in ambiguous findings. Here the temperature-dependence of intrinsic kinetic isotope effects of the enzyme formate dehydrogenase is used to examine the distribution of H-donor to H-acceptor distance as a function of the protein's mass. The protein dynamics are altered in the heavy enzyme to diminish motions that determine the transition state sampling in the native enzyme, in accordance with a Born-Oppenheimer-like effect on bond activation. Findings of this work suggest components related to fast frequencies that can be explained by Born-Oppenheimer enzyme hypothesis (vibrational) and also slower time scale events that are non-Born-Oppenheimer in nature (electrostatic), based on evaluations of protein mass dependence of donor-acceptor distance and forward commitment to catalysis along with steady state and single turnover measurements. Together, the findings suggest that the mass modulation affected both local, fast, protein vibrations associated with the catalyzed chemistry and the protein's macromolecular electrostatics at slower time scales; that is, both Born-Oppenheimer and non-Born-Oppenheimer effects are observed. Comparison to previous studies leads to the conclusion that isotopic labeling of the protein may have different effects on different systems, however, making heavy enzyme studies a very exciting technique for exploring the dynamics link to catalysis in proteins.

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 12C, 14N, and nonexchangeable 1H by 13C, 15N, and 2H, 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.

Flavin-dependent thymidylate synthase: N5 of flavin as a Methylene carrier.

Thymidylate is synthesized de novo in all living organisms for replication of genomes. The chemical transformation is reductive methylation of deoxyuridylate at C5 to form deoxythymidylate. All eukaryotes including humans complete this well-understood transformation with thymidylate synthase utilizing 6R-N5-N10-methylene-5,6,7,8-tetrahydrofolate as both a source of methylene and a reducing hydride. In 2002, flavin-dependent thymidylate synthase was discovered as a new pathway for de novo thymidylate synthesis. The flavin-dependent catalytic mechanism is different than thymidylate synthase because it requires flavin as a reducing agent and methylene transporter. This catalytic mechanism is not well-understood, but since it is known to be very different from thymidylate synthase, there is potential for mechanism-based inhibitors that can selectively inhibit the flavin-dependent enzyme to target many human pathogens with low host toxicity.

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 Met20 loop is dynamically coupled to the chemical step during catalysis. A more basic, yet unanswered question is whether the Met20 loop adopts a closed conformation during the chemical hydride transfer step. To examine the most likely conformation of the Met20 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 Met20 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 Met20 loop states. We found that the conformation of the Met20 loop has a dramatic effect on the ordering of active site hydration, although the Met20 loop conformation only has a moderate effect on the hydride transfer rate and donor-acceptor distance dynamics. Finally, we evaluated the pKa of the substrate N5 position in closed and disordered Met20 loop states and found a strong correlation between N5 basicity and the conformation of the Met20 loop.

Benchmarking Quantum Mechanics/Molecular Mechanics (QM/MM) Methods on the Thymidylate Synthase-Catalyzed Hydride Transfer.

Given the ubiquity of hydride-transfer reactions in enzyme-catalyzed processes, identifying the appropriate computational method for evaluating such biological reactions is crucial to perform theoretical studies of these processes. In this paper, the hydride-transfer step catalyzed by thymidylate synthase (TSase) is studied by examining hybrid quantum mechanics/molecular mechanics (QM/MM) potentials via multiple semiempirical methods and the M06-2X hybrid density functional. Calculations of protium and tritium transfer in these reactions across a range of temperatures allowed calculation of the temperature dependence of kinetic isotope effects (KIE). Dynamics and quantum-tunneling effects are revealed to have little effect on the reaction rate, but are significant in determining the KIEs and their temperature dependence. A good agreement with experiments is found, especially when computed for RM1/MM simulations. The small temperature dependence of quantum tunneling corrections and the quasiclassical contribution term cancel each other, while the recrossing transmission coefficient seems to be temperature-independent over the interval of 5-40 °C.

Preserved hydride transfer mechanism in evolutionarily divergent thymidylate synthases.

Thymidylate synthase (TSase) catalyzes a hydride transfer in the last step of the de novo biosynthesis of the DNA nucleotide thymine. We compared two isozymes, namely, TSase from Escherichia coli (ecTSase) and TSase from Bacillus subtilis (bsTSase) that represent a case of divergent evolution. Interestingly, a highly conserved histidine (H147 of ecTSase) was proposed to serve a critical role in catalysis, but in bsTSase it is naturally substituted by valine (Val). Yet, bsTSase is more active than ecTSase, and the intrinsic kinetic isotope effects (KIEs) of both are temperature-independent, suggesting a similarly well-organized transition state (TS) for the catalyzed hydride transfer. To examine the role of that histidine (His) in TSase catalysis, we examined the kinetics of H147V ecTSase, which "bridges" between these two TSases. In contrast to both wild-type TSases, the single mutation results in deficient catalysis. The mutation leads to intrinsic KIEs that are temperature-dependent, indicating a substantial imperfection in its TS. The findings reveal two important features: a direct role of H147 in the hydride transfer step catalyzed by the ecTSase and the evolutionary compensation for its deficiency in bsTSase via extensive polymorphism across the protein. Very different active site residues are observed for these evolutionarily divergent isozymes, which result in a well-organized TS for both. It is suggested that evolutionary pressure compensated for the H to V substitution at the active site of bsTSase by polymorphism leading to a well-organized TS in both enzymes.

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.

Noncovalent Intermediate of Thymidylate Synthase: Fact or Fiction?

Thymidylate synthase is an attractive target for antibiotic and anticancer drugs due to its essential role in the de novo biosynthesis of the DNA nucleotide thymine. The enzymatic reaction is initiated by a nucleophilic activation of the substrate via formation of a covalent bond to an active site cysteine. The traditionally accepted mechanism is then followed by a series of covalently bound intermediates, where that bond is only cleaved upon product release. Recent computational and experimental studies suggest that the covalent bond between the protein and substrate is actually quite labile. Importantly, these findings predict the existence of a noncovalently bound bisubstrate intermediate, not previously anticipated, which could be the target of a novel class of drugs inhibiting DNA biosynthesis. Here we report the synthesis of the proposed intermediate and findings supporting its chemical and kinetic competence. These findings substantiate the predicted nontraditional mechanism and the potential of this intermediate as a new drug lead.