Asparagine 229 in thymidylate synthase contributes to, but is not essential for, catalysis

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.

Flavin-dependent thymidylate synthase X limits chromosomal DNA replication

We have investigated the hitherto unexplored possibility that differences in the catalytic efficiencies of thymidylate synthases ThyX and ThyA, enzymes that produce the essential DNA precursor dTMP, have influenced prokaryotic genome evolution. We demonstrate that DNA replication speed in bacteria and archaea that contain the low-activity ThyX enzyme is up to 10-fold decreased compared with species that contain the catalytically more efficient ThyA. Our statistical studies of >400 genomes indicated that ThyA proteins are preferred for the replication of large genomes, providing further evidence that the thymidylate metabolism is limiting expansion of prokaryotic genomes. Because both ThyX and ThyA participate in frequent reciprocal gene replacement events, our observations indicate that the bacterial metabolism continues to modulate the size and composition of prokaryotic genomes. We also propose that the increased kinetic efficiency of thymidylate synthesis has contributed to extending the prokaryotic evolutionary potential.

AdoMet-dependent Methyl-transfer: Glu119 Is Essential for DNA C5-Cytosine Methyltransferase M.HhaI

The role of Glu119 in S-adenosyl-L-methionine-dependent DNA methyltransferase M.HhaI-catalyzed DNA methylation was studied. Glu119 belongs to the highly conserved Glu/Asn/Val motif found in all DNA C5-cytosine methyltransferases, and its importance for M.HhaI function remains untested. We show that formation of the covalent intermediate between Cys81 and the target cytosine requires Glu119, since conversion to Ala, Asp or Gln lowers the rate of methyl transfer 10(2)-10(6) fold. Further, unlike the wild-type M.HhaI, these mutants are not trapped by the substrate in which the target cytosine is replaced with the mechanism-based inhibitor 5-fluorocytosine. The DNA binding affinity for the Glu119Asp mutant is decreased 10(3)-fold. Thus, the ability of the enzyme to stabilize the extrahelical cytosine is coupled directly to tight DNA binding. The structures of the ternary protein/DNA/AdoHcy complexes for both the Glu119Ala and Glu119Gln mutants (2.70 A and 2.75 A, respectively) show that the flipped base is positioned nearly identically with that observed in the wild-type M.HhaI complex. A single water molecule in the Glu119Ala structure between Ala119 and the extrahelical cytosine N3 is lacking in the Glu119Gln and wild-type M.HhaI structures, and most likely accounts for this mutant's partial activity. Glu119 has essential roles in activating the target cytosine for nucleophilic attack and contributes to tight DNA binding.

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.

Experimental Verification of the Crucial Roles of Glu73 in the Catalytic Activity and Structural Stability of Goose Type Lysozyme

The roles of Glu(73), which has been proposed to be a catalytic residue of goose type (G-type) lysozyme based on X-ray structural studies, were investigated by means of its replacement with Gln, Asp, and Ala using ostrich egg-white lysozyme (OEL) as a model. No remarkable differences in secondary structure or substrate binding ability were observed between the wild type and Glu(73)-mutated proteins, as evaluated by circular dichroism (CD) spectroscopy and chitin-coated celite chromatography. Substitution of Glu(73) with Gln or Ala abolished the enzymatic activity toward both the bacterial cell substrate and N-acetylglucosamine pentamer, (GlcNAc)(5), while substitution with Asp did not abolish but drastically reduced the activity of OEL. These results demonstrate that the carboxyl group of Glu(73) is directly involved in the catalytic action of G-type lysozyme. Furthermore, the stabilities of all three mutants, which were determined from the thermal and guanidine hydrochloride (GdnHCl) unfolding curves, respectively, were significantly decreased relative to those of the wild type. The results obtained clearly indicate the crucially important roles of Glu(73) in the structural stability as well as in the catalytic activity of G-type lysozyme.

Relative free energies of binding to thymidylate synthase of 2- and/or 4-thio and/or 5-fluoro analogues of dUMP

