Flavin-Dependent Thymidylate Synthase as a Drug Target for Deadly Microbes: Mutational Study and a Strategy for Inhibitor Design

Irimpan I Mathews^*

Stanford Synchrotron Radiation Lightsource, Stanford University, Menlo Park, CA 94025, USA

*Corresponding Author:

Irimpan I Mathews
Stanford Synchrotron Radiation Lightsource
Stanford University Menlo Park
CA 94025, USA
Tel: (650) 926 5105
Fax: (650) 926 2258/3292
E-mail: iimathew@slac.stanford.edu

Received Date: February 14, 2013; Accepted Date: April 18, 2013; Published Date: April 22, 2013

Citation: Mathews II (2013) Flavin-Dependent Thymidylate Synthase as a Drug Target for Deadly Microbes: Mutational Study and a Strategy for Inhibitor Design. J Bioterr Biodef S12:004.. doi: 10.4172/2157-2526.S12-004

Copyright: © 2013 Mathews II. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are credited.

Visit for more related articles at Journal of Bioterrorism & Biodefense

Abstract

The identification of flavin-dependent thymidylate synthase (FDTS) as an essential enzyme and its occurrence in several pathogenic microbes opens opportunities for using FDTS enzyme as an excellent target for new antimicrobial drug discovery. In contrast to the human thymidylate synthase enzyme that utilizes methylene-tetrahydrofolate (CH2H4 folate) for the conversion of dUMP to dTMP, the microbial enzymes utilize an additional non-covalently bound FAD molecule for the hydride transfer from NAD(P)H. The structural and mechanistic differences between the human and microbial enzymes present an attractive opportunity for the design of antimicrobial compounds specific for the pathogens. We have determined the crystal structure of FDTS enzyme in complex with the methyl donor, CH2H4 folate. We describe here the structure of a FDTS mutant and compare it with other FDTS complex structures, including a FDTS-CH2H4 folate complex. We identified a conformational change essential for substrate binding and propose a strategy for the design of FDTS specific inhibitors.

Keywords

FDTS enzyme; Structure/function studies; Drug design; Crystallography

Introduction

The de novo biosynthesis of thymidylate (2’-deoxythymine-5’- monophosphate; dTMP), one of the four bases of DNA, requires the enzyme thymidylate synthase [1]. Two types of thymidylate synthases have been described and both of them use 2’-deoxyuridine-5’- monophosphate (dUMP) as the substrate [1,2]. The classical thymidylate synthases (TS) use N⁵,N¹⁰-methylene-5,6,7,8-tetrahydrofolate (CH₂H₄ folate) to reductively methylate dUMP producing dTMP, while the recently identified flavin-dependent thymidylate synthase (FDTS) uses a non-covalently bound flavin adenine nucleotide (FAD) for the reduction [2]. FDTS is found in ~30% of microbial genome. The two families of thymidylate synthases are mechanistically and structurally different [1-4]. Our recent studies have shown that, unlike the classical enzyme which uses a cysteine residue to form a covalent bond with dUMP, the flavin-dependent enzyme doesn’t use an enzymatic nucleophile for the reaction [3]. The uniqueness of the FDTS enzyme is also revealed by a novel fold of its structure [4]. The structures of FDTS from various organisms share similar fold, and the high level of sequence similarity of FDTS from other organisms indicates very similar structures for all of them [5-7].

The rise in bacterial resistance has stimulated new interest in finding novel targets for the development of effective antimicrobial agents. The presence of FDTS in many pathogenic organisms (Figure 1) and its absence in human make FDTS as an attractive target for antimicrobials [2] and a number of studies are in progress to develop specific inhibitors for the FDTS enzymes [8,9]. The catalytic mechanism of classical enzyme is well understood and has facilitated the development of several inhibitors, some of which are in clinical use as anticancer drugs (e.g., 5-flouro-uracil, tomudex (Raltitrexed)) [1,10]. Several structures of the classical enzyme, including ternary complexes with various combinations of substrate and folate cofactor, along with their analogs are available [1,11]. Unfortunately, the inhibitors for the classical thymidylate synthase are not specific to the FDTS enzymes [12]. The complexity of the FDTS reaction mechanism and the conformational flexibility of the active site region make it difficult to perform rational drug design with the currently available information. There are opposing views regarding the most important methylene-transfer step, with some studies proposing an indirect methylene-transfer through an arginine residue [13] while other studies indicating a direct methylene transfer from CH₂H₄ folate to dUMP [3,6,12,14]. Therefore, it is important to understand the details of the FDTS mechanism and determine its structures in various complexes and intermediates.

