Escherichia coli thymidylate synthase: amino acid substitutions by suppression of amber nonsense mutations

The genetic landscape of a metabolic interaction

While much prior work has explored the constraints on protein sequence and evolution induced by physical protein-protein interactions, the sequence-level constraints emerging from non-binding functional interactions in metabolism remain unclear. To quantify how variation in the activity of one enzyme constrains the biochemical parameters and sequence of another, we focus on dihydrofolate reductase (DHFR) and thymidylate synthase (TYMS), a pair of enzymes catalyzing consecutive reactions in folate metabolism. We use deep mutational scanning to quantify the growth rate effect of 2696 DHFR single mutations in 3 TYMS backgrounds under conditions selected to emphasize biochemical epistasis. Our data are well-described by a relatively simple enzyme velocity to growth rate model that quantifies how metabolic context tunes enzyme mutational tolerance. Together our results reveal the structural distribution of epistasis in a metabolic enzyme and establish a foundation for the design of multi-enzyme systems.

Functional Prokaryotic-Like Deoxycytidine Triphosphate Deaminases and Thymidylate Synthase in Eukaryotic Social Amoebae: Vertical, Endosymbiotic, or Horizontal Gene Transfer?

The de novo synthesis of deoxythymidine triphosphate uses several pathways: gram-negative bacteria use deoxycytidine triphosphate deaminase to convert deoxycytidine triphosphate into deoxyuridine triphosphate, whereas eukaryotes and gram-positive bacteria instead use deoxycytidine monophosphate deaminase to transform deoxycytidine monophosphate to deoxyuridine monophosphate. It is then unusual that in addition to deoxycytidine monophosphate deaminases, the eukaryote Dictyostelium discoideum has 2 deoxycytidine triphosphate deaminases (Dcd1Dicty and Dcd2Dicty). Expression of either DcdDicty can fully rescue the slow growth of an Escherichia coli dcd knockout. Both DcdDicty mitigate the hydroxyurea sensitivity of a Schizosaccharomyces pombe deoxycytidine monophosphate deaminase knockout. Phylogenies show that Dcd1Dicty homologs may have entered the common ancestor of the eukaryotic groups of Amoebozoa, Obazoa, Metamonada, and Discoba through an ancient horizontal gene transfer from a prokaryote or an ancient endosymbiotic gene transfer from a mitochondrion, followed by horizontal gene transfer from Amoebozoa to several other unrelated groups of eukaryotes. In contrast, the Dcd2Dicty homologs were a separate horizontal gene transfer from a prokaryote or a virus into either Amoebozoa or Rhizaria, followed by a horizontal gene transfer between them. ThyXDicty, the D. discoideum thymidylate synthase, another enzyme of the deoxythymidine triphosphate biosynthesis pathway, was suggested previously to be acquired from the ancestral mitochondria or by horizontal gene transfer from alpha-proteobacteria. ThyXDicty can fully rescue the E. coli thymidylate synthase knockout, and we establish that it was obtained by the common ancestor of social amoebae not from mitochondria but from a bacterium. We propose horizontal gene transfer and endosymbiotic gene transfer contributed to the enzyme diversity of the deoxythymidine triphosphate synthesis pathway in most social amoebae, many Amoebozoa, and other eukaryotes.

The Genetic Landscape of a Metabolic Interaction

Enzyme abundance, catalytic activity, and ultimately sequence are all shaped by the need of growing cells to maintain metabolic flux while minimizing accumulation of deleterious intermediates. While much prior work has explored the constraints on protein sequence and evolution induced by physical protein-protein interactions, the sequence-level constraints emerging from non-binding functional interactions in metabolism remain unclear. To quantify how variation in the activity of one enzyme constrains the biochemical parameters and sequence of another, we focused on dihydrofolate reductase (DHFR) and thymidylate synthase (TYMS), a pair of enzymes catalyzing consecutive reactions in folate metabolism. We used deep mutational scanning to quantify the growth rate effect of 2,696 DHFR single mutations in 3 TYMS backgrounds under conditions selected to emphasize biochemical epistasis. Our data are well-described by a relatively simple enzyme velocity to growth rate model that quantifies how metabolic context tunes enzyme mutational tolerance. Together our results reveal the structural distribution of epistasis in a metabolic enzyme and establish a foundation for the design of multi-enzyme systems.

Trans-editing by aminoacyl-tRNA synthetase-like editing domains

Aminoacyl-tRNA synthetases (aaRS) are ubiquitous enzymes responsible for aminoacyl-tRNA (aa-tRNA) synthesis. Correctly formed aa-tRNAs are necessary for proper decoding of mRNA and accurate protein synthesis. tRNAs possess specific nucleobases that promote selective recognition by cognate aaRSs. Selecting the cognate amino acid can be more challenging because all amino acids share the same peptide backbone and several are isosteric or have similar side chains. Thus, aaRSs can misactivate non-cognate amino acids and produce mischarged aa-tRNAs. If left uncorrected, mischarged aa-tRNAs deliver their non-cognate amino acid to the ribosome resulting in misincorporation into the nascent polypeptide chain. This changes the primary protein sequence and potentially causes misfolding or formation of non-functional proteins that impair cell survival. A variety of proofreading or editing pathways exist to prevent and correct mistakes in aa-tRNA formation. Editing may occur before the amino acid transfer step of aminoacylation via hydrolysis of the aminoacyl-adenylate. Alternatively, post-transfer editing, which occurs after the mischarged aa-tRNA is formed, may be carried out via a distinct editing site on the aaRS where the mischarged aa-tRNA is deacylated. In recent years, it has become clear that most organisms also encode factors that lack aminoacylation activity but resemble aaRS editing domains and function to clear mischarged aa-tRNAs in trans. This review focuses on these trans-editing factors, which are encoded in all three domains of life and function together with editing domains present within aaRSs to ensure that the accuracy of protein synthesis is sufficient for cell survival.

