An essential role for water in an enzyme reaction mechanism: the crystal structure of the thymidylate synthase mutant E58Q

Synthetic Organic "Aquachemistry" that Relies on Neither Cosolvents nor Surfactants

There is a growing awareness of the underlying power of catalytic reactions in water that is not limited to innate sustainability alone. Some Type III reactions are catalytically accelerated without dissolution of reactants and are occasionally highly selective, as shown by comparison with the corresponding reactions run in organic solvents or under solvent-free conditions. Such catalysts are highly diversified, including hydrophilic, lipophilic, and even solid catalysts. In this Outlook, we highlight the impressive characteristics of illustrative catalysis that is exerted despite the immiscibility of the substrates and reveal the intrinsic benefits of these enigmatic reactions for synthetic organic chemistry, albeit with many details remaining unclear. We hope that this brief introduction to the expanding field of synthetic organic "aquachemistry" will inspire organic chemists to use the platform to invent new transformations.

Role of water in cyclooxygenase catalysis and design of anti-inflammatory agents targeting two sites of the enzyme

While designing the anti-inflammatory agents targeting cyclooxygenase-2 (COX-2), we first identified a water loop around the heme playing critical role in the enzyme catalysis. The results of molecular dynamic studies supported by the strong hydrogen-bonding equilibria of the participating atoms, radical stabilization energies, the pK_a of the H-donor/acceptor sites and the cyclooxygenase activity of pertinent muCOX-2 ravelled the working of the water–peptide channel for coordinating the flow of H·/electron between the heme and Y385. Based on the working of H·/electron transfer channel between the 12.5 Å distant radical generation and the radical disposal sites, a series of molecules was designed and synthesized. Among this category of compounds, an appreciably potent anti-inflammatory agent exhibiting IC₅₀ 0.06 μM against COX-2 and reversing the formalin induced analgesia and carageenan induced inflammation in mice by 90% was identified. Further it was revealed that, justifying its bidentate design, the compound targets water loop (heme bound site) and the arachidonic acid binding pockets of COX-2.

A Single Mutation Traps a Half-Sites Reactive Enzyme in Midstream, Explaining Asymmetry in Hydride Transfer

In Escherichia coli thymidylate synthase (EcTS), rate-determining hydride transfer from the cofactor 5,10-methylene-5,6,7,8-tetrahydrofolate to the intermediate 5-methylene-2'-deoxyuridine 5'-monophosphate occurs by hydrogen tunneling, requiring precise alignment of reactants and a closed binding cavity, sealed by the C-terminal carboxyl group. Mutations that destabilize the closed conformation of the binding cavity allow small molecules such as β-mercaptoethanol (β-ME) to enter the active site and compete with hydride for addition to the 5-methylene group of the intermediate. The C-terminal deletion mutant of EcTS produced the β-ME adduct in proportions that varied dramatically with cofactor concentration, from 50% at low cofactor concentrations to 0% at saturating cofactor conditions, suggesting communication between active sites. We report the 2.4 Å X-ray structure of the C-terminal deletion mutant of E. coli TS in complex with a substrate and a cofactor analogue, CB3717. The structure is asymmetric, with reactants aligned in a manner consistent with hydride transfer in only one active site. In the second site, CB3717 has shifted to a site where the normal cofactor would be unlikely to form 5-methylene-2'-deoxyuridine 5'-monophosphate, consistent with no formation of the β-ME adduct. The structure shows how the binding of the cofactor at one site triggers hydride transfer and borrows needed stabilization from substrate binding at the second site. It indicates pathways through the dimer interface that contribute to allostery relevant to half-sites reactivity.

The general base in the thymidylate synthase catalyzed proton abstraction

