Hydride transfer versus hydrogen radical transfer in thymidylate synthase

Enhancing Protein Crystallization under a Magnetic Field

High-quality crystals are essential to ensure high-resolution structural information. Protein crystals are controlled by many factors, such as pH, temperature, and the ion concentration of crystalline solutions. We previously reported the development of a device dedicated to protein crystallization. In the current study, we have further modified and improved our device. Exposure to external magnetic field leads to alignment of the crystal toward a preferred direction depending on the magnetization energy. Each material has different magnetic susceptibilities depending on the individual direction of their unit crystal cells. One of the strategies to acquire a large crystal entails controlling the nucleation rate. Furthermore, exposure of a crystal to a magnetic field may lead to new morphologies by affecting the crystal volume, shape, and quality.

Thermodynamic Parameters of Elementary Steps for 3,5-Disubstituted 1,4-Dihydropyridines To Release Hydride Anions in Acetonitrile

A series of 3,5-disubstituted 1,4-dihydropyridine derivatives including the derivative with two chiral centers, 6H (R² = CH₃, CH₂Ph), as a new type of organic hydride source were synthesized and characterized. The thermodynamic driving forces (defined as enthalpy changes or standard redox potentials) of the 6 elementary steps for the organic hydrides to release hydride ions in acetonitrile were measured by isothermal titration calorimetry and electrochemical methods. The impacts of the substituents and functional groups bearing the N1 and C3/C5 positions on the thermodynamic driving forces of the 6 elementary steps were examined and analyzed. Moreover, the results showed that the reaction mechanism between the chiral organic hydride and activated ketone (ethyl benzoylformate) was identified as the concerted hydride transfer pathway based on the thermodynamic analysis platform. These valuable and crucial thermodynamic parameters will provide a broadly beneficial impact on the applications of 3,5-disubstituted 1,4-dihydropyridine derivatives in organic synthesis and pharmaceutical chemistry.

Where and How Does an Organic Molecule Having a C-X Bond Release X- Anion Like an Inorganic Compound? A Theoretical Study

In this research, at first, a comparative DFT study on hydride or fluoride release from a number of known orthoamides, their fluorine derivative having a central C-F bond and some simple organic compounds, is reported. The obtained data show that orthoamides release hydride or fluoride anions much easier than do other known organic compounds studied here. Interestingly, three simulated orthoamides having a central C-F bond, spontaneously and like an inorganic compound, release fluoride anion upon dissolving in polar solvents. The calculations confirmed that hyperconjugation interactions in the initial orthoamides facilitate the anion release from these compounds. However, the data clearly show that the proper overlap of an empty p orbital of the central carbon atom with adjacent lone-pair orbitals of nitrogen atoms ( lpN → lp*C or, in other words, the pπ-pπ interaction) in the resulting carbocations is a more important factor that must be taken into consideration when designing the hydride, fluoride, or other anion releasing agents. In addition, for the first time, a mechanism of hydride and fluoride removal from the orthoamides either through their reaction with BH₃ and BF₃ molecules, respectively, or upon their protonation is provided. To continue and for generalization of results to other groups attached to the central carbon atom of orthoamides, the reactions of H⁺ with orthoamides having the CH₃, OH, CN, NH₂, or C₅H₅ substituents were studied. The results showed that in all cases and in both gas and solution phases the above substituents easily leave the carbon atom as an anion and bond to H⁺. Thus, we have to conclude that upon the reaction of orthoamides with usual Lewis acids, many types of substituents on their central carbon atom can leave these compounds as anions.

Widespread Perturbation of Function, Structure, and Dynamics by a Conservative Single-Atom Substitution in Thymidylate Synthase

