Plastic adaptation toward mutations in proteins: structural comparison of thymidylate synthases

Effects of Osmolytes on Ligand Binding to Dihydropteroate Synthase from Bacillus anthracis

Osmolyte interactions with ligands can affect their affinity for proteins and are dependent upon the cosolute and the functional groups of the ligand. Here, we explored ligand binding to Bacillus anthracis dihydropteroate synthase (BaDHPS) under osmotic stress conditions. Osmolyte effects were specific to the cosolute and ligand, suggesting interaction of the osmolytes with the free ligands in solution. The association rates of pterin pyrophosphate were mostly unaffected by the osmolytes, except for a 2-fold decrease in the presence of 1 M trehalose, while the dissociation rates decreased in most osmolyte solutions. The viscosity and dielectric constant of the solution did not correlate with the effects of the osmolytes. Experimental results were compared with predicted preferential interaction coefficients (Δμ₂₃/RT) between the osmolytes and ligands. The Δμ₂₃/RT were able to predict the experimental data for most of the osmolytes. Trehalose and proline effects did not correlate with the predicted values, indicating that these two osmolytes may affect binding in more complex ways than simple preferential interactions. Additionally, osmolytes weakly interacted with the sulfa drug sulfathiazole, which altered its affinity for BaDHPS, suggesting that these types of weak interactions can also impact drug binding. As osmolytes affect ligands binding to two different folate cycle enzymes (DHFRs and DHPS), we predicted how ligand binding to other folate cycle enzymes will be altered by the presence of osmolytes.

Thermodynamic and NMR Assessment of Ligand Cooperativity and Intersubunit Communication in Symmetric Dimers: Application to Thymidylate Synthase

Thymidylate synthase (TS) is a homodimeric enzyme with evidence for negative regulation of one protomer while the other protomer acts on substrate, so called half-the-sites reactivity. The mechanisms by which multisubunit allosteric proteins communicate between protomers is not well understood, and the simplicity of dimeric systems has advantages for observing conformational and dynamic processes that functionally connect distance-separated active sites. This review considers progress in overcoming the inherent challenges of accurate thermodynamic and atomic-resolution characterization of interprotomer communication mechanisms in symmetric protein dimers, with TS used as an example. Isothermal titration calorimetry (ITC) is used to measure ligand binding cooperativity, even in cases where the two binding enthalpies are similar, and NMR spectroscopy is used to detect site-specific changes occurring in the two protomers. The NMR approach makes use of mixed-labeled dimers, enabling protomer-specific detection of signals in the singly ligated state. The rich informational content of the NMR signals from the singly ligated state, relative to the apo and saturated states, requires new considerations that do not arise in simple cases of 1:1 protein-ligand interactions.

Isolation and characterization of a N4-like lytic bacteriophage infecting Vibrio splendidus, a pathogen of fish and bivalves

A novel virulent bacteriophage, vB_VspP_pVa5, infecting a strain of Vibrio splendidus was isolated from a sea-cage aquaculture farm in Greece, and characterized using microbiological methods and genomic analysis. Bacteriophage vB_VspP_pVa5 is a N4-like podovirus with an icosahedral head measuring 85 nm in length and a short non-contractile tail. The phage had a narrow host range infecting only the bacterial host, a latent period of 30 min and a burst size of 24 virions per infected bacterium. Its genome size was 78,145 bp and genomic analysis identified 107 densely-packed genes, 40 of which could be annotated. In addition to the very large virion encapsulated DNA-dependent RNA polymerase which is the signature of the N4-like genus, an interesting feature of the novel phage is the presence of a self-splicing group I intron in the thymidylate synthase gene. A tRNA^Stop interrupted by a ~2.5kb open reading frame–containing area was also identified. The absence of genes related to lysogeny along with the high efficacy observed during in vitro cell lysis trials, indicate that the vB_VspP_pVa5 is a potential candidate component in a bacteriophage cocktail suitable for the biological control of V. splendidus in aquaculture.

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.

The influence of active site conformations on the hydride transfer step of the thymidylate synthase reaction mechanism

The hydride transfer from C6 of tetrahydrofolate to the reaction's exocyclic methylene-dUMP intermediate is the rate limiting step in thymidylate synthase (TSase) catalysis. This step has been studied by means of QM/MM molecular dynamics simulations to generate the corresponding free energy surfaces. The use of two different initial X-ray structures has allowed exploring different conformational spaces and the existence of chemical paths with not only different reactivities but also different reaction mechanisms. The results confirm that this chemical conversion takes place preferentially via a concerted mechanism where the hydride transfer is conjugated to thiol-elimination from the product. The findings also confirm the labile character of the substrate-enzyme covalent bond established between the C6 of the nucleotide substrate and a conserved cysteine residue. The calculations also reproduce and rationalize a normal H/T 2° kinetic isotope effect measured for that step. From a computational point of view, the results demonstrate that the use of an incomplete number of coordinates to describe the real reaction coordinate can render biased results.

Bacterial Thymidylate Synthase Binds Two Molecules of Substrate and Cofactor without Cooperativity

Thymidylate synthase (TSase) is a clinically important enzyme because it catalyzes synthesis of the sole de novo source of deoxy-thymidylate. Without this enzyme, cells die a "thymineless death" since they are starved of a crucial DNA synthesis precursor. As a drug target, TSase is well studied in terms of its structure and reaction mechanism. An interesting mechanistic feature of dimeric TSase is that it is "half-the-sites reactive", which is a form of negative cooperativity. Yet, the basis for this is not well-understood. Some experiments point to cooperativity at the binding steps of the reaction cycle as being responsible for the phenomenon, but the literature contains conflicting reports. Here we use ITC and NMR to resolve these inconsistencies. This first detailed thermodynamic dissection of multisite binding of dUMP to E. coli TSase shows the nucleotide binds to the free and singly bound forms of the enzyme with nearly equal affinity over a broad range of temperatures and in multiple buffers. While small but significant differences in ΔC°P for the two binding events show that the active sites are not formally equivalent, there is little-to-no allostery at the level of ΔG°bind. In addition NMR titration data reveal that there is minor intersubunit cooperativity in formation of a ternary complex with the mechanism based inhibitor, 5F-dUMP, and cofactor. Taken together, the data show that functional communication between subunits is minimal for both binding steps of the reaction coordinate.

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.

2'-Deoxyuridine 5'-monophosphate substrate displacement in thymidylate synthase through 6-hydroxy-2H-naphtho[1,8-bc]furan-2-one derivatives

Thymidylate synthase (TS) is a target for antifolate-based chemotherapies of microbial and human diseases. Here, ligand-based, synthetic, and X-ray crystallography studies led to the discovery of 6-(3-cyanobenzoyloxy)-2-oxo-2H-naphto[1,8-bc]furan, a novel inhibitor with a Ki of 310 nM against Pneumocystis carinii TS. The X-ray ternary complex with Escherichia coli TS revealed, for the first time, displacement of the substrate toward the dimeric protein interface, thus providing new opportunities for further design of specific inhibitors of microbial pathogens.