Experimental and Computational Studies Delineate the Role of Asparagine 177 in Hydride Transfer for E. coli Thymidylate Synthase

Thymidylate synthase (TSase), an enzyme responsible for the de novo biosynthesis of 2'-deoxythymidine 5'-monophosphate (thymidylate, dTMP) necessary for DNA synthesis, has been a drug target for decades. TSase is a highly conserved enzyme across species ranging from very primitive organisms to mammals. Among the many conserved active site residues, an asparagine (N177, using Escherichia coli residues numbering) appears to make direct hydrogen bonds with both the C4=O4 carbonyl of the 2'-deoxyuridine 5'-monophosphate (uridylate, dUMP) substrate and its pyrimidine ring's N3. Recent studies have reassessed the TSase catalytic mechanism, focusing on the degree of negative charge accumulation at the O4 carbonyl of the substrate during two critical H-transfers - a proton abstraction and a hydride transfer. To obtain insights into the role of this conserved N177 on the hydride transfer, we examined its aspartic acid (D) and serine (S) mutants - each of which is expected to alter hydrogen bonding and charge stabilization around the C4=O4 carbonyl of the 2'-deoxyuridine 5'-monophosphate (uridylate, dUMP) substrate. Steady-state kinetics, substrate binding order studies and temperature-dependency analysis of intrinsic KIEs for the hydride transfer step of the TSase catalytic cycle suggest the active site of N177D is not precisely organized for that step. A smaller disruption was observed for N177S, which could be rationalized by partial compensation by water molecules and rearrangement of other residues toward preparation of the system for the hydride transfer under study. These experimental findings are qualitatively mirrored by QM/MM computational simulations, thereby shedding light on the sequence and synchronicity of steps in the TSase-catalyzed reaction. This information could potentially inform the design of mechanism-based drugs targeting this enzyme.

Implementation of a loss-of-function system to determine growth and stress-associated mutagenesis in Bacillus subtilis

A forward mutagenesis system based on the acquisition of mutations that inactivate the thymidylate synthase gene (TMS) and confer a trimethoprim resistant (Tmpr) phenotype was developed and utilized to study transcription-mediated mutagenesis (TMM). In addition to thyA, Bacillus subtilis possesses thyB, whose expression occurs under conditions of cell stress; therefore, we generated a thyB- thyA+ mutant strain. Tmpr colonies of this strain were produced with a spontaneous mutation frequency of ~1.4 × 10-9. Genetic disruption of the canonical mismatch (MMR) and guanine oxidized (GO) repair pathways increased the Tmpr frequency of mutation by ~2-3 orders of magnitude. A wide spectrum of base substitutions as well as insertion and deletions in the ORF of thyA were found to confer a Tmpr phenotype. Stationary-phase-associated mutagenesis (SPM) assays revealed that colonies with a Tmpr phenotype, accumulated over a period of ten days with a frequency of ~ 60 ×10-7. The Tmpr system was further modified to study TMM by constructing a ΔthyA ΔthyB strain carrying an IPTG-inducible Pspac-thyA cassette. In conditions of transcriptional induction of thyA, the generation of Tmpr colonies increased ~3-fold compared to conditions of transcriptional repression. Further, the Mfd and GreA factors were necessary for the generation of Tmpr colonies in the presence of IPTG in B. subtilis. Because GreA and Mfd facilitate transcription-coupled repair, our results suggest that TMM is a mechanim to produce genetic diversity in highly transcribed regions in growth-limited B. subtilis cells.

Hotspots in an obligate homodimeric anticancer target. Structural and functional effects of interfacial mutations in human thymidylate synthase

Human thymidylate synthase (hTS), a target for antiproliferative drugs, is an obligate homodimer. Single-point mutations to alanine at the monomer-monomer interface may enable the identification of specific residues that delineate sites for drugs aimed at perturbing the protein-protein interactions critical for activity. We computationally identified putative hotspot residues at the interface and designed mutants to perturb the intersubunit interaction. Dimer dissociation constants measured by a FRET-based assay range from 60 nM for wild-type hTS up to about 1 mM for single-point mutants and agree with computational predictions of the effects of these mutations. Mutations that are remote from the active site retain full or partial activity, although the substrate KM values were generally higher and the dimer was less stable. The lower dimer stability of the mutants can facilitate access to the dimer interface by small molecules and thereby aid the design of inhibitors that bind at the dimer interface.

Modeling functional changes to Escherichia coli thymidylate synthase upon single residue replacements: a structure-based approach