We have recently reported the first structures of the quaternary complexes of FDTS from Thermotoga maritima (TmFDTS) with FAD, dUMP and CH₂H₄ folate and CH₂H₄ folate mimics. Since several of the inhibitors of classical thymidylate synthase are based on the folate binding site and not selective for FDTS enzymes, it is expected that novel compounds utilizing the unique folate binding modes may provide new avenues for FDTS specific inhibitor design [15]. This emphasizes the importance of a proper understanding of the binding interactions near the folate binding site.

One of the residues implicated in the folate binding interaction in FDTS is histidine 53 (T. maritima numbering). This residue is fully conserved among the FDTS from various organisms and previous studies showed the essential role of this residue in NAD(P)H oxidation or methyl transfer [6]. The methylene transfer step is one of the least understood processes in the FDTS catalysis. The recent structures of the ternary complexes of TmFDTS with FAD, dUMP and CH₂H₄ folate and identified the folate binding site and proposed it as a binding site for NADPH [16]. One of the residues implicated in the folate binding interaction is histidine 53. We mutated this residue to aspartic acid (H53D) and present the structures of the H53D-FAD and H53D-FADdUMP complexes and a comparison with native enzyme structures. Earlier we reported the crystal structure of the H53A mutant and its complex with FAD, dUMP and CH₂H₄ folate [16]. We also reported that both the H53A and H53D mutants showed dTMP formation with dramatically reduced activity (Table S2 of reference 16).

Results and Discussion

We have crystallized and solved the structures of H53D mutant of the Thermotoga maritima FDTS with FAD and in complex with FAD and dUMP (Table 1). The structures of the H53D mutant complexes are very similar to the native enzyme, which forms a biologically active tetramer. An extensive array of hydrogen bonding and hydrophobic interactions stabilize the tetrameric structure with ~2000Ų surface area buried per monomer. Previous crystallographic and activity studies have confirmed the presence of each active site at the interface of the three subunits [4,17]. The two interacting active sites in each side of the enzyme form a large active site grove spanning around 50Å. A tightly bound FAD molecule is observed in the all of the reported structures. However, a structure of the apoenzyme obtained by removing the FAD using high amounts of NaCl showed that FAD molecule is not essential for the stabilization of the tetramer [4]. The structures of the complexes presented here show that the substrate-binding loop can be stabilized in two conformations and this affects the binding of the molecules at the substrate binding site.

Table 1: Crystallographic parameters, data collection and refinement statistics.

FAD binding site

FAD acts as the reducing agent in the FDTS reaction. The ribityl and the AMP groups are strongly bound in the active site with the catalytically important flavin ring exposed to the solvent [4]. In the dUMP complexes, flavin ring of the FAD molecule stacks with the pyrimidine ring of the dUMP. It has been reported that the flavin ring of the FAD molecule is generally disordered in structures without the dUMP [4]. This is true for the current H53D complex with FAD. However, in the viral enzyme and the coryne enzyme flavin ring is stabilized by stacking interaction with the histidine 53 side chain [6,18]. Interestingly, the flavin ring uses the si-face and re-face for the stacking interaction in the viral and coryne enzymes, respectively. In the reported structure of the quaternary complex with FAD, dUMP and CH₂H₄ folate, the flavin ring uses the re-face to stack with the histidine side chain. It is also interesting to note that during the folate stacking histidine 53 side chain flips to the opposite side (torsion angle N-Cα- Cβ-Cγ=~-172° for viral and coryne enzymes and ~-56° for the folate bound complex). It is important to note that flavin ring uses the si-face to stack with dUMP [4] as well as the CH₂H₄ folate [16]. The folate/FADdependent tRNA T54 methyltransferase (TrmFO), which catalyzes the same net reaction as the FDTS enzyme, the re-face of the flavin is stacked with the folate [19]. Our earlier studies with two mutants of FDTS (E144R and R174K (ref 17) (R174K+FAD+dUMP work is not published)) with FAD and in complex with FAD and dUMP indicated that the flavin is able to rotate in the active site during the formation of the dUMP complex [16].

