Catalytic cysteine of thymidylate synthase is activated upon substrate binding

L-Glyceraldehyde Inhibits Neuroblastoma Cell Growth via a Multi-Modal Mechanism on Metabolism and Signaling

Glyceraldehyde (GA) is a three-carbon monosaccharide that can be present in cells as a by-product of fructose metabolism. Bruno Mendel and Otto Warburg showed that the application of GA to cancer cells inhibits glycolysis and their growth. However, the molecular mechanism by which this occurred was not clarified. We describe a novel multi-modal mechanism by which the L-isomer of GA (L-GA) inhibits neuroblastoma cell growth. L-GA induces significant changes in the metabolic profile, promotes oxidative stress and hinders nucleotide biosynthesis. GC-MS and 13C-labeling was employed to measure the flow of carbon through glycolytic intermediates under L-GA treatment. It was found that L-GA is a potent inhibitor of glycolysis due to its proposed targeting of NAD(H)-dependent reactions. This results in growth inhibition, apoptosis and a redox crisis in neuroblastoma cells. It was confirmed that the redox mechanisms were modulated via L-GA by proteomic analysis. Analysis of nucleotide pools in L-GA-treated cells depicted a previously unreported observation, in which nucleotide biosynthesis is significantly inhibited. The inhibitory action of L-GA was partially relieved with the co-application of the antioxidant N-acetyl-cysteine. We present novel evidence for a simple sugar that inhibits cancer cell proliferation via dysregulating its fragile homeostatic environment.

Covalent Modification of the Flavin in Proline Dehydrogenase by Thiazolidine-2-Carboxylate

Proline dehydrogenase (PRODH) catalyzes the first step of proline catabolism, the FAD-dependent 2-electron oxidation of l-proline to Δ¹-pyrroline-5-carboxylate. PRODH has emerged as a possible cancer therapy target, and thus the inhibition of PRODH is of interest. Here we show that the proline analogue thiazolidine-2-carboxylate (T2C) is a mechanism-based inactivator of PRODH. Structures of the bifunctional proline catabolic enzyme proline utilization A (PutA) determined from crystals grown in the presence of T2C feature strong electron density for a 5-membered ring species resembling l-T2C covalently bound to the N5 of the FAD in the PRODH domain. The modified FAD exhibits a large butterfly bend angle, indicating that the FAD is locked into the 2-electron reduced state. Reduction of the FAD is consistent with the crystals lacking the distinctive yellow color of the oxidized enzyme and stopped-flow kinetic data showing that T2C is a substrate for the PRODH domain of PutA. A mechanism is proposed in which PRODH catalyzes the oxidation of T2C at the C atom adjacent to the S atom of the thiazolidine ring (C5). Then, the N5 atom of the reduced FAD attacks the C5 of the oxidized T2C species, resulting in the covalent adduct observed in the crystal structure. To our knowledge, this is the first report of T2C inactivating (or inhibiting) PRODH or any other flavoenzyme. These results may inform the design of new mechanism-based inactivators of PRODH for use as chemical probes to study the roles of proline metabolism in cancer.

Parallel reaction pathways and noncovalent intermediates in thymidylate synthase revealed by experimental and computational tools

Thymidylate synthase was one of the most studied enzymes due to its critical role in molecular pathogenesis of cancer. Nevertheless, many atomistic details of its chemical mechanism remain unknown or debated, thereby imposing limits on design of novel mechanism-based anticancer therapeutics. Here, we report unprecedented isolation and characterization of a previously proposed intact noncovalent bisubstrate intermediate formed in the reaction catalyzed by thymidylate synthase. Free-energy surfaces of the bisubstrate intermediates interconversions computed with quantum mechanics/molecular mechanics (QM/MM) methods and experimental assessment of the corresponding kinetics indicate that the species is the most abundant productive intermediate along the reaction coordinate, whereas accumulation of the covalent bisubstrate species largely occurs in a parallel nonproductive pathway. Our findings not only substantiate relevance of the previously proposed noncovalent intermediate but also support potential implications of the overstabilized covalent intermediate in drug design targeting DNA biosynthesis.

New Insight into the Octamer of TYMS Stabilized by Intermolecular Cys43-Disulfide

