Ensemble-based virtual screening of African natural products to target human thymidylate synthase

Human thymidylate synthase (hTS) is a validated drug target for chemotherapy. A virtual screening experiment was used to prioritize a list of compounds from African Natural Products Databases docked against the orthosteric binding pocket of hTS. Consensus scores of binding affinities from ensemble-based virtual screening, hydrated docking and MM-PBSA calculations ranked compounds NEA4433 and NEA4434 as the best candidates owing to binding affinity scores in the picomolar order, their excellent ADMET profiles and the good stability of the protein-ligand complexes formed. The current study demonstrates the role of water in small molecule binding to hTS in mediating protein-ligand interactions. Similarly, the robust ensemble docking (relaxed scheme complex) ranked NEA4433 and NEA4434 as the best candidates. Furthermore, the best candidates prioritized were shown to strongly interact with the same residues that interacted with hTS substrate and cofactor.

Investigation of the Anticancer Effect of α-Aminophosphonates and Arylidine Derivatives of 3-Acetyl-1-aminoquinolin-2(1H)-one on the DMBA Model of Breast Cancer in Albino Rats with In Silico Prediction of Their Thymidylate Synthase Inhibitory Effect

Breast cancer is a major cause of death in women worldwide. In this study, 60 female rats were classified into 6 groups; negative control, α-aminophosphonates, arylidine derivatives of 3-acetyl-1-aminoquinolin-2(1H)-one, DMBA, DMBA & α-aminophosphonates, and DMBA & arylidine derivatives of 3-acetyl-1-aminoquinolin-2(1H)-one. New α-aminophosphonates and arylidine derivatives of 3-acetyl-1-aminoquinolin-2(1H)-one were synthesized and elucidated by different spectroscopic and elemental analysis. Histopathological examination showed marked proliferation of cancer cells in the DMBA group. Treatment with α-aminophosphonates mainly decreased tumor mass. Bcl2 expression increased in DMBA-administered rats and then declined in the treated groups, mostly with α-aminophosphonates. The level of CA15-3 markedly declined in DMBA groups treated with α-aminophosphonates and arylidine derivatives of 3-acetyl-1-aminoquinolin-2(1H)-one. Gene expression of GST-P, PCNA, PDK, and PIK3CA decreased in the DMBA group treated with α-aminophosphonates and arylidine derivatives of 3-acetyl-1-aminoquinolin-2(1H)-one, whereas PIK3R1 and BAX increased in the DMBA group treated with α-aminophosphonates and arylidine derivatives of 3-acetyl-1-aminoquinolin-2(1H)-one. The molecular docking postulated that the investigated compounds can inhibt the Thymidylate synthase TM due to high hydrophobicity charachter.

Caught in Action: X-ray Structure of Thymidylate Synthase with Noncovalent Intermediate Analog

Methylation of 2-deoxyuridine-5'-monophosphate (dUMP) at the C5 position by the obligate dimeric thymidylate synthase (TSase) in the sole de novo biosynthetic pathway to thymidine 5'-monophosphate (dTMP) proceeds by forming a covalent ternary complex with dUMP and cosubstrate 5,10-methylenetetrahydrofolate. The crystal structure of an analog of this intermediate gives important mechanistic insights but does not explain the half-of-the-sites activity of the enzyme. Recent experiments showed that the C5 proton and the catalytic Cys are eliminated in a concerted manner from the covalent ternary complex to produce a noncovalent bisubstrate intermediate. Here, we report the crystal structure of TSase with a close synthetic analog of this intermediate in which it has partially reacted with the enzyme but in only one protomer, consistent with the half-of-the-sites activity of this enzyme. Quantum mechanics/molecular mechanics simulations confirmed that the analog could undergo catalysis. The crystal structure shows a new water 2.9 Å from the critical C5 of the dUMP moiety, which in conjunction with other residues in the network, may be the elusive general base that abstracts the C5 proton of dUMP during the reaction.

Advanced Spectroscopy and APBS Modeling for Determination of the Role of His190 and Trp103 in Mouse Thymidylate Synthase Interaction with Selected dUMP Analogues

A homo-dimeric enzyme, thymidylate synthase (TS), has been a long-standing molecular target in chemotherapy. To further elucidate properties and interactions with ligands of wild-type mouse thymidylate synthase (mTS) and its two single mutants, H190A and W103G, spectroscopic and theoretical investigations have been employed. In these mutants, histidine at position 190 and tryptophan at position 103 are substituted with alanine and glycine, respectively. Several emission-based spectroscopy methods used in the paper demonstrate an especially important role for Trp 103 in TS ligands binding. In addition, the Advanced Poisson-Boltzmann Solver (APBS) results show considerable differences in the distribution of electrostatic potential around Trp 103, as compared to distributions observed for all remaining Trp residues in the mTS family of structures. Together, spectroscopic and APBS results reveal a possible interplay between Trp 103 and His190, which contributes to a reduction in enzymatic activity in the case of H190A mutation. Comparison of electrostatic potential for mTS complexes, and their mutants, with the substrate, dUMP, and inhibitors, FdUMP and N4-OH-dCMP, suggests its weaker influence on the enzyme-ligand interactions in N4OH-dCMP-mTS compared to dUMP-mTS and FdUMP-mTS complexes. This difference may be crucial for the explanation of the "abortive reaction" inhibitory mechanism of N4OH-dCMP towards TS. In addition, based on structural analyses and the H190A mutant capacity to form a denaturation-resistant complex with N4-OH-dCMP in the mTHF-dependent reaction, His190 is apparently responsible for a strong preference of the enzyme active center for the anti rotamer of the imino inhibitor form.

Structure based designing of benzimidazole/benzoxazole derivatives as anti-leishmanial agents

Folates are essential biomolecules required to carry out many crucial processes in leishmania parasite. Dihydrofolate reductase-thymidylate synthase (DHFR-TS) and pteridine reductase 1 (PTR1) involved in folate biosynthesis in leishmania have been established as suitable targets for development of chemotherapy against leishmaniasis. In the present study, various computational tools such as homology modelling, pharmacophore modelling, docking, molecular dynamics and molecular mechanics have been employed to design dual DHFR-TS and PTR1 inhibitors. Two designed molecules, i.e. 2-(4-((4-nitrobenzyl)oxy)phenyl)-1H-benzo[d]imidazole and 2-(4-((2,4-dichlorobenzyl)oxy)phenyl)-1H-benzo[d]oxazolemolecules were synthesized. MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide) assay was performed to evaluate in vitro activity of molecules against promastigote form of Leishmania donovani using Miltefosine as standard. 2-(4-((4-nitrobenzyl)oxy)phenyl)-1H-benzo[d]imidazole and 2-(4-((2,4-dichlorobenzyl)oxy)phenyl)-1H-benzo[d]oxazolemolecules were found to be moderately active with showed IC₅₀ = 68 ± 2.8 µM and 57 ± 4.2 µM, respectively.

