Structures of thymidylate synthase with a C-terminal deletion: role of the C-terminus in alignment of 2'-deoxyuridine 5'-monophosphate and 5,10-methylenetetrahydrofolate

Mutational and structural studies of (βα)8-barrel fold methylene-tetrahydropterin reductases utilizing a common catalytic mechanism

Methylene-tetrahydropterin reductases catalyze the reduction of a methylene to a methyl group bound to a reduced pterin as C1 carrier in various one-carbon (C1) metabolisms. F420-dependent methylene-tetrahydromethanopterin (methylene-H4MPT) reductase (Mer) and the flavin-independent methylene-tetrahydrofolate (methylene-H4F) reductase (Mfr) use a ternary complex mechanism for the direct transfer of a hydride from F420H2 and NAD(P)H to the respective methylene group, whereas FAD-dependent methylene-H4F reductase (MTHFR) uses FAD as prosthetic group and a ping-pong mechanism to catalyze the reduction of methylene-H4F. A ternary complex structure and a thereof derived catalytic mechanism of MTHFR is available, while no ternary complex structures of Mfr or Mer are reported. Here, Mer from Methanocaldococcus jannaschii (jMer) was heterologously produced and the crystal structures of the enzyme with and without F420 were determined. A ternary complex of jMer was modeled on the basis of the jMer-F420 structure and the ternary complex structure of MTHFR by superimposing the polypeptide after fixing hydride-transferring atoms of the flavins on each other, and by the subsequent transfer of the methyl-tetrahydropterin from MTHFR to jMer. Mutational analysis of four functional amino acids, which are similarly positioned in the three reductase structures, indicated despite the insignificant sequence identity, a common catalytic mechanism with a 5-iminium cation of methylene-tetrahydropterin as intermediate protonated by a shared glutamate. According to structural, mutational and phylogenetic analysis, the evolution of the three reductases most likely proceeds via a convergent development although a divergent scenario requiring drastic structural changes of the common ancestor cannot be completely ruled out.

An enzymatic activation of formaldehyde for nucleotide methylation

Folate enzyme cofactors and their derivatives have the unique ability to provide a single carbon unit at different oxidation levels for the de novo synthesis of amino-acids, purines, or thymidylate, an essential DNA nucleotide. How these cofactors mediate methylene transfer is not fully settled yet, particularly with regard to how the methylene is transferred to the methylene acceptor. Here, we uncovered that the bacterial thymidylate synthase ThyX, which relies on both folate and flavin for activity, can also use a formaldehyde-shunt to directly synthesize thymidylate. Combining biochemical, spectroscopic and anaerobic crystallographic analyses, we showed that formaldehyde reacts with the reduced flavin coenzyme to form a carbinolamine intermediate used by ThyX for dUMP methylation. The crystallographic structure of this intermediate reveals how ThyX activates formaldehyde and uses it, with the assistance of active site residues, to methylate dUMP. Our results reveal that carbinolamine species promote methylene transfer and suggest that the use of a CH₂O-shunt may be relevant in several other important folate-dependent reactions.

A cytosine modification mechanism revealed by the structure of a ternary complex of deoxycytidylate hydroxymethylase from bacteriophage T4 with its cofactor and substrate

To protect viral DNA against the host bacterial restriction system, bacterio-phages utilize a special modification system - hydroxymethylation - in which dCMP hydroxymethylase (dCH) converts dCMP to 5-hydroxymethyl-dCMP (5hm-dCMP) using N5,N10-methylenetetrahydrofolate as a cofactor. Despite shared similarity with thymidylate synthase (TS), dCH catalyzes hydroxylation through an exocyclic methylene intermediate during the last step, which is different from the hydride transfer that occurs with TS. In contrast to the extensively studied TS, the hydroxymethylation mechanism of a cytosine base is not well understood due to the lack of a ternary complex structure of dCH in the presence of both its substrate and cofactor. This paper reports the crystal structure of the ternary complex of dCH from bacteriophage T4 (T4dCH) with dCMP and tetrahydrofolate at 1.9 Å resolution. The authors found key residues of T4dCH for accommodating the cofactor without a C-terminal tail, an optimized network of ordered water molecules and a hydrophobic gating mechanism for cofactor regulation. In combination with biochemical data on structure-based mutants, key residues within T4dCH and a substrate water molecule for hydroxymethylation were identified. Based on these results, a complete enzyme mechanism of dCH and signature residues that can identify dCH enzymes within the TS family have been proposed. These findings provide a fundamental basis for understanding the pyrimidine modification system.

In Vivo Titration of Folate Pathway Enzymes

How enzymes behave in cells is likely different from how they behave in the test tube. Previous in vitro studies find that osmolytes interact weakly with folate. Removal of the osmolyte from the solvation shell of folate is more difficult than removal of water, which weakens binding of folate to its enzyme partners. To examine if this phenomenon occurs in vivo, osmotic stress titrations were performed with Escherichia coli Two strategies were employed: resistance to an antibacterial drug and complementation of a knockout strain by the appropriate gene cloned into a plasmid that allows tight control of expression levels as well as labeling by a degradation tag. The abilities of the knockout and complemented strains to grow under osmotic stress were compared. Typically, the knockout strain could grow to high osmolalities on supplemented medium, while the complemented strain stopped growing at lower osmolalities on minimal medium. This pattern was observed for an R67 dihydrofolate reductase clone rescuing a ΔfolA strain, for a methylenetetrahydrofolate reductase clone rescuing a ΔmetF strain, and for a serine hydroxymethyltransferase clone rescuing a ΔglyA strain. Additionally, an R67 dihydrofolate reductase clone allowed E. coli DH5α to grow in the presence of trimethoprim until an osmolality of ∼0.81 is reached, while cells in a control titration lacking antibiotic could grow to 1.90 osmol.IMPORTANCEE. coli can survive in drought and flooding conditions and can tolerate large changes in osmolality. However, the cell processes that limit bacterial growth under high osmotic stress conditions are not known. In this study, the dose of four different enzymes in E. coli was decreased by using deletion strains complemented by the gene carried in a tunable plasmid. Under conditions of limiting enzyme concentration (lower than that achieved by chromosomal gene expression), cell growth can be blocked by osmotic stress conditions that are normally tolerated. These observations indicate that E. coli has evolved to deal with variations in its osmotic environment and that normal protein levels are sufficient to buffer the cell from environmental changes. Additional factors involved in the osmotic pressure response may include altered protein concentration/activity levels, weak solute interactions with ligands which can make it more difficult for proteins to bind their substrates/inhibitors/cofactors in vivo, and/or viscosity effects.