Free energy perturbation calculations have been applied to evaluate the relative free energies of binding of 2'-deoxyuridine-5'-monophosphate (dUMP) and its 2- and/or 4-thio and/or 5-fluoro analogues to the wild-type E. coli thymidylate synthase (ecTS). The results accurately reproduce experimentally measured differences in the free energy of binding of dUMP versus 5-fluoro-dUMP to thymidylate synthase. They indicate that preferred binding of dUMP compared to 5-fluoro-dUMP in the binary complex is equally related to (i) more favorable electrostatic interactions of the dUMP molecule in the enzyme active site, and (ii) its less favorable solvation in the aqueous solution. The relative free energies of binding in the binary complex show moderate and qualitatively indistinguishable discrimination among the studied fluorinated and non-fluorinated 2- and/or 4-thio analogues of dUMP. The binding free energies of monothio analogues of dUMP and 5-fluoro-dUMP correspond quite well with experimentally measured activities of these nucleotides in the thymidylate synthase reaction. On the other hand, the binding free energies of both dithio analogues, 2,4-dithio-dUMP and 2,4-dithio-FdUMP, show lack of such correlation. The latter suggests that very low activities of the dithio analogues of dUMP and 5-fluoro-dUMP may relate more to the covalent reaction of these nucleotides within the ternary complex with TS and 5,10-methylenetetrahydrofolate, than to their pre-covalent binding. We speculate that a lack of substrate activity of 2,4-dithio-dUMP is related to the high aromaticity of its pyrimidine ring that prevents the Michael addition of the active site cysteine thiol to the pyrimidine C6 atom. A stronger affinity of the fluorinated analogues of dUMP to thymidylate synthase, compared to the non-fluorinated congeners, results from the fluorine substituent producing a local strain in the C6 region in the pyrimidine ring, thus sensitizing C6 to the Michael addition of the cysteine thiol.

Enzyme-substrate interactions in the purine-specific nucleoside hydrolase from Trypanosoma vivax

Nucleoside hydrolases are key enzymes in the purine salvage pathway of Trypanosomatidae and are considered as targets for drug design. We previously reported the first x-ray structure of an inosine-adenosine-guanosine preferring nucleoside hydrolase (IAG-NH) from Trypanosoma vivax (). Here we report the 2.0-A crystal structure of the slow D10A mutant in complex with the inhibitor 3-deaza-adenosine and the 1.6-A crystal structure of the same enzyme in complex with a genuine substrate inosine. The enzyme-substrate complex shows the substrate bound to the enzyme in a different conformation from 3-deaza-adenosine and provides a snapshot along the reaction coordinate of the enzyme-catalyzed reaction. The chemical groups on the substrate important for binding and catalysis are mapped. The 2'-OH, 3'-OH, and 5'-OH contribute 4.6, 7.5, and 5.4 kcal/mol to k(cat)/K(m), respectively. Specific interactions with the exocyclic groups on the purine ring are not required for catalysis. Site-directed mutagenesis indicates that the purine specificity of the IAG-NHs is imposed by a parallel aromatic stacking interaction involving Trp(83) and Trp(260). The pH profiles of k(cat) and k(cat)/K(m) indicate the existence of one or more proton donors, possibly involved in leaving group activation. However, mutagenesis of the active site residues around the nucleoside base and an alanine scan of a flexible loop near the active site fail to identify this general acid. The parallel aromatic stacking seems to provide the most likely alternative mechanism for leaving group activation.

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

Crystal structures of 5-fluoro-dUrd and its 2 and/or 4-thio analogues: models of substituted dUMP pyrimidine ring interacting with thymidylate synthase

In order to understand the influence on thymidylate synthase interactions with dUMP analogues of the pyrimidine ring 2- and/or 4-thio, and 5-fluoro substitutions, X-ray diffractions by crystals of 5-fluoro-dUrd and its 2- and 4-thio, and 2,4-dithio analogues were measured, the four structures solved and refined. The following conclusions were suggested by results of comparative analyses of structural parameters (bond lengths, valence angles), followed by theoretical considerations based on calculated resonance structure distributions and aromaticity indices of the uracil, thiouracil, fluorouracil and fluorothiouracil rings. The effect of 4-thio substitution of FdUMP, altering specificity of inactivation of thymidylate synthases from various sources, is probably due to weaker proton acceptor power of the 4-thio substituent and increasing acidity (enhanced proton-donor power) of the N(3)-H moiety, resulting in an impaired fitness into the network of hydrogen bonds in the enzyme active center cleft. 2,4-Dithio substitution results in (i) impaired pyrimidine ring recognition by the enzyme active center, due to the 4-thio substituent (ii) increased pyrimidine ring aromaticity in dUMP, leading to resistance of C(6) to nucleophilic attack by the enzyme active center cysteine and (iii) altered planarity of the pyrimidine ring and deflections, with respect to the ring plane, of substituents at C(2), C(4) and C(5). 5-Fluoro substitution apparently activates the pyrimidine ring towards the interaction with thymidylate synthase by producing local strain, which results in an increased reactivity as predicted by the Walsh-Bent rule.