The enzyme thymidylate synthase (TSase), an important chemotherapeutic drug target, catalyzes the formation of 2'-deoxythymidine-5'-monophosphate (dTMP), a precursor of one of the DNA building blocks. TSase catalyzes a multi-step mechanism that includes the abstraction of a proton from the C5 of the substrate 2'-deoxyuridine-5'-monophosphate (dUMP). Previous studies on ecTSase proposed that an active-site residue, Y94 serves the role of the general base abstracting this proton. However, since Y94 is neither very basic, nor connected to basic residues, nor located close enough to the pyrimidine proton to be abstracted, the actual identity of this base remains enigmatic. Based on crystal structures, an alternative hypothesis is that the nearest potential proton-acceptor of C5 of dUMP is a water molecule that is part of a hydrogen bond (H-bond) network comprised of several water molecules and several protein residues including H147, E58, N177, and Y94. Here, we examine the role of the residue Y94 in the proton abstraction step by removing its hydroxyl group (Y94F mutant). We investigated the effect of the mutation on the temperature dependence of intrinsic kinetic isotope effects (KIEs) and found that these KIEs are more temperature dependent than those of the wild-type enzyme (WT). These results suggest that the phenolic -OH of Y94 is a component of the transition state for the proton abstraction step. The findings further support the hypothesis that no single functional group is the general base, but a network of bases and hydroxyls (from water molecules and tyrosine) sharing H-bonds across the active site can serve the role of the general base to remove the pyrimidine proton.

Infrared Pump-Probe Study of Nanoconfined Water Structure in Reverse Micelle

The influence of nanoconfinement on water structure is studied with time- and frequency-resolved vibrational spectroscopy of hydrazoic acid (HN3) encapsulated in reverse micelle. The azido stretch mode of HN3 is found to be a promising infrared probe for studying the structure and local hydrogen-bond environment of confined and interfacial water in reverse micelle due to its narrow spectral bandwidth and large transition dipole moment. The results show a clear separation between the core and shell spectral components, making it advantageous over the previously studied infrared probes. The measured vibrational lifetimes appear to be substantially different for the interfacial and bulk-like environments but show no remarkable size dependency, which indicates that water structures around this IR probe are distinctively different in the core and shell regions. The influence of local hydrogen bond network in the first and higher solvation shells on the vibrational dynamics of HN3 is further discussed.

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.

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.

The structure of Enterococcus faecalis thymidylate synthase provides clues about folate bacterial metabolism

Drug resistance to therapeutic antibiotics poses a challenge to the identification of novel targets and drugs for the treatment of infectious diseases. Infections caused by Enterococcus faecalis are a major health problem. Thymidylate synthase (TS) from E. faecalis is a potential target for antibacterial therapy. The X-ray crystallographic structure of E. faecalis thymidylate synthase (EfTS), which was obtained as a native binary complex composed of EfTS and 5-formyltetrahydrofolate (5-FTHF), has been determined. The structure provides evidence that EfTS is a half-of-the-sites reactive enzyme, as 5-FTHF is bound to two of the four independent subunits present in the crystal asymmetric unit. 5-FTHF is a metabolite of the one-carbon transfer reaction catalysed by 5-formyltetrahydrofolate cyclo-ligase. Kinetic studies show that 5-FTHF is a weak inhibitor of EfTS, suggesting that the EfTS-5-FTHF complex may function as a source of folates and/or may regulate one-carbon metabolism. The structure represents the first example of endogenous 5-FTHF bound to a protein involved in folate metabolism.

Explaining an unusually fast parasitic enzyme: folate tail-binding residues dictate substrate positioning and catalysis in Cryptosporidium hominis thymidylate synthase