The details mentioned above show that isoalloxazine (flavin) ring of FAD binds in a big pocket that tolerates large movements of the isoalloxazine ring. Importantly, the isoalloxazine ring is able to rotate in the binding pocket and utilize same face of the ring to bind to substrate and cofactors. This is in contrast to the relatively rigid binding mode observed for the isoalloxazine ring in most of the enzymes that use FAD as the cofactor [20-23]. The presence of the large active site cavity in FDTS that tolerates major conformational movements of the ligands makes the design of specific inhibitors very challenging.

The FAD molecules in the H53D+FAD complex show very weak density for the whole FAD molecule and no density for the flavin ring (Table 2, Figure 2a). The FAD molecules in the H53D+FAD+dUMP complex also showed weak electron density indicating poor binding (Table 2, Figure 2b). This is in contrast to the flavin ring only disorder observed for the native enzyme with FAD complex and the very good electron density observed for FAD and dUMP in the FAD-dUMP complex (Table 2) [4].

Table 2: Comparison of the temperature factors of various protein chains, FAD molecules and dUMP molecules in the mutant and native enzyme^a.

Substrate binding site

In general, dUMP and analogs are strongly bound in the enzyme with several direct and water mediated hydrogen bonds to the protein. Furthermore, the pyrimidine ring of dUMP is stacked to the flavin ring of FAD in complexes with FAD. It has also been reported that substrate induced conformational changes near the active site is important in the stabilization of the substrate binding site [4]. A main difference between the current and the reported structures is the very weak electron density observed for the dUMP (Table 2, Figure 2b). Only two of the active sites showed good electron density for dUMP, while the third active site showed weak density for dUMP, the fourth one showed very weak density only for the phosphate group. It is not clear whether differences in electron density between the four active sites indicate any allosteric interaction amongst the active sites.

Open and closed confirmations

There are several mechanisms proposed for the FDTS catalysis with various suggestions for the binding and release of the substrate and other cofactors [3]. Unfortunately, the large conformational flexibility of the FDTS active site makes it difficult to give a structural perspective to the biochemical results. It has been reported that the conformational changes during FAD and dUMP binding brings various conserved residues into close proximity to these molecules.

We compared the native enzyme structure with the FAD complex, with FAD and dUMP complex, and FAD, dUMP and CH₂H₄ folate complex and identified two major conformational changes during various binding processes (Figure 3). Various combinations of these conformational changes take place during the binding of the substrate and/or cofactors. The close to open conformational change of the 90-loop/substrate-binding loop is very important because this conformational change brings important residues to the substrate binding site [4]. In the open conformation of the substrate-binding loop, residues from Ser88 to Arg90 make hydrogen-bonding interactions with the substrate. While the Ser88 Oγ and Gly89 N atoms H-bonds to the phosphate group of the substrate, the Arg90 side chain H-bonds to one of the oxygen atoms of the pyrimidine base. The Ser88 and Arg90 are highly conserved residues [16].

A comparison of the active sites of the H53D+dUMP complex shows that the substrate-binding loop conformational change plays an important role in the stabilization of the dUMP binding (Table 2, Figure 4). The active sites that show good electron density for dUMP (chains A and B) showed closed conformation for the substrate-binding loop. The dUMP molecule in chain C showed weaker density and the substratebinding loop showed double conformation. The open confirmation observed in chain D showed very weak density for dUMP with density for the phosphate group only. This shows that the open conformation of the substrate-binding loop doesn’t favor the substrate binding. These conformational changes may also be important for the binding and release of the substrate and product.