Effect of halogen substitutions on dUMP to stability of thymidylate synthase/dUMP/mTHF ternary complex using molecular dynamics simulation

The stability of the thymidylate synthase (TS)/2-deoxyuridine-5-monophosphate (dUMP)/5,10-methylene-5,6,7,8-tetrahydrofolate (mTHF) ternary complex formation and Michael addition are considered as important steps that are involved in the inhibition mechanism of the anticancer prodrug 5-fluorouracil (5-FU). Here, the effect of three different halogen substitutions on the C-5 position of the dUMP (XdUMPs = FdUMP, CldUMP, and BrdUMP), the normal substrate, on the stability of the TS/dUMP and TS/dUMP/mTHF binary and ternary complexes, respectively, was investigated via molecular dynamics simulation. The simulated results revealed that the stability of all the systems was substantially increased by mTHF binding to the catalytic pocket. In the ternary complex, a much greater stabilization of the dUMP and XdUMPs through electrostatic interactions, including charge-charge and hydrogen bond interactions, was found compared to mTHF. An additional unique hydrogen bond between the substituted fluorine of FdUMP and the hydroxyl group of the TS Y94 residue was observed in both the binary and ternary complexes. The distance between the S(-) atom of the TS C146 residue and the C6 atom of dUMP, at <4 Å in all systems, suggested that a Michael addition with the formation of a S-C6 covalent bond potentially occurred, although the hydrogen atom on C6 of dUMP is substituted by a halogen atom. The MM/PBSA binding free energy revealed the significant role of the bridging waters around the ligands in the increased binding affinity (∼10 kcal/mol) of dUMP/XdUMP, either alone or together with mTHF, toward TS. The order of the averaged binding affinity in the ternary systems was found to be CldUMP ≈ FdUMP > dUMP > BrdUMP, suggesting that CldUMP could be a potent candidate TS inhibitor, the same as FdUMP (the metabolite form of 5-FU).

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.

The structure of Enterococcus faecalis thymidylate synthase provides clues about folate bacterial metabolism

Drug resistance to therapeutic antibiotics poses a challenge to the identification of novel targets and drugs for the treatment of infectious diseases. Infections caused by Enterococcus faecalis are a major health problem. Thymidylate synthase (TS) from E. faecalis is a potential target for antibacterial therapy. The X-ray crystallographic structure of E. faecalis thymidylate synthase (EfTS), which was obtained as a native binary complex composed of EfTS and 5-formyltetrahydrofolate (5-FTHF), has been determined. The structure provides evidence that EfTS is a half-of-the-sites reactive enzyme, as 5-FTHF is bound to two of the four independent subunits present in the crystal asymmetric unit. 5-FTHF is a metabolite of the one-carbon transfer reaction catalysed by 5-formyltetrahydrofolate cyclo-ligase. Kinetic studies show that 5-FTHF is a weak inhibitor of EfTS, suggesting that the EfTS-5-FTHF complex may function as a source of folates and/or may regulate one-carbon metabolism. The structure represents the first example of endogenous 5-FTHF bound to a protein involved in folate metabolism.

Synthesis and evaluation of 6-aza-2'-deoxyuridine monophosphate analogs as inhibitors of thymidylate synthases, and as substrates or inhibitors of thymidine monophosphate kinase in Mycobacterium tuberculosis

A series of 5-substituted analogs of 6-aza-2'-deoxyuridine 5'-monophosphate, 6-aza-dUMP, has been synthesized and evaluated as potential inhibitors of the two mycobacterial thymidylate synthases (i.e., a flavin-dependent thymidylate synthase, ThyX, and a classical thymidylate synthase, ThyA). Replacement of C(6) of the natural substrate dUMP by a N-atom in 6-aza-dUMP 1a led to a derivative with weak ThyX inhibitory activity (33% inhibition at 50 μM). Introduction of alkyl and aryl groups at C(5) of 1a resulted in complete loss of inhibitory activity, whereas the attachment of a 3-(octanamido)prop-1-ynyl side chain in derivative 3 retained the weak level of mycobacterial ThyX inhibition (40% inhibition at 50 μM). None of the synthesized derivatives displayed any significant inhibitory activity against mycobacterial ThyA. The compounds have also been evaluated as potential inhibitors of mycobacterial thymidine monophosphate kinase (TMPKmt). None of the derivatives showed any significant TMPKmt inhibition. However, replacement of C(6) of the natural substrate (dTMP) by a N-atom furnished 6-aza-dTMP (1b), which still was recognized as a substrate by TMPKmt.

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.

Evolution of metamorphism in thymidylate synthases within the primate lineages

Crystal structures of human thymidylate synthase (hTS) revealed that the protein exists in active and inactive conformations, defined by the position of a loop containing the active site nucleophile. TS is highly homologous among diverse species; however, the residue at position 163 (hTS) differs among species. Arginine at this position is predicted by structural modeling to enable conformational switching. Arginine or lysine is reported at this position in all mammals in the GenBank and Ensembl databases, with arginine reported in only primates. Sequence analysis of the TS gene of representative primates revealed that arginine occurs at this relative position in all primates except a representative of prosimians. Mutant human proteins were created with residues at position 163 that occur in TSs from prokaryotes and eukaryotes. Catalytic constants (k(cat)) of mutant enzymes were 45-149% of hTS, with the lysine mutant (R163K) exhibiting the highest k(cat). The effect of lysine substitution on solution structure and on ligand binding was investigated. R163K exhibited higher intrinsic fluorescence, a more negative molar ellipticity, and higher dissociation constants (K(d)) for ligands that modulate protein conformation than hTS. Temperature effects on intrinsic fluorescence and catalytic activity of hTS and R163K are consistent with proteins populating different conformational states. The data indicate that the enzyme with arginine at the position corresponding to 163 (hTS) evolved after the divergence of prosimians and simians and that substitution of lysine by arginine confers unique structural and functional properties to the enzyme expressed in simian primates.

Application of Grote-Hynes theory to the reaction catalyzed by thymidylate synthase