Thymidylate synthase (TSase) is responsible for synthesizing the sole de novo source of dTMP in all organisms. TSase is a drug target, and as such, it has been well studied in terms of both structure and reaction mechanism. Cysteine 146 in Escherichia coli TSase is universally conserved because it serves as the nucleophile in the enzyme mechanism. Here we use the C146S mutation to probe the role of the sulfur atom in early events in the catalytic cycle beyond serving as the nucleophile. Surprisingly, the single-atom substitution severely decreases substrate binding affinity, and the unfavorable ΔΔG°_bind is comprised of roughly equal enthalpic and entropic components at 25 °C. Chemical shifts in the free and dUMP-bound states show the mutation causes perturbations throughout TSase, including regions important for complex stability, in agreement with a less favorable enthalpy change. We measured the nuclear magnetic resonance methyl symmetry axis order parameter (S²_axis), a proxy for conformational entropy, for TSase at all vertices of the dUMP binding/C146S mutation thermodynamic cycle and found that the calculated TΔΔS°_conf is similar in sign and magnitude to the calorimetric TΔΔS°. Further, we ascribed minor resonances in wild-type-dUMP spectra to a state with a covalent bond between Sγ of C146 and C6 of dUMP and find S²_axis values are unaffected by covalent bond formation, indicating this reaction step is neutral with respect to ΔS°_conf. Lastly, the C146S mutation allowed us to measure cofactor analog binding by isothermal titration calorimetry without the confounding heat signature of covalent bond formation. Raltitrexed binds free and singly bound TSase with similar affinities, yet the two binding events have different enthalpy changes, providing further evidence of communication between the two active sites.

The bioinspired design of a reagent allows the functionalization of Cα-H of α,β-unsaturated carbonyl compounds via the Baylis-Hillman chemistry under ambient conditions

A rationally designed reagent capable of affecting alkylation at Cα of α,β-unsaturated carbonyl compounds is reported. The reaction proceeded at room temperature without any additives. The pH and H-bond formation during the reaction play a key role in the working of the reagent.

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.

Activation of Two Sequential H-transfers in the Thymidylate Synthase Catalyzed Reaction

Thymidylate synthase (TSase) catalyzes the de novo biosynthesis of thymidylate, a precursor for DNA, and is thus an important target for chemotherapeutics and antibiotics. Two sequential C-H bond cleavages catalyzed by TSase are of particular interest: a reversible proton abstraction from the 2'-deoxy-uridylate substrate, followed by an irreversible hydride transfer forming the thymidylate product. QM/MM calculations of the former predicted a mechanism where the abstraction of the proton leads to formation of a novel nucleotide-folate intermediate that is not covalently bound to the enzyme (Wang, Z.; Ferrer, S.; Moliner, V.; Kohen, A. Biochemistry2013, 52, 2348-2358). Existence of such intermediate would hold promise as a target for a new class of drugs. Calculations of the subsequent hydride transfer predicted a concerted H-transfer and elimination of the enzymatic cysteine (Kanaan, N.; Ferrer, S.; Marti, S.; Garcia-Viloca, M.; Kohen, A.; Moliner, V. J. Am. Chem. Soc.2011, 133, 6692-6702). A key to both C-H activations is a highly conserved arginine (R166) that stabilizes the transition state of both H-transfers. Here we test these predictions by studying the R166 to lysine mutant of E. coli TSase (R166K) using intrinsic kinetic isotope effects (KIEs) and their temperature dependence to assess effects of the mutation on both chemical steps. The findings confirmed the predictions made by the QM/MM calculations, implicate R166 as an integral component of both reaction coordinates, and thus provide critical support to the nucleotide-folate intermediate as a new target for rational drug design.

Thymidylate synthase inspired biomodel reagent for the conversion of uracil to thymine

Inspired by TSase catalysis for dUMP conversion to dTMP, a biomodel reagent is developed. The presence of NH2, Gly-(S)-Cys and (S)-oxiran methyl, at C5, C4 and N-10 of acridine, respectively, in addition to the pH of the reaction mixture, allows for good complementary inter- and intra-molecular interactions and chiral discrimination for the reagent to achieve conversion of uracil to thymine.

Concerted versus stepwise mechanism in thymidylate synthase

Thymidylate synthase (TSase) catalyzes the intracellular de novo formation of thymidylate (a DNA building block) in most living organisms, making it a common target for chemotherapeutic and antibiotic drugs. Two mechanisms have been proposed for the rate-limiting hydride transfer step in TSase catalysis: a stepwise mechanism in which the hydride transfer precedes the cleavage of the covalent bond between the enzymatic cysteine and the product and a mechanism where both happen concertedly. Striking similarities between the enzyme-bound enolate intermediates formed in the initial and final step of the reaction supported the first mechanism, while QM/MM calculations favored the concerted mechanism. Here, we experimentally test these two possibilities using secondary kinetic isotope effect (KIE), mutagenesis study, and primary KIEs. The findings support the concerted mechanism and demonstrate the critical role of an active site arginine in substrate binding, activation of enzymatic nucleophile, and the hydride transfer studied here. The elucidation of this reduction/substitution sheds light on the critical catalytic step in TSase and may aid future drug or biomimetic catalyst design.