A closer examination of the open and closed conformation of the substrate-binding loop shows that the open conformation is stabilized by hydrogen bonding interaction of the tyrosine 91 hydroxyl group to the mutated aspartic acid (Figure 5). Similar hydrogen bonding interaction of the tyrosine 91 from the open loop with histidine 53 is observed in the native enzyme FAD complex (PDB code: 1O2A). This hydrogen bonding interaction is absent in the closed conformation and the distance between the corresponding atoms in the closed conformation is around 8 Å. The structural changes accompanying the open conformation also brings the conserved arginine 90 to the vicinity of tyrosine 47. In the closed conformation of the substratebinding loop, arginine 90 side chain is involved in hydrogen bonding interactions with the substrate and protein atoms from the neighboring protein chain. These interactions stabilize the substrate binding site. The tyrosine 47 and 91 residues generally show good conservation among the FDTS enzymes [16].

The observed stabilization of the closed conformation substratebinding loop in the mutated protein suggests the possibility of using chemical compounds to lock the open conformation of the substratebinding loop. Since closed conformation of the substrate-binding loop is very important for substrate binding, design of chemicals to lock the open conformation may be a good strategy to develop inhibitors specific for the FDTS enzymes. The recently discovered 150-cavity in group-1 influenza A neuraminidase provided a target for rational structurebased drug development and novel techniques have been developed to lock open the 150-loop as a strategy for the inhibition [24,25]. An analysis of the reported structures of various FDTS enzymes shows that FDTS tolerates large movements of the ligands in the binding pocket, thus making the design of specific inhibitors very challenging.

Conclusions

FDTS is an essential enzyme found in several pathogenic microbes. Because of the structural and mechanistic differences between FDTS and the human enzyme and the crucial role of FDTS enzyme in bacterial cells, the FDTS enzymes have been proposed as a priority target for developing new anti-microbial compounds [2,26]. Unfortunately, because of the complex nature of the FDTS reaction catalysis and the non-specificity of the known TS inhibitors for FDTS enzyme, it has been difficult to develop FDTS specific inhibitors. We have shown that conformational changes of active site are important for the binding of the substrate and various cofactors. Our data shows that the closed conformation of the substrate-binding loop is essential for substrate binding. We propose the development of compounds that can lock the open conformation of the substrate-binding loop as a strategy for FDTS specific inhibitor design.

Materials and Methods

Chemicals

All chemicals were reagent grade and used as purchased without further purification, unless specified.

Protein expression and purification

The H53D mutant of FDTS from T. maritima (TM0449, GenBank accession number NP228259) was expressed and purified as previously described [27].

Crystallization and structure determination

The crystals of the H53D mutant with FAD and with FAD and dUMP were crystallized at 22°C in 50-60% (w/v) PEG 200 and 100 mM Tris buffer, pH 8.0. The FAD molecule stays bound during purification and no further FAD was included in the crystallization trials. The dUMP complex was prepared by treating the FAD complex with 10 mM dUMP. The crystals were flash cooled directly from the drop. Diffraction data were collected at the Stanford Synchrotron Radiation Lightsource (SSRL) beamline 9-2 using Q315 detector. The wavelengths used for the data collection of the H53D with FAD and the dUMP complexes were 0.9795 Å and 1.0 Å, respectively. All data were integrated using the XDS package [28]. These crystals belonged to the P2₁2₁2₁ space group.

Structures of the complexes were solved by molecular replacement (MOLREP [29]) or rigid body refinement using the T. maritima tetramer (PDB code: 1O26) as the search template. Model building and refinement were performed by Coot [30] and REFMAC [31]. The Ramachandran statistics for the final structures showed no outliers Structures of the complexes were solved by molecular replacement (MOLREP [29]) or rigid body refinement using the T. maritima tetramer (PDB code: 1O26) as the search template. Model building and refinement were performed by Coot [30] and REFMAC [31]. The Ramachandran statistics for the final structures showed no outliers