A theoretical study of dynamic effects on the rate-limiting step of the thymidylate synthase catalyzed reaction has been carried out by means of Grote-Hynes theory, successfully predicting the values of the recrossing effects for a chemical reaction that involves the transfer of a classical light particle. The transmission coefficients, obtained at 278, 293, 303, and 313 K, are almost invariant and in all cases far from unity, revealing a significant coupling of the environment motions and the reaction coordinate. Nevertheless, their energetic contribution to the activation free energy represents less than 0.50 kcal/mol for each of the four tested temperatures. Calculation of the transmission coefficient for the isotopically labeled hydride transfer has rendered almost the same values, in agreement with the experimentally observed temperature-independent KIEs. Fourier transform of the time-dependent friction kernel at these four temperatures has allowed obtaining the transition-state friction spectra, which present very small dependence with temperature. Their analysis has led to the identification of some key vibrational modes governing the coupling between the reaction coordinate and the protein environment, thus identifying the relevant motions in the active site and obtaining a full picture of the role of each amino acid.

High resolution mapping and positional cloning of ENU-induced mutations in the Rw region of mouse chromosome 5

Background

Forward genetic screens in mice provide an unbiased means to identify genes and other functional genetic elements in the genome. Previously, a large scale ENU mutagenesis screen was conducted to query the functional content of a ~50 Mb region of the mouse genome on proximal Chr 5. The majority of phenotypic mutants recovered were embryonic lethals.

Results

We report the high resolution genetic mapping, complementation analyses, and positional cloning of mutations in the target region. The collection of identified alleles include several with known or presumed functions for which no mutant models have been reported (Tbc1d14, Nol14, Tyms, Cad, Fbxl5, Haus3), and mutations in genes we or others previously reported (Tapt1, Rest, Ugdh, Paxip1, Hmx1, Otoe, Nsun7). We also confirmed the causative nature of a homeotic mutation with a targeted allele, mapped a lethal mutation to a large gene desert, and localized a spermiogenesis mutation to a region in which no annotated genes have coding mutations. The mutation in Tbc1d14 provides the first implication of a critical developmental role for RAB-GAP-mediated protein transport in early embryogenesis.

Conclusion

This collection of alleles contributes to the goal of assigning biological functions to all known genes, as well as identifying novel functional elements that would be missed by reverse genetic approaches.

QM/MM study of thymidylate synthase: enzymatic motions and the temperature dependence of the rate limiting step

Thymidylate synthase (TS) is an enzyme that catalyzes a complex cascade of reactions. A theoretical study of the reduction of an exocyclic methylene intermediate by hydride transfer from the 6S position of 5,6,7,8-tetrahydrofolate (H(4)folate), to form 2'-deoxyuridine 5'-monophosphate (dTMP) and 7,8-dihydrofolate (H(2)folate), has been carried out using hybrid quantum mechanics/molecular mechanics methods. This step is of special interest because it is the rate-limiting step of the reaction catalyzed by TS. The acceptor of this hydride is an intermediate that is covalently bound to the enzyme via a thioether bond to an overall conserved active site cysteine residue (Cys146 in Escherichia coli). Heretofore, whether the hydride transfer precedes the thiol abstraction that releases the product from the enzyme or whether these two processes are concerted has been an open question. We have examined this step in terms of free energy surfaces obtained at the same temperatures we previously used in experimental studies of this mechanistic step (273-313 K). Analysis of the results reveals that substantial features of the reaction and the nature of the H-transfer seem to be temperature independent, in agreement with our experimental data. The findings also indicate that the hydride transfer and the scission of Cys146 take place in a concerted but asynchronous fashion. This 1,3-S(N)2 substitution is assisted by arginine 166 and several other arginine residues in the active site that polarize the carbon-sulfur bond and stabilize the charge transferred from cofactor to substrate. Finally, the simulation elucidates the molecular details of the enzyme's motion that brings the system to its transition state and, in accordance with the experimental data, indicates that this "tunneling ready" conformation is temperature independent.

Formation of catalytically active cross-species heterodimers of thymidylate synthase from Plasmodium falciparum and Plasmodium vivax

Thymidylate synthase (TS) of Plasmodium dihydrofolate reductase-thymidylate synthase (DHFR-TS) functions as a homodimeric enzyme with two active sites located near the subunit interface. The dimerization is essential for catalysis, since the active site of each subunit contains amino acid residues contributed from the other TS domain. In P. falciparum DHFR-TS, it has been shown that the active sites require Cys-490 from one domain and Arg-470 donated from the other domain. Mutants of these two series can complement one another giving rise to active enzyme. Here, the potential to form cross-species heterodimers between P. falciparum and P. vivax TS has been explored. Formation of cross-species heterodimer was tested by co-transformation of TS-inactive Cys-490 mutants of P. falciparum or P. vivax with corresponding TS-inactive Arg-486 mutants of P. vivax or P. falciparum into thymidine-requiring Escherichia coli. Active heterodimers were detected by subunit complementation and 6-[(3)H]-FdUMP binding assays. All combinations of the mutants tested, except for (Pf)R470A+(Pv)C506Y, were able to form catalytically active cross-species heterodimers. The single active site formed by (Pf)R470D+(Pv)C506Y and (Pv)R486D+(Pf)C490A pairs of cross-species heterodimers has k(cat) and K(m) values similar to those of intra-species heterodimers of P. falciparum and P. vivax. This is the first report to demonstrate that the TS subunit interface between Plasmodium species is sufficiently conserved to allow formation of fully active cross-species heterodimer.

Variants of human thymidylate synthase with loop 181-197 stabilized in the inactive conformation

Loop 181-197 of human thymidylate synthase (hTS) populates two major conformations, essentially corresponding to the loop flipped by 180 degrees . In one of the conformations, the catalytic Cys195 residue lies distant from the active site making the enzyme inactive. Ligands stabilizing this inactive conformation may function as allosteric inhibitors. To facilitate the search for such inhibitors, we have expressed and characterized several mutants designed to shift the equilibrium toward the inactive conformer. In most cases, the catalytic efficiency of the mutants was only somewhat impaired with values of k(cat)/K(m) reduced by factors in a 2-12 range. One of the mutants, M190K, is however unique in having the value of k(cat)/K(m) smaller by a factor of approximately 7500 than the wild type. The crystal structure of this mutant is similar to that of the wt hTS with loop 181-197 in the inactive conformation. However, the direct vicinity of the mutation, residues 188-194 of this loop, assumes a different conformation with the positions of C(alpha) shifted up to 7.2 A. This affects region 116-128, which became ordered in M190K while it is disordered in wt. The conformation of 116-128 is however different than that observed in hTS in the active conformation. The side chain of Lys190 does not form contacts and is in solvent region. The very low activity of M190K as compared to another mutant with a charged residue in this position, M190E, suggests that the protein is trapped in an inactive state that does not equilibrate easily with the active conformer.

Crystal structure of an RluF-RNA complex: a base-pair rearrangement is the key to selectivity of RluF for U2604 of the ribosome