Escherichia coli thymidylate synthase (TS) is an enzyme that is indispensable to DNA synthesis and cell division, as it provides the only de novo source of dTMP by catalyzing the reductive methylation of dUMP, thus making it a key target for chemotherapeutic agents. High resolution X-ray crystallographic structures are available for TS and, owing to its relatively small size, successful experimental mutagenesis studies have been conducted on the enzyme. In this study, an in silico mutagenesis technique is used to investigate the effects of single amino acid substitutions in TS on enzymatic activity, one that employs the TS protein structure as well as a knowledge-based, four-body statistical potential. For every single residue TS variant, this approach yields both a global structural perturbation score and a set of local environmental perturbation scores that characterize the mutated position as well as all structurally neighboring residues. Global scores for the TS variants are capable of uniquely characterizing groups of residue positions in the enzyme according to their physicochemical, functional, or structural properties. Additionally, these global scores elucidate a statistically significant structure–function relationship among a collection of 372 single residue TS variants whose activity levels have been experimentally determined. Predictive models of TS variant activity are subsequently trained on this dataset of experimental mutants, whose respective feature vectors encode information regarding the mutated position as well as its six nearest residue neighbors in the TS structure, including their environmental perturbation scores.

Genetic determinants involved in p-aminosalicylic acid resistance in clinical isolates from tuberculosis patients in northern China from 2006 to 2012

p-Aminosalicylic acid (PAS) is an important compound for treating multidrug-resistant tuberculosis (TB). Previous studies showed that thyA mutations are often related to PAS resistance in clinical isolates. We performed a systematic analysis of isolate genotypes and detected mutations in three folate pathway genes (folC, thyA, and ribD) in 61.1% (127/208) of PAS-resistant isolates, including 11 double mutants. This result expands our knowledge about the distribution and frequency of mutations related to PAS resistance in mycobacterial clinical isolates.

Inhibitor-bound complexes of dihydrofolate reductase-thymidylate synthase from Babesia bovis

Babesiosis is a tick-borne disease caused by eukaryotic Babesia parasites which are morphologically similar to Plasmodium falciparum, the causative agent of malaria in humans. Like Plasmodium, different species of Babesia are tuned to infect different mammalian hosts, including rats, dogs, horses and cattle. Most species of Plasmodium and Babesia possess an essential bifunctional enzyme for nucleotide synthesis and folate metabolism: dihydrofolate reductase-thymidylate synthase. Although thymidylate synthase is highly conserved across organisms, the bifunctional form of this enzyme is relatively uncommon in nature. The structural characterization of dihydrofolate reductase-thymidylate synthase in Babesia bovis, the causative agent of babesiosis in livestock cattle, is reported here. The apo state is compared with structures that contain dUMP, NADP and two different antifolate inhibitors: pemetrexed and raltitrexed. The complexes reveal modes of binding similar to that seen in drug-resistant malaria strains and point to the utility of applying structural studies with proven cancer chemotherapies towards infectious disease research.

pFARs, plasmids free of antibiotic resistance markers, display high-level transgene expression in muscle, skin and tumour cells

BACKGROUND

Nonviral gene therapy requires a high yield and a low cost production of eukaryotic expression vectors that meet defined criteria such as biosafety and quality of pharmaceutical grade. To fulfil these objectives, we designed a novel antibiotic-free selection system.

METHODS

The proposed strategy relies on the suppression of a chromosomal amber mutation by a plasmid-borne function. We first introduced a nonsense mutation into the essential Escherichia coli thyA gene, resulting in thymidine auxotrophy. The bacterial strain was optimized for the production of small and novel plasmids free of antibiotic resistance markers (pFARs) and encoding an amber suppressor t-RNA. Finally, the potentiality of pFARs as eukaryotic expression vectors was assessed by monitoring luciferase activities after electrotransfer of LUC-encoding plasmids into various tissues.

RESULTS

The introduction of pFARs into the optimized bacterial strain restored normal growth to the auxotrophic mutant and allowed an efficient production of monomeric supercoiled plasmids. The electrotransfer of LUC-encoding pFAR into muscle led to high luciferase activities, demonstrating an efficient gene delivery. In transplanted tumours, transgene expression levels were superior after electrotransfer of the pFAR derivative compared to a plasmid carrying a kanamycin resistance gene. Finally, in skin, whereas luciferase activities decreased within 3 weeks after intradermal electrotransfer of a conventional expression vector, sustained luciferase expression was observed with the pFAR plasmid.

CONCLUSIONS

Thus, we have designed a novel strategy for the efficient production of biosafe plasmids and demonstrated their potentiality for nonviral gene delivery and high-level transgene expression in several tissues.

Quality control despite mistranslation caused by an ambiguous genetic code

A high level of accuracy during protein synthesis is considered essential for life. Aminoacyl-tRNA synthetases (aaRSs) translate the genetic code by ensuring the correct pairing of amino acids with their cognate tRNAs. Because some aaRSs also produce misacylated aminoacyl-tRNA (aa-tRNA) in vivo, we addressed the question of protein quality within the context of missense suppression by Cys-tRNA(Pro), Ser-tRNA(Thr), Glu-tRNA(Gln), and Asp-tRNA(Asn). Suppression of an active-site missense mutation leads to a mixture of inactive mutant protein (from translation with correctly acylated aa-tRNA) and active enzyme indistinguishable from the wild-type protein (from translation with misacylated aa-tRNA). Here, we provide genetic and biochemical evidence that under selective pressure, Escherichia coli not only tolerates the presence of misacylated aa-tRNA, but can even require it for growth. Furthermore, by using mass spectrometry of a reporter protein not subject to selection, we show that E. coli can survive the ambiguous genetic code imposed by misacylated aa-tRNA tolerating up to 10% of mismade protein. The editing function of aaRSs to hydrolyze misacylated aa-tRNA is not essential for survival, and the EF-Tu barrier against misacylated aa-tRNA is not absolute. Rather, E. coli copes with mistranslation by triggering the heat shock response that stimulates nonoptimized polypeptides to achieve a native conformation or to be degraded. In this way, E. coli ensures the presence of sufficient functional protein albeit at a considerable energetic cost.

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.

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.