Isotope effects as probes for enzyme catalyzed hydrogen-transfer reactions

Kinetic Isotope effects (KIEs) have long served as a probe for the mechanisms of both enzymatic and solution reactions. Here, we discuss various models for the physical sources of KIEs, how experimentalists can use those models to interpret their data, and how the focus of traditional models has grown to a model that includes motion of the enzyme and quantum mechanical nuclear tunneling. We then present two case studies of enzymes, thymidylate synthase and alcohol dehydrogenase, and discuss how KIEs have shed light on the C-H bond cleavages those enzymes catalyze. We will show how the combination of both experimental and computational studies has changed our notion of how these enzymes exert their catalytic powers.

Relationship of femtosecond-picosecond dynamics to enzyme-catalyzed H-transfer

At physiological temperatures, enzymes exhibit a broad spectrum of conformations, which interchange via thermally activated dynamics. These conformations are sampled differently in different complexes of the protein and its ligands, and the dynamics of exchange between these conformers depends on the mass of the group that is moving and the length scale of the motion, as well as restrictions imposed by the globular fold of the enzymatic complex. Many of these motions have been examined and their role in the enzyme function illuminated, yet most experimental tools applied so far have identified dynamics at time scales of seconds to nanoseconds, which are much slower than the time scale for H-transfer between two heavy atoms. This chemical conversion and other processes involving cleavage of covalent bonds occur on picosecond to femtosecond time scales, where slower processes mask both the kinetics and dynamics. Here we present a combination of kinetic and spectroscopic methods that may enable closer examination of the relationship between enzymatic C-H → C transfer and the dynamics of the active site environment at the chemically relevant time scale. These methods include kinetic isotope effects and their temperature dependence, which are used to study the kinetic nature of the H-transfer, and 2D IR spectroscopy, which is used to study the dynamics of transition-state- and ground-state-analog complexes. The combination of these tools is likely to provide a new approach to examine the protein dynamics that directly influence the chemical conversion catalyzed by enzymes.

Experimental and theoretical studies of enzyme-catalyzed hydrogen-transfer reactions

The mechanisms of enzyme-catalyzed reactions are medicinally important and present a fascinating intellectual challenge. Many experimental and theoretical techniques can shed light on these mechanisms, and here, we shall focus on the utility of kinetic isotope effects (KIEs) to study enzymatic reactions that involve hydrogen transfers. We will provide a short background on the prevailing models to interpret KIEs and then present more detailed reviews of two model enzymes: alcohol dehydrogenase and thymidylate synthase. These two examples provide a context to describe the types of experiments and theoretical calculations that drive this field forward and the kind of information each can furnish. We emphasize the importance of cooperation between experimentalists and theoreticians to continue the progress toward a comprehensive theory of enzyme catalysis.

A remote mutation affects the hydride transfer by disrupting concerted protein motions in thymidylate synthase

The role of protein flexibility in enzyme-catalyzed activation of chemical bonds is an evolving perspective in enzymology. Here we examine the role of protein motions in the hydride transfer reaction catalyzed by thymidylate synthase (TSase). Being remote from the chemical reaction site, the Y209W mutation of Escherichia coli TSase significantly reduces the protein activity, despite the remarkable similarity between the crystal structures of the wild-type and mutant enzymes with ligands representing their Michaelis complexes. The most conspicuous difference between these two crystal structures is in the anisotropic B-factors, which indicate disruption of the correlated atomic vibrations of protein residues in the mutant. This dynamically altered mutant allows a variety of small thiols to compete for the reaction intermediate that precedes the hydride transfer, indicating disruption of motions that preorganize the protein environment for this chemical step. Although the mutation causes higher enthalpy of activation of the hydride transfer, it only shows a small effect on the temperature dependence of the intrinsic KIE, suggesting marginal changes in the geometry and dynamics of the H-donor and -acceptor at the tunneling ready state. These observations suggest that the mutation disrupts the concerted motions that bring the H-donor and -acceptor together during the pre- and re-organization of the protein environment. The integrated structural and kinetic data allow us to probe the impact of protein motions on different time scales of the hydride transfer reaction within a complex enzymatic mechanism.