Contributions of orientation and hydrogen bonding to catalysis in Asn229 mutants of thymidylate synthase

We have determined structures of binary and ternary complexes of five Asn229 variants of thymidylate synthase (TS) and related their structures to the kinetic constants measured previously. Asn229 forms two hydrogen bonds to the pyrimidine ring of the substrate 2'-deoxyuridine-5'-monophosphate (dUMP). These hydrogen bonds constrain the orientation of dUMP in binary complexes with dUMP, and in ternary complexes with dUMP and the TS cofactor, 5,10-methylene-5,6,7,8-tetrahydrofolate. In N229 mutants, where these hydrogen bonds cannot be made, dUMP binds in a misoriented or more disordered fashion. Most N229 mutants exhibit no activity for the dehalogenation of 5-bromo-dUMP, which requires correct orientation of dUMP against Cys198. Since bound dUMP forms the binding surface against which the pterin ring of cofactor binds, misorientation of dUMP results in higher Km values for cofactor. At the same time, binding of the cofactor aids in ordering and positioning dUMP for catalysis. Hydrophobic mutants, such as N229I, favor an arrangement of solvent molecules and side-chains around the ligands similar to that in a proposed transition state for ternary complex formation in wild-type TS, and kcat values are similar to the wild-type value. Smaller, more hydrophilic mutants favor arrangements of the solvent and side-chains surrounding the ligands that do not resemble the proposed transition state. These changes correspond to decreases in kcat of up to 2000-fold, with only modest increases in Km or Kd. These results are consistent with the proposal that the hydrogen-bonding network between water, dUMP and side-chains in the active-site cavity contributes to catalysis in TS. Asn229 has the unique ability to maintain this critical network, without sterically interfering with dUMP binding.

Modeling of reaction steps relevant to deoxyuridylate (dUMP) enzymatic methylation and thymidylate synthase mechanism-based inhibition

Theoretical quantum mechanical ab initio Hartree-Fock calculations on molecular systems, modeling processes related to the specificity of thymidylate synthase inactivation are reported. We considered several steps of the methylation of the substrate dUMP and 4- or 5-mono- and 4,5-bisubstituted dUMP analogs, as well. The following reactions were modeled: the cysteine residue (Cys198 in the L.casei enzyme) nucleophilic attack on the substrate and the substrate C(5)-H proton abstraction. The substrate was modeled by the 1-methyluracil molecule and its structural analogs. The cysteine Cys198 residue was modeled by the methylmercaptane molecule. The substrate-enzyme binary complex was modeled by the 1-methyl-5,6-dihydro-6-thiomethyl-uracil (P1) molecule. The present theoretical calculations suggest that the cysteine nucleophilic attack on the substrate may result in the SH-group addition to the pyrimidine C(5)=C(6) bond in the course of a weakly exothermic reaction. The formerly presumed enolate carbanion appeared to be weakly stable or unstable and it can readily split into the thiol and pyrimidine residues. The s2-thio- (P2) and s2,4-dithio- (P3) substrate analogs should form stable thiolate anions after cysteine residue attachment to the C(6) position of the pyrimidine ring. Studies of the deformed P1 molecule interacting with a water molecule bound to the pyrimidine C(4)=O carbonyl residue allow a suggestion that this water molecule may be directly involved in the C(5)-H proton abstraction and may serve as a proton transmitter between the substrate and the proton acceptor residue, possibly located on the cofactor N10-nitrogen. Interaction of the pyrimidine C(4)=O group, or its modification, with the N5,10-methylenetetrahydrofolate N(10) nitrogen atom is suggested as an additional factor influencing the inhibition process.

Thymidylate synthase inhibition: a structure-based rationale for drug design