A role for glutamine 183 in the folate oxidative half-reaction of methylenetetrahydrofolate reductase from Escherichia coli

The flavoprotein methylenetetrahydrofolate reductase (MTHFR) from Escherichia coli catalyzes a ping-pong reaction with NADH and 5,10-methylenetetrahydrofolate (CH₂-H₄folate) to produce NAD⁺ and 5-methyltetrahydrofolate (CH₃-H₄folate). This work focuses on the function of the invariant, active-site aminoacyl residue Gln183. X-ray structures of the enzyme complexes E_red(wild-type)•NADH and E_ox(Glu28Gln)•CH₃-H₄folate indicate that Gln183 makes key hydrogen-bonding interactions with both NADH and folate in their respective half-reactions, suggesting roles in binding each substrate. We propose that the polarity of Gln183 may also aid in stabilizing the proposed 5-iminium cation intermediate during catalysis in the oxidative half-reaction with folate. We have prepared mutants Gln183Ala and Gln183Glu, which we hypothesize to have altered charge/polarity and hydrogen bonding properties. We have examined the enzymes by steady-state and stopped-flow kinetics and by measurement of the flavin redox potentials. In the reductive half-reaction, NADH binding affinity and the rate of flavin reduction have not been hindered by either mutation. By contrast, our results support a minor role for Gln183 in the oxidative half-reaction. The Gln183Ala variant exhibited a 6-10 fold lower rate of folate reduction and bound CH₂-H₄folate with 7-fold lower affinity, whereas the Gln183Glu mutant displayed catalytic constants within 3-fold of the wild-type enzyme.

Classification of ligand molecules in PDB with graph match-based structural superposition

The fast heuristic graph match algorithm for small molecules, COMPLIG, was improved by adding a structural superposition process to verify the atom-atom matching. The modified method was used to classify the small molecule ligands in the Protein Data Bank (PDB) by their three-dimensional structures, and 16,660 types of ligands in the PDB were classified into 7561 clusters. In contrast, a classification by a previous method (without structure superposition) generated 3371 clusters from the same ligand set. The characteristic feature in the current classification system is the increased number of singleton clusters, which contained only one ligand molecule in a cluster. Inspections of the singletons in the current classification system but not in the previous one implied that the major factors for the isolation were differences in chirality, cyclic conformations, separation of substructures, and bond length. Comparisons between current and previous classification systems revealed that the superposition-based classification was effective in clustering functionally related ligands, such as drugs targeted to specific biological processes, owing to the strictness of the atom-atom matching.

Enzyme Tunnels and Gates As Relevant Targets in Drug Design

Many enzymes contain tunnels and gates that are essential to their function. Gates reversibly switch between open and closed conformations and thereby control the traffic of small molecules-substrates, products, ions, and solvent molecules-into and out of the enzyme's structure via molecular tunnels. Many transient tunnels and gates undoubtedly remain to be identified, and their functional roles and utility as potential drug targets have received comparatively little attention. Here, we describe a set of general concepts relating to the structural properties, function, and classification of these interesting structural features. In addition, we highlight the potential of enzyme tunnels and gates as targets for the binding of small molecules. The different types of binding that are possible and the potential pharmacological benefits of such targeting are discussed. Twelve examples of ligands bound to the tunnels and/or gates of clinically relevant enzymes are used to illustrate the different binding modes and to explain some new strategies for drug design. Such strategies could potentially help to overcome some of the problems facing medicinal chemists and lead to the discovery of more effective drugs.

Structural Analysis of Thymidylate Synthase from Kaposi's Sarcoma-Associated Herpesvirus with the Anticancer Drug Raltitrexed

Kaposi's sarcoma-associated herpesvirus (KSHV) is a highly infectious human herpesvirus that causes Kaposi's sarcoma. KSHV encodes functional thymidylate synthase, which is a target for anticancer drugs such as raltitrexed or 5-fluorouracil. Thymidylate synthase catalyzes the conversion of 2'-deoxyuridine-5'-monophosphate (dUMP) to thymidine-5'-monophosphate (dTMP) using 5,10-methylenetetrahydrofolate (mTHF) as a co-substrate. The crystal structures of thymidylate synthase from KSHV (apo), complexes with dUMP (binary), and complexes with both dUMP and raltitrexed (ternary) were determined at 1.7 Å, 2.0 Å, and 2.4 Å, respectively. While the ternary complex structures of human thymidylate synthase and E. coli thymidylate synthase had a closed conformation, the ternary complex structure of KSHV thymidylate synthase was observed in an open conformation, similar to that of rat thymidylate synthase. The complex structures of KSHV thymidylate synthase did not have a covalent bond between the sulfhydryl group of Cys219 and C6 atom of dUMP, unlike the human thymidylate synthase. The catalytic Cys residue demonstrated a dual conformation in the apo structure, and its sulfhydryl group was oriented toward the C6 atom of dUMP with no covalent bond upon ligand binding in the complex structures. These structural data provide the potential use of antifolates such as raltitrexed as a viral induced anticancer drug and structural basis to design drugs for targeting the thymidylate synthase of KSHV.