Structural and Functional Characterization of the Human Thymidylate Synthase (hTS) Interface Variant R175C, New Perspectives for the Development of hTS Inhibitors

Human thymidylate synthase (hTS) is pivotal for cell survival and proliferation, indeed it provides the only synthetic source of dTMP, required for DNA biosynthesis. hTS represents a validated target for anticancer chemotherapy. However, active site-targeting drugs towards hTS have limitations connected to the onset of resistance. Thus, new strategies have to be applied to effectively target hTS without inducing resistance in cancer cells. Here, we report the generation and the functional and structural characterization of a new hTS interface variant in which Arg175 is replaced by a cysteine. Arg175 is located at the interface of the hTS obligate homodimer and protrudes inside the active site of the partner subunit, in which it provides a fundamental contribution for substrate binding. Indeed, the R175C variant results catalytically inactive. The introduction of a cysteine at the dimer interface is functional for development of new hTS inhibitors through innovative strategies, such as the tethering approach. Structural analysis, performed through X-ray crystallography, has revealed that a cofactor derivative is entrapped inside the catalytic cavity of the hTS R175C variant. The peculiar binding mode of the cofactor analogue suggests new clues exploitable for the design of new hTS inhibitors.

Evidence of Destabilization of the Human Thymidylate Synthase (hTS) Dimeric Structure Induced by the Interface Mutation Q62R

In human cells, thymidylate synthase (TS) provides the only source of 2′-deoxythymidyne-5′-monophosphate (dTMP), which is required for DNA biosynthesis. Because of its pivotal role, human TS (hTS) represents a validated target for anticancer chemotherapy. Nonetheless, the efficacy of drugs blocking the hTS active site has limitations due to the onset of resistance in cancer cells, requiring the identification of new strategies to effectively inhibit this enzyme. Human TS works as an obligate homodimer, making the inter-subunit interface an attractive targetable area. Here, we report the design and investigation of a new hTS variant, in which Gln62, located at the dimer interface, has been replaced by arginine in order to destabilize the enzyme quaternary assembly. The hTS Q62R variant has been characterized though kinetic assay, thermal denaturation analysis and X-ray crystallography. Our results provide evidence that hTS Q62R has a reduced melting temperature. The effective destabilization of the TS quaternary structure is also confirmed by structural analysis, showing that the introduced mutation induces a slight aperture of the hTS dimer. The generation of hTS variants having a more accessible interface area can facilitate the screening of interface-targeting molecules, providing key information for the rational design of innovative hTS interface inhibitors.

Human Thymidylate Synthase Inhibitors Halting Ovarian Cancer Growth

Human thymidylate synthase (hTS) has an important role in DNA biosynthesis, thus it is essential for cell survival. TS is involved in the folate pathways, specifically in the de novo pyrimidine biosynthesis. Structure and functions are intimately correlated, account for cellular activity and, in a broader view, with in vivo mechanisms. hTS is a target for anticancer agents, some of which are clinical drugs. The understanding of the detailed mechanism of TS inhibition by currently used drugs and of the interaction with the mechanism of action of other anticancer agents can suggest new perspective of TS inhibition able to improve the anticancer effect and to overcome drug resistance. TS-targeting drugs in therapy today are inhibitors that bind at the active site and that mostly resemble the substrates. Nonsubstrate analogs offer an opportunity for allosteric binding and novel mode of inhibition in the cancer cells. This chapter illustrates the relationship among the large number of hTS actions at molecular and clinical levels, its role as a target for ovarian cancer therapy, in particular in cases of overexpression of hTS and other folate proteins such as those induced by platinum drug treatments, and address the potential combination of TS inhibitors with other suitable anticancer agents.

Structural analyses of human thymidylate synthase reveal a site that may control conformational switching between active and inactive states

Thymidylate synthase (TS) is the sole enzyme responsible for de novo biosynthesis of thymidylate (TMP) and is essential for cell proliferation and survival. Inhibition of human TS (hTS) has been extensively investigated for cancer chemotherapy, but several aspects of its activity and regulation are still uncertain. In this study, we performed comprehensive structural and biophysical studies of hTS using crystallography and thermal shift assay and provided the first detailed structural information on the conformational changes induced by ligand binding to the hTS active site. We found that upon binding of the antifolate agents raltitrexed and nolatrexed, the two insert regions in hTS, the functions of which are unclear, undergo positional shifts toward the catalytic center. We investigated the inactive conformation of hTS and found that the two insert regions are also involved in the conformational transition between the active and inactive state of hTS. Moreover, we identified a ligand-binding site in the dimer interface, suggesting that the cavity in the dimer interface could serve as an allosteric site of hTS to regulate the conformational switching between the active and inactive states. On the basis of these findings, we propose a regulatory mechanism of hTS activity that involves allosteric regulation of interactions of hTS with its own mRNA depending on cellular demands for TMP.

DNA and RNA Pyrimidine Nucleobase Alkylation at the Carbon-5 Position

The carbon 5 of pyrimidine nucleobases is a privileged position in terms of nucleoside modification in both DNA and RNA. The simplest modification of uridine at this position is methylation leading to thymine. Thymine is an integral part of the standard nucleobase repertoire of DNA that is synthesized at the nucleotide level. However, it also occurs in RNA, where it is synthesized posttranscriptionally at the polynucleotide level. The cytidine analogue 5-methylcytidine also occurs in both DNA and RNA, but is introduced at the polynucleotide level in both cases. The same applies to a plethora of additional derivatives found in nature, resulting either from a direct modification of the 5-position by electrophiles or by further derivatization of the 5-methylpyrimidines. Here, we review the structural diversity of these modified bases, the variety of cofactors that serve as carbon donors, and the common principles shared by enzymatic mechanisms generating them.

Understanding the molecular mechanism of substrate channeling and domain communication in protozoal bifunctional TS-DHFR

Most species, such as humans, have monofunctional forms of thymidylate synthase (TS) and dihydrofolate reductase (DHFR) that are key folate metabolism enzymes making critical folate components required for DNA synthesis. In contrast, several parasitic protozoa, including Leishmania major (Lm), Plasmodium falciparum (Pf), Toxoplasma gondii (Tg) and Cryptosporidium hominis (Ch), contain a unique bifunctional thymidylate synthase-dihydrofolate reductase (TS-DHFR) having the two sequential catalytic activities contained on a single polypeptide chain. It has been suggested that the bifunctional nature of the two catalytic activities may enable substrate channeling. The 3D structures for each of these enzymes reveals distinct features for each species. While three of the four species (Pf, Tg and Ch) contain a junctional region linking the two domains, this is lacking in Lm. The Lm and Pf contain N-terminal amino acid extensions. A multidisciplinary approach using structural studies and transient kinetic analyses combined with mutational analysis has investigated the roles of these unique structural features for each enzyme. Additionally, the possibility of substrate channeling behavior was explored. These studies have identified unique, functional regions in both the TS and DHFR domains that govern efficient catalysis for each species. Surprisingly, even though there are structural similarities among the species, each is regulated in a distinct manner. This structural and mechanistic information was also used to exploit species-specific inhibitor design.