Thymidylate synthase (TS) is a very interesting target in antiproliferative diseases. Its inhibition causes thimineless death of the cells and compounds inhibiting TS are widely used in anticancer therapy. The classical antifolate TS inhibitors are structural analogs of the folate cofactor; they often share the same metabolic pathways and this causes the development of resistance inside the cells. A detailed analysis of the available x-ray crystal structures of the complexes of the enzyme with different substrates and inhibitors support the finding of a structural basis of their biological activity. TS inhibitors nonstructural analog of folate, non-analog antifolate inhibitors (NAAI), are welcome as a new interesting research topic. Among the most recent and interesting ones, compounds from Agouron related to the indole structure, are independent on the folate metabolism, highly active and specific for human TS. Other compounds, phthalein derivatives, can inhibit TS enzymes from various sources and show an interesting biological activity profile: they inhibit better bacterial and fungal TS than human TS. The x-ray crystal structures of some of these inhibitors with TS show that they bind in a different binding site from that of the classical folate TS inhibitors. This indicates a potential allosteric binding site useful for future drug discovery studies.

Structure of pvu II DNA-(cytosine N4) methyltransferase, an example of domain permutation and protein fold assignment

We have determined the structure of Pvu II methyltransferase (M. Pvu II) complexed with S -adenosyl-L-methionine (AdoMet) by multiwavelength anomalous diffraction, using a crystal of the selenomethionine-substituted protein. M. Pvu II catalyzes transfer of the methyl group from AdoMet to the exocyclic amino (N4) nitrogen of the central cytosine in its recognition sequence 5'-CAGCTG-3'. The protein is dominated by an open alpha/beta-sheet structure with a prominent V-shaped cleft: AdoMet and catalytic amino acids are located at the bottom of this cleft. The size and the basic nature of the cleft are consistent with duplex DNA binding. The target (methylatable) cytosine, if flipped out of the double helical DNA as seen for DNA methyltransferases that generate 5-methylcytosine, would fit into the concave active site next to the AdoMet. This M. Pvu IIalpha/beta-sheet structure is very similar to those of M. Hha I (a cytosine C5 methyltransferase) and M. Taq I (an adenine N6 methyltransferase), consistent with a model predicting that DNA methyltransferases share a common structural fold while having the major functional regions permuted into three distinct linear orders. The main feature of the common fold is a seven-stranded beta-sheet (6 7 5 4 1 2 3) formed by five parallel beta-strands and an antiparallel beta-hairpin. The beta-sheet is flanked by six parallel alpha-helices, three on each side. The AdoMet binding site is located at the C-terminal ends of strands beta1 and beta2 and the active site is at the C-terminal ends of strands beta4 and beta5 and the N-terminal end of strand beta7. The AdoMet-protein interactions are almost identical among M. Pvu II, M. Hha I and M. Taq I, as well as in an RNA methyltransferase and at least one small molecule methyltransferase. The structural similarity among the active sites of M. Pvu II, M. Taq I and M. Hha I reveals that catalytic amino acids essential for cytosine N4 and adenine N6 methylation coincide spatially with those for cytosine C5 methylation, suggesting a mechanism for amino methylation.

A tetrad of ionizable amino acids is important for catalysis in barley beta-glucanases

Determination of the crystal structures of a 1,3-beta-D-glucanase (E.C. 3.2.1.39) and a 1,3-1,4-beta-D-glucanase (E.C. 3.2.1.73) from barley (Hordeum vulgare) (Varghese, J.N, Garrett, T. P. J., Colman, P. M., Chen, L., Høj, P. B., and Fincher, G. B. (1994) Proc. Natl. Acad. Sci. U.S.A. 91, 2785-2789) showed the spatial positions of the catalytic residues in the substrate-binding clefts of the enzymes and also identified highly conserved neighboring amino acid residues. Site-directed mutagenesis of the 1,3-beta-glucanase has now been used to investigate the importance of these residues. Substitution of glutamine for the catalytic nucleophile Glu231 (mutant E231Q) reduced the specific activity about 20,000-fold. In contrast, substitution of glutamine for the catalytic acid Glu288 (mutant E288Q) had less severe consequences, reducing kcat approximately 350-fold with little effect on Km. Substitution of two neighboring and strictly conserved active site-located residues Glu279 (mutant E279Q) and Lys282 (mutant K282M) led to 240- and 2500-fold reductions of Kcat, respectively, with small increases in Km. Thus, a tetrad of ionizable amino acids is required for efficient catalysis in barley beta-glucanases. The active site-directed inhibitor 2,3-epoxypropyl beta-laminaribioside was soaked into native crystals. Crystallographic refinement revealed all four residues (Glu231, Glu279, Lys282, and Glu288) to be in contact with the bound inhibitor, and the orientation of bound substrate in the active site of the glucanase was deduced.