Structure of the Varicella Zoster Virus Thymidylate Synthase Establishes Functional and Structural Similarities as the Human Enzyme and Potentiates Itself as a Target of Brivudine

Varicella zoster virus (VZV) is a highly infectious human herpesvirus that is the causative agent for chicken pox and shingles. VZV encodes a functional thymidylate synthase (TS), which is the sole enzyme that produces dTMP from dUMP de novo. To study substrate binding, the complex structure of TS_VZV with dUMP was determined to a resolution of 2.9 Å. In the absence of a folate co-substrate, dUMP binds in the conserved TS active site and is coordinated similarly as in the human encoded TS (TS_HS) in an open conformation. The interactions between TS_VZV with dUMP and a cofactor analog, raltitrexed, were also studied using differential scanning fluorimetry (DSF), suggesting that TS_VZV binds dUMP and raltitrexed in a sequential binding mode like other TS. The DSF also revealed interactions between TS_VZV and in vitro phosphorylated brivudine (BVDU_P), a highly potent anti-herpesvirus drug against VZV infections. The binding of BVDU_P to TS_VZV was further confirmed by the complex structure of TS_VZV and BVDU_P solved at a resolution of 2.9 Å. BVDU_P binds similarly as dUMP in the TS_HS but it induces a closed conformation of the active site. The structure supports that the 5-bromovinyl substituent on BVDU_P is likely to inhibit TS_VZV by preventing the transfer of a methylene group from its cofactor and the subsequent formation of dTMP. The interactions between TS_VZV and BVDU_P are consistent with that TS_VZV is indeed a target of brivudine in vivo. The work also provided the structural basis for rational design of more specific TS_VZV inhibitors.

Conservation and Role of Electrostatics in Thymidylate Synthase

Conservation of function across families of orthologous enzymes is generally accompanied by conservation of their active site electrostatic potentials. To study the electrostatic conservation in the highly conserved essential enzyme, thymidylate synthase (TS), we conducted a systematic species-based comparison of the electrostatic potential in the vicinity of its active site. Whereas the electrostatics of the active site of TS are generally well conserved, the TSs from minimal organisms do not conform to the overall trend. Since the genomes of minimal organisms have a high thymidine content compared to other organisms, the observation of non-conserved electrostatics was surprising. Analysis of the symbiotic relationship between minimal organisms and their hosts, and the genetic completeness of the thymidine synthesis pathway suggested that TS from the minimal organism Wigglesworthia glossinidia (W.g.b.) must be active. Four residues in the vicinity of the active site of Escherichia coli TS were mutated individually and simultaneously to mimic the electrostatics of W.g.b TS. The measured activities of the E. coli TS mutants imply that conservation of electrostatics in the region of the active site is important for the activity of TS, and suggest that the W.g.b. TS has the minimal activity necessary to support replication of its reduced genome.

Impact of Mutating the Key Residues of a Bifunctional 5,10-Methylenetetrahydrofolate Dehydrogenase-Cyclohydrolase from Escherichia coli on Its Activities

Methylenetetrahydrofolate dehydrogenase-cyclohydrolase (FolD) catalyzes interconversion of 5,10-methylene-tetrahydrofolate and 10-formyl-tetrahydrofolate in the one-carbon metabolic pathway. In some organisms, the essential requirement of 10-formyl-tetrahydrofolate may also be fulfilled by formyltetrahydrofolate synthetase (Fhs). Recently, we developed an Escherichia coli strain in which the folD gene was deleted in the presence of Clostridium perfringens fhs (E. coli ΔfolD/p-fhs) and used it to purify FolD mutants (free from the host-encoded FolD) and determine their biological activities. Mutations in the key residues of E. coli FolD, as identified from three-dimensional structures (D121A, Q98K, K54S, Y50S, and R191E), and a genetic screen (G122D and C58Y) were generated, and the mutant proteins were purified to determine their kinetic constants. Except for the R191E and K54S mutants, others were highly compromised in terms of both dehydrogenase and cyclohydrolase activities. While the R191E mutant showed high cyclohydrolase activity, it retained only a residual dehydrogenase activity. On the other hand, the K54S mutant lacked the cyclohydrolase activity but possessed high dehydrogenase activity. The D121A and G122D (in a loop between two helices) mutants were highly compromised in terms of both dehydrogenase and cyclohydrolase activities. In vivo and in vitro characterization of wild-type and mutant (R191E, G122D, D121A, Q98K, C58Y, K54S, and Y50S) FolD together with three-dimensional modeling has allowed us to develop a better understanding of the mechanism for substrate binding and catalysis by E. coli FolD.

Activation of a unique flavin-dependent tRNA-methylating agent

TrmFO is a tRNA methyltransferase that uses methylenetetrahydrofolate (CH2THF) and flavin adenine dinucleotide hydroquinone as cofactors. We have recently shown that TrmFO from Bacillus subtilis stabilizes a TrmFO-CH2-FADH adduct and an ill-defined neutral flavin radical. The adduct contains a unique N-CH2-S moiety, with a methylene group bridging N5 of the isoalloxazine ring and the sulfur of an active-site cysteine (Cys53). In the absence of tRNA substrate, this species is remarkably stable but becomes catalytically competent for tRNA methylation following tRNA addition using the methylene group as the source of methyl. Here, we demonstrate that this dormant methylating agent can be activated at low pH, and we propose that this process is triggered upon tRNA addition. The reaction proceeds via protonation of Cys53, cleavage of the C-S bond, and generation of a highly reactive [FADH(N5)═CH2]+ iminium intermediate, which is proposed to be the actual tRNA-methylating agent. This mechanism is fully supported by DFT calculations. The radical present in TrmFO is characterized here by optical and EPR/ENDOR spectroscopy approaches together with DFT calculations and is shown to be the one-electron oxidized product of the TrmFO-CH2-FADH adduct. It is also relatively stable, and its decomposition is facilitated by high pH. These results provide new insights into the structure and reactivity of the unique flavin-dependent methylating agent used by this class of enzymes.

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.