Flavin-Dependent Thymidylate Synthase as a New Antibiotic Target

In humans de novo synthesis of 2'-deoxythymidine-5'-monophosphate (dTMP), an essential building block of DNA, utilizes an enzymatic pathway requiring thymidylate synthase (TSase) and dihydrofolate reductase (DHFR). The enzyme flavin-dependent thymidylate synthase (FDTS) represents an alternative enzymatic pathway to synthesize dTMP, which is not present in human cells. A number of pathogenic bacteria, however, depend on this enzyme in lieu of or in conjunction with the analogous human pathway. Thus, inhibitors of this enzyme may serve as antibiotics. Here, we review the similarities and differences of FDTS vs. TSase including aspects of their structure and chemical mechanism. In addition, we review current progress in the search for inhibitors of flavin dependent thymidylate synthase as potential novel therapeutics.

The general base in the thymidylate synthase catalyzed proton abstraction

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

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.

Role of long-range protein dynamics in different thymidylate synthase catalyzed reactions

Recent studies of Escherichia coli thymidylate synthase (ecTSase) showed that a highly conserved residue, Y209, that is located 8 Å away from the reaction site, plays a key role in the protein’s dynamics. Those crystallographic studies indicated that Y209W mutant is a structurally identical but dynamically altered relative to the wild type (WT) enzyme, and that its turnover catalytic rate governed by a slow hydride-transfer has been affected. The most challenging test of an examination of a fast chemical conversion that precedes the rate-limiting step has been achieved here. The physical nature of both fast and slow C-H bond activations have been compared between the WT and mutant by means of observed and intrinsic kinetic isotope effects (KIEs) and their temperature dependence. The findings indicate that the proton abstraction step has not been altered as much as the hydride transfer step. Additionally, the comparison indicated that other kinetic steps in the TSase catalyzed reaction were substantially affected, including the order of the substrate binding. Enigmatically, although Y209 is H-bonded to 3'-OH of 2'-deoxyuridine-5'-mono­phosphate (dUMP), its altered dynamics is more pronounced on the binding of the remote cofactor, (6R)-N⁵,N¹⁰-methylene-5,6,7,8-tetrahydrofolate (CH₂H₄folate), revealing the importance of long-range dynamics of the enzymatic complex and its catalytic function.

High-level QM/MM calculations support the concerted mechanism for Michael addition and covalent complex formation in thymidylate synthase

Thymidylate synthase (TS) is a promising cancer target, due to its crucial function in thymine synthesis. It performs the reductive methylation of 2'-deoxyuridine-5'-phosphate (dUMP) to thymidine-5'-phosphate (dTMP), using N-5,10-methylene-5,6,7,8-tetrahydrofolate (mTHF) as a cofactor. After the formation of the dUMP/mTHF/TS noncovalent complex, and subsequent conformational activation, this complex has been proposed to react via nucleophilic attack (Michael addition) by Cys146, followed by methylene-bridge formation to generate the ternary covalent intermediate. Herein, QM/MM (B3LYP-D/6-31+G(d)-CHARMM27) methods are used to model the formation of the ternary covalent intermediate. A two-dimensional potential energy surface reveals that the methylene-bridged intermediate is formed via a concerted mechanism, as indicated by a single transition state on the minimum energy pathway and the absence of a stable enolate intermediate. A range of different QM methods (B3LYP, MP2 and SCS-MP2, and different basis sets) are tested for the calculation of the activation energy barrier for the formation of the methylene-bridged intermediate. We test convergence of the QM/MM results with respect to size of the QM region. Inclusion of Arg166, which interacts with the nucleophilic thiolate, in the QM region is important for reliable results; the MM model apparently does not reproduce energies for distortion of the guanidinium side chain correctly. The spin component scaled-Møller-Plessett perturbation theory (SCS-MP2) approach was shown to be in best agreement (within 1.1 kcal/mol) while the results obtained with MP2 and B3LYP also yielded acceptable values (deviating by less than 3 kcal/mol) compared with the barrier derived from experiment. Our results indicate that using a dispersion-corrected DFT method, or a QM method with an accurate treatment of electron correlation, increases the agreement between the calculated and experimental activation energy barriers, compared with the semiempirical AM1 method. These calculations provide important insight into the reaction mechanism of TS and may be useful in the design of new TS inhibitors.

A simple extension to the CMASA method for the prediction of catalytic residues in the presence of single point mutations

The automatic identification of catalytic residues still remains an important challenge in structural bioinformatics. Sequence-based methods are good alternatives when the query shares a high percentage of identity with a well-annotated enzyme. However, when the homology is not apparent, which occurs with many structures from the structural genome initiative, structural information should be exploited. A local structural comparison is preferred to a global structural comparison when predicting functional residues. CMASA is a recently proposed method for predicting catalytic residues based on a local structure comparison. The method achieves high accuracy and a high value for the Matthews correlation coefficient. However, point substitutions or a lack of relevant data strongly affect the performance of the method. In the present study, we propose a simple extension to the CMASA method to overcome this difficulty. Extensive computational experiments are shown as proof of concept instances, as well as for a few real cases. The results show that the extension performs well when the catalytic site contains mutated residues or when some residues are missing. The proposed modification could correctly predict the catalytic residues of a mutant thymidylate synthase, 1EVF. It also successfully predicted the catalytic residues for 3HRC despite the lack of information for a relevant side chain atom in the PDB file.

Backbone and ILV methyl resonance assignments of E. coli thymidylate synthase bound to cofactor and a nucleotide analogue

Thymidylate synthase (TSase) is a 62 kDa homodimeric enzyme required for de novo synthesis of thymidine monophosphate in most organisms. This makes the enzyme an excellent target for anticancer and microbial antibiotic drugs. In addition, TSase has been shown to exhibit negative cooperativity and half-the-sites reactivity. For these collective reasons, TSase is widely studied, and much is known about its kinetics and structure as it progresses through a multi-step catalytic cycle. Recently, nuclear magnetic resonance spin relaxation has been instrumental in demonstrating the critical role of dynamics in enzyme function in small model systems. These studies raise questions about how dynamics affect function in larger enzymes with more complex reaction coordinates. TSase is an ideal candidate given its size, oligomeric state, cooperativity, and status as a drug target. Here, as a pre-requisite to spin relaxation studies, we present the backbone and ILV methyl resonance assignments of TSase from Escherichia coli bound to a substrate analogue and cofactor.

Ultrafast real-time visualization of active site flexibility of flavoenzyme thymidylate synthase ThyX

