Conservation and Role of Electrostatics in Thymidylate Synthase

Molecular Mechanism of Thymidylate Synthase Inhibition by N4-Hydroxy-dCMP in View of Spectrophotometric and Crystallographic Studies

Novel evidence is presented allowing further clarification of the mechanism of the slow-binding thymidylate synthase (TS) inhibition by N⁴-hydroxy-dCMP (N⁴-OH-dCMP). Spectrophotometric monitoring documented time- and temperature-, and N⁴-OH-dCMP-dependent TS-catalyzed dihydrofolate production, accompanying the mouse enzyme incubation with N⁴-OH-dCMP and N^5,10-methylenetetrahydrofolate, known to inactivate the enzyme by the covalent binding of the inhibitor, suggesting the demonstrated reaction to be uncoupled from the pyrimidine C(5) methylation. The latter was in accord with the hypothesis based on the previously presented structure of mouse TS (cf. PDB ID: 4EZ8), and with conclusions based on the present structure of the parasitic nematode Trichinella spiralis, both co-crystallized with N⁴-OH-dCMP and N^5,10-methylenetetrahdrofolate. The crystal structure of the mouse TS-N⁴-OH-dCMP complex soaked with N^5,10-methylenetetrahydrofolate revealed the reaction to run via a unique imidazolidine ring opening, leaving the one-carbon group bound to the N(10) atom, thus too distant from the pyrimidine C(5) atom to enable the electrophilic attack and methylene group transfer.

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.

Substrate Channeling via a Transient Protein-Protein Complex: The case of D-Glyceraldehyde-3-Phosphate Dehydrogenase and L-Lactate Dehydrogenase

Substrate channeling studies have frequently failed to provide conclusive results due to poor understanding of this subtle phenomenon. We analyzed the mechanism of NADH-channeling from D-glyceraldehyde-3-phosphate dehydrogenase (GAPDH) to L-lactate Dehydrogenase (LDH) using enzymes from different cells. Enzyme kinetics studies showed that LDH activity with free NADH and GAPDH-NADH complex always take place in parallel. The channeling is observed only in assays that mimic cytosolic conditions where free NADH concentration is negligible and the GAPDH-NADH complex is dominant. Molecular dynamics and protein-protein interaction studies showed that LDH and GAPDH can form a leaky channeling complex only at the limiting NADH concentrations. Surface calculations showed that positive electric field between the NAD(H) binding sites on LDH and GAPDH tetramers can merge in the LDH-GAPDH complex. NAD(H)-channeling within the LDH-GAPDH complex can be an extension of NAD(H)-channeling within each tetramer. In the case of a transient LDH-(GAPDH-NADH) complex, the relative contribution from the channeled and the diffusive paths depends on the overlap between the off-rates for the LDH-(GAPDH-NADH) complex and the GAPDH-NADH complex. Molecular evolution or metabolic engineering protocols can exploit substrate channeling for metabolic flux control by fine-tuning substrate-binding affinity for the key enzymes in the competing reaction paths.

High basal heat-shock protein expression in bats confers resistance to cellular heat/oxidative stress

Bats, unique among mammals with powered flight, have many species with the longest size-proportionate lifespan of all mammals. Evolutionary adaptations would have been required to survive the elevated body temperatures during flight. Heat shock protein (HSP), highly conserved master regulators of cell stress, expression was examined across tissues and various cell lines in bats. Basal expression level of major HSPs (HSP70 and HSP90) is significantly higher in two different bat species compared to other mammals. This HSP expression could be a bat-unique, key factor to modulate cellular stress and death. Consequently, bat cells survive prolonged heat treatment, along with other stress stimuli, in a HSP-dependent manner, whereas other mammalian cells succumbed. This suggests HSP expression in bats could be an important adaption to intrinsic metabolic stresses like flight and therefore an important model to study stress resilience and longevity in general.

A novel thymidylate synthase from the Vibrionales, Alteromonadales, Aeromonadales, and Pasteurellales (VAAP) clade with altered nucleotide and folate binding sites

Thymidylate synthase (TS, E.C. 2.1.1.45) is a crucial enzyme for de novo deoxythymidine monophosphate (dTMP) biosynthesis. The gene for this enzyme is thyA, which encodes the folate-dependent TS that converts deoxyuridine monophosphate group (dUMP) into (dTMP) using the cofactor 5,10-methylenetetrahydrofolate (mTHF) as a carbon donor. We identified the thyA gene in the genome of the Vibrio parahaemolyticus strain FIM-S1708+ that is innocuous to humans but pathogenic to crustaceans. Surprisingly, we found changes in the residues that bind the substrate dUMP and mTHF, previously postulated as invariant among all TSs known (Finer-Moore, Santi & Stroud, 2003). Interestingly, those amino acid changes were also found in a clade of microorganisms that contains Vibrionales, Alteromonadales, Aeromonadales, and Pasteurellales (VAAP) from the Gammaproteobacteria class. In this work, we studied the biochemical properties of recombinant TS from V. parahemolyticus FIM-S1708+ (VpTS) to address the natural changes in the TS amino acid sequence of the VAAP clade. Interestingly, the K_m for dUMP was 27.3 ± 4.3 µM, about one-fold larger compared to other TSs. The K_m for mTHF was 96.3 ± 18 µM, about three- to five-fold larger compared to other species, suggesting also loss of affinity. Thus, the catalytic efficiency was between one or two orders of magnitude smaller for both substrates. We used trimethoprim, a common antibiotic that targets both TS and DHFR for inhibition studies. The IC₅₀ values obtained were high compared to other results in the literature. Nonetheless, this molecule could be a lead for the design antibiotics towards pathogens from the VAAP clade. Overall, the experimental results also suggest that in the VAAP clade the nucleotide salvage pathway is important and should be investigated, since the de novo dTMP synthesis appears to be compromised by a less efficient thymidylate synthase.

Crystal structures of nematode (parasitic T. spiralis and free living C. elegans), compared to mammalian, thymidylate synthases (TS). Molecular docking and molecular dynamics simulations in search for nematode-specific inhibitors of TS

Three crystal structures are presented of nematode thymidylate synthases (TS), including Caenorhabditis elegans (Ce) enzyme without ligands and its ternary complex with dUMP and Raltitrexed, and binary complex of Trichinella spiralis (Ts) enzyme with dUMP. In search of differences potentially relevant for the development of species-specific inhibitors of the nematode enzyme, a comparison was made of the present Ce and Ts enzyme structures, as well as binary complex of Ce enzyme with dUMP, with the corresponding mammalian (human, mouse and rat) enzyme crystal structures. To complement the comparison, tCONCOORD computations were performed to evaluate dynamic behaviors of mammalian and nematode TS structures. Finally, comparative molecular docking combined with molecular dynamics and free energy of binding calculations were carried out to search for ligands showing selective affinity to T. spiralis TS. Despite an overall strong similarity in structure and dynamics of nematode vs mammalian TSs, a pool of ligands demonstrating predictively a strong and selective binding to TsTS has been delimited. These compounds, the E63 family, locate in the dimerization interface of TsTS where they exert species-specific interactions with certain non-conserved residues, including hydrogen bonds with Thr174 and hydrophobic contacts with Phe192, Cys191 and Tyr152. The E63 family of ligands opens the possibility of future development of selective inhibitors of TsTS and effective agents against trichinellosis.