A catalytic intermediate and several flavin redox states stabilized by folate-dependent tRNA methyltransferase from Bacillus subtilis

The flavoprotein TrmFO catalyzes the C5 methylation of uridine 54 in the TΨC loop of tRNAs using 5,10-methylenetetrahydrofolate (CH(2)THF) as a methylene donor and FAD as a reducing agent. Here, we report biochemical and spectroscopic studies that unravel the remarkable capability of Bacillus subtilis TrmFO to stabilize, in the presence of oxygen, several flavin-reduced forms, including an FADH(•) radical, and a catalytic intermediate endowed with methylating activity. The FADH(•) radical was characterized by high-field electron paramagnetic resonance and electron nuclear double-resonance spectroscopies. Interestingly, the enzyme exhibited tRNA methylation activity in the absence of both an added carbon donor and an external reducing agent, indicating that a reaction intermediate, containing presumably CH(2)THF and FAD hydroquinone, is present in the freshly purified enzyme. Isolation by acid treatment, under anaerobic conditions, of noncovalently bound molecules, followed by mass spectrometry analysis, confirmed the presence in TrmFO of nonmodified FAD. Addition of formaldehyde to the purified enzyme protects the reduced flavins from decay by probably preventing degradation of CH(2)THF. The absence of air-stable reduced FAD species during anaerobic titration of oxidized TrmFO, performed in the absence or presence of added CH(2)THF, argues against their thermodynamic stabilization but rather implicates their kinetic trapping by the enzyme. Altogether, the unexpected isolation of a stable catalytic intermediate suggests that the flavin-binding pocket of TrmFO is a highly insulated environment, diverting the reduced FAD present in this intermediate from uncoupled reactions.

Mechanism of influence of phosphorylation on serine 124 on a decrease of catalytic activity of human thymidylate synthase

Regulation by phosphorylation is a well-established mechanism for controlling biological activity of proteins. Recently, phosphorylation of serine 124 in human thymidylate synthase (hTS) has been shown to lower the catalytic activity of the enzyme. To clarify a possible mechanism of the observed influence, molecular dynamics (MD), essential dynamics (ED) and MM-GBSA studies were undertaken. Structures derived from the MD trajectories reveal incorrect binding alignment between the pyrimidine ring of the substrate, dUMP, and the pterine ring of the cofactor analogue, THF, in the active site of the phosphorylated enzyme. The ED analysis indicates changes in the behavior of collective motions in the phosphorylated enzyme, suggesting that the formation of the closed ternary complex is hindered. Computed free energies, in agreement with structural analysis, predict that the binding of dUMP and THF to hTS is favored in the native compared to phosphorylated state of the enzyme. The paper describes at the structural level how phosphorylation at the distant site influences the ligand binding. We propose that the 'phosphorylation effect' is transmitted from the outside loop of Ser 124 into the active site via a subtle mechanism initiated by the long-range electrostatic repulsion between the phosphate groups of dUMP and Ser124. The mechanism can be described in terms of the interplay between the two groups of amino acids: the link (residues 125-134) and the patch (residues 189-192), resulting in the change of orientation of the pyrimidine ring of dUMP, which, in turn, prevents the correct alignment between the latter ring and the pterin ring of THF.

Functional role for the conformationally mobile phenylalanine 223 in the reaction of methylenetetrahydrofolate reductase from Escherichia coli

The flavoprotein methylenetetrahydrofolate reductase from Escherichia coli catalyzes the reduction of 5,10-methylenetetrahydrofolate (CH(2)-H(4)folate) by NADH via a ping-pong reaction mechanism. Structures of the reduced enzyme in complex with NADH and of the oxidized Glu28Gln enzyme in complex with CH(3)-H(4)folate [Pejchal, R., Sargeant, R., and Ludwig, M. L. (2005) Biochemistry 44, 11447-11457] have revealed Phe223 as a conformationally mobile active site residue. In the NADH complex, the NADH adopts an unusual hairpin conformation and is wedged between the isoalloxazine ring of the FAD and the side chain of Phe223. In the folate complex, Phe223 swings out from its position in the NADH complex to stack against the p-aminobenzoate ring of the folate. Although Phe223 contacts each substrate in E. coli MTHFR, this residue is not invariant; for example, a leucine occurs at this site in the human enzyme. To examine the role of Phe223 in substrate binding and catalysis, we have constructed mutants Phe223Ala and Phe223Leu. As predicted, our results indicate that Phe223 participates in the binding of both substrates. The Phe223Ala mutation impairs NADH and CH(2)-H(4)folate binding each 40-fold yet slows catalysis of both half-reactions less than 2-fold. Affinity for CH(2)-H(4)folate is unaffected by the Phe223Leu mutation, and the variant catalyzes the oxidative half-reaction 3-fold faster than the wild-type enzyme. Structures of ligand-free Phe223Leu and Phe223Leu/Glu28Gln MTHFR in complex with CH(3)-H(4)folate have been determined at 1.65 and 1.70 A resolution, respectively. The structures show that the folate is bound in a catalytically competent conformation, and Leu223 undergoes a conformational change similar to that observed for Phe223 in the Glu28Gln-CH(3)-H(4)folate structure. Taken together, our results suggest that Leu may be a suitable replacement for Phe223 in the oxidative half-reaction of E. coli MTHFR.