The Bacterial YbaK Protein Is a Cys-tRNAPro andCys-tRNACysDeacylase

Bacterial prolyl-tRNA synthetases and some smaller paralogs, YbaK and ProX, can hydrolyze misacylated Cys-tRNA Pro or Ala-tRNA Pro. To assess the significance of this quality control editing reaction in vivo, we tested Escherichia coli ybaK for its ability to suppress the E. coli thymidylate synthase thyA:146CCA missense mutant strain, which requires Cys-tRNA(Pro) for growth in the absence of thymine. Missense suppression was observed in a ybaK deletion background, suggesting that YbaK functions as a Cys-tRNA Pro deacylase in vivo. In vitro studies with the full set of 20 E. coli aminoacyl-tRNAs revealed that the Haemophilus influenzae and E. coli YbaK proteins are moderately general aminoacyl-tRNA deacylases that preferentially hydrolyze Cys-tRNA Pro and Cys-tRNA Cys and are also weak deacylases that cleave Gly-tRNA, Ala-tRNA, Ser-tRNA, Pro-tRNA, and Met-tRNA. The ProX protein acted as an aminoacyl-tRNA deacylase that cleaves preferentially Ala-tRNA and Gly-tRNA. The potential of H. influenzae YbaK to hydrolyze in vivo correctly charged Cys-tRNA Cys was tested in E. coli strain X2913 (ybaK+). Overexpression of H. influenzae ybaK decreased the in vivo ratio of Cys-tRNA Cys to tRNA Cys from 65 to 35% and reduced the growth rate of strain X2913 by 30% in LB medium. These data suggest that YbaK-mediated hydrolysis of aminoacyl-tRNA can influence cell growth.

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.

An alternative flavin-dependent mechanism for thymidylate synthesis

Although deoxythymidylate cannot be provided directly by ribonucleotide reductase, the gene encoding thymidylate synthase ThyA is absent from the genomes of a large number of nonsymbiotic microbes. We show that ThyX (Thy1) proteins of previously unknown function form a large and distinct class of thymidylate synthases. ThyX has a wide but sporadic phylogenetic distribution, almost exclusively limited to microbial genomes lacking thyA. ThyX and ThyA use different reductive mechanisms, because ThyX activity is dependent on reduced flavin nucleotides. Our findings reveal complexity in the evolution of thymidine in present-day DNA. Because ThyX proteins are found in many pathogenic microbes, they present a previously uncharacterized target for antimicrobial compounds.

Expanding the genetic code of Escherichia coli

A unique transfer RNA (tRNA)/aminoacyl-tRNA synthetase pair has been generated that expands the number of genetically encoded amino acids in Escherichia coli. When introduced into E. coli, this pair leads to the in vivo incorporation of the synthetic amino acid O-methyl-l-tyrosine into protein in response to an amber nonsense codon. The fidelity of translation is greater than 99%, as determined by analysis of dihydrofolate reductase containing the unnatural amino acid. This approach should provide a general method for increasing the genetic repertoire of living cells to include a variety of amino acids with novel structural, chemical, and physical properties not found in the common 20 amino acids.

Genetic analysis of the Escherichia coli FtsZ.ZipA interaction in the yeast two-hybrid system. Characterization of FtsZ residues essential for the interactions with ZipA and with FtsA

The recruitment of ZipA to the septum by FtsZ is an early, essential step in cell division in Escherichia coli. We have used polymerase chain reaction-mediated random mutagenesis in the yeast two-hybrid system to analyze this interaction and have identified residues within a highly conserved sequence at the C terminus of FtsZ as the ZipA binding site. A search for suppressors of a mutation that causes a loss of interaction (ftsZ(D373G)) identified eight different changes at two residues within this sequence. In vitro, wild type FtsZ interacted with ZipA with a high affinity in an enzyme-linked immunosorbent assay, whereas FtsZ(D373G) failed to interact. Two mutant proteins examined restored this interaction significantly. In vivo, the alleles tested are significantly more toxic than the wild type ftsZ and cannot complement a deletion. We have shown that a fusion, which encodes the last 70 residues of FtsZ in the two-hybrid system, is sufficient for the interaction with FtsA and ZipA. However, when the wild type sequence is compared with one that encodes FtsZ(D373G), no interaction was seen with either protein. Mutations surrounding Asp-373 differentially affected the interactions of FtsZ with ZipA and FtsA, indicating that these proteins bind the C terminus of FtsZ differently.

Replacement set mutagenesis of the four phosphate-binding arginine residues of thymidylate synthase