The essential enzyme TS-DHFR from Cryptosporidium hominis undergoes an unusually rapid rate of catalysis at the conserved TS domain, facilitated by two nonconserved residues, Ala287 and Ser290, in the folate tail-binding region. Mutation of these two residues to their conserved counterparts drastically affects multiple steps of the TS catalytic cycle. We have determined the crystal structures of all three mutants (A287F, S290G, and A287F/S290G) in complex with active site ligands dUMP and CB3717. The structural data show two effects of the mutations: an increased distance between the ligands in the active site and increased flexibility of the folate ligand in the partially open enzyme state that precedes conformational change to the active catalytic state. The latter effect is able to be rescued by the mutants containing the A287F mutation. In addition, the conserved water network of TS is altered in each of the mutants. The structural results point to a role of the folate tail-binding residues in closely positioning ChTS ligands and restricting ligand flexibility in the partially open state to allow for a rapid transition to the active closed state and enhanced rate of catalysis. These results provide an explanation on how folate tail-binding residues at one end of the active site affect long-range interactions throughout the TS active site and validate these residues as targets for species-specific drug design.

The role of protein dynamics in thymidylate synthase catalysis: variants of conserved 2'-deoxyuridine 5'-monophosphate (dUMP)-binding Tyr-261

The enzyme thymidylate synthase (TS) catalyzes the reductive methylation of 2'-deoxyuridine 5'-monophosphate (dUMP) to 2'-deoxythymidine 5'-monophosphate. Using kinetic and X-ray crystallography experiments, we have examined the role of the highly conserved Tyr-261 in the catalytic mechanism of TS. While Tyr-261 is distant from the site of methyl transfer, mutants at this position show a marked decrease in enzymatic activity. Given that Tyr-261 forms a hydrogen bond with the dUMP 3'-O, we hypothesized that this interaction would be important for substrate binding, orientation, and specificity. Our results, surprisingly, show that Tyr-261 contributes little to these features of the mechanism of TS. However, the residue is part of the structural core of closed ternary complexes of TS, and conservation of the size and shape of the Tyr side chain is essential for maintaining wild-type values of kcat/Km. Moderate increases in Km values for both the substrate and cofactor upon mutation of Tyr-261 arise mainly from destabilization of the active conformation of a loop containing a dUMP-binding arginine. Besides binding dUMP, this loop has a key role in stabilizing the closed conformation of the enzyme and in shielding the active site from the bulk solvent during catalysis. Changes to atomic vibrations in crystals of a ternary complex of Escherichia coli Tyr261Trp are associated with a greater than 2000-fold drop in kcat/Km. These results underline the important contribution of dynamics to catalysis in TS.

Exploring structurally conserved solvent sites in protein families

Protein-bound water molecules are important components of protein structure, and therefore, protein function and energetics. Although structural conservation of solvent has been studied in a few protein families, a lack of suitable computational tools has hindered more comprehensive analyses. Herein we present a semiautomated computational approach for identifying solvent sites that are conserved among proteins sharing a common three-dimensional structure. This method is tested on six protein families: (1) monodomain cytochrome c, (2) fatty-acid binding protein, (3) lactate/malate dehydrogenase, (4) parvalbumin, (5) phospholipase A2, and (6) serine protease. For each family, the method successfully identified previously known conserved solvent sites. Moreover, the method discovered 22 novel conserved solvent sites, some of which have higher degrees of conservation than the previously known sites. All six families studied had solvent sites with more than 90% conservation and these sites were invariably located in regions of the protein with very high sequence conservation. These results suggest that highly conserved solvent sites, by virtue of their proximity to conserved residues, should be considered as one of the defining three-dimensional structural characteristics of protein families and folds.

The structure and thermal motion of the B800-850 LH2 complex from Rps.acidophila at 2.0A resolution and 100K: new structural features and functionally relevant motions

