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

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. *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 February 14, 2013; Accepted April 18, 2013; Published 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. J o u r na l o f B iote rrorism&iodefe n s e ISSN: 2157-2526 Journal of Bioterrorism & Biodefense 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


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 5 ,N 10 -methylene-5,6,7,8-tetrahydrofolate (CH2H4 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][2][3][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][6][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 CH2H4 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 CH2H4 folate and CH2H4 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 2 H 4 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 complex with FAD, dUMP and CH2H4 folate [16]. We also reported that both the H53A and H53D mutants showed dTMP formation with dramatically reduced activity (Table S2 of reference 17).

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Å 2 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.

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 CH2H4 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 CH2H4 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][21][22][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  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].

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 2 H 4 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].

Crystallographic parameters
Space group Unit-cell dimensions (Å)   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 substrate-binding 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.   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.

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 1 2 1 2 1 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 (Table 1). The figures were generated using PyMOL graphic program [32].

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