Thymidylate synthase (TYMS) is an essential enzyme for the de novo synthesis of deoxythymidine monophosphate (dTMP) and has been a primary target for cancer chemotherapy. Although the physical structure of TYMS and the molecular mechanisms of TYMS catalyzing the conversion of deoxyuridine monophosphate (dUMP) to dTMP have been the subject of thorough studies, its oligomeric structure remains unclear. Here, we show that human TYMS not only exists in dimer form but also as an octamer by intermolecular Cys43-disulfide formation. We optimized the expression conditions of recombinant human TYMS using the Escherichia coli system. Using high-performance liquid chromatography⁻tandem mass spectrometry (HPLC⁻MS/MS), we have shown that purified TYMS has catalytic activity for producing dTMP. In the absence of reductant β-mercaptoethanol, SDS-PAGE and size exclusion chromatography (SEC) showed that the size of the TYMS protein is about 35 kDa, 70 kDa, and 280 kDa. When the Cys43 was mutated to Gly, the band of ~280 kDa and the peak of the octamer disappeared. Therefore, TYMS was determined to form an octamer, depending on the presence of Cys43-disulfide. By measuring steady-state parameters for the monomer, dimer, and octamer, we found the k_cat of the octamer was increased slightly more than the monomer. On the basis of these findings, we suggest that the octamer in the active state might have a potential influence on the design of new drug targets.

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.

The cysteine proteome

The cysteine (Cys) proteome is a major component of the adaptive interface between the genome and the exposome. The thiol moiety of Cys undergoes a range of biologic modifications enabling biological switching of structure and reactivity. These biological modifications include sulfenylation and disulfide formation, formation of higher oxidation states, S-nitrosylation, persulfidation, metalation, and other modifications. Extensive knowledge about these systems and their compartmentalization now provides a foundation to develop advanced integrative models of Cys proteome regulation. In particular, detailed understanding of redox signaling pathways and sensing networks is becoming available to allow the discrimination of network structures. This research focuses attention on the need for atlases of Cys modifications to develop systems biology models. Such atlases will be especially useful for integrative studies linking the Cys proteome to imaging and other omics platforms, providing a basis for improved redox-based therapeutics. Thus, a framework is emerging to place the Cys proteome as a complement to the quantitative proteome in the omics continuum connecting the genome to the exposome.

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.

The R163K mutant of human thymidylate synthase is stabilized in an active conformation: structural asymmetry and reactivity of cysteine 195

Loop 181-197 of human thymidylate synthase (hTS) populates two conformational states. In the first state, Cys195, a residue crucial for catalytic activity, is in the active site (active conformer); in the other conformation, it is about 10 A away, outside the active site (inactive conformer). We have designed and expressed an hTS variant, R163K, in which the inactive conformation is destabilized. The activity of this mutant is 33% higher than that of wt hTS, suggesting that at least one-third of hTS populates the inactive conformer. Crystal structures of R163K in two different crystal forms, with six and two subunits per asymmetric part of the unit cells, have been determined. All subunits of this mutant are in the active conformation while wt hTS crystallizes as the inactive conformer in similar mother liquors. The structures show differences in the environment of catalytic Cys195, which correlate with Cys195 thiol reactivity, as judged by its oxidation state. Calculations show that the molecular electrostatic potential at Cys195 differs between the subunits of the dimer. One of the dimers is asymmetric with a phosphate ion bound in only one of the subunits. In the absence of the phosphate ion, that is in the inhibitor-free enzyme, the tip of loop 47-53 is about 11 A away from the active site.

Functional analysis of the Mycobacterium tuberculosis FAD-dependent thymidylate synthase, ThyX, reveals new amino acid residues contributing to an extended ThyX motif