Escherichia coli pseudouridine synthase RluF is dedicated to modifying U2604 in a stem-loop of 23S RNA, while a homologue, RluB, modifies the adjacent base, U2605. Both uridines are in the same RNA stem, separated by approximately 4 A. The 3.0 A X-ray crystal structure of RluF bound to the isolated stem-loop, in which U2604 is substituted by 5-fluorouridine to prevent catalytic turnover, shows RluF distinguishes closely spaced bases in similar environments by a selectivity mechanism based on a frameshift in base pairing. The RNA stem-loop is bound to a conserved binding groove in the catalytic domain. A base from a bulge in the stem, A2602, has folded into the stem, forcing one strand of the RNA stem to translate by one position and thus positioning U2604 to flip into the active site. RluF does not modify U2604 in mutant stem-loops that lack the A2602 bulge and shows dramatically higher activity for a stem-loop with a mutation designed to facilitate A2602 refolding into the stem with concomitant RNA strand translation. Residues whose side chains contact rearranged bases in the bound stem-loop, while conserved among RluFs, are not conserved between RluFs and RluBs, suggesting that RluB does not bind to the rearranged stem loop.

Identification and characterization of a novel thymidylate synthase from deep-sea thermophilic bacteriophage Geobacillus virus E2

Thymidylate synthase (TS) is essential for de novo synthesis of dTMP and is a key enzyme involved in DNA synthesis and transcriptional regulation of organisms. Due to their biologic importance, TSs have been intensively studied. In this investigation, a thermostable TS was identified from a deep-sea thermophilic bacteriophage Geobacillus virus E2 (GVE2). It was demonstrated that GVE2-TS was highly homologous to known TSs and contained five characteristic conserved domains. The temporal analyses by Northern and Western blots revealed that the GVE2-TS was transcribed and expressed early after Geobacillus virus E2 infection, identifying it as a viral early gene. As shown by gel mobility shift assays, the recombinant GVE2-TS protein had the capacity to bind its own mRNA. Our study presented the first report on thymidylate synthase from deep-sea thermophilic bacteriophage.

The crystal structure of E. coli rRNA pseudouridine synthase RluE

Pseudouridine synthase RluE modifies U2457 in a stem of 23 S RNA in Escherichia coli. This modification is located in the peptidyl transferase center of the ribosome. We determined the crystal structures of the C-terminal, catalytic domain of E. coli RluE at 1.2 A resolution and of full-length RluE at 1.6 A resolution. The crystals of the full-length enzyme contain two molecules in the asymmetric unit and in both molecules the N-terminal domain is disordered. The protein has an active site cleft, conserved in all other pseudouridine synthases, that contains invariant Asp and Tyr residues implicated in catalysis. An electropositive surface patch that covers the active site cleft is just wide enough to accommodate an RNA stem. The RNA substrate stem can be docked to this surface such that the catalytic Asp is adjacent to the target base, and a conserved Arg is positioned to help flip the target base out of the stem into the enzyme active site. A flexible RluE specific loop lies close to the conserved region of the stem in the model, and may contribute to substrate specificity. The stem alone is not a good RluE substrate, suggesting RluE makes additional interactions with other regions in the ribosome.

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.

A Unique RNA Fold in the RumA-RNA-Cofactor Ternary Complex Contributes to Substrate Selectivity and Enzymatic Function

A single base (U1939) within E. coli 23S ribosomal RNA is methylated by its dedicated enzyme, RumA. The structure of RumA/RNA/S-adenosylhomocysteine uncovers the mechanism for achieving unique selectivity. The single-stranded substrate is "refolded" on the enzyme into a compact conformation with six key intra-RNA interactions. The RNA substrate contributes directly to catalysis. In addition to the target base, a second base is "flipped out" from the core loop to stack against the adenine of the cofactor S-adenosylhomocysteine. Nucleotides in permuted sequence order are stacked into the site vacated by the everted target U1939 and compensate for the energetic penalty of base eversion. The 3' hairpin segment of the RNA binds distal to the active site and provides binding energy that contributes to enhanced catalytic efficiency. Active collaboration of RNA in catalysis leads us to conclude that RumA and its substrate RNA may reflect features from the earliest RNA-protein era.

Identification of the thymidylate synthase within the genome of white spot syndrome virus

Thymidylate synthase (TS) (EC 2.1.1.45) is essential for thede novosynthesis of dTMP in prokaryotic and eukaryotic organisms. Within the white spot syndrome virus (WSSV) genome, an open reading frame (WSV067) that encodes a 289 amino acid polypeptide showed significant homology to all known TSs from species including mammals, plants, fungi, protozoa, bacteria and DNA viruses. In this study, WSV067 was expressed inEscherichia coli, and the purified recombinant protein showed TS activity in dUMP−folate-binding assays using ultraviolet difference spectroscopy. RT-PCR and Western blot analyses showed that WSV067 was a genuine and early gene. Phylogenetic analysis revealed that WSSV-TS was more closely related to the TSs of eukaryotes than to those from prokaryotes.

Inhibitor Specificity via Protein Dynamics

Structure-based drug design of species-specific inhibitors generally exploits structural differences in proteins from different organisms. Here, we demonstrate how achieving specificity can be aided by targeting differences in the dynamics of proteins. Thymidylate synthase (TS) is a good target for anticancer agents and a potential target for antibacterial agents. Most inhibitors are folate-analogs that bind at the folate binding site and are not species specific. In contrast, alpha156 is not a folate-analog and is specific for bacterial TS; it has been shown crystallographically to bind in a nonconserved binding site. Docking calculations and crystal structure-based estimation of the essential dynamics of TSs from five different species show that differences in the dynamics of TSs make the active site more accessible to alpha156 in the prokaryotic than in the eukaryotic TSs and thereby enhance the specificity of alpha156.

The only active mutant of thymidylate synthase D169, a residue far from the site of methyl transfer, demonstrates the exquisite nature of enzyme specificity

Cysteine is the only variant of D169, a cofactor-binding residue in thymidylate synthase, that shows in vivo activity. The 2.4 A crystal structure of Escherichia coli thymidylate synthase D169C in a complex with dUMP and the antifolate CB3717 shows it to be an asymmetric dimer, with only one active site covalently bonded to dUMP. At the active site with covalently bound substrate, C169 S gamma adopts the roles of both carboxyl oxygens of D169, making a 3.6 A S...H[bond]N hydrogen bond to 3-NH of CB3717 and a 3.4 A water-mediated hydrogen bond to H212. Analogous hydrogen bonds formed during the enzyme reaction are important for cofactor binding and are postulated to contribute to catalysis. The C169 side chain is likely to be ionized, making it a better hydrogen bond acceptor than a neutral sulfhydryl group. At the second active site, C169 S gamma makes a shorter (3 A) hydrogen bond to the 3-NH of CB3717, CB3717 is approximately 1.5 A out of its binding site and there is no covalent bond between dUMP and the catalytic cysteine. Changes to partitioning among productive and non-productive conformations of reaction intermediates may contribute as much, if not more, to the diminished activity of this mutant than reduced stabilization of transition states.