In many bacteria the flavoenzyme thymidylate synthase ThyX produces the DNA nucleotide deoxythymidine monophosphate from dUMP, using methylenetetrahydrofolate as carbon donor and NADPH as hydride donor. Because all three substrates bind in close proximity to the catalytic flavin adenine dinucleotide group, substantial flexibility of the ThyX active site has been hypothesized. Using femtosecond time-resolved fluorescence spectroscopy, we have studied the conformational heterogeneity and the conformational interconversion dynamics in real time in ThyX from the hyperthermophilic bacterium Thermotoga maritima. The dynamics of electron transfer to excited flavin adenine dinucleotide from a neighboring tyrosine residue are used as a sensitive probe of the functional dynamics of the active site. The fluorescence decay spanned a full three orders of magnitude, demonstrating a very wide range of conformations. In particular, at physiological temperatures, multiple angstrom cofactor-residue displacements occur on the picoseconds timescale. These experimental findings are supported by molecular dynamics simulations. Binding of the dUMP substrate abolishes this flexibility and stabilizes the active site in a configuration where dUMP closely interacts with the flavin cofactor and very efficiently quenches fluorescence itself. Our results indicate a dynamic selected-fit mechanism where binding of the first substrate dUMP at high temperature stabilizes the enzyme in a configuration favorable for interaction with the second substrate NADPH, and more generally have important implications for the role of active site flexibility in enzymes interacting with multiple poly-atom substrates and products. Moreover, our data provide the basis for exploring the effect of inhibitor molecules on the active site dynamics of ThyX and other multisubstrate flavoenzymes.

Mg2+ binds to the surface of thymidylate synthase and affects hydride transfer at the interior active site

Thymidylate synthase (TSase) produces the sole intracellular de novo source of thymidine (i.e., the DNA base T) and thus is a common target for antibiotic and anticancer drugs. Mg(2+) has been reported to affect TSase activity, but the mechanism of this interaction has not been investigated. Here we show that Mg(2+) binds to the surface of Escherichia coli TSase and affects the kinetics of hydride transfer at the interior active site (16 Å away). Examination of the crystal structures identifies a Mg(2+) near the glutamyl moiety of the folate cofactor, providing the first structural evidence for Mg(2+) binding to TSase. The kinetics and NMR relaxation experiments suggest that the weak binding of Mg(2+) to the protein surface stabilizes the closed conformation of the ternary enzyme complex and reduces the entropy of activation on the hydride transfer step. Mg(2+) accelerates the hydride transfer by ~7-fold but does not affect the magnitude or temperature dependence of the intrinsic kinetic isotope effect. These results suggest that Mg(2+) facilitates the protein motions that bring the hydride donor and acceptor together, but it does not change the tunneling ready state of the hydride transfer. These findings highlight how variations in cellular Mg(2+) concentration can modulate enzyme activity through long-range interactions in the protein, rather than binding at the active site. The interaction of Mg(2+) with the glutamyl tail of the folate cofactor and nonconserved residues of bacterial TSase may assist in designing antifolates with polyglutamyl substitutes as species-specific antibiotic drugs.

An innovative strategy for dual inhibitor design and its application in dual inhibition of human thymidylate synthase and dihydrofolate reductase enzymes

Due to the diligence of inherent redundancy and robustness in many biological networks and pathways, multitarget inhibitors present a new prospect in the pharmaceutical industry for treatment of complex diseases. Nevertheless, to design multitarget inhibitors is concurrently a great challenge for medicinal chemists. We have developed a novel computational approach by integrating the affinity predictions from structure-based virtual screening with dual ligand-based pharmacophore to discover potential dual inhibitors of human Thymidylate synthase (hTS) and human dihydrofolate reductase (hDHFR). These are the key enzymes in folate metabolic pathway that is necessary for the biosynthesis of RNA, DNA, and protein. Their inhibition has found clinical utility as antitumor, antimicrobial, and antiprotozoal agents. A druglike database was utilized to perform dual-target docking studies. Hits identified through docking experiments were mapped over a dual pharmacophore which was developed from experimentally known dual inhibitors of hTS and hDHFR. Pharmacophore mapping procedure helped us in eliminating the compounds which do not possess basic chemical features necessary for dual inhibition. Finally, three structurally diverse hit compounds that showed key interactions at both active sites, mapped well upon the dual pharmacophore, and exhibited lowest binding energies were regarded as possible dual inhibitors of hTS and hDHFR. Furthermore, optimization studies were performed for final dual hit compound and eight optimized dual hits demonstrating excellent binding features at target systems were also regarded as possible dual inhibitors of hTS and hDHFR. In general, the strategy used in the current study could be a promising computational approach and may be generally applicable to other dual target drug designs.

Crystal structure and enzymatic characterization of thymidylate synthase X from Helicobacter pylori strain SS1

Thymidylate synthase X (ThyX) catalyzes the methylation of dUMP to form dTMP in bacterial life cycle and is regarded as a promising target for antibiotics discovery. Helicobacter pylori is a human pathogen associated with a number of human diseases. Here, we cloned and purified the ThyX enzyme from H. pylori SS1 strain (HpThyX). The recombinant HpThyX was discovered to exhibit the maximum activity at pH 8.5, and K(m) values of the two substrates dUMP and CH(2) H(4) folate were determined to be 15.3 ± 1.25 μM and 0.35 ± 0.18 mM, respectively. The analyzed crystal structure of HpThyX with the cofactor FAD and the substrate dUMP (at 2.31 Å) revealed that the enzyme was a tetramer bound to four dUMP and four FAD molecules. Different from the catalytic feature of the classical thymidylate synthase (ThyA), N5 atom of the FAD functioned as a nucleophile in the catalytic reaction instead of Ser84 and Ser85 residues. Our current work is expected to help better understand the structural and enzymatic features of HpThyX thus further providing valuable information for anti-H. pylori inhibitor discovery.

Temperature dependence of the kinetic isotope effects in thymidylate synthase. A theoretical study

In recent years, the temperature dependence of primary kinetic isotope effects (KIE) has been used as indicator for the physical nature of enzyme-catalyzed H-transfer reactions. An interactive study where experimental data and calculations examine the same chemical transformation is a critical means to interpret more properly temperature dependence of KIEs. Here, the rate-limiting step of the thymidylate synthase-catalyzed reaction has been studied by means of hybrid quantum mechanics/molecular mechanics (QM/MM) simulations in the theoretical framework of the ensemble-averaged variational transition-state theory with multidimensional tunneling (EA-VTST/MT) combined with Grote-Hynes theory. The KIEs were calculated across the same temperature range examined experimentally, revealing a temperature independent behavior, in agreement with experimental findings. The calculations show that the H-transfer proceeds with ∼91% by tunneling in the case of protium and ∼80% when the transferred protium is replaced by tritium. Dynamic recrossing coefficients are almost invariant with temperature and in all cases far from unity, showing significant coupling between protein motions and the reaction coordinate. In particular, the relative movement of a conserved arginine (Arg166 in Escherichia coli ) promotes the departure of a conserved cysteine (Cys146 in E. coli ) from the dUMP by polarizing the thioether bond thus facilitating this bond breaking that takes place concomitantly with the hydride transfer. These promoting vibrations of the enzyme, which represent some of the dimensions of the real reaction coordinate, would limit the search through configurational space to efficiently find those decreasing both barrier height and width, thereby enhancing the probability of H-transfer by either tunneling (through barrier) or classical (over-the-barrier) mechanisms. In other words, the thermal fluctuations that are coupled to the reaction coordinate, together with transition-state geometries and tunneling, are the same in different bath temperatures (within the limited experimental range examined). All these terms contribute to the observed temperature independent KIEs in thymidylate synthase.