A novel FAD-dependent thymidylate synthase, ThyX, is present in a variety of eubacteria and archaea, including the mycobacteria. A short motif found in all thyX genes, RHRX(7-8)S, has been identified. The three-dimensional structure of the Mycobacterium tuberculosis ThyX enzyme has been solved. Building upon this information, we used directed mutagenesis to produce 67 mutants of the M. tuberculosis thyX gene. Each enzyme was assayed to determine its ability to complement the defect in thymidine biosynthesis in a delta thyA strain of Escherichia coli. Enzymes from selected strains were then tested in vitro for their ability to catalyze the oxidation of NADPH and the release of a proton from position 5 of the pyrimidine ring of dUMP. The results defined an extended motif of amino acids essential to enzyme activity in M. tuberculosis (Y44X(24)H69X(25)R95HRX(7)S105XRYX(90)R199 [with the underlined histidine acting as the catalytic residue and the underlined serine as the nucleophile]) and provided insight into the ThyX reaction mechanism. ThyX is found in a variety of bacterial pathogens but is absent in humans, which depend upon an unrelated thymidylate synthase, ThyA. Therefore, ThyX is a potential target for development of antibacterial drugs.

Role of Y94 in proton and hydride transfers catalyzed by thymidylate synthase

Thymidylate synthase (TS) catalyzes the substitution of a carbon-bound proton in a uracil base by a methyl group to yield thymine in the de novo biosynthesis of this DNA base. The enzymatic mechanism involves making and breaking several covalent bonds. Traditionally, a conserved tyrosine (Y94 in Escherichia coli, Y146 in Lactobacillus casei, and Y135 in humans) was assumed to serve as the general base catalyzing the proton abstraction. That assumption was examined here by comparing the nature of the proton abstraction using wild-type (wt) E. coli TS (ecTS) and its Y94F mutant (with a turnover rate reduced by 2 orders of magnitude). A subsequent hydride transfer was also studied using the wt and Y94F. The physical nature of both H-transfer steps was examined by determining intrinsic kinetic isotope effects (KIEs). Surprisingly, the findings did not suggest a direct role for Y94 in the proton abstraction step. The effect of this mutation on the subsequent hydride transfer was examined by a comparison of the temperature dependency of the intrinsic KIE on both the wt and the mutant. The intrinsic KIEs for Y94F at physiological temperatures were slightly smaller than those for wt but, otherwise, were as temperature-independent, suggesting a perfectly preorganized reaction coordinate for both enzymes. At reduced temperatures, however, the KIE for the mutant increased with a decrease in temperature, indicating a poorly preorganized reaction coordinate. Other kinetic and structural properties were also compared, and the findings suggested that Y94 is part of a H-bond network that plays a critical role at a step between the proton and the hydride transfers, presumably the dissociation of H4folate from the covalently bound intermediate. The possibility that no single residue serves as the general base in question but, rather, that the whole network of H-bonds at the active site catalyzes proton abstraction is discussed.

Molecular analysis of the thymidine-auxotrophic small colony variant phenotype of Staphylococcus aureus

Thymidine-auxotrophic small colony variants (SCVs) of Staphylococcus aureus are frequently isolated from the chronically infected airways of patients suffering from cystic fibrosis. To date, little is known regarding the molecular mechanisms leading to the formation of this special phenotype, but the auxotrophism for thymidine suggests that impaired thymidine metabolism might play a major role. Sequence analysis of the thymidylate synthase-encoding thyA gene of six clinical thymidine-auxotrophic S. aureus SCVs revealed that all isolates had mutations within thyA. In five isolates the function of the thymidylate synthase was definitely impaired: three of them showed a truncation of the thyA coding sequence by nonsense or frame-shift mutations, in one further isolate the active site of the enzyme was affected by an internal 12-bp deletion, and another isolate had a 173-bp deletion spanning the 5'-terminal region of thyA and the preceding DNA sequence. The sixth isolate showed two amino acid substitutions within the thyA gene product. To confirm the importance of impaired thymidylate synthase synthesis or activity for the formation of the thymidine-auxotrophic SCV phenotype, we constructed a thyA knock-out mutant of a wild-type S. aureus strain. This mutant showed all characteristics of clinical SCVs, such as slow growth, decreased pigment production, reduced hemolytic activity, auxotrophism for thymidine, resistance to trimethoprim/sulfamethoxazol, and reduced plasma coagulase activity. Complementation of the thyA knock-out mutant with intact thyA in trans nearly restored the normal phenotype. In conclusion, these data confirm at the molecular level that impaired thymidylate synthase function is causative for the formation of the thymidine-auxotrophic SCV phenotype in S. aureus.

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.