Cofactor chemogenomics

Cofactors are organic molecules, most of them originating from vitamins, that bind to enzymes making them able to catalyze defined reactions. A cofactor-based chemogenomics approach exploits the presence of a cofactor-binding domain to develop compound scaffolds tailored to mimic the cofactor and to replace it within target enzyme classes. As a result, a loss of function is observed. An expansion of the cofactor scaffold to include structural/chemical features derived from the substrate, that usually binds at cofactor adjacent sites, increases the specificity of the enzyme fishing. This approach has been so far applied only to NAD(P)(+)-dependent enzymes. However, it is suitable for all other cofactors, with difficulties, for some of them, originated by very tight binding. In the case of cofactors covalently bound to the enzyme, the competition between the natural cofactor and the cofactor scaffold mimic can only occur during enzyme folding.

Biomolecular approach to the design of potential drugs

The approach to drug design on the basis of molecular-level information on biological processes is being driven by the expanding knowledge of the details of molecular events in biological systems. We have directed attention to the design of potentially active compounds based on the aforementioned "biomolecular" concepts. Selected examples from our studies are discussed. This paper presents three case studies of approaches to the development of potential medicinal agents whose design has evolved from considerations of molecular mechanisms of processes in selected biological systems.

Role of N-terminal residues in the ubiquitin-independent degradation of human thymidylate synthase