Synthesis, antitumor activity and molecular docking study of novel benzofuran-2-yl pyrazole pyrimidine derivatives

A new series of (benzofuran-2-yl)-1-phenyl-1H-pyrazol-4-yl) pyrimidine derivatives were synthesized from 3-(benzofuran-2-yl)-1-phenyl-1H-pyrazol-4-carbaldehyde (1) through different routes of cyclocondensation reactions. Condensation of 1 with active methylene compounds afforded compounds 2-8. The cyclization of 2 with chloroacetic acid, ortho substituted benzoic acid and/or ethanolamine gave compounds 9-12. Also condensation of 2 with hydrazine hydrate followed by cyclocondensation afforded corresponding triazines and pyrazole derivatives 18-27. Some docking studies of the newly prepared compounds as thymidylate synthase inhibitors have been done. Also the cytotoxic activity of some of the prepared compounds as a representative examples was evaluated against HEPG2 (human liver carcinoma cell line) in comparison with 5-fluorouracil (5-Fu).

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.

Crystallization and preliminary crystallographic studies of a flavin-dependent thymidylate synthase from Helicobacter pylori

The ThyX enzymes that have recently been identified in various bacteria, including some important human pathogens such as Helicobacter pylori and Mycobacterium tuberculosis, are flavin-dependent thymidylate synthases that function in the place of classic thymidylate synthase enzymes in the biosynthesis of dTMP, which is one of the building blocks of DNA. They are promising targets for the development of novel antibiotics because they utilize catalytic mechanisms that are distinct from those of the classic thymidylate synthases found in most organisms, including humans. In this study, H. pylori ThyX was purified and crystallized in complex with flavin adenine dinucleotide (FAD) and a diffraction data set was collected to 2.5 A resolution. The crystals belonged to space group C2, with unit-cell parameters a = 221.92, b = 49.43, c = 143.02 A, beta = 98.84 degrees .

Synthesis, biological evaluation and molecular docking of aryl hydrazines and hydrazides for anticancer activity

Aryl hydrazine and hydrazide analogues were synthesized based on p-tolyl hydrazine, isolated as a breakdown product of a secondary metabolite from the mushroom, Agaricus bisporus, and tested to be highly active molecule than 5-fluorouracil in in vitro anticancer studies. The synthesized analogues were tested for anticancer activity using NCI protocol. Anolgues 12 and 15 emerged as molecules with significant in vitro anticancer activity. Molecular docking study revealed the binding orientations of aryl hydrazines and hydrazides analogues in the active sites of thymidylate synthase.

Phosphorylation of thymidylate synthase from various sources by human protein kinase CK2 and its catalytic subunits

Thymidylate synthase (TS) was found to be a substrate for both catalytic subunits of human CK2, with phosphorylation by CK2alpha and CK2alpha' characterized by similar K(m) values, 4.6microM and 4.2microM, respectively, but different efficiencies, the apparent turnover number with CK2alpha being 10-fold higher. With both catalytic subunits, phosphorylation of human TS, like calmodulin and BID, was strongly inhibited in the presence of the regulatory subunit CK2beta, the holoenzyme being activated by polylysine. Phosphorylation of recombinant human, rat, mouse and Trichinella spiralis TSs proteins was compared, with the human enzyme being apparently a much better substrate than the others. Following hydrolysis and TLC, phosphoserine was detected in human and rat, and phosphotyrosine in T. spiralis, TS, used as substrates for CK2alpha. MALDI-TOF MS analysis led to identification of phosphorylated Ser(124) in human TS, within a sequence LGFS(124)TREEGD, atypical for a CK2 substrate recognition site. The phosphorylation site is located in a region considered important for the catalytic mechanism or regulation of human TS, corresponding to the loop 107-128. Following phosphorylation by CK2alpha, resulting in incorporation of 0.4mol of phosphate per mol of dimeric TS, human TS exhibits unaltered K(m) values for dUMP and N(5,10)-methylenetetrahydrofolate, but a 50% lower turnover number, pointing to a strong influence of Ser(124) phosphorylation on its catalytic efficiency.

Multiparameter screening on SlipChip used for nanoliter protein crystallization combining free interface diffusion and microbatch methods

This paper describes two SlipChip-based approaches to protein crystallization: a SlipChip-based free interface diffusion (FID) method and a SlipChip-based composite method that simultaneously performs microbatch and FID crystallization methods in a single device. The FID SlipChip was designed to screen multiple reagents, each at multiple diffusion equilibration times, and was validated by screening conditions for crystallization of two proteins, enoyl-CoA hydratase from Mycobacterium tuberculosis and dihydrofolate reductase/thymidylate synthase from Babesia bovis, against 48 different reagents at five different equilibration times each, consuming 12 microL of each protein for a total of 480 experiments using three SlipChips. The composite SlipChip was designed to screen multiple reagents, each at multiple mixing ratios and multiple equilibration times, and was validated by screening conditions for crystallization of two proteins, enoyl-CoA hydratase from Mycobacterium tuberculosis and dihydrofolate reductase/thymidylate synthase from Babesia bovis. To prevent cross-contamination while keeping the solution in the neck channels for FID stable, the plates of the SlipChip were etched with a pattern of nanowells. This nanopattern was used to increase the contact angle of aqueous solutions on the surface of the silanized glass. The composite SlipChip increased the number of successful crystallization conditions and identified more conditions for crystallization than separate FID and microbatch screenings. Crystallization experiments were scaled up in well plates using conditions identified during the SlipChip screenings, and X-ray diffraction data were obtained to yield the protein structure of dihydrofolate reductase/thymidylate synthase at 1.95 A resolution. This free-interface diffusion approach provides a convenient and high-throughput method of setting up gradients in microfluidic devices and may find additional applications in cell-based assays.

Crystallization and preliminary crystallographic studies of dihydrofolate reductase-thymidylate synthase from Trypanosoma cruzi, the Chagas disease pathogen

Trypanosoma cruzi dihydrofolate reductase-thymidylate synthase (TcDHFR-TS) was crystallized in complexes with the dihydrotriazine-based or quinazoline-based antifolates C-448, cycloguanil (CYC) and Q-8 in order to gain insight into the interactions of this DHFR enzyme with classical and novel inhibitors. The TcDHFR-TS-C-448-NDP-dUMP crystal belonged to space group C222(1) with two molecules per asymmetric unit and diffracted to 2.37 angstrom resolution. The TcDHFR-TS-CYC, TcDHFR-TS-CYC-NDP and TcDHFR-TS-Q-8-NDP crystals belonged to space group P2(1) with four molecules per asymmetric unit and diffracted to 2.1, 2.6 and 2.8 angstrom resolution, respectively. Crystals belonging to the two different space groups were suitable for structure determination.