Macropodid herpesvirus 1 encodes genes for both thymidylate synthase and ICP34.5

Macropodid herpesvirus 1 (MaHV-1) is an unclassified alphaherpesvirus linked with the fatal infections of kangaroos and other marsupials. During the characterisation of the internal repeat region of MaHV-1, an open reading frame (ORF) encoding for thymidylate synthase (TS) gene was identified and completely sequenced. Southern blot analysis confirmed the presence of two copies of the TS gene in the MaHV-1 genome as expected. Computer analysis of the TS ORF showed it was 948 nucleotides in length. A putative polyadenylation signal was identified 17-22 bp inside the ORF implying a minimal or absent 3' untranslated region. The predicted polypeptide was 316 amino acid residues in length and contained the highly conserved motifs for folate binding and F-dUMP binding, typical of all TS enzymes. Interestingly, MaHV-1 TS polypeptide had highest similarity to the human TS polypeptide (81%) compared to the TS polypeptides of other herpesviruses (72-75%). Immediately upstream of the TS gene, a second ORF of 510 bp, encoding a polypeptide with 170 amino acid residues, was identified. The carboxyl domain of this MaHV-1 polypeptide shared 68% similarity to a 59 amino acid motif of human herpesvirus 1 ICP34.5, identifying it as the MaHV-1 ICP34.5 homologue. This is the first report of a herpesvirus that encodes for both TS and ICP34.5.

Predicting and harnessing protein flexibility in the design of species-specific inhibitors of thymidylate synthase

BACKGROUND

Protein plasticity in response to ligand binding abrogates the notion of a rigid receptor site. Thus, computational docking alone misses important prospective drug design leads. Bacterial-specific inhibitors of an essential enzyme, thymidylate synthase (TS), were developed using a combination of computer-based screening followed by in-parallel synthetic elaboration and enzyme assay [Tondi et al. (1999) Chem. Biol. 6, 319-331]. Specificity was achieved through protein plasticity and despite the very high sequence conservation of the enzyme between species.

RESULTS

The most potent of the inhibitors synthesized, N,O-didansyl-L-tyrosine (DDT), binds to Lactobacillus casei TS (LcTS) with 35-fold higher affinity and to Escherichia coli TS (EcTS) with 24-fold higher affinity than to human TS (hTS). To reveal the molecular basis for this specificity, we have determined the crystal structure of EcTS complexed with DDT and 2'-deoxyuridine-5'-monophosphate (dUMP). The 2.0 A structure shows that DDT binds to EcTS in a conformation not predicted by molecular docking studies and substantially differently than other TS inhibitors. Binding of DDT is accompanied by large rearrangements of the protein both near and distal to the enzyme's active site with movement of C alpha carbons up to 6 A relative to other ternary complexes. This protein plasticity results in novel interactions with DDT including the formation of hydrogen bonds and van der Waals interactions to residues conserved in bacterial TS but not hTS and which are hypothesized to account for DDT's specificity. The conformation DDT adopts when bound to EcTS explains the activity of several other LcTS inhibitors synthesized in-parallel with DDT suggesting that DDT binds to the two enzymes in similar orientations.

CONCLUSIONS

Dramatic protein rearrangements involving both main and side chain atoms play an important role in the recognition of DDT by EcTS and highlight the importance of incorporating protein plasticity in drug design. The crystal structure of the EcTS/dUMP/DDT complex is a model system to develop more selective TS inhibitors aimed at pathogenic bacterial species. The crystal structure also suggests a general formula for identifying regions of TS and other enzymes that may be treated as flexible to aid in computational methods of drug discovery.

Crystal structure of a deletion mutant of human thymidylate synthase Delta (7-29) and its ternary complex with Tomudex and dUMP

The crystal structures of a deletion mutant of human thymidylate synthase (TS) and its ternary complex with dUMP and Tomudex have been determined at 2.0 A and 2.5 A resolution, respectively. The mutant TS, which lacks 23 residues near the amino terminus, is as active as the wild-type enzyme. The ternary complex is observed in the open conformation, similar to that of the free enzyme and to that of the ternary complex of rat TS with the same ligands. This is in contrast to Escherichia coli TS, where the ternary complex with Tomudex and dUMP is observed in the closed conformation. While the ligands interact with each other in identical fashion regardless of the enzyme conformation, they are displaced by about 1.0 A away from the catalytic cysteine in the open conformation. As a result, the covalent bond between the catalytic cysteine sulfhydryl and the base of dUMP, which is the first step in the reaction mechanism of TS and is observed in all ternary complexes of the E. coli enzyme, is not formed. This displacement results from differences in the interactions between Tomudex and the protein that are caused by differences in the environment of the glutamyl tail of the Tomudex molecule. Despite the absence of the closed conformation, Tomudex inhibits human TS ten-fold more strongly than E. coli TS. These results suggest that formation of a covalent bond between the catalytic cysteine and the substrate dUMP is not required for effective inhibition of human TS by cofactor analogs and could have implications for drug design by eliminating this as a condition for lead compounds.

Site-directed ligand discovery

We report a strategy (called "tethering") to discover low molecular weight ligands ( approximately 250 Da) that bind weakly to targeted sites on proteins through an intermediary disulfide tether. A native or engineered cysteine in a protein is allowed to react reversibly with a small library of disulfide-containing molecules ( approximately 1,200 compounds) at concentrations typically used in drug screening (10 to 200 microM). The cysteine-captured ligands, which are readily identified by MS, are among the most stable complexes, even though in the absence of the covalent tether the ligands may bind very weakly. This method was applied to generate a potent inhibitor for thymidylate synthase, an essential enzyme in pyrimidine metabolism with therapeutic applications in cancer and infectious diseases. The affinity of the untethered ligand (K(i) approximately 1 mM) was improved 3,000-fold by synthesis of a small set of analogs with the aid of crystallographic structures of the tethered complex. Such site-directed ligand discovery allows one to nucleate drug design from a spatially targeted lead fragment.

The crystal structure of thymidylate synthase from Pneumocystis carinii reveals a fungal insert important for drug design

Thymidylate synthase from Pneumocystis carinii (PcTS) is an especially important drug target, since P. carinii is a fungus that causes opportunistic pneumonia infections in immune-compromised patients and is among the leading causes of death of AIDS patients. Thymidylate synthase (TS) is the sole enzyme responsible for the de novo production of deoxythymidine monophosphate and hence is crucial for DNA replication in every organism. Inhibitors selective for P. carinii TS over human TS would be greatly beneficial in combating this disease. The crystal structure of TS from P. carinii bound to its substrate, dUMP, and a cofactor mimic, CB3717, was determined to 2.6 A resolution. A comparison with other species of TS shows that the volume of the closed PcTS active-site is 20 % larger than that of five other TS closed active-sites. A two-residue proline insert that is strictly conserved among all fungal species of TS, and a novel C-terminal closing interaction involving a P. carinii-specific tyrosine residue are primarily responsible for this increase in volume. The structure suggests several options for designing an inhibitor specific to PcTS and avoiding interactions with human TS. Taking advantage of the residue substitutions of P. carinii TS over human TS enables the design of a selective inhibitor. Additionally, the larger volume of the active-site of PcTS is an important advantage for designing de novo inhibitors that will exclude the human TS active-site through steric hindrance.