Functional analysis of substrate and cofactor complex structures of a thymidylate synthase-complementing protein

Like thymidylate synthase (TS) in eukaryotes, the thymidylate synthase-complementing proteins (TSCPs) are mandatory for cell survival of many prokaryotes in the absence of external sources of thymidylate. Details of the mechanism of this novel family of enzymes are unknown. Here, we report the structural and functional analysis of a TSCP from Thermotoga maritima and its complexes with substrate, analogs, and cofactor. The structures presented here provide a basis for rationalizing the TSCP catalysis and reveal the possibility of the design of an inhibitor. We have identified a new helix-loop-strand FAD binding motif characteristic of the enzymes in the TSCP family. The presence of a hydrophobic core with residues conserved among the TSCP family suggests a common overall fold.

Interaction of thymidylate synthase with the 5'-thiophosphates, 5'-dithiophosphates, 5'-H-phosphonates and 5'-S-thiosulfates of 2'-deoxyuridine, thymidine and 5-fluoro-2'-deoxyuridine

New analogs of dUMP, dTMP and 5-fluoro-dUMP, including the corresponding 5'-thiophosphates (dUMPS, dTMPS and FdUMPS), 5'-dithiophosphates (dUMPS2, dTMPS2 and FdUMPS2), 5'-H-phosphonates (dUMP-H, dTMP-H and FdUMP-H) and 5'-S-thiosulfates (dUSSO3, dTSSO3 and FdUSSO3), have been synthesized and their interactions studied with highly purified mammalian thymidylate synthase. dUMPS and dUMPS2 proved to be good substrates, and dTMPS and dTMPS2 classic competitive inhibitors, only slightly weaker than dTMP. Their 5-fluoro congeners behaved as potent, slow-binding inhibitors. By contrast, the corresponding 5'-H-phosphonates and 5'-S-thiosulfates displayed weak activities, only FdUMP-H and FdUSSO3 exhibiting significant interactions with the enzyme, as weak competitive slow-binding inhibitors versus dUMR The pH-dependence of enzyme time-independent inhibition by FdUMP and FdUMPS was found to correlate with the difference in pKa values of the phosphate and thiophosphate groups, the profile of FdUMPS being shifted (approximately 1 pH unit) toward lower pH values, so that binding of dUMP and its analogs is limited by the phosphate secondary hydroxyl ionization. Hence, together with the effects of 5'-H-phosphonate and 5'-S-thiosulfate substituents, the much weaker interactions of the nucleotide analogs (3-5 orders of magnitude lower than for the parent 5'-phosphates) with the enzyme is further evidence that the enzyme's active center prefers the dianionic phosphate group for optimum binding.

Structure of human thymidylate synthase suggests advantages of chemotherapy with noncompetitive inhibitors

Thymidylate synthase (TS) is a major target in the chemotherapy of colorectal cancer and some other neoplasms. The emergence of resistance to the treatment is often related to the increased levels of TS in cancer cells, which have been linked to the elimination of TS binding to its own mRNA upon drug binding, a feedback regulatory mechanism, and/or to the increased stability to intracellular degradation of TS.drug complexes (versus unliganded TS). The active site loop of human TS (hTS) has a unique conformation resulted from a rotation by 180 degrees relative to its orientation in bacterial TSs. In this conformation, the enzyme must be inactive, because the catalytic cysteine is no longer positioned in the ligand-binding pocket. The ordered solvent structure obtained from high resolution crystallographic data (2.0 A) suggests that the inactive loop conformation promotes mRNA binding and intracellular degradation of the enzyme. This hypothesis is supported by fluorescence studies, which indicate that in solution both active and inactive forms of hTS are present. The binding of phosphate ion shifts the equilibrium toward the inactive conformation; subsequent dUMP binding reverses the equilibrium toward the active form. Thus, TS inhibition via stabilization of the inactive conformation should lead to less resistance than is observed with presently used drugs, which are analogs of its substrates, dUMP and CH(2)H(4)folate, and bind in the active site, promoting the active conformation. The presence of an extension at the N terminus of native hTS has no significant effect on kinetic properties or crystal structure.