Arginines R23, R178, R179 and R218 in thymidylate synthase (TS, EC 2. 1.1.45) are hydrogen bond donors to the phosphate moiety of the substrate, dUMP. In order to investigate how these arginines contribute to enzyme function, we prepared complete replacement sets of mutants at each of the four sites in Lactobacillus casei TS. Mutations of R23 increase K:(m) for dUMP 2-20-fold, increase K:(m) for cofactor 8-40-fold and decrease k(cat) 9-20-fold, reflecting the direct role of the R23 side chain in binding and orienting the cofactor in ternary complexes of the enzyme. Mutations of R178 increase K:(m) for dUMP 40-2000-fold, increase K:(m) for cofactor 3-20-fold and do not significantly affect k(cat). These results are consistent with the fact that this residue is an integral part of the dUMP-binding wall and contributes to the orientation and ordering of several other dUMP binding residues. Kinetic parameters for all R179 mutations except R179P were not significantly different from wild-type values, reflecting the fact that this external arginine does not directly contact the cofactor or other ligand-binding residues. R218 is essential for the structure of the catalytic site and all mutations of this arginine except R218K were inactive.

Isolation and characterization of a thymidylate synthase-deficient human colon tumor cell line

Following mutagenesis of the human colorectal tumor cell line HCT C with ethyl methanesulfonate, clonal sublines were isolated that survived on medium toxic to cells expressing thymidylate synthase (TS). The subline exhibiting the lowest TS activity, designated as C18, was characterized. Extracts from C18 cells were mixed with extracts from parental C cells to determine whether the TS-deficient phenotype is trans-acting. No effect was observed on the activity of TS in parental extracts. The levels of functional TS in C18 cells were analyzed by the binding of the mechanism-based inhibitor 5-fluoro-2'-deoxyuridylate (FdUMP) under conditions that allowed for the detection of 10 fmol of TS. Only a low level of FdUMP-TS complexes was detected in C18 extracts. The level of TS expression in C18 cells was similar to that in parental C cells, as indicated by immunoblot and RNA analyses. DNA sequence analysis of TS cDNA from C18 cells revealed the existence of a point mutation (C-->T) at nucleotide 647 that predicts the replacement of Ser216 by a leucine residue. That the C18 cell line was homozygous for this mutation was indicated by restriction fragment-length polymorphism analysis and by primer extension analysis. To provide additional evidence that substitution of Ser216 by a leucine residue created a defective protein, a TS-deficient bacterial strain was transformed with an expression vector containing the mutated human TS cDNA. The transformed strain exhibited thymidine auxotrophy, indicating that the mutant TS (Leu216) is nonfunctional.

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.

Reassigning cysteine in the genetic code of Escherichia coli

We investigated directed deviations from the universal genetic code. Mutant tRNAs that incorporate cysteine at positions corresponding to the isoleucine AUU, AUC, and AUA and methionine AUG codons were introduced in Escherichia coli K12. Missense mutations at the cysteine catalytic site of thymidylate synthase were systematically crossed with synthetic suppressor tRNACys genes coexpressed from compatible plasmids. Strains harboring complementary codon/anticodon associations could be stably propagated as thymidine prototrophs. A plasmid-encoded tRNACys reading the codon AUA persisted for more than 500 generations in a strain requiring its suppressor activity for thymidylate biosynthesis, but was eliminated from a strain not requiring it. Cysteine miscoding at the codon AUA was also enforced in the active site of amidase, an enzyme found in Helicobacter pylori and not present in wild-type E. coli. Propagating the amidase missense mutation in E. coli with an aliphatic amide as nitrogen source required the overproduction of Cys-tRNA synthetase together with the complementary suppressor tRNACys. The toxicity of cysteine miscoding was low in all our strains. The small size and amphiphilic character of this amino acid may render it acceptable as a replacement at most protein positions and thus apt to overcome the steric and polar constraints that limit evolution of the genetic code.

Isolation and characterization of thymitaq (AG337) and 5-fluoro-2-deoxyuridylate-resistant mutants of human thymidylate synthase from ethyl methanesulfonate-exposed human sarcoma HT1080 cells

Thymidylate synthase plays an essential role in the synthesis of DNA. Recently, several new and specific thymidylate synthase inhibitors that occupy the folate binding site, including Tomudex(R), BW1843U89, and Thymitaq, have demonstrated therapeutic activity in patients with advanced cancer. In order to find drug-resistant forms of human thymidylate synthase for gene therapy applications, human sarcoma HT1080 cells were exposed to ethyl methanesulfonate and Thymitaq selection. Thymitaq-resistant clonal derived sublines were established, and analysis indicated that both gene amplification and point mutations contributed to drug resistance. Eight mutant cDNAs that were identified from Thymitaq-resistant sublines were generated by site-directed mutagenesis and transfected into thymidylate synthase-negative cells. Only K47E, D49G, or G52S mutants retain enzyme activity. Moreover, cytotoxicity studies demonstrated that D49G and G52S transfected cells, besides displaying resistance to Thymitaq with IC50 values 40- and 12-fold greater than wild-type enzyme transfected cells, respectively, also lead to fluorodeoxyuridine resistance (26- and 97-fold in IC50 values, respectively) but not to Tomudex or BW1843U89. Characterization of the purified altered enzymes obtained from expression in Escherichia coli is consistent with the cell growth inhibition results. We postulate that the D49G or G52S mutation leads to the structural perturbation of the highly conserved Arg50 loop, decreasing the binding of thymidylate synthase to the inhibitors, Thymitaq and fluorodeoxyuridylate.