Temperature dependence of the kinetic isotope effects in thymidylate synthase. A theoretical study

In recent years, the temperature dependence of primary kinetic isotope effects (KIE) has been used as indicator for the physical nature of enzyme-catalyzed H-transfer reactions. An interactive study where experimental data and calculations examine the same chemical transformation is a critical means to interpret more properly temperature dependence of KIEs. Here, the rate-limiting step of the thymidylate synthase-catalyzed reaction has been studied by means of hybrid quantum mechanics/molecular mechanics (QM/MM) simulations in the theoretical framework of the ensemble-averaged variational transition-state theory with multidimensional tunneling (EA-VTST/MT) combined with Grote-Hynes theory. The KIEs were calculated across the same temperature range examined experimentally, revealing a temperature independent behavior, in agreement with experimental findings. The calculations show that the H-transfer proceeds with ∼91% by tunneling in the case of protium and ∼80% when the transferred protium is replaced by tritium. Dynamic recrossing coefficients are almost invariant with temperature and in all cases far from unity, showing significant coupling between protein motions and the reaction coordinate. In particular, the relative movement of a conserved arginine (Arg166 in Escherichia coli ) promotes the departure of a conserved cysteine (Cys146 in E. coli ) from the dUMP by polarizing the thioether bond thus facilitating this bond breaking that takes place concomitantly with the hydride transfer. These promoting vibrations of the enzyme, which represent some of the dimensions of the real reaction coordinate, would limit the search through configurational space to efficiently find those decreasing both barrier height and width, thereby enhancing the probability of H-transfer by either tunneling (through barrier) or classical (over-the-barrier) mechanisms. In other words, the thermal fluctuations that are coupled to the reaction coordinate, together with transition-state geometries and tunneling, are the same in different bath temperatures (within the limited experimental range examined). All these terms contribute to the observed temperature independent KIEs in thymidylate synthase.

Thymidylate synthase catalyzed H-transfers: two chapters in one tale

Examination of the nature of different bond activations along the same catalytic path is of general interest in chemistry and biology. In this report, we compare the physical nature of two sequential H-transfers in the same enzymatic reaction. Thymidylate synthase (TSase) catalyzes a complex reaction that involves many chemical transformations including two different C-H bond cleavages, a rate-limiting C-H-C hydride transfer and a non-rate-limiting C-H-O proton transfer. Although the large kinetic complexity imposes difficulties in studying the proton transfer catalyzed by TSase, we are able to experimentally extract the intrinsic kinetic isotope effects (KIEs) on both steps. In contrast with the hydride transfer, the intrinsic KIEs of the proton transfer are temperature dependent. The results are interpreted within the framework of the Marcus-like model. This interpretation suggests that TSase optimizes the donor-acceptor geometries for the slower and overall rate-limiting hydride transfer but not for the faster proton transfer.

RNA-assisted catalysis in a protein enzyme: The 2'-hydroxyl of tRNA(Thr) A76 promotes aminoacylation by threonyl-tRNA synthetase

Aminoacyl-tRNA synthetases (aaRSs) join amino acids to 1 of 2 terminal hydroxyl groups of their cognate tRNAs, thereby contributing to the overall fidelity of protein synthesis. In class II histidyl-tRNA synthetase (HisRS) the nonbridging S(p)-oxygen of the adenylate is a potential general base for aminoacyl transfer. To test for conservation of this mechanism in other aaRSs and the role of terminal hydroxyls of tRNA in aminoacyl transfer, we investigated the class II Escherichia coli threonyl-tRNA synthetase (ThrRS). As with other class II aaRSs, the rate-determining step for ThrRS is amino acid activation. In ThrRS, however, the 2'-OH of A76 of tRNA(Thr) and a conserved active-site histidine (His-309) collaborate to catalyze aminoacyl transfer by a mechanism distinct from HisRS. Conserved residues in the ThrRS active site were replaced with alanine, and then the resulting mutant proteins were analyzed by steady-state and rapid kinetics. Nearly all mutants preferentially affected the amino acid activation step, with only a modest effect on aminoacyl transfer. By contrast, H309A ThrRS decreased transfer 242-fold and imposed a kinetic block to CCA accommodation. His-309 hydrogen bonds to the 2'-OH of A76, and substitution of the latter by hydrogen or fluorine decreased aminoacyl transfer by 763- and 94-fold, respectively. The proton relay mechanism suggested by these data to promote aminoacylation is reminiscent of the NAD(+)-dependent mechanisms of alcohol dehydrogenases and sirtuins and the RNA-mediated catalysis of the ribosomal peptidyl transferase center.