Genetic complementation and resistance to 5-fluoro-2'-deoxyuridine in thymidine auxotrophs expressing a highly defective mutant of human thymidylate synthase

A mutant human thymidylate synthase (TS) has been created in which a glutamine residue at position 214 has been replaced by glutamate. Glutamine at position 214 is postulated to be involved in maintaining the enzyme in a conformation that facilitates the binding of the substrate dUMP. Although the kcat/Km of the mutant protein for the substrate, dUMP, is 10(3) lower than that of wild-type TS, the mutant TS confers thymidine prototrophy on a TS-deficient bacterial strain when expressed at high levels. In the present investigation, a TS-deficient Chinese hamster lung cell line was transfected with DNA encoding the defective protein. Thymidine prototrophs were isolated that expressed the defective protein at levels that were physiologically relevant. The activities of the enzymes expressed endogenously in representative prototrophs were consistent with the activities observed for the purified proteins. At similar levels of TS expression, thymidine prototrophs expressing Glu214 TS were 8-fold more resistant to 5-fluoro-2'-deoxyuridine (FdUrd) cytotoxicity than are prototrophs expressing Gln214 TS. FdUrd is a prodrug of the tight-binding TS inhibitor, 5-fluoro-2'-deoxyuridine-5'-monophosphate (FdUMP). The resistance to FdUrd was associated with a significant decrease in the binding of FdUMP to the purified mutant enzyme. The data are consistent with the interpretation that TSs that are highly defective are capable of sufficient dTMP production for cell survival and optimal growth, yet may confer resistance to TS-directed inhibitors.

Kinetic studies on drug-resistant variants of Escherichia coli thymidylate synthase: functional effects of amino acid substitutions at residue 4

A naturally occurring mutant of human thymidylate synthase (hTS) that contains a Tyr to His mutation at residue 33 was found to confer 4-fold resistance to 5-fluoro-2'-deoxyuridine (FdUrd), a prodrug of 5-fluoro-2'-deoxyuridine 5'-monophosphate (FdUMP). The crystal structure of hTS implicated this Tyr residue in a drug resistance mechanistic role that may include both substrate binding and catalysis (Schiffer et al., Biochemistry, 34, 16279-16287, 1995). Because of the existence of a defined kinetic scheme and the development of a bacterial expression vector for the overproduction of Escherichia coli TS (ecTS), we chose to initially study the corresponding residue in the bacterial enzyme, Tyr 4 of ecTS. Nine mutant ecTS enzymes that differed in sequence at position 4 were generated. Mutants with a charged or polar side chain (Ser, Cys, Asp, and Arg) and Gly precipitated in the cell paste, resulting in no catalytic activity in cell-free extracts. Although most of the His 4 mutant precipitated, sufficient amounts remained in the cell-free extract to permit isolation to near homogeneity. Wild-type ecTS and mutants with a hydrophobic side chain (Phe, Ile, and Val) were expressed at nearly 30% of the total cellular protein. The k(cat) values for the isolatable mutants were 2- to 10-fold lower than that of the wild-type enzyme, while the K(m) values for 2'-deoxyuridylate (dUMP) and 5,10-methylenetetrahydrofolate (CH(2)H(4)folate) were similar for all the mutants. Dissociation constants for binary complex formation determined by stopped-flow spectroscopy were similar for the wild-type and mutant enzymes for both dUMP and 2'-deoxythymidylate, indicating that this mutation does not significantly alter the binding of the natural nucleotide ligands. However, each mutant enzyme had three- to 5-fold lower affinity for FdUMP in the binary complex compared with the wild-type enzyme, and only His 4 showed a lower affinity for FdUMP in the ternary complex. Analysis of k(burst) showed that the initial binding of CH(2)H(4)folate is weaker for each mutant compared to the wild-type enzyme and that lower k(cat) values were due to compromised rates that govern the chemical transformation of bound substrates to bound products.

Regulation of p53 expression by thymidylate synthase

Previous studies showed that thymidylate synthase (TS), as an RNA binding protein, regulates its own synthesis by impairing the translation of TS mRNA. In this report, we present evidence that p53 expression is affected in a similar manner by TS. For these studies, we used a TS-depleted human colon cancer HCT-C cell that had been transfected with either the human TS cDNA or the Escherichia coli TS gene. The level of p53 protein in transfected cells overexpressing human TS was significantly reduced when compared with its corresponding parent HCT-C cells. This suppression of p53 expression was the direct result of decreased translational efficiency of p53 mRNA. Similar results were obtained upon transfection of HCT-C cells with pcDNA 3.1 (+) containing the E. coli TS gene. These findings provide evidence that TS, from diverse species, specifically regulates p53 expression at the translational level. In addition, TS-overexpressing cells with suppressed levels of p53 are significantly impaired in their ability to arrest in G1 phase in response to exposure to a DNA-damaging agent such as gamma-irradiation. These studies provide support for the in vivo biological relevance of the interaction between TS and p53 mRNA and identify a molecular pathway for controlling p53 expression.

Functional assessment of surface loops: deletion of eukaryote-specific peptide inserts in thymidylate synthase of Saccharomyces cerevisiae

The primary structure is known for at least 29 thymidylate synthases and the crystal structure is known for several from both prokaryotes and eukaryotes. All these are markedly similar making thymidylate synthase one of the most highly conserved enzymes known. There are, however, two surface loops, one near the active site and the other near the dimer interface, which exist in distinctly prokaryotic and eukaryotic versions. Specifically, in eukaryotes these two surface loops have small peptide inserts conserved in size and partly conserved in sequence, that are not present in the prokaryotic thymidylate synthases. To address the possibility that these inserts provide eukaryote-specific functions the Saccharomyces cerevisiae loops were individually modified to mimic their prokaryotic counterparts. Altering the surface loop near the active site increased Km for the nucleotide substrate and decreased apparent Vmax. Mutant variants with alterations in the other surface loop were unable to dimerize. Therefore these surface loops have acquired, perhaps by way of the eukaryotic inserts, characteristics that are important for catalytic activity and quaternary structure respectively.

Identification of a thymidylate synthase gene within the genome of Chilo iridescent virus