The structure at 100K of integral membrane light-harvesting complex II (LH2) from Rhodopseudomonas acidophila strain 10050 has been refined to 2.0A resolution. The electron density has been significantly improved, compared to the 2.5A resolution map, by high resolution data, cryo-cooling and translation, libration, screw (TLS) refinement. The electron density reveals a second carotenoid molecule, the last five C-terminal residues of the alpha-chain and a carboxy modified alpha-Met1 which forms the ligand of the B800 bacteriochlorophyll. TLS refinement has enabled the characterisation of displacements between molecules in the complex. B850 bacteriochlorophyll molecules are arranged in a ring of 18 pigments composed of nine approximate dimers. These pigments are strongly coupled and at their equilibrium positions the excited state dipole interaction energies, within and between dimers, are approximately 370cm(-1) and 280cm(-1), respectively. This difference in coupling energy is similar in magnitude to changes in interaction energies arising from the pigment displacements described by TLS tensors. The displacements appear to be non-random in nature and appear to be designed to optimise the modulation of pigment energy interactions. This is the first time that LH2 pigment displacements have been quantified experimentally. The calculated energy changes indicate that there may be significant contributions to inter-pigment energy interactions from molecular displacements and these may be of importance to photosynthetic energy transfer.

Multi-targeted antifolates aimed at avoiding drug resistance form covalent closed inhibitory complexes with human and Escherichia coli thymidylate synthases

Crystal structures of four pyrrolo(2,3-d)pyrimidine-based antifolate compounds, developed as inhibitors of thymidylate synthase (TS) in a strategy to circumvent drug-resistance, have been determined in complexes with their in vivo target, human thymidylate synthase, and with the structurally best-characterized Escherichia coli enzyme, to resolutions of 2.2-3.0 A. The 2.9 A crystal structure of a complex of human TS with one of the inhibitors, the multi-targeted antifolate LY231514, demonstrates that this compound induces a "closed" enzyme conformation and leads to formation of a covalent bond between enzyme and substrate. This structure is one of the first liganded human TS structures, and its solution was aided by mutation to facilitate crystallization. Structures of three other pyrrolo(2,3-d)pyrimidine-based antifolates in complex with Escherichia coli TS confirm the orientation of this class of inhibitors in the active site. Specific interactions between the polyglutamyl moiety and a positively charged groove on the enzyme surface explain the marked increase in affinity of the pyrrolo(2,3-d)pyrimidine inhibitors once they are polyglutamylated, as mediated in vivo by the cellular enzyme folyl polyglutamate synthetase.

Crystal structure of a deletion mutant of human thymidylate synthase Delta (7-29) and its ternary complex with Tomudex and dUMP

The crystal structures of a deletion mutant of human thymidylate synthase (TS) and its ternary complex with dUMP and Tomudex have been determined at 2.0 A and 2.5 A resolution, respectively. The mutant TS, which lacks 23 residues near the amino terminus, is as active as the wild-type enzyme. The ternary complex is observed in the open conformation, similar to that of the free enzyme and to that of the ternary complex of rat TS with the same ligands. This is in contrast to Escherichia coli TS, where the ternary complex with Tomudex and dUMP is observed in the closed conformation. While the ligands interact with each other in identical fashion regardless of the enzyme conformation, they are displaced by about 1.0 A away from the catalytic cysteine in the open conformation. As a result, the covalent bond between the catalytic cysteine sulfhydryl and the base of dUMP, which is the first step in the reaction mechanism of TS and is observed in all ternary complexes of the E. coli enzyme, is not formed. This displacement results from differences in the interactions between Tomudex and the protein that are caused by differences in the environment of the glutamyl tail of the Tomudex molecule. Despite the absence of the closed conformation, Tomudex inhibits human TS ten-fold more strongly than E. coli TS. These results suggest that formation of a covalent bond between the catalytic cysteine and the substrate dUMP is not required for effective inhibition of human TS by cofactor analogs and could have implications for drug design by eliminating this as a condition for lead compounds.

Crystal structure of deoxycytidylate hydroxymethylase from bacteriophage T4, a component of the deoxyribonucleoside triphosphate-synthesizing complex