Thymidylate synthase (TS) catalyses the reductive methylation of dUMP to form dTMP, a reaction that is essential for maintenance of nucleotide pools during cell growth. Because the enzyme is indispensable for DNA replication in actively dividing cells, it is an important target for cytotoxic drugs used in cancer chemotherapy, including fluoropyrimidines (e.g. 5-fluorouracil and 5-fluoro-2'-deoxyuridine) and anti-folates (e.g. raltitrexed, LY231514, ZD9331 and BW1843U89). These drugs generate metabolites that bind to the enzyme's active site and inhibit catalytic activity, leading to thymidylate deprivation and cellular apoptosis. Ligand binding to TS results in stabilization of the enzyme and an increase in its intracellular concentration. Previously, we showed that degradation of the TS polypeptide is carried out by the 26 S proteasome in a ubiquitin-independent manner. Such degradation is directed by the disordered N-terminal region of the TS polypeptide, and is abrogated by ligand binding. In the present study, we have verified the ubiquitin-independent nature of TS proteolysis by showing that a 'lysine-less' polypeptide, in which all lysine residues were replaced by arginine, is still subject to proteasome-mediated degradation. In addition, we have mapped the structural determinants of intracellular TS degradation in more detail and show that residues at the N-terminal end of the molecule, particularly the penultimate amino acid Pro2, play an important role in governing the half-life of the enzyme. This region is capable on its own of destabilizing an evolutionarily distinct TS molecule that normally lacks this domain, indicating that it functions as a degradation signal. Interestingly, degradation of an intrinsically unstable mutant form of TS, containing a Pro-->Leu substitution at residue 303, is directed by C-terminal, rather than N-terminal, sequences. The implications of these findings for the control of TS expression, and for the regulation of protein degradation in general, are discussed.

Kinetic Characterization of Bifunctional Thymidylate Synthase-Dihydrofolate Reductase (TS-DHFR) from Cryptosporidium hominis

This study presents a kinetic characterization of the recently crystallized bifunctional thymidylate synthasedihydrofolate reductase (TS-DHFR) enzyme from the apicomplexa parasite, Cryptosporidium hominis. Our study focuses on determination of the C. hominis TS-DHFR kinetic mechanism, substrate channeling behavior, and domain-domain communication. Unexpectedly, the unique mechanistic features of C. hominis TS-DHFR involve the highly conserved TS domain. At 45 s(-) (1), C. hominis TS activity is 10-40-fold faster than other TS enzymes studied and a new kinetic mechanism was required to simulate C. hominis TS behavior. A large accumulation of dihydrofolate produced at TS and a lag in product formation at DHFR were observed. These observations make C. hominis TS-DHFR the first bifunctional TS-DHFR enzyme studied for which there is clear evidence against dihydrofolate substrate channeling. Furthermore, whereas with Leishmania major TS-DHFR there are multiple lines of evidence for domain-domain communication (ligand binding at one active site affecting activity of the other enzyme), no such effects were observed with C. hominis TS-DHFR.

GEMDOCK: A generic evolutionary method for molecular docking

We have developed an evolutionary approach for flexible ligand docking. This approval, GEMDOCK, uses a Generic Evolutionary Method for molecular DOCKing and an empirical scoring function. The former combines both discrete and continuous global search strategies with local search strategies to speed up convergence, whereas the latter results in rapid recognition of potential ligands. GEMDOCK was tested on a diverse data set of 100 protein-ligand complexes from the Protein Data Bank. In 79% of these complexes, the docked lowest energy ligand structures had root-mean-square derivations (RMSDs) below 2.0 A with respect to the corresponding crystal structures. The success rate increased to 85% if the structure water molecules were retained. We evaluated GEMDOCK on two cross-docking experiments in which each ligand of a protein ensemble was docked into each protein of the ensemble. Seventy-six percent of the docked structures had RMSDs below 2.0 A when the ligands were docked into foreign structures. We analyzed and validated GEMDOCK with respect to various search spaces and scoring functions, and found that if the scoring function was perfect, then the predicted accuracy was also essentially perfect. This study suggests that GEMDOCK is a useful tool for molecular recognition and may be used to systematically evaluate and thus improve scoring functions.

Probing electrostatic channeling in protozoal bifunctional thymidylate synthase-dihydrofolate reductase using site-directed mutagenesis

In this study we used site-directed mutagenesis to test the hypothesis that substrate channeling in the bifunctional thymidylate synthase-dihydrofolate reductase enzyme from Leishmania major occurs via electrostatic interactions between the negatively charged dihydrofolate produced at thymidylate synthase and a series of lysine and arginine residues on the surface of the protein. Accordingly, 12 charge reversal or charge neutralization mutants were made, with up to 6 putative channel residues changed at once. The mutants were assessed for impaired channeling using two criteria: a lag in product formation at dihydrofolate reductase and an increase in dihydrofolate accumulation. Surprisingly, none of the mutations produced changes consistent with impaired channeling, so our findings do not support the electrostatic channeling hypothesis. Burst experiments confirmed that the mutants also did not interfere with intermediate formation at thymidylate synthase. One mutant, K282E/R283E, was found to be thymidylate synthase-dead because of an impaired ability to form the covalent enzyme-methylene tetrahydrofolate-deoxyuridate complex prerequisite for chemical catalysis.

Methylenetetrahydrofolate reductase: a common human polymorphism and its biochemical implications

Methlenetetrahydrofolate (CH2-H4folate) is required for the conversion of homocysteine to methionine and of dUMP to dTMP in support of DNA synthesis, and also serves as a major source of one carbon unit for purine biosynthesis. This review presents biochemical studies of a human polymorphism in methylenetetrahydrofolate reductase, which catalyzes the reaction shown below. The mutation decreases the flux of CH2-H4folate into CH3-H4folate, and is associated with both beneficial and deleterious effects that can be traced to the molecular effect of the substitution of alanine 222 by valine.

The only active mutant of thymidylate synthase D169, a residue far from the site of methyl transfer, demonstrates the exquisite nature of enzyme specificity

Cysteine is the only variant of D169, a cofactor-binding residue in thymidylate synthase, that shows in vivo activity. The 2.4 A crystal structure of Escherichia coli thymidylate synthase D169C in a complex with dUMP and the antifolate CB3717 shows it to be an asymmetric dimer, with only one active site covalently bonded to dUMP. At the active site with covalently bound substrate, C169 S gamma adopts the roles of both carboxyl oxygens of D169, making a 3.6 A S...H[bond]N hydrogen bond to 3-NH of CB3717 and a 3.4 A water-mediated hydrogen bond to H212. Analogous hydrogen bonds formed during the enzyme reaction are important for cofactor binding and are postulated to contribute to catalysis. The C169 side chain is likely to be ionized, making it a better hydrogen bond acceptor than a neutral sulfhydryl group. At the second active site, C169 S gamma makes a shorter (3 A) hydrogen bond to the 3-NH of CB3717, CB3717 is approximately 1.5 A out of its binding site and there is no covalent bond between dUMP and the catalytic cysteine. Changes to partitioning among productive and non-productive conformations of reaction intermediates may contribute as much, if not more, to the diminished activity of this mutant than reduced stabilization of transition states.

Mechanistic characterization of Toxoplasma gondii thymidylate synthase (TS-DHFR)-dihydrofolate reductase. Evidence for a TS intermediate and TS half-sites reactivity

This study describes the use of rapid transient kinetic methods to characterize the bifunctional thymidylate synthase-dihydrofolate reductase (TS-DHFR) enzyme from Toxoplasma gondii. In addition to elucidating the detailed kinetic scheme for this enzyme, this work provides the first direct kinetic evidence for the formation of a TS intermediate and for half-sites TS reactivity in human and Escherichia coli monofunctional TS and in T. gondii and Leishmania major bifunctional TS-DHFR. Comparison of the T. gondii TS-DHFR catalytic mechanism to that of the L. major enzyme reveals the mechanistic differences to be predominantly in DHFR activity. Specifically, TS ligand induced domain-domain communication involving DHFR activation is observed only in the L. major enzyme and, whereas both DHFR activities involve a rate-limiting conformational change, the change occurs at different positions along the kinetic pathway.

Ligand-receptor docking with the Mining Minima optimizer

The optimizer developed for the Mining Minima algorithm, which uses ideas from Genetic Algorithms, the Global Underestimator Method, and Poling, has been adapted for use in ligand-receptor docking. The present study describes the resulting methodology and evaluates its accuracy and speed for 27 test systems. The performance of the new docking algorithm appears to be competitive with that of previously published methods. The energy model, an empirical force field with a distance-dependent dielectric treatment of solvation, is adequate for a number of test cases, although incorrect low-energy conformations begin to compete with the correct conformation for larger sampling volumes and for highly solvent-exposed binding sites that impose little steric constraint on the ligand.

Structures of trihydroxynaphthalene reductase-fungicide complexes: implications for structure-based design and catalysis

BACKGROUND

Trihydroxynaphthalene reductase catalyzes two intermediate steps in the fungal melanin biosynthetic pathway. The enzyme, a typical short-chain dehydrogenase, is the biochemical target of three commercial fungicides. The fungicides bind preferentially to the NADPH form of the enzyme.

RESULTS

Three X-ray structures of the Magnaporthe grisea enzyme complexed with NADPH and two commercial and one experimental fungicide were determined at 1.7 A (pyroquilon), 2.0 A (2,3-dihydro-4-nitro-1H-inden-1-one, 1), and 2.1 A (phthalide) resolutions. The chemically distinct inhibitors occupy similar space within the enzyme's active site. The three inhibitors share hydrogen bonds with the side chain hydroxyls of Ser-164 and Tyr-178 via a carbonyl oxygen (pyroquilon and 1) or via a carbonyl oxygen and a ring oxygen (phthalide). Active site residues occupy similar positions among the three structures. A buried water molecule that is hydrogen bonded to the NZ nitrogen of Lys-182 in each of the three structures likely serves to stabilize the cationic form of the residue for participation in catalysis.

CONCLUSIONS

The pro S hydrogen of NADPH (which is transferred as a hydride to the enzyme's naphthol substrates) is directed toward the carbonyl carbon of the inhibitors that mimic an intermediate along the reaction coordinate. Modeling tetrahydroxynaphthalene and trihydroxynaphthalene in the active site shows steric and electrostatic repulsion between the extra hydroxyl oxygen of the former substrate and the sulfur atom of Met-283 (the C-terminal residue), which accounts, in part, for the 4-fold greater substrate specificity for trihydroxynaphthalene over tetrahydroxynaphthalene.

The second naphthol reductase of fungal melanin biosynthesis in Magnaporthe grisea: tetrahydroxynaphthalene reductase

Mutants of Magnaporthe grisea harboring a defective gene for 1,3, 8-trihydroxynaphthalene reductase retain the capability to produce scytalone, thus suggesting the existence of a second naphthol reductase that can catalyze the reduction of 1,3,6, 8-tetrahydroxynaphthalene to scytalone within the fungal melanin biosynthetic pathway. The second naphthol reductase gene was cloned from M. grisea by identification of cDNA fragments with weak homology to the cDNA of trihydroxynaphthalene reductase. The amino acid sequence for the second naphthol reductase is 46% identical to that of trihydroxynaphthalene reductase. The second naphthol reductase was produced in Esherichia coli and purified to homogeneity. Substrate competition experiments indicate that the second reductase prefers tetrahydroxynaphthalene over trihydroxynaphthalene by a factor of 310; trihydroxynaphthalene reductase prefers trihydroxynaphthalene over tetrahydroxynaphthalene by a factor of 4.2. On the basis of the 1300-fold difference in substrate specificities between the two reductases, the second reductase is designated tetrahydroxynaphthalene reductase. Tetrahydroxynaphthalene reductase has a 200-fold larger K(i) for the fungicide tricyclazole than that of trihydroxynaphthalene reductase, and this accounts for the latter enzyme being the primary physiological target of the fungicide. M. grisea mutants lacking activities for both trihydroxynaphthalene and tetrahydroxynaphthalene reductases do not produce scytalone, indicating that there are no other metabolic routes to scytalone.

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.

Random sequence mutagenesis and resistance to 5-fluorouridine in human thymidylate synthases

Thymidylate synthase (TS) catalyzes the methylation of dUMP to dTMP and is the target for the widely used chemotherapeutic agent 5-fluorouracil. We used random sequence mutagenesis to replace 13 codons within the active site of TS and obtain variants that are resistant to 5-fluorodeoxyuridine (5-FdUR). The resulting random library was selected for its ability to complement a TS-deficient Escherichia coli strain, and sequence analysis of survivors found multiple substitutions to be tolerable within the targeted region. An independent selection of the library was carried out in the presence of 5-FdUR, resulting in a more limited spectrum of mutations. One specific mutation, C199L, was observed in more than 46% of 5-FdUR-resistant clones. A 5-FdUR-resistant triple mutant, A197V/L198I/C199F, was purified to apparent homogeneity. Kinetic studies with the substrate dUMP indicate that this mutant is similar to the wild type in regards to kcat and Km values for dUMP and the cosubstrate CH2H4-folate. In contrast, equilibrium binding studies with the inhibitor, FdUMP, demonstrate that the dissociation constant (Kd) for FdUMP binding into the ternary complex was 20-fold higher than values obtained for the wild-type enzyme. This 5-FdUMP-resistant mutant, or others similarly selected, is a candidate for use in gene therapy to render susceptible normal cells resistant to the toxic effects of systemic 5-fluorouracil.

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.

Development and validation of a genetic algorithm for flexible docking

Prediction of small molecule binding modes to macromolecules of known three-dimensional structure is a problem of paramount importance in rational drug design (the "docking" problem). We report the development and validation of the program GOLD (Genetic Optimisation for Ligand Docking). GOLD is an automated ligand docking program that uses a genetic algorithm to explore the full range of ligand conformational flexibility with partial flexibility of the protein, and satisfies the fundamental requirement that the ligand must displace loosely bound water on binding. Numerous enhancements and modifications have been applied to the original technique resulting in a substantial increase in the reliability and the applicability of the algorithm. The advanced algorithm has been tested on a dataset of 100 complexes extracted from the Brookhaven Protein DataBank. When used to dock the ligand back into the binding site, GOLD achieved a 71% success rate in identifying the experimental binding mode.

The complex of the anti-cancer therapeutic, BW1843U89, with thymidylate synthase at 2.0 A resolution: implications for a new mode of inhibition

BACKGROUND

Thymidylate synthase (TS) is critical to DNA synthesis as it catalyzes the rate limiting step in the only biosynthetic pathway for deoxythymidine monophosphate (dTMP) production. TS is therefore an important target for anti-proliferative and anti-cancer drug design. The TS enzymatic mechanism involves the reductive methylation of the substrate, deoxyuridine monophosphate (dUMP), by transfer of a methylene group from the co-factor, methylenetetrahydrofolate (CH2H4folate), resulting in the production of deoxythymidine monophosphate (dTMP) and dihydrofolate (H2folate). Previous drug design efforts based on co-factor analogues have produced good inhibitors of TS, but poor bioavailability and toxicity have limited their usefulness. BW1843U89, a folate analogue, is a recently developed compound which is an exceptionally strong inhibitor (Ki = 0.09 nM), has good bioavailability and in clinical trials thus far has not demonstrated significant toxicity.

RESULTS

We report the crystal structure of E. coli TS in ternary complex with dUMP and BW1843U89 at 2.0 A resolution. Although the benzoquinazoline ring system of the inhibitor binds to TS in much the same manner as previously determined for H2folate and CB3717, the larger size of the ligand is accommodated by the enzyme through a local distortion of the active site, that is not strictly conserved in both monomers in the asymmetric unit. Several conserved waters that had been previously implicated in mechanistic roles have been displaced.

CONCLUSIONS

BW1843U89 forms a ternary complex with dUMP and completes with CH2H4 folate at the active site. Inhibition of TS by BW1843U89 shows four unique aspects in its mechanism of action. BW1843U89 prevents the Michael addition of dUMP to Cys146, in contrast to the mechanisms implicated from crystallography of other quinazoline based inhibitors; displaces a catalytic water from the active site; reorders a peptide loop (Leu72-Trp83) in the active site; and is unique amongst the antifolates in inactivating TS at a stoichiometric ratio of one molecule per TS dimer. Thus, it exploits the principles of negative cooperativity that are increasingly being recognized in the catalytic mechanism of the enzyme per se. The structure suggests that this 'half-the-sites' effect is catalytic and not related to ligand binding. Therefore BW1843U89 is both a competitive inhibitor (at the binding site) and a non-competitive inhibitor at the other site.

Electrostatic guidance of catalysis by a conserved glutamic acid in Escherichia coli dTMP synthase and bacteriophage T4 dCMP hydroxymethylase

Thymidylate synthase (TS) and dCMP hydroxymethylase (CH) are homologous enzymes which catalyze the alkylation of C5 of pyrimidine nucleotides. One of the first catalytic steps is isomerization of the alkyl donor, methylenetetrahydrofolate, from its N5,N10 bridged form to the N5 iminium ion upon enzyme binding. Glu58 in TS has been postulated [Matthews et al. (1990) J. Mol. Biol. 214, 937-948] to be involved in this isomerization and the deprotonation of C5 of the nucleotide. Substitution by Asp or Gln of Glu58 in Escherichia coli TS, or of the corresponding Glu60 in CH from phage T4, decreases the activity of either enzyme. Alkylation is slowed much more than deprotonation, indicating uncoupling of steps which are tightly coupled for the wild-type enzymes. The data support minor roles for Glu58/60 in nucleotide binding and in isomerization of methylenetetrahydrofolate, but no major roles in nucleotide deprotonation, product dissociation, or hydration catalyzed by CH. The primary role of Glu58/60 is to accelerate bond cleavage between N5 of tetrahydrofolate and the methylene being transferred. The influence of Glu58/60 on the rate of bond cleavage is proposed to arise from electrostatic destabilization due to the proximity of the glutamyl carboxylate, of the anionic species formed when C5 of the nucleotide is deprotonated. The proposal explains the uncoupling of deprotonation and alkylation with the Glu58/60 variants and the reduced kinetic isotope effect on hydride transfer for TS(Glu58Gln). The inability of 5-deazatetrahydrofolate to stimulate enzyme-catalyzed tritium exchange from [5-(3H)]nucleotides into solvent suggests that N5 of tetrahydrofolate is the base which deprotonates the nucleotide.

Evidence from 18O exchange studies for an exocyclic methylene intermediate in the reaction catalyzed by T4 deoxycytidylate hydroxymethylase

18O exchange experiments were designed to identify the final intermediate in the catalytic mechanism of bacteriophage T4 deoxycytidylate (dCMP) hydroxymethylase (CH). CH catalyzes the formation of 5-(hydroxymethyl)-dCMP (HmdCMP) from dCMP and methylenetetrahydrofolate (CH2-THF). CH resembles thymidylate synthase (TS), an enzyme of known three-dimensional structure, in both amino acid sequence and the reaction catalyzed. The final intermediate in the reaction catalyzed by TS or CH has been proposed to be the nucleotide with an exocyclic 5-methylene group covalently linked to the enzyme. This intermediate is then hydrated to HmdCMP (by CH) or reduced to deoxythymidylate (by TS). We report here that CH catalyzes the incorporation of 18O from solvent water into the product, HmdCMP, in the presence of tetrahydrofolate (THF). The cause of this exchange is a reverse reaction followed by a resynthesis. CH also catalyzes the exchange of 18O from solvent water into HmdCMP in the absence of exogenous THF and in the presence of THF analogues that lack N-5. N-5 is the nitrogen that is likely to be bound to the methylene as it is transferred to dCMP. A CH variant that lacks the nucleophilic Cys 148 is incapable of promoting these 18O exchange reactions. The THF analogues lacking N-5 do not promote a CH-catalyzed reverse reaction. Rather, we propose that the CH-catalyzed 18O exchange reaction promoted by these THF analogues occurs via 5-methylene-dCMP linked to the enzyme through Cys 148. We conclude here that enzyme-bound 5-methylene-dCMP is the final intermediate during catalysis by CH, as has also been proposed for TS and dUMP.

Structure of and kinetic channelling in bifunctional dihydrofolate reductase-thymidylate synthase

The bifunctional enzyme dihydrofolate reductase-thymidylate synthase catalyses both the reductive methylation of 2'-deoxyuridylate and the subsequent reduction of dihydrofolate to yield 2'-deoxythymidylate and tetrahydrofolate at two spacially discrete sites situated on different protein domains. The X-ray structure of dihydrofolate reductase-thymidylate synthase from Leishmania major indicates that transfer of dihydrofolate between these sites does not occur by transient binding at both sites but rather by movement of dihydrofolate across the surface of the protein. The enzyme has an unusual surface charge distribution that could account for this channelling of dihydrofolate between active sites.