The thymidylate synthase (TS, EC 2.1.1.45) is essential for the de novo synthesis of dTMP in pro- and eucaryotic organisms. Consequently it plays a major role in the replication of the DNA genome of a cell or a DNA virus. The gene encoding the TS of Chilo iridescent virus (CIV) was identified by nucleotide sequence analysis of the viral genome and was mapped within the EcoRI CIV DNA fragments G and R. Computer assisted analysis of the DNA nucleotide sequence between the genome coordinates 0.482 and 0.489 revealed an open reading frame (ORF) of 885 nucleotides. This ORF was found to encode a polypeptide of 295 amino acid residues (33.9 kDa) that showed significant homologies to known TS of different species including mammals, plants, fungi, protozoa, bacteria, and DNA viruses. The highest amino acid homologies were found between the CIV-TS and the TS of herpesvirus ateles (54.0%), Saccharomyces cerevisiae (51.8%), herpesvirus saimiri (51.0%), rhesus monkey rhadinovirus (50.7%), mouse (50.5%), rat (50.2%), varicella-zoster virus (50.2%), equine herpesvirus 2 (50.0%), and the human TS (48.4%). The CIV-TS contains six amino acid domains that are highly conserved in the TS of other species. Within these domains the major amino acid residues are present for which a functional role has been reported. The CIV-TS was found to be more closely related to the TS of eucaryotes than to the TS of procaryotes indicating the phylogenetic origin of the CIV-TS gene. The identification of a TS gene in the genome of CIV is the first report of a viral TS that is not encoded by a herpesvirus or a bacteriophage.

Inactivity of N229A thymidylate synthase due to water-mediated effects: isolating a late stage in methyl transfer

Mutation of thymidylate synthase N229(177) to alanine results in an essentially inactive enzyme, yet it leads to formation of a stable ternary complex. The kinetics of N229(177)A show that kcat for Escherichia coli is reduced by 200-fold while the Km for dUMP is increased 200-fold and the Km for folate increased by tenfold versus the wild-type enzyme. The crystal structures of N229(177)A in complex with dUMP and CB3717, and in complex with dUMP alone are determined at 2.4 A, and 2.5 A resolution. These structures identify the covalently bound ternary complex and show how N229(177)A traps an intermediate, and so becomes inactive in a later step of the reaction. Since the smaller alanine side-chain at N229(177)A does not directly sterically impair binding of ligands, the structures implicate, and place quantitative limits on the involvement of the structured water network in the active site of thymidylate synthase in both catalysis and in determining the binding affinity for dUMP (in contrast, the N229(177)V mutation in Lactobacillus casei has minimal effect on activity).

Metal ion-mediated substrate-assisted catalysis in type II restriction endonucleases

The 2.15-A resolution cocrystal structure of EcoRV endonuclease mutant T93A complexed with DNA and Ca2+ ions reveals two divalent metals bound in one of the active sites. One of these metals is ligated through an inner-sphere water molecule to the phosphate group located 3' to the scissile phosphate. A second inner-sphere water on this metal is positioned approximately in-line for attack on the scissile phosphate. This structure corroborates the observation that the pro-SP phosphoryl oxygen on the adjacent 3' phosphate cannot be modified without severe loss of catalytic efficiency. The structural equivalence of key groups, conserved in the active sites of EcoRV, EcoRI, PvuII, and BamHI endonucleases, suggests that ligation of a catalytic divalent metal ion to this phosphate may occur in many type II restriction enzymes. Together with previous cocrystal structures, these data allow construction of a detailed model for the pretransition state configuration in EcoRV. This model features three divalent metal ions per active site and invokes assistance in the bond-making step by a conserved lysine, which stabilizes the attacking hydroxide ion nucleophile.

Random sequence mutagenesis and resistance to 5-fluorouridine in human thymidylate synthases

Thymidylate synthase (TS) catalyzes the methylation of dUMP to dTMP and is the target for the widely used chemotherapeutic agent 5-fluorouracil. We used random sequence mutagenesis to replace 13 codons within the active site of TS and obtain variants that are resistant to 5-fluorodeoxyuridine (5-FdUR). The resulting random library was selected for its ability to complement a TS-deficient Escherichia coli strain, and sequence analysis of survivors found multiple substitutions to be tolerable within the targeted region. An independent selection of the library was carried out in the presence of 5-FdUR, resulting in a more limited spectrum of mutations. One specific mutation, C199L, was observed in more than 46% of 5-FdUR-resistant clones. A 5-FdUR-resistant triple mutant, A197V/L198I/C199F, was purified to apparent homogeneity. Kinetic studies with the substrate dUMP indicate that this mutant is similar to the wild type in regards to kcat and Km values for dUMP and the cosubstrate CH2H4-folate. In contrast, equilibrium binding studies with the inhibitor, FdUMP, demonstrate that the dissociation constant (Kd) for FdUMP binding into the ternary complex was 20-fold higher than values obtained for the wild-type enzyme. This 5-FdUMP-resistant mutant, or others similarly selected, is a candidate for use in gene therapy to render susceptible normal cells resistant to the toxic effects of systemic 5-fluorouracil.

Recognition of flanking DNA sequences by EcoRV endonuclease involves alternative patterns of water-mediated contacts

The 2.1-A cocrystal structure of EcoRV endonuclease bound to 5'-CGGGATATCCC, in a crystal lattice isomorphous with the cocrystallized undecamer 5'-AAAGATATCTT previously determined, shows novel base recognition in the major groove of the DNA flanking the GATATC target site. Lys104 of the enzyme interacts through water molecules with the exocyclic N-4 amino groups of flanking cytosines. Steric exclusion of water molecule-binding sites by the 5-methyl group of thymine drives the adoption of alternative water-mediated contacts with AT versus GC flanks. This structure provides a rare example of structural adaptability in the recognition of different DNA sequences by a protein and suggests preferred strategies for the expansion of target site specificity by EcoRV.

The additivity of substrate fragments in enzyme-ligand binding

BACKGROUND

Enzymes have evolved to recognise their target substrates with exquisite selectivity and specificity. Whether fragments of the substrate--perhaps never available to the evolving enzyme--are bound in the same manner as the parent substrate addresses the fundamental basis of specificity. An understanding of the relative contributions of individual portions of ligand molecules to the enzyme-binding interaction may offer considerable insight into the principles of substrate recognition.

RESULTS

We report 12 crystal structures of Escherichia coli thymidylate synthase in complexes with available fragments of the substrate (dUMP), both with and without the presence of a cofactor analogue. The structures display considerable fidelity of binding mode and interactions. These complexes reveal several interesting features: the cofactor analogue enhances the localisation of substrate and substrate fragments near the reactive thiol; the ribose moiety reduces local disorder through additional specific enzyme-ligand interactions; the pyrimidine has multiple roles, ranging from stereospecificity to mechanistic competence; and the glycosidic linkage has an important role in the formation of a covalent attachment between substrate and enzyme.