Bacteriophage T4 deoxycytidylate hydroxymethylase (EC 2.1.2.8), a homodimer of 246-residue subunits, catalyzes hydroxymethylation of the cytosine base in deoxycytidylate (dCMP) to produce 5-hydroxymethyl-dCMP. It forms part of a phage DNA protection system and appears to function in vivo as a component of a multienzyme complex called deoxyribonucleoside triphosphate (dNTP) synthetase. We have determined its crystal structure in the presence of the substrate dCMP at 1.6 A resolution. The structure reveals a subunit fold and a dimerization pattern in common with thymidylate synthases, despite low (approximately 20%) sequence identity. Among the residues that form the dCMP binding site, those interacting with the sugar and phosphate are arranged in a configuration similar to the deoxyuridylate binding site of thymidylate synthases. However, the residues interacting directly or indirectly with the cytosine base show a more divergent structure and the presumed folate cofactor binding site is more open. Our structure reveals a water molecule properly positioned near C-6 of cytosine to add to the C-7 methylene intermediate during the last step of hydroxymethylation. On the basis of sequence comparison and crystal packing analysis, a hypothetical model for the interaction between T4 deoxycytidylate hydroxymethylase and T4 thymidylate synthase in the dNTP-synthesizing complex has been built.

Crystal structure of thymidylate synthase A from Bacillus subtilis

Thymidylate synthase (TS) converts dUMP to dTMP by reductive methylation, where 5,10-methylenetetrahydrofolate is the source of both the methylene group and reducing equivalents. The mechanism of this reaction has been extensively studied, mainly using the enzyme from Escherichia coli. Bacillus subtilis contains two genes for TSs, ThyA and ThyB. The ThyB enzyme is very similar to other bacterial TSs, but the ThyA enzyme is quite different, both in sequence and activity. In ThyA TS, the active site histidine is replaced by valine. In addition, the B. subtilis enzyme has a 2.4-fold greater k(cat) than the E. coli enzyme. The structure of B. subtilis thymidylate synthase in a ternary complex with 5-fluoro-dUMP and 5,10-methylenetetrahydrofolate has been determined to 2.5 A resolution. Overall, the structure of B. subtilis TS (ThyA) is similar to that of the E. coli enzyme. However, there are significant differences in the structures of two loops, the dimer interface and the details of the active site. The effects of the replacement of histidine by valine and a serine to glutamine substitution in the active site area, and the addition of a loop over the carboxy terminus may account for the differences in k(cat) found between the two enzymes.

Inactivity of N229A thymidylate synthase due to water-mediated effects: isolating a late stage in methyl transfer

Mutation of thymidylate synthase N229(177) to alanine results in an essentially inactive enzyme, yet it leads to formation of a stable ternary complex. The kinetics of N229(177)A show that kcat for Escherichia coli is reduced by 200-fold while the Km for dUMP is increased 200-fold and the Km for folate increased by tenfold versus the wild-type enzyme. The crystal structures of N229(177)A in complex with dUMP and CB3717, and in complex with dUMP alone are determined at 2.4 A, and 2.5 A resolution. These structures identify the covalently bound ternary complex and show how N229(177)A traps an intermediate, and so becomes inactive in a later step of the reaction. Since the smaller alanine side-chain at N229(177)A does not directly sterically impair binding of ligands, the structures implicate, and place quantitative limits on the involvement of the structured water network in the active site of thymidylate synthase in both catalysis and in determining the binding affinity for dUMP (in contrast, the N229(177)V mutation in Lactobacillus casei has minimal effect on activity).

The additivity of substrate fragments in enzyme-ligand binding

BACKGROUND

Enzymes have evolved to recognise their target substrates with exquisite selectivity and specificity. Whether fragments of the substrate--perhaps never available to the evolving enzyme--are bound in the same manner as the parent substrate addresses the fundamental basis of specificity. An understanding of the relative contributions of individual portions of ligand molecules to the enzyme-binding interaction may offer considerable insight into the principles of substrate recognition.

RESULTS