A molecular modeling study of the interaction of 2'-fluoro-substituted analogues of dUMP/FdUMP with thymidylate synthase

Molecular dynamics simulations and free energy calculations are presented, exploring previously described experimentally studied interactions of a series of 2'-fluoro-substituted dUMP/FdUMP analogues with thymidylate synthase (TS). The results show the inhibitory behaviors of 2'-F-ara-UMP, 2',2''-diF-dUMP and 2',5-diF-ara-UMP to be dependent upon the binding positions and orientations adopted by the molecules of these compounds in the active site of TS. The binding mode of 2',5-diF-ara-UMP suggests a novel role of the active site residue Trp 80, stabilizing through hydrophobic stacking the binding position of the pyrimidine ring in 2',5-diF-ara-UMP.

Role of an invariant lysine residue in folate binding on Escherichia coli thymidylate synthase: calorimetric and crystallographic analysis of the K48Q mutant

Thymidylate synthase (TS) catalyzes the reductive methylation of deoxyuridine monophosphate (dUMP) using methylene tetrahydrofolate (CH(2)THF) as cofactor, the glutamate tail of which forms a water-mediated hydrogen bond with an invariant lysine residue of this enzyme. To understand the role of this interaction, we studied the K48Q mutant of Escherichia coli TS using structural and biophysical methods. The k(cat) of the K48Q mutant was 430-fold lower than wild-type TS in activity, while the K(m) for the (R)-stereoisomer of CH(2)THF was 300 microM, about 30-fold larger than K(m) from the wild-type TS. Affinity constants were determined using isothermal titration calorimetry, which showed that binding was reduced by one order of magnitude for folate-like TS inhibitors, such as propargyl-dideazafolate (PDDF) or compounds that distort the TS active site like BW1843U89 (U89). The crystal structure of the K48Q-dUMP complex revealed that dUMP binding is not impaired in the mutant, and that U89 in a ternary complex of K48Q-nucleotide-U89 was bound in the active site with subtle differences relative to comparable wild-type complexes. PDDF failed to form ternary complexes with K48Q and dUMP. Thermodynamic data correlated with the structural determinations, since PDDF binding was dominated by enthalpic effects while U89 had an important entropic component. In conclusion, K48 is critical for catalysis since it leads to a productive CH(2)THF binding, while mutation at this residue does not affect much the binding of inhibitors that do not make contact with this group.

Role of Y94 in proton and hydride transfers catalyzed by thymidylate synthase

Thymidylate synthase (TS) catalyzes the substitution of a carbon-bound proton in a uracil base by a methyl group to yield thymine in the de novo biosynthesis of this DNA base. The enzymatic mechanism involves making and breaking several covalent bonds. Traditionally, a conserved tyrosine (Y94 in Escherichia coli, Y146 in Lactobacillus casei, and Y135 in humans) was assumed to serve as the general base catalyzing the proton abstraction. That assumption was examined here by comparing the nature of the proton abstraction using wild-type (wt) E. coli TS (ecTS) and its Y94F mutant (with a turnover rate reduced by 2 orders of magnitude). A subsequent hydride transfer was also studied using the wt and Y94F. The physical nature of both H-transfer steps was examined by determining intrinsic kinetic isotope effects (KIEs). Surprisingly, the findings did not suggest a direct role for Y94 in the proton abstraction step. The effect of this mutation on the subsequent hydride transfer was examined by a comparison of the temperature dependency of the intrinsic KIE on both the wt and the mutant. The intrinsic KIEs for Y94F at physiological temperatures were slightly smaller than those for wt but, otherwise, were as temperature-independent, suggesting a perfectly preorganized reaction coordinate for both enzymes. At reduced temperatures, however, the KIE for the mutant increased with a decrease in temperature, indicating a poorly preorganized reaction coordinate. Other kinetic and structural properties were also compared, and the findings suggested that Y94 is part of a H-bond network that plays a critical role at a step between the proton and the hydride transfers, presumably the dissociation of H4folate from the covalently bound intermediate. The possibility that no single residue serves as the general base in question but, rather, that the whole network of H-bonds at the active site catalyzes proton abstraction is discussed.