CONCLUSIONS

The requirements of ligand-protein binding can be understood in terms of the binding of separate fragments of the ligand. Fragments which are subsystems of the natural substrate for the enzyme confer specific contributions to the binding affinity, orientation or electrostatics of the enzymatic mechanism. This ligand-binding analysis provides a complementary method to the more prevalent approaches utilising site-directed mutagenesis. In addition, these observations suggest a modular approach for rational drug design utilising chemical fragments.

Flexible ligand docking using conformational ensembles

Molecular docking algorithms suggest possible structures for molecular complexes. They are used to model biological function and to discover potential ligands. A present challenge for docking algorithms is the treatment of molecular flexibility. Here, the rigid body program, DOCK, is modified to allow it to rapidly fit multiple conformations of ligands. Conformations of a given molecule are pre-calculated in the same frame of reference, so that each conformer shares a common rigid fragment with all other conformations. The ligand conformers are then docked together, as an ensemble, into a receptor binding site. This takes advantage of the redundancy present in differing conformers of the same molecule. The algorithm was tested using three organic ligand protein systems and two protein-protein systems. Both the bound and unbound conformations of the receptors were used. The ligand ensemble method found conformations that resembled those determined in X-ray crystal structures (RMS values typically less than 1.5 A). To test the method's usefulness for inhibitor discovery, multi-compound and multi-conformer databases were screened for compounds known to bind to dihydrofolate reductase and compounds known to bind to thymidylate synthase. In both cases, known inhibitors and substrates were identified in conformations resembling those observed experimentally. The ligand ensemble method was 100-fold faster than docking a single conformation at a time and was able to screen a database of over 34 million conformations from 117,000 molecules in one to four CPU days on a workstation.

Contributions of orientation and hydrogen bonding to catalysis in Asn229 mutants of thymidylate synthase

We have determined structures of binary and ternary complexes of five Asn229 variants of thymidylate synthase (TS) and related their structures to the kinetic constants measured previously. Asn229 forms two hydrogen bonds to the pyrimidine ring of the substrate 2'-deoxyuridine-5'-monophosphate (dUMP). These hydrogen bonds constrain the orientation of dUMP in binary complexes with dUMP, and in ternary complexes with dUMP and the TS cofactor, 5,10-methylene-5,6,7,8-tetrahydrofolate. In N229 mutants, where these hydrogen bonds cannot be made, dUMP binds in a misoriented or more disordered fashion. Most N229 mutants exhibit no activity for the dehalogenation of 5-bromo-dUMP, which requires correct orientation of dUMP against Cys198. Since bound dUMP forms the binding surface against which the pterin ring of cofactor binds, misorientation of dUMP results in higher Km values for cofactor. At the same time, binding of the cofactor aids in ordering and positioning dUMP for catalysis. Hydrophobic mutants, such as N229I, favor an arrangement of solvent molecules and side-chains around the ligands similar to that in a proposed transition state for ternary complex formation in wild-type TS, and kcat values are similar to the wild-type value. Smaller, more hydrophilic mutants favor arrangements of the solvent and side-chains surrounding the ligands that do not resemble the proposed transition state. These changes correspond to decreases in kcat of up to 2000-fold, with only modest increases in Km or Kd. These results are consistent with the proposal that the hydrogen-bonding network between water, dUMP and side-chains in the active-site cavity contributes to catalysis in TS. Asn229 has the unique ability to maintain this critical network, without sterically interfering with dUMP binding.

Thymidylate synthase inhibition: a structure-based rationale for drug design

Thymidylate synthase (TS) is a very interesting target in antiproliferative diseases. Its inhibition causes thimineless death of the cells and compounds inhibiting TS are widely used in anticancer therapy. The classical antifolate TS inhibitors are structural analogs of the folate cofactor; they often share the same metabolic pathways and this causes the development of resistance inside the cells. A detailed analysis of the available x-ray crystal structures of the complexes of the enzyme with different substrates and inhibitors support the finding of a structural basis of their biological activity. TS inhibitors nonstructural analog of folate, non-analog antifolate inhibitors (NAAI), are welcome as a new interesting research topic. Among the most recent and interesting ones, compounds from Agouron related to the indole structure, are independent on the folate metabolism, highly active and specific for human TS. Other compounds, phthalein derivatives, can inhibit TS enzymes from various sources and show an interesting biological activity profile: they inhibit better bacterial and fungal TS than human TS. The x-ray crystal structures of some of these inhibitors with TS show that they bind in a different binding site from that of the classical folate TS inhibitors. This indicates a potential allosteric binding site useful for future drug discovery studies.

Binding of biotin to streptavidin stabilizes intersubunit salt bridges between Asp61 and His87 at low pH

The remarkable stability of the streptavidin tetramer towards subunit dissociation becomes even greater upon binding of biotin. At two equivalent extensive monomer-monomer interfaces, monomers tightly associate into dimers that in turn associate into the tetramer at a less extensive dimer-dimer interface. To probe the structural basis for the enhancement of the stability of streptavidin by biotin, the crystal structures of apostreptavidin and its complexes with biotin and other small molecule and cyclic peptide ligands were determined and compared at resolutions as high as 1.36 A over a range of pH values from as low as 1.39. At low pH dramatic changes occur in the conformation and intersubunit hydrogen bonds involving the loop comprising Asp61 to Ser69. The hydrogen-bonded salt bridge between Asp61 Odelta2 and His87 Ndelta1, observed at higher pH, is replaced with a strong hydrogen bond between Asp61 Odelta1 and Asn85 Odelta1. Through crystallography at multiple pH values, the pH where this conformational change occurs, and thus the pKa of Asp61, was determined in crystals of space group I222 and/or I4122 of apostreptavidin and complexes. A range in pKa values for Asp61 was observed in these structures, the lowest being 1.78+/-0.19 for I222 streptavidin-biotin in 2.9 M (NH4)2SO4. At low pH the decrease in pKa of Asp61 and preservation of the intersubunit Asp61 Odelta2-Ndelta1 His87 hydrogen-bonded salt bridge in streptavidin-biotin versus apostreptavidin or streptavidin-peptide complexes is associated with an ordering of the flexible flap comprising residues Ala46 to Glu51, that in turn orders the Arg84 side-chain of a neighboring loop through resulting hydrogen bonds. Ordering of Arg84 in close proximity to the strong intersubunit interface appears to stabilize the conformation associated with the Asp61 Odelta2-Ndelta1 His87 hydrogen-bonded salt bridge. Thus, in addition to the established role of biotin in tetramer stabilization by direct mediation of intersubunit interactions at the weak interface through contact with Trp120, biotin may enhance tetramer stability at the strong interface more indirectly by ordering loop residues.