We report 12 crystal structures of Escherichia coli thymidylate synthase in complexes with available fragments of the substrate (dUMP), both with and without the presence of a cofactor analogue. The structures display considerable fidelity of binding mode and interactions. These complexes reveal several interesting features: the cofactor analogue enhances the localisation of substrate and substrate fragments near the reactive thiol; the ribose moiety reduces local disorder through additional specific enzyme-ligand interactions; the pyrimidine has multiple roles, ranging from stereospecificity to mechanistic competence; and the glycosidic linkage has an important role in the formation of a covalent attachment between substrate and enzyme.

CONCLUSIONS

The requirements of ligand-protein binding can be understood in terms of the binding of separate fragments of the ligand. Fragments which are subsystems of the natural substrate for the enzyme confer specific contributions to the binding affinity, orientation or electrostatics of the enzymatic mechanism. This ligand-binding analysis provides a complementary method to the more prevalent approaches utilising site-directed mutagenesis. In addition, these observations suggest a modular approach for rational drug design utilising chemical fragments.

Crystal structures of a marginally active thymidylate synthase mutant, Arg 126-->Glu

Thymidylate synthase (TS) is a long-standing target for anticancer drugs and is of interest for its rich mechanistic features. The enzyme catalyzes the conversion of dUMP to dTMP using the co-enzyme methylenetetrahydrofolate, and is perhaps the best studied of enzymes that catalyze carbon-carbon bond formation. Arg 126 is found in all TSs but forms only 1 of 13 hydrogen bonds to dUMP during catalysis, and just one of seven to the phosphate group alone. Despite this, when Arg 126 of TS from Escherichia coli was changed to glutamate (R126E), the resulting protein had kcat reduced 2000-fold and Km reduced 600-fold. The crystal structure of R126E was determined under two conditions--in the absence of bound ligand (2.4 A resolution), and with dUMP and the antifolate CB3717 (2.2 A resolution). The first crystals, which did not contain dUMP despite its presence in the crystallization drop, displayed Glu 126 in a position to sterically and electrostatically interfere with binding of the dUMP phosphate. The second crystals contained both dUMP and CB3717 in the active site, but Glu 126 formed three hydrogen bonds to nearby residues (two through water) and was in a position that partially overlapped with the normal phosphate binding site, resulting in a approximately 1 A shift in the phosphate group. Interestingly, the protein displayed the typical ligand-induced conformational change, and the covalent bond to Cys 146 was present in one of the protein's two active sites.

Thymidylate synthase: structure, inhibition, and strained conformations during catalysis

Thymidylate synthase (TS) is a long-standing target for chemotherapeutic agents because of its central role in DNA synthesis, and it is also of interest because of its rich mechanistic features. The reaction catalyzed by TS is the methylation of dUMP, with the transferred methyl group provided by the cofactor methylenetetrahydrofolate (CH2THF). Recently, several crystal structure determinations and mechanistic studies have led to a deeper understanding of the TS reaction mechanism, and address the role of conformational change in TS catalysis and inhibition. Included among these structures are complexes of TS bound to substrate dUMP; cofactor CH2THF; the nucleotide analogs 5-fluoro-dUMP, 5-nitro-dUMP and dGMP; and the promising antifolates BW1843, ZD1694, and AG337. From these studies, a picture of TS emerges where ligand-induced conformational changes play key roles in catalysis by straining the thiol adduct that occurs during the reaction; by protecting the highly reactive reaction intermediates; and by providing a means to stabilize a high-energy conformer of the cofactor after initial binding of a low-energy conformer. The best inhibitors of TS also induce and stabilize a conformational change in TS. One inhibitor, BW1843, distorts the active site on binding, and intercalates into a hydrophobic patch between two mobile subdomains in the protein. Also discussed are recent developments in the cell biology and regulation of eukaryotic TS and the use of structure-based drug design in the development of the antifolates currently in clinical trial for the treatment of cancer.