Coordinates

Coordinates for the complexes have been deposited in the Protein Data Bank (accession codes: 4KAR (H53D+FAD complex) and 4KAS (H53D+FAD+dUMP complex).

Acknowledgements

I thank S. A. Lesley, H. Klock, and E. Ambing (The Genomics Institute of the Novartis Research Foundation) for the protein samples and Q. Xu and A. Kumar for critical reading of the manuscript. I thank members of the SMB group at SSRL for helpful discussions and support. Portions of this research were carried out at the Stanford Synchrotron Radiation Lightsource, a Directorate of SLAC National Accelerator Laboratory and an Office of Science User Facility operated for the U.S. Department of Energy Office of Science by Stanford University. The SSRL Structural Molecular Biology Program is supported by the DOE Office of Biological and Environmental Research, and by the National Institutes of Health, National Center for Research Resources, Biomedical Technology Program (P41RR001209), and the National Institute of General Medical Sciences.

References

1. Finer-Moore JS, Santi DV, Stroud RM (2003) Lessons and conclusions from dissecting the mechanism of a bisubstrate enzyme: thymidylate synthase mutagenesis, function, and structure. Biochemistry 42: 248-256.
2. Myllykallio H, Lipowski G, Leduc D, Filee J, Forterre P, et al. (2002) An alternative flavin-dependent mechanism for thymidylate synthesis. Science 297: 105-107.
3. Koehn EM, Fleischmann T, Conrad JA, Palfey BA, Lesley SA, et al. (2009) An unusual mechanism of thymidylate biosynthesis in organisms containing the thyX gene. Nature 458: 919-923.
4. Mathews II, Deacon AM, Canaves JM, McMullan D, Lesley SA, et al. (2003) Functional analysis of substrate and cofactor complex structures of a thymidylate synthase-complementing protein. Structure 11: 677-690.
5. Wang K, Wang Q, Chen J, Chen L, Jiang H, et al. (2011) Crystal structure and enzymatic characterization of thymidylate synthase X from Helicobacter pylori strain SS1. Protein Sci 20: 1398-1410.
6. Graziani S, Xia Y, Gurnon JR, Van Etten JL, Leduc D, et al. (2004) Functional analysis of FAD-dependent thymidylate synthase ThyX from Paramecium bursaria Chlorella virus-1. J Biol Chem 279: 54340-54347.
7. Sampathkumar P, Turley S, Ulmer JE, Rhie HG, Sibley CH, et al. (2005) Structure of the Mycobacterium tuberculosis flavin dependent thymidylate synthase (MtbThyX) at 2.0A resolution. J Mol Biol 352: 1091-1104.
8. Basta T, Boum Y, Briffotaux J, Becker HF, Lamarre-Jouenne I, et al. (2012) Mechanistic and structural basis for inhibition of thymidylate synthase ThyX. Open Biol 2: 120120.
9. Esra Onen F, Boum Y, Jacquement C, Spanedda MV, Jaber N, et al. (2008) Design, synthesis and evaluation of potent thymidylate synthase X inhibitors. Bioorg Med Chem Lett 18: 3628-3631.
10. Cunningham D, Zalcberg J, Smith I, Gore M, Pazdur R, et al. (1996) 'Tomudex' (ZD1694): a novel thymidylate synthase inhibitor with clinical antitumour activity in a range of solid tumours. 'Tomudex' International Study Group. Ann Oncol 7: 179-182.
11. Carreras CW, Santi DV (1995) The catalytic mechanism and structure of thymidylate synthase. Annu Rev Biochem 64: 721-762.
12. Mishanina TV, Koehn EM, Conrad JA, Palfey BA, Lesley SA, et al. (2012) Trapping of an intermediate in the reaction catalyzed by flavin-dependent thymidylate synthase. J Am Chem Soc 134: 4442-4448.
13. Griffin J, Roshick C, Iliffe-Lee E, McClarty G (2005) Catalytic mechanism of Chlamydia trachomatis flavin-dependent thymidylate synthase. J Biol Chem 280: 5456-5467.
14. Hunter JH, Gujjar R, Pang CK, Rathod PK (2008) Kinetics and ligand-binding preferences of Mycobacterium tuberculosis thymidylate synthases, ThyA and ThyX. PLoS One 3: e2237.
15. Chernyshev A, Fleischmann T, Kohen A (2007) Thymidyl biosynthesis enzymes as antibiotic targets. Appl Microbiol Biotechnol 74: 282-289.
16. Koehn EM, Perissinotti LL, Moghram S, Prabhakar A, Lesley SA, et al. (2012) Folate binding site of flavin-dependent thymidylate synthase. Proc Natl Acad Sci U S A 109: 15722-15727.
17. Leduc D, Graziani S, Lipowski G, Marchand C, Le Maréchal P, et al. (2004) Functional evidence for active site location of tetrameric thymidylate synthase X at the interphase of three monomers. Proc Natl Acad Sci U S A 101: 7252-7257.
18. FDTS from Corynebacterium glutamicum-PDB code: 3FNN.
19. Nishimasu H, Ishitani R, Yamashita K, Iwashita C, Hirata A, et al. (2009) Atomic structure of a folate/FAD-dependent tRNA T54 methyltransferase. Proc Natl Acad Sci U S A 106: 8180-8185.
20. Park HW, Kim ST, Sancar A, Deisenhofer J (1995) Crystal structure of DNA photolyase from Escherichia coli. Science 268: 1866-1872.
21. Sridhar Prasad G, Kresge N, Muhlberg AB, Shaw A, Jung YS, et al. (1998) The crystal structure of NADPH:ferredoxin reductase from Azotobacter vinelandii. Protein Sci 7: 2541-2549.
22. Fraaije MW, van den Heuvel RH, van Berkel WJ, Mattevi A (1999) Covalent flavinylation is essential for efficient redox catalysis in vanillyl-alcohol oxidase. J Biol Chem 274: 35514-35520.
23. Cunane LM, Chen ZW, Shamala N, Mathews FS, Cronin CN, et al. (2000) Structures of the flavocytochrome p-cresol methylhydroxylase and its enzyme-substrate complex: gated substrate entry and proton relays support the proposed catalytic mechanism. J Mol Biol 295: 357-374.
24. Amaro RE, Swift RV, Votapka L, Li WW, Walker RC, et al. (2011) Mechanism of 150-cavity formation in influenza neuraminidase. Nat Commun 2: 388.
25. Rudrawar S, Dyason JC, Rameix-Welti MA, Rose FJ, Kerry PS, et al. (2010) Novel sialic acid derivatives lock open the 150-loop of an influenza A virus group-1 sialidase. Nat Commun 1: 113.
26. Chernyshev A, Fleischmann T, Kohen A (2007) Thymidyl biosynthesis enzymes as antibiotic targets. Appl Microbiol Biotechnol 74: 282-289.
27. Lesley SA, Kuhn P, Godzik A, Deacon AM, Mathews I, et al. (2002) Structural genomics of the Thermotoga maritima proteome implemented in a high-throughput structure determination pipeline. Proc Natl Acad Sci U S A 99: 11664-11669.
28. Kabsch W (2010) XDS. Acta Crystallogr D Biol Crystallogr 66: 125-132.
29. Vagin A, Teplyakov A (1997) MOLREP: an automated program for molecular replacement. J Appl Cryst 30: 1022-1025.
30. Emsley P, Cowtan K (2004) Coot: model-building tools for molecular graphics. Acta Crystallogr D Biol Crystallogr 60: 2126-2132.
31. Murshudov GN, Vagin AA, Dodson EJ (1997) Refinement of macromolecular structures by the maximum-likelihood method. Acta Crystallogr D Biol Crystallogr 53: 240-255.
32. The PyMOL Molecular Graphics System, Version 1.5.0.4 Schrödinger, LLC.

--