Crystal structures of a marginally active thymidylate synthase mutant, Arg 126-->Glu

Thymidylate synthase (TS) is a long-standing target for anticancer drugs and is of interest for its rich mechanistic features. The enzyme catalyzes the conversion of dUMP to dTMP using the co-enzyme methylenetetrahydrofolate, and is perhaps the best studied of enzymes that catalyze carbon-carbon bond formation. Arg 126 is found in all TSs but forms only 1 of 13 hydrogen bonds to dUMP during catalysis, and just one of seven to the phosphate group alone. Despite this, when Arg 126 of TS from Escherichia coli was changed to glutamate (R126E), the resulting protein had kcat reduced 2000-fold and Km reduced 600-fold. The crystal structure of R126E was determined under two conditions--in the absence of bound ligand (2.4 A resolution), and with dUMP and the antifolate CB3717 (2.2 A resolution). The first crystals, which did not contain dUMP despite its presence in the crystallization drop, displayed Glu 126 in a position to sterically and electrostatically interfere with binding of the dUMP phosphate. The second crystals contained both dUMP and CB3717 in the active site, but Glu 126 formed three hydrogen bonds to nearby residues (two through water) and was in a position that partially overlapped with the normal phosphate binding site, resulting in a approximately 1 A shift in the phosphate group. Interestingly, the protein displayed the typical ligand-induced conformational change, and the covalent bond to Cys 146 was present in one of the protein's two active sites.

Conversion of thymidylate synthase into an HIV protease substrate

Thymidylate synthase (TS) is an essential enzyme of DNA metabolism. We have carried out an extensive insertional mutagenesis of the Escherichia coli TS gene (thyA) using three different methods. Insertion of exogenous sequences at unique restriction sites or at random positions produced defective mutants, whereas comparison of TS sequences from different species allowed us to identify six zones permissive for insertions of exogenous sequences. The insertion of Human immunodeficiency virus type 1 (HIV-1) protease substrate sequences into the permissive sites converted TS to an HIV-1 protease substrate, and the in vivo cleavage of these insertions by the cloned HIV-1 protease conferred a thymidylate synthase-deficient phenotype in some of our E. coli mutant strains. In agreement with crystallographic data, these results show that the permissive sites are located in regions of the TS protein not essential for enzyme activity and accessible to cleavage by HIV protease. These results also show that it is possible to control a growth phenotype in E. coli through the protease-mediated destruction of an essential metabolic enzyme. Because both wild type and thymidylate synthase-deficient phenotypes are selectable on the appropriate growth medium, these thyA mutants could be used for genetic selections of protease inhibitors and analysis of protease specificities.

Site-directed mutagenesis of varicella-zoster virus thymidylate synthase. Analysis of two highly conserved regions of the enzyme

We have constructed a series of mutants to study the role of two structurally and functionally important regions of thymidylate synthase (TS) from varicella-zoster virus (VZV). The first centres on a conserved glycine residue in the beta-kink of beta-strand i, a partially buried region of the protein that is important for dimer interactions and the formation of the active site. We show that the glycine residue located in beta-strand i is not essential for enzyme activity and that beta-strand i can readily accommodate several amino acid substitutions and also an insertion. A covariant residue that accommodates these changes was also identified. The second region of interest was the solvent-exposed and highly mobile C-terminal residue which is an essential component of the active site in TS from Lactobacillus casei and Escherichia coli. We demonstrate that removal of the C-terminal residue from VZV TS does not completely inactivate the enzyme, implying that there are significant structural differences between the virus and bacterial enzymes. By combining site-directed mutagenesis and molecular modelling we have identified these differences and propose a model that explains the contrasting activities.

Nicotinic receptor binding site probed with unnatural amino acid incorporation in intact cells

The nonsense codon suppression method for unnatural amino acid incorporation has been applied to intact cells and combined with electrophysiological analysis to probe structure-function relations in the nicotinic acetylcholine receptor. Functional receptors were expressed in Xenopus oocytes when tyrosine and phenylalanine derivatives were incorporated at positions 93, 190, and 198 in the binding site of the alpha subunit. Subtle changes in the structure of an individual side chain produced readily detectable changes in the function of this large channel protein. At each position, distinct features of side chain structure dominated the dose-response relation, probably by governing the agonist-receptor binding.

Towards engineering proteins by site-directed incorporation in vivo of non-natural amino acids

Altering protein structure via the techniques of protein engineering has already allowed the development of proteins displaying both modified and novel activities. The only limitation of conventional site-directed mutagenesis, the cornerstone of protein engineering, is that substitutions are restricted to the 20 naturally occurring, proteinogenic amino acids. However, the discovery of a 21st amino acid, selenocysteine, and the development of novel in vitro translation systems have demonstrated that considerably more substitutions are possible. To this end, a number of experimental approaches have been developed that allow the incorporation of synthetic amino acids into proteins. Some of these have already been successfully applied in vitro and efforts to transfer this technology to in vivo systems are now underway.

Expression, purification, and characterization of thymidylate synthase from Lactococcus lactis