An analysis of hydrophobic interactions of thymidylate synthase with methotrexate: free energy calculations involving mutant and native structures bound to methotrexate

Since the human body for many reasons can adapt and become resistant to drugs, it is important to develop and validate computer aided drug design (CADD) methods that could help predict binding affinity changes that can result from these resistant enzymes. The free energy perturbation (FEP) methodology is the most accurate means of estimating relative binding affinities between inhibitors and protein variants. In this paper, we describe the role played by hydrophobic residues lining the active site region, particularly (79)Ile and (176)Phe, in the binding of methotrexate to the Escherichia coli (E. coli) thymidylate synthase (TS) enzyme, using the thermodynamic cycle perturbation (TCP) approach. The computed binding free energy differences on the binding of methotrexate to the native and some mutant E. coli TS structures have been compared with experimental results. Computationally, four different 'mutations' have been simulated on the TS enzyme with methotrexate (MTX): (79)Ile --> (79)Val; (79)Ile --> (79)Ala; (79)Ile --> (79)Leu; and (176)Phe --> (176)Ile. The calculated results indicate that in each of these cases, the native residues ((79) Ile and (176) Phe) interact more favorably with methotrexate than the mutant residues and these results are corroborated by experimental measurements. Binding preference to wild type residues can be rationalized in terms of their better hydrophobic contacts with the phenyl ring of methotrexate.

Thymidylate synthase and dihydropyrimidine dehydrogenase levels are associated with response to 5-fluorouracil in Caenorhabditis elegans

5-Fluorouracil (5-FU), a pyrimidine antagonist, has a long history in cancer treatment. The targeted pyrimidine biosynthesis pathway includes dihydropyrimidine dehydrogenase (DPD), which converts 5-FU to an inactive metabolite, and thymidylate synthase (TS), which is a major target of 5-FU. Using Caenorhabditis elegans as a model system to study the functional and resistance mechanisms of anti-cancer drugs, we examined these two genes in order to determine the extent of molecular conservation between C. elegans and humans. Overexpression of the worm DPD and TS homologs (DPYD-1 and Y110A7A.4, respectively) suppressed germ cell death following 5-FU exposure. In addition, DPYD-1 depletion by RNAi resulted in 5-FU sensitivity, while treatment with Y110A7A.4 RNAi and 5-FU resulted in similar patterns of embryonic death. Thus, the pathway of 5-FU function appears to be highly conserved between C. elegans and humans at the molecular level.

Role of an invariant lysine residue in folate binding on Escherichia coli thymidylate synthase: calorimetric and crystallographic analysis of the K48Q mutant

Thymidylate synthase (TS) catalyzes the reductive methylation of deoxyuridine monophosphate (dUMP) using methylene tetrahydrofolate (CH(2)THF) as cofactor, the glutamate tail of which forms a water-mediated hydrogen bond with an invariant lysine residue of this enzyme. To understand the role of this interaction, we studied the K48Q mutant of Escherichia coli TS using structural and biophysical methods. The k(cat) of the K48Q mutant was 430-fold lower than wild-type TS in activity, while the K(m) for the (R)-stereoisomer of CH(2)THF was 300 microM, about 30-fold larger than K(m) from the wild-type TS. Affinity constants were determined using isothermal titration calorimetry, which showed that binding was reduced by one order of magnitude for folate-like TS inhibitors, such as propargyl-dideazafolate (PDDF) or compounds that distort the TS active site like BW1843U89 (U89). The crystal structure of the K48Q-dUMP complex revealed that dUMP binding is not impaired in the mutant, and that U89 in a ternary complex of K48Q-nucleotide-U89 was bound in the active site with subtle differences relative to comparable wild-type complexes. PDDF failed to form ternary complexes with K48Q and dUMP. Thermodynamic data correlated with the structural determinations, since PDDF binding was dominated by enthalpic effects while U89 had an important entropic component. In conclusion, K48 is critical for catalysis since it leads to a productive CH(2)THF binding, while mutation at this residue does not affect much the binding of inhibitors that do not make contact with this group.

Hydride transfer versus hydrogen radical transfer in thymidylate synthase

The nature of a H-transfer in the thymidylate synthase catalyzed reaction was investigated by comparison of the wild-type enzyme with the W80M mutant. The nature of the H-transfer was not affected, as indicated by intrinsic isotope effects and their temperature dependence. These findings support a single-step hydride transfer instead of a two-step radical transfer.

SAR index: quantifying the nature of structure-activity relationships

Structure-activity relationships (SARs) can display very different features. Small chemical modifications of active molecules often dramatically alter biological responses. By contrast, structurally diverse molecules can have similar activity. SARs can also be heterogeneous in nature. For example, for structurally diverse molecules with similar activity, closely related analogs might have significant differences in potency. Given the inherent complexity of SARs, it has been very difficult to estimate SAR characteristics from molecular structure. On the basis of systematic correlation of 2D structural similarity and compound potency, we have developed a function termed "SAR Index" that quantitatively describes the nature of SARs and establishes different SAR categories: continuous, discontinuous, heterogeneous-relaxed, and heterogeneous-constrained. These heterogeneous SAR categories are described for the first time. Given a set of active compounds and their potency values, SAR Index calculations can estimate how likely it is to identify structurally distinct molecules having similar activity.

Significance of mutations on the structural perturbation of thymidylate synthase: implications for their involvement in subunit exchange

Wild-type thymidylate synthase (WT-TS) from Escherichia coli and several of its mutants showed varying degrees of susceptibility to trypsin. While WT-TS was resistant to trypsin as were the mutants C146S, K48E, and R126K, others such as Y94A, Y94F, C146W, and R126E were digested but at different rates from one another. The peptides released from the mutants were identified by mass spectrometry and Edman sequence analysis. The known crystal structures for WT-TS, Y94F, and R126E, surprisingly, showed no structural differences that could explain the difference in their susceptibility to trypsin. One explanation is that the mutations could perturb the dynamic equilibrium of the dimeric state of the mutants as to increase their dissociation to monomers, which being less structured than the dimer, would be hydrolyzed more readily by trypsin. Earlier studies appear to support this proposal since conditions that promote subunit dissociation in solutions of R126E with other inactive mutants, such as dilution, low concentrations of urea, and elevated pH, greatly enhance the rate of restoration of TS activity. Analytic ultracentrifuge studies with various TSs in urea, or at pH 9.0, or that have been highly diluted are, for the most part, in agreement with this thesis, since these conditions are associated with an increase in dissociation to monomers, particularly with the mutant TSs. However, these studies do not rule out the possibility that conformation differences among the various TS dimers are responsible for the differences in susceptibility to trypsin, particularly at high concentrations of protein where the WT-TS and mutants are mainly dimers.

