The influence of active site conformations on the hydride transfer step of the thymidylate synthase reaction mechanism

The hydride transfer from C6 of tetrahydrofolate to the reaction's exocyclic methylene–dUMP intermediate is the rate limiting step in thymidylate synthase (TSase) catalysis. This step has been studied by means of QM/MM molecular dynamics simulations to generate the corresponding free energy surfaces. The use of two different initial X-ray structures has allowed exploring different conformational spaces and the existence of chemical paths with not only different reactivities but also different reaction mechanisms. The results confirm that this chemical conversion takes place preferentially via a concerted mechanism where the hydride transfer is conjugated to thiol-elimination from the product. The findings also confirm the labile character of the substrate–enzyme covalent bond established between the C6 of the nucleotide substrate and a conserved cysteine residue. The calculations also reproduce and rationalize a normal H/T 2° kinetic isotope effect measured for that step. From a computational point of view, the results demonstrate that the use of an incomplete number of coordinates to describe the real reaction coordinate can render biased results

QM/MM Calculations Suggest a Novel Intermediate Following the Proton Abstraction Catalyzed by Thymidylate Synthase

The cleavage of covalent C–H bonds is one of the most energetically demanding, yet biologically essential, chemical transformations. Two C–H bond cleavages are involved in the reaction catalyzed by thymidylate synthase (TSase), which provides the sole de novo source of thymidylate (i.e., the DNA base T) for most organisms. Our QM/MM free energy calculations show that the C–H → O proton transfer has three transition states that are energetically similar but structurally diverse. These characteristics are different from our previous calculation results on the C–H → C hydride transfer, providing an explanation for differences in temperature dependences of KIEs on these two C–H bond activation steps. The calculations also suggest that the traditionally proposed covalent bond between the protein and substrate (the C6–S bond) is very labile during the multistep catalytic reaction. Collective protein motions not only assist cleavage of the C6–S bond to stabilize the transition state of the proton transfer step but also rearrange the H-bond network at the end of this step to prepare the active site for subsequent chemical steps. These computational results illustrate functionalities of specific protein residues that reconcile many previous experimental observations and provide guidance for future experiments to examine the proposed mechanisms. The synchronized conformational changes in the protein and ligands observed in our simulations demonstrate participation of protein motions in the reaction coordinate of enzymatic reactions. Our computational findings suggest the existence of new reaction intermediates not covalently bound to TSase, which may lead to a new class of drugs targeting DNA biosynthesis

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) that are sensitive to the physical nature of the chemical step, and protein mass-modulation that slows down fast dynamics (femto- to picosecond timescale) 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–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 13C, 15N, and 2H at nonexchangeable hydrogens yield 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

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

Protein mass-modulated effects in the catalytic mechanism of dihydrofolate reductase : beyond promoting vibrations

This work was supported by NIH research grants GM068036 (V.L.S.) and GM65368 (A.K.), and NSF grant CHE-0133117 (A.K.).The role of fast protein dynamics in enzyme catalysis has been of great interest in the past decade. Recent “heavy enzyme” studies demonstrate that protein mass-modulated vibrations are linked to the energy barrier for the chemical step of catalyzed reactions. However, the role of fast dynamics in the overall catalytic mechanism of an enzyme has not been addressed. Protein mass-modulated effects in the catalytic mechanism of Escherichia coli dihydrofolate reductase (ecDHFR) are explored by isotopic substitution (13C, 15N, and non-exchangeable 2H) of the wild-type ecDHFR (l-DHFR) to generate a vibrationally perturbed “heavy ecDHFR” (h-DHFR). Steady-state, pre-steady-state, and ligand binding kinetics, intrinsic kinetic isotope effects (KIEint) on the chemical step, and thermal unfolding experiments of both l- and h-DHFR show that the altered protein mass affects the conformational ensembles and protein–ligand interactions, but does not affect the hydride transfer at physiological temperatures (25–45 °C). Below 25 °C, h-DHFR shows altered transition state (TS) structure and increased barrier-crossing probability of the chemical step compared with l-DHFR, indicating temperature-dependent protein vibrational coupling to the chemical step. Protein mass-modulated vibrations in ecDHFR are involved in TS interactions at cold temperatures and are linked to dynamic motions involved in ligand binding at physiological temperatures. Thus, mass effects can affect enzymatic catalysis beyond alterations in promoting vibrations linked to chemistry.Publisher PDFPeer reviewe

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.This work was supported by National Institutes of Health (NIH) Grant R01 GM024485 to R.M.S. and NIH Grant R01 GM65368 to A.K., the Spanish Ministerio de Ciencia, Innovación y Universidades (Grants PGC2018–094852–B–C21 and PID2019–107098RJ–I00), Generalitat Valenciana (Grants AICO/2019/195, SEJI/2020/007, and APOSTD/2020/015), and Universitat Jaume I (UJI–A2019–04 and UJI–B2020–03)

Mg 2+ Binds to the Surface of Thymidylate Synthase and Affects Hydride Transfer at the Interior Active Site

Thymidylate synthase (TSase) produces the sole intracellular de novo source of thymidine (i.e. the DNA base T) and thus is a common target for antibiotic and anticancer drugs. Mg2+ has been reported to affect TSase activity, but the mechanism of this interaction has not been investigated. Here we show that Mg2+ binds to the surface of Escherichia coli TSase and affects the kinetics of hydride transfer at the interior active site (16 Å away). Examination of the crystal structures identifies a Mg2+ near the glutamyl moiety of the folate cofactor, providing the first structural evidence for Mg2+ binding to TSase. The kinetics and NMR relaxation experiments suggest that the weak binding of Mg2+ to the protein surface stabilizes the closed conformation of the ternary enzyme complex and reduces the entropy of activation on the hydride transfer step. Mg2+ accelerates the hydride transfer by ca. 7-fold but does not affect the magnitude or temperature-dependence of the intrinsic kinetic isotope effect. These results suggest that Mg2+ facilitates the protein motions that bring the hydride donor and acceptor together, but it does not change the tunneling ready state of the hydride transfer. These findings highlight how variations in cellular Mg2+ concentration can modulate enzyme activity through long-range interactions in the protein, rather than binding at the active site. The interaction of Mg2+ with the glutamyl-tail of the folate cofactor and nonconserved residues of bacterial TSase may assist in designing antifolates with poly-glutamyl substitutes as species-specific antibiotic drugs

Dihydrofolate reductase as a model for studies of enzyme dynamics and catalysis [version 1; referees: 2 approved]

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