The thymidylate synthase (TS) gene from Lactococcus lactis has been highly expressed in Escherichia coli. The TS protein was purified by sequential chromatography on Q-Sepharose and phenyl-Sepharose. Six grams of cell pellet yielded 140 mg of homogeneous TS. TS is a highly conserved enzyme, and several of the conserved amino acid residues that have been implicated in catalytic function are altered in L. lactis TS. By use of a 3-dimensional homology model, we have predicted covariant changes that might compensate for these differences. With the large amounts of L. lactis TS now available, studies can be pursued to understand the structure-function relationships of this enzyme compared to other TSs and to confirm the presumed roles of the compensatory changes predicted in the homology model.

Functions of the gene products of Escherichia coli

A list of currently identified gene products of Escherichia coli is given, together with a bibliography that provides pointers to the literature on each gene product. A scheme to categorize cellular functions is used to classify the gene products of E. coli so far identified. A count shows that the numbers of genes concerned with small-molecule metabolism are on the same order as the numbers concerned with macromolecule biosynthesis and degradation. One large category is the category of tRNAs and their synthetases. Another is the category of transport elements. The categories of cell structure and cellular processes other than metabolism are smaller. Other subjects discussed are the occurrence in the E. coli genome of redundant pairs and groups of genes of identical or closely similar function, as well as variation in the degree of density of genetic information in different parts of the genome.

Asparagine 229 in thymidylate synthase contributes to, but is not essential for, catalysis

The conserved Asn-229 (N229) of thymidylate synthase (TS, EC 2.1.1.45) provides the only side chain that directly hydrogen bonds with the pyrimidine ring of the substrate dUMP. The carboxamide moiety forms a cyclic hydrogen bond network with the NH-3 and O-4 of the base and is a prime candidate for assisting proton-transfer reactions that occur at O-4 of the pyrimidine ring of dUMP. A complete replacement set of mutants at position 229 of Lactobacillus casei TS (N229 mutants) has been prepared, purified, and characterized. Fifteen of the 19 TS mutants were catalytically active. Steady-state kinetic parameters of N229 mutants varied 17- and 115-fold in the Km values for 5,10-methylene-5,6,7,8-tetrahydrofolate and dUMP, respectively, 1000-fold in kcat values, and 10,000-fold in kcat/Km values. Wild-type TS possesses lower Km and higher kcat and kcat/Km values than any of the TS N229 mutants. We conclude that N229 contributes to, but is not essential for, binding and catalysis. When the wild-type enzyme was not considered, there were excellent correlations between log kcat and the hydrophobicity of the side chains at position 229, in which the more hydrophobic side chains showed higher values. Our results suggest a unique interaction between N229 and the substrates that seems important in appropriately positioning the uracil heterocycle for catalysis. We propose that in the absence of N229, the electrophilic catalyst that transfers protons to the O-4 and stabilizes enol intermediates is a highly conserved molecule of water.

Recognition of tRNA(Cys) by Escherichia coli cysteinyl-tRNA synthetase

A study of the recognition of tRNA(Cys) by Escherichia coli cysteinyl-tRNA synthetase using in vivo and in vitro methods was performed. All three anticodon nucleotides, the discriminator nucleotide (73), and some elements within the tertiary domain (the D stem/loop, the T psi C stem/loop, and the variable loop) are important for recognition; the anticodon stem and acceptor stem appear to contain no essential elements. A T7 RNA polymerase transcript corresponding to tRNA(Cys) is only a 5.5-fold worse substrate than native tRNA(Cys) (in terms of the specificity constant, kcat/Km), mainly due to an increase in the value of Km for the transcript. The greatest loss of specificity caused by mutation of a single nucleotide occurs when the discriminator U73 is changed; kcat/Km declines 3-4 orders of magnitude depending on the substitution. Mutations in the wobble nucleotide of the anticodon also cause reductions in the specificity constant of 3 orders of magnitude, while mutations in the other anticodon nucleotides caused lesser effects. Interestingly, a C35A mutation (with the phenylalanine anticodon GAA) had no effect on aminoacylation by the cysteinyl-tRNA synthetase. Several amber suppressor tRNAs were constructed whose in vivo identity did not correlate with their in vitro specificity, indicating the need for both types of experiments to understand the factors which maintain tRNA specificity.

Asn177 in Escherichia coli thymidylate synthase is a major determinant of pyrimidine specificity

The substrate preference of recombinant Escherichia coli thymidylate synthase (TS) has been altered from 2'-deoxyuridylate (dUMP) to 2'-deoxycytidylate (dCMP) by site-directed mutagenesis of the codon for Asn177, which was changed to aspartic acid. The side-chain amide of Asn177 forms hydrogen bonds with O4 and N3 of dUMP bound to the crystalline enzyme [Montfort, W. R., Perry, K. M., Fauman, E. B., Finer-Moore, J. S., Maley, G. F., Hardy, L., Maley, F. & Stroud, R. M. (1990) Biochemistry 29, 6964-6977]. This Asn is invariant in all natural sequences for TS known. The values of kcat for the mutant enzyme, TS(N177D), with dCMP and dUMP are, respectively, 0.09 and 0.002 times the value of kcat of wild-type TS with dUMP as substrate. TS(N177D) turns over dCMP at 35 times its rate of dUMP turnover, whereas wild-type TS turns over dCMP at < 10(-5) of its rate of dUMP turnover. Thus Asn177 is a major determinant of the pyrimidine nucleotide specificity of TS. The mutant enzyme, like wild-type TS, forms a covalent complex with 5-fluoro-dUMP in the presence of 5,10-methylenetetrahydrofolate. TS(N177D) also has a newly acquired ability to be transiently inactivated by dUMP. This time-dependent inactivation requires the presence of methylenetetrahydrofolate and may be due to the accumulation of the enzyme in the form of a catalytic intermediate. The likely mechanistic basis for discrimination by TS between dUMP and dCMP is their differing requirements for charge stabilization during covalent catalysis.