Interactions of 2'-fluoro-substituted dUMP analogues with thymidylate synthase

A series of 2'-fluoro-substituted dUMP/FdUMP analogues were synthesized, their interaction with human recombinant thymidylate synthase investigated, and structural (1)H and (19)F NMR study of the corresponding nucleosides performed. While 2'-F-dUMP (fluorine in the "down" configuration), in striking contrast to 2'-F-ara-UMP (fluorine in the "up" configuration) and 2',2''-diF-dUMP, showed substrate activity, 2'-F-ara-UMP and 2',2''-diF-dUMP were classic inhibitors, and 2',5-diF-ara-UMP behaved as a strong slow-binding inhibitor, suggesting the 2'-F substituent in the "up" position to interfere with the active center cysteine thiol addition to the pyrimidine C(6) and the pyrimidine C(5)-F to prevent this interference. In support, the direct through space heteronuclear coupling J(HF) was observed for the fluorine "up" derivatives, 2'-F-ara-U and 2',5-diF-ara-U, causing the splitting of the H(6) resonance lines. The absence of such splitting in 2',2''-diF-dUrd, indicating an unusual orientation of the base in relation to the furanose, was associated with an exceptionally weak interaction with the enzyme.

Nonconserved residues Ala287 and Ser290 of the Cryptosporidium hominis thymidylate synthase domain facilitate its rapid rate of catalysis

Cryptosporidium hominis TS-DHFR exhibits an unusually high rate of catalysis at the TS domain, at least 10-fold greater than those of other TS enzymes. Using site-directed mutagenesis, we have mutated residues Ala287 and Ser290 in the folate-binding helix to phenylalanine and glycine, respectively, the corresponding residues in human and most other TS enzymes. Our results show that the mutant A287F, the mutant S290G, and the double mutant all have reduced affinities for methylene tetrahydrofolate and reduced rates of reaction at the TS domain. Interestingly, the S290G mutant enzyme had the lowest TS activity, with a catalytic efficiency approximately 200-fold lower than that of the wild type (WT). The rate of conformational change of the S290G mutant is approximately 80 times slower than that of WT, resulting in a change in the rate-limiting step from hydride transfer to covalent ternary complex formation. We have determined the crystal structure of ligand-bound S290G mutant enzyme, which shows that the primary effect of the mutation is an increase in the distance between the TS ligands. The kinetic and crystal structure data presented here provide the first evidence explaining the unusually fast TS rate in C. hominis.

Cooperative inhibition of human thymidylate synthase by mixtures of active site binding and allosteric inhibitors

Thymidylate synthase (TS) is a target in the chemotherapy of colorectal cancer and some other neoplasms. It catalyzes the transfer of a methyl group from methylenetetrahydrofolate to dUMP to form dTMP. On the basis of structural considerations, we have introduced 1,3-propanediphosphonic acid (PDPA) as an allosteric inhibitor of human TS (hTS); it is proposed that PDPA acts by stabilizing an inactive conformer of loop 181-197. Kinetic studies showed that PDPA is a mixed (noncompetitive) inhibitor versus dUMP. In contrast, versus methylenetrahydrofolate at concentrations lower than 0.25 microM, PDPA is an uncompetitive inhibitor, while at PDPA concentrations higher than 1 microM the inhibiton is noncompetive, as expected. At the concentrations corresponding to uncompetitive inhibition, PDPA shows positive cooperativity with an antifolate inhibitor, ZD9331, which binds to the active conformer. PDPA binding leads to the formation of hTS tetramers, but not higher oligomers. These data are consistent with a model in which hTS exists preferably as an asymmetric dimer with one subunit in the active conformation of loop 181-197 and the other in the inactive conformation.

A quantum mechanics/molecular mechanics study of the catalytic mechanism of the thymidylate synthase

A theoretical study of the molecular mechanism of the thymidylate synthase-catalyzed reaction has been carried out using hybrid quantum mechanics/molecular mechanics methods. We have examined all of the stationary points (reactants, intermediates, transition structures, and products) on the multidimensional potential energy surfaces for the multistep enzymatic process. The characterization of these relevant structures facilitates the gaining of insight into the role of the different residues in the active site. Furthermore, analysis of the full energy profile has revealed that the step corresponding to the reduction of the exocyclic methylene intermediate by hydride transfer from the 6S position of 5,6,7,8-tetrahydrofolate (H4folate), forming dTMP and 7,8-dihydrofolate (H2folate), is the rate-limiting step, in accordance with the experimental data. In this step, the hydride transfer and the scission of an overall conserved active site cysteine residue (Cys146 in Escherichia coli) take place in a concerted but very asynchronous way. These findings have also been tested with primary and secondary deuterium, tritium, and sulfur kinetic isotope effects, and the calculations have been compared to experimental data. Finally, the incorporation of high-level quantum mechanical corrections to the semiempirical AM1 Hamiltonian into our hybrid scheme has allowed us to obtain reasonable values of the energy barrier for the rate-limiting step. The resulting picture of the complete multistep enzyme mechanism that is obtained reveals several new features of substantial mechanistic interest.

The effect of 5-substitution in the pyrimidine ring of dUMP on the interaction with thymidylate synthase: molecular modeling and QSAR

Thymidylate synthase (TS) is a target enzyme for a number of anticancer agents including the 5-fluorouracil metabolite, FdUMP. The present paper reports on molecular modeling studies of the effect of substitution at C(5) position in the pyrimidine ring of the TS substrate, dUMP, on the binding affinity for the enzyme. The results of molecular dynamics simulations show that the binding of C(5) analogues of dUMP to TS in the binary complexes does not undergo changes, unless a substituent with a large steric effect, such as the propyl group, is involved. On the other hand, apparent differences in the binding of the TS cofactor, resulting from varying substitution at dUMP C(5), are observed in the modeled structures of the ternary complexes of TS. These binding characteristics are supplemented with a classical QSAR model quantifying the relation between the affinity for TS and the substituent electronic and steric effects of C(5) analogues of dUMP. Based on the findings from the present work, the perspectives for finding promising new C(5) analogues of dUMP as potential agents targeted against TS are discussed.

Thymidyl biosynthesis enzymes as antibiotic targets

The two long-known "classical" enzymes of uridyl-5-methylation, thymidylate synthase and ribothymidyl synthase, have been joined by two alternative methylation enzymes, flavin-dependent thymidylate synthase and folate-dependent ribothymidyl synthase. These two newly discovered enzymes have much in common: both contain flavin cofactors, utilize methylenetetrahydrofolate as a source of methyl group, and perform thymidylate synthesis via chemical pathways distinct from those of their classic counterparts. Several severe human pathogens (e.g., typhus, anthrax, tuberculosis, and more) depend on these "alternative" enzymes for reproduction. These and other distinctive properties make the alternative enzymes and their corresponding genes appealing targets for new antibiotics.

Sequence analysis and functional study of thymidylate synthase from zebrafish, Danio rerio

The thymidylate synthase (TS), an important target for many anticancer drugs, has been cloned from different species. But the cDNA property and function of TS in zebrafish are not well documented. In order to use zebrafish as an animal model for screening novel anticancer agents, we isolated TS cDNA from zebrafish and compared its sequence with those from other species. The open reading frame (ORF) of zebrafish TS cDNA sequence was 954 nucleotides, encoding a 318-amino acid protein with a calculated molecular mass of 36.15 kDa. The deduced amino acid sequence of zebrafish TS was similar to those from other organisms, including rat, mouse and humans. The zebrafish TS protein was expressed in Escherichia coli and purified to homogeneity. The purified zebrafish TS showed maximal activity at 28 degrees C with similar K(m) value to human TS. Western immunoblot assay confirmed that TS was expressed in all the developmental stages of zebrafish with a high level of expression at the 1-4 cell stages. To study the function of TS in zebrafish embryo development, a short hairpin RNA (shRNA) expression vector, pSilencer 4.1-CMV/TS, was constructed which targeted the protein-coding region of zebrafish TS mRNA. Significant change in the development of tail and epiboly was found in zebrafish embryos microinjected pSilencer4.1-CMV/TS siRNA expression vector.

Antibacterial Agent Discovery Using Thymidylate Synthase Biolibrary Screening

Thymidylate synthase (TS, ThyA) catalyzes the reductive methylation of 2'-deoxyuridine 5'-monophosphate to 2'-deoxythymidine 5'-monophosphate, an essential precursor for DNA synthesis. A specific inhibition of this enzyme induces bacterial cell death. As a second round lead optimization design, new 1,2-naphthalein derivatives have been synthesized and tested against a TS-based biolibrary, including human thymidylate synthase (hTS). Docking studies have been performed to rationalize the experimentally observed affinity profiles of 1,2-naphthalein compounds toward Lactobacillus casei TS and hTS. The best TS inhibitors have been tested against a number of clinical isolates of Gram-positive-resistant bacterial strains. Compound 3,3-bis(3,5-dibromo-4-hydroxyphenyl)-1H,3H-naphtho[1,2-c]furan-1-one (5) showed significant antibacterial activity, no in vitro toxicity, and dose-response effects against Staphylococcus epidermidis (MIC=0.5-2.5 microg/mL) clinical isolate strains, which are resistant to at least 17 of the best known antibacterial agents, including vancomycin. So far this compound can be regarded as a leading antibacterial agent.

Structure of the Y94F mutant of Escherichia coli thymidylate synthase

Tyr94 of Escherichia coli thymidylate synthase is thought to be involved, either directly or by activation of a water molecule, in the abstraction of a proton from C5 of the 2'-deoxyuridine 5'-monophosphate (dUMP) substrate. Mutation of Tyr94 leads to a 400-fold loss in catalytic activity. The structure of the Y94F mutant has been determined in the native state and as a ternary complex with thymidine 5'-monophosphate (dTMP) and 10-propargyl 5,8-dideazafolate (PDDF). There are no structural changes ascribable to the mutation other than loss of a water molecule hydrogen bonded to the tyrosine OH, which is consistent with a catalytic role for the phenolic OH.

A proposed clinical test for monitoring fluoropyrimidine therapy: detection and stability of thymidylate synthase ternary complexes

5-fluorouracil forms classic (covalent, ternary) complexes consisting of thymidylate synthase, fluoro-deoxyuridine monophosphate, and 5,10-methylene tetrahydrofolate. Despite a high pharmacologic interest in the classic complexes formed in cells treated with fluorouracil anticancer agents, the in vivo stability of the complexes and the possible interference in complex formation by other coadministered compounds have not been adequately described. We visualized classic complexes unaccompanied by unbound thymidylate synthase, inferring complete enzymatic inhibition, in 5-fluorouracil-treated S. cerevisiae and cancer cells in vitro and in murine tumors in vivo treated with 5-fluorouracil. Classic complexes persisted 13 days in cancer cells after a pulse of 5-fluorouracil. Classic complexes were reduced to absent in cancer cells in which the older antifolates methotrexate and aminopterin, or the modern antifolates pemetrexed and tomudex, were coadministered with 5-fluorouracil. Classic complexes were, however, detected when an alternate drug, 5-fluorodeoxyuridine, was administered with methotrexate. We visualized classic complexes at fifteen minutes to seven days after an acute single dose of 5-fluorouracil in mouse tumor models, in tumors and normal tissues. Using the same assay, we detected unbound thymidylate synthase in untreated human tissues, supporting the future use of this assay in evaluating the most appropriate dose of fluoropyrimidine and coadministered agents in clinical settings.

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

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

Crystal structures of thymidylate synthase mutant R166Q: structural basis for the nearly complete loss of catalytic activity

Thymidylate synthase (TS) catalyzes the folate-dependent methylation of deoxyuridine monophosphate (dUMP) to form thymidine monophosphate (dTMP). We have investigated the role of invariant arginine 166, one of four arginines that contact the dUMP phosphate, using site-directed mutagenesis, X-ray crystallography, and TS from Escherichia coli. The R166Q mutant was crystallized in the presence of dUMP and a structure determined to 2.9 A resolution, but neither the ligand nor the sulfate from the crystallization buffer was found in the active site. A second structure determined with crystals prepared in the presence of dUMP and the antifolate 10-propargyl-5,8-dideazafolate revealed that the inhibitor was bound in an extended, nonproductive conformation, partially occupying the nucleotide-binding site. A sulfate ion, rather than dUMP, was found in the nucleotide phosphate-binding site. Previous studies have shown that the substitution at three of the four arginines of the dUMP phosphate-binding site is permissive; however; for Arg166, all the mutations lead to a near-inactive mutant. The present structures of TS R166Q reveal that the phosphate-binding site is largely intact, but with a substantially reduced affinity for phosphate, despite the presence of the three remaining arginines. The position of Cys146, which initiates catalysis, is shifted in the mutant and resides in a position that interferes with the binding of the dUMP pyrimidine moiety.

Catalytic mechanism and structure of viral flavin-dependent thymidylate synthase ThyX

By using biochemical and structural analyses, we have investigated the catalytic mechanism of the recently discovered flavin-dependent thymidylate synthase ThyX from Paramecium bursaria chlorella virus-1 (PBCV-1). Site-directed mutagenesis experiments have identified several residues implicated in either NADPH oxidation or deprotonation activity of PBCV-1 ThyX. Chemical modification by diethyl pyrocarbonate and mass spectroscopic analyses identified a histidine residue (His53) crucial for NADPH oxidation and located in the vicinity of the redox active N-5 atom of the FAD ring system. Moreover, we observed that the conformation of active site key residues of PBCV-1 ThyX differs from earlier reported ThyX structures, suggesting structural changes during catalysis. Steady-state kinetic analyses support a reaction mechanism where ThyX catalysis proceeds via formation of distinct ternary complexes without formation of a methyl enzyme intermediate.