Amino acid substitution analysis of E. coli thymidylate synthase: the study of a highly conserved region at the N-terminus

Amino acid substitution analysis within a highly conserved region of Escherichia coli thymidylate synthase (TS), using suppression of amber mutations by tRNA suppressors, has yielded a bank of 124 new mutationally altered TS proteins. These mutant proteins have been used to study the structure-function relationship of the Escherichia coli TS protein at the N-terminus corresponding to residues 20 through 35. This region contains a block of amino acids whose sequence has been well conserved among other known TS proteins from various organisms. Positions 20 through 25 contain a surface loop structure and positions 26 through 35 encompass a beta-strand. We find that residues surrounding a beta-bulge structure within the beta-strand are particularly sensitive to amino acid substitution, suggesting that this structure is maintained by a highly ordered packing arrangement. Three residues in the surface loop that are present at the base of the substrate binding pocket are also sensitive to amino acid substitution. The remainder of the conserved sites, including those at the dimer interface, are tolerant to most, if not all, of the substitutions tested.

Mutation of asparagine 229 to aspartate in thymidylate synthase converts the enzyme to a deoxycytidylate methylase

The conserved Asn 229 of thymidylate synthase (TS) forms a cyclic hydrogen bond network with the 3-NH and 4-O of the nucleotide substrate dUMP. The Asn 229 to Asp mutant of Lactobacillus casei thymidylate synthase (TS N229D) has been prepared, purified, and investigated. Steady-state kinetic parameters of TS N229D show 3.5- and 10-fold increases in the Km values of CH2H4folate and dUMP, respectively, and a 1000-fold decrease in kcat. Most important, the Asp 229 mutation changes the substrate specificity of TS to an enzyme which recognizes and methylates dCMP in preference to dUMP. With TS N229D the Km for dCMP is bout 3-fold higher than for dUMP, and the Km for CH2H4folate is increased about 5-fold; however, the kcat for dCMP methylation is 120-fold higher than that for dUMP methylation. Specificity for dCMP versus dUMP, as measured by kcat/Km, changes from negligible with wild-type TS to about a 40-fold increase with TS N229D. TS N229D reacts with CH2H4folate and FdUMP or FdCMP to form ternary complexes which are analogous to the TS-FdUMP-CH2H4folate complex. From what is known of the mechanism and structure of TS, the dramatic change in substrate specificity of TS N229D is proposed to involve a hydrogen bond network between Asp 229 and the 3-N and 4-NH2 of the cytosine heterocycle, causing protonation of the 3-N and stabilization of a reactive imino tautomer. A similar mechanism is proposed for related enzymes which catalyze one-carbon transfers to cytosine heterocycles.

Substitutions of a cysteine conserved among DNA cytosine methylases result in a variety of phenotypes

The proposed mechanism for DNA (cytosine-5)-methyltransferases envisions a key role for a cysteine residue. It is expected to form a covalent link with carbon 6 of the target cytosine, activating the normally inactive carbon 5 for methyl transfer. There is a single conserved cysteine among all DNA (cytosine-5)-methyltransferases making it the candidate nucleophile. We have changed this cysteine to other amino acids for the EcoRII methylase; which methylates the second cytosine in the sequence 5'-CCWGG-3'. Mutants were tested for their methyl transferring ability and for their ability to form covalent complexes with DNA. The latter property was tested indirectly with the use of a genetic assay involving sensitivity of cells to 5-azacytidine. Replacement of the conserved cysteine with glycine, valine, tryptophan or serine led to an apparent loss of methyl transferring ability. Interestingly, cells carrying the mutant with serine did show sensitivity to 5-azacytidine, suggesting the ability to link to DNA. Unexpectedly, substitution of the cysteine with glycine results in the inhibition of cell growth and the mutant allele can be maintained in the cells only when it is poorly expressed. These results suggest that the conserved cysteine in the EcoRII methylase is essential for methylase action and it may play more than one role in it.

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

The structure of thymidylate synthase (TS) from Escherichia coli was solved from cubic crystals with a = 133 A grown under reducing conditions at pH 7.0, and refined to R = 22% at 2.1 A resolution. The structure is compared with that from Lactobacillus casei solved to R = 21% at 2.3 A resolution. The structures are compared using a difference distance matrix, which identifies a common core of residues that retains the same relationship to one another in both species. After subtraction of the effects of a 50 amino acid insert present in Lactobacillus casei, differences in position of atoms correlate with temperature factors and with distance from the nearest substituted residue. The dependence of structural difference on thermal factor is parameterized and reflects both errors in coordinates that correlate with thermal factor, and the increased width of the energy well in which atoms of high thermal factor lie. The dependence of structural difference on distance from the nearest substitution also depends on thermal factors and shows an exponential dependence with half maximal effect at 3.0 A from the substitution. This represents the plastic accommodation of the protein which is parameterized in terms of thermal B factor and distance from a mutational change.