Computer-aided gene design.

A computer program, which runs on MS-DOS personal computers, is described that assists in the design of synthetic genes coding for proteins. The goal of the program is the design of a gene which (i) contains as many unique restriction sites as possible and (ii) uses a specific codon usage. The gene designed according to the criteria above is (i) suitable for 'modular mutagenesis' experiments and (ii) optimized for expression. The program 'reverse-translates' protein sequences into degenerated DNA sequences, generates a map of potential restriction sites and locates sequence positions where unique restriction sites can be accommodated. The nucleic acid sequence is then 'refined' according to a specific codon usage to remove any degeneration. Unique restriction sites, if potentially present, can be 'forced' into the degenerated nucleic acid sequence by using 'priority codes' assigned to different restriction sequences.     

Creating new restriction sites by silent changes in coding sequences

We present methods for identifying a useful type of DNA site--one that can be mutated to create a new restriction site within a coding region without changing the amino acid sequence. These "latent sites" are abundant--silent mutations creating one of 44 different 6-bp or 8-bp recognition sites were found at relatively high density, roughly one latent site per 9 bp, in the eleven genes tested. Our analysis suggests that site-directed mutagenesis can be used to refashion coding sequences at will for flexible analysis.     

A set of Macintosh computer programs for the design and analysis of synthetic genes

Computer programs that can be used for the design of syntheticgenes and that are run on an Apple Macintosh computer are described.These programs determine nucleic acid sequences encoding aminoacid sequences. They select DNA sequences based on codon usageas specified by the user, and determine the placement of basechanges that can be used to create restriction enzyme siteswithout altering the amino acid sequence. A new algorithm forfinding restriction sites by translating the restriction endonucleasetarget sequence in all three reading frames and then searchingthe given peptide or protein amino acid sequence with theseshort restriction enzyme peptide sequences is described. Examplesare given for the creation of synthetic DNA sequences for thebovine prethrombin-2 and ribonuclease A genes Received on October 18, 1988; accepted on December 9, 1988     

Chemical synthesis of a mitochondrial gene designed for expression in the yeast nucleus

A novel DNA sequence coding for subunit 8 of the mitochondrial ATPase of Saccharomyces cerevisiae has been constructed by chemical synthesis. The synthetic gene, termed NAP1, is designed for expression in the yeast nucleus and codes for a 48 amino acid polypeptide identical to that encoded by the mitochondrial aap1 gene of S. cerevisiae. The codons chosen for the NAP1 sequence correspond almost exclusively to those most frequently occurring in highly expressed yeast genes. The NAP1 coding region differs in 31 codons from that of aap1, and is flanked by sequences carrying restriction enzyme sites useful for cloning and for gene expression. A 170 bp double stranded DNA molecule was constructed by assembling 12 oligonucleotides (12 to 45 bases in length) in a single annealing/ligation mixture. This synthetic gene will provide a route for the systematic manipulation, through in vitro mutagenesis, of the structure of a protein normally encoded by mitochondrial DNA.     

Directed reconstruction of the influenza virus hemagglutinin gene

A new approach to create chimeric genes by directed exchange of oligonucleotide fragments was developed. By oligonucleotide-directed mutagenesis a few deletion mutants of the influenza virus hemagglutinin (HA) gene were obtained. These variants of HA gene contain unique restriction sites in DNA regions coding for the A and B epitopes of the HA molecule. The obtained special vectors may be used for cloning DNA fragments coding for new amino acid sequences in internal sites of the HA gene.     

Polymorphic Sites and the Mechanism of Evolution in Human Mitochondrial DNA

Twelve restriction enzymes were used to screen for the presence or absence of cleavage sites at 441 locations in the mitochondrial DNA of 112 humans from four continents. Cleavage maps were constructed by comparison of DNA fragment sizes with those expected from the published sequence for one human mtDNA. One hundred and sixty-three of the sites were polymorphic, i.e., present in some individuals but absent from others, 278 sites being invariant. These polymorphisms probably result from single base substitutions and occur in all functional regions of the genome.--In 77 cases, it was possible to specify the exact nature and location (within a restriction site) of the mutation responsible for the absence of a restriction site in a known human mtDNA sequence and its presence in another human mtDNA. Fifty-two of these 77 gain mutations occur in genes coding for proteins, 34 being silent and 18 causing amino acid replacements; moreover, nine of the replacements are radical.--Notable also is the anomalous ratio of transitions to transversions required to account for these 77 restriction site differences between the known human mtDNA sequences and other human mtDNAs. This ratio is lower for most groups of restriction sites than has been reported from sequence comparisons of limited parts of the mtDNA genome in closely related mammals, perhaps indicating a special functional role or sensitivity to mutagenesis for palindromic regions containing high levels of guanine and cytosine.--From the genomic distribution of the 163 polymorphic sites, it is inferred that the level of point mutational variability in tRNA and rRNA genes is nearly as high as in protein-coding genes but lower than in noncoding mtDNA. Thus, the functional constraints operating on components of the protein-synthetic apparatus may be lower for mitochondria than for other systems. Furthermore, the mitochondrial genes for tRNAs that recognize four codons are more variable than those recognizing only two codons.--Among the more variable of the human mitochondrial genes coding for proteins is that for subunit 2 of cytochrome oxidase; this polypeptide appears to have been evolving about five times faster in primates than in other mammals. Cytochrome c, a nuclearly encoded protein that interacts directly with the oxidase 2 subunit in electron transport, has also evolved faster in primates than in rodents or ungulates. This example, along with that for the mitochondrial rRNA genes and the nuclear genes coding for mitochondrial ribosomal proteins, provides evidence for coevolution between specific nuclear and mitochondrial genes.     

Heterotachy and functional shift in protein evolution

Study of structure/function relationships constitutes an important field of research, especially for modification of protein function and drug design. However, the fact that rational design (i.e. the modification of amino acid sequences by means of directed mutagenesis, based on knowledge of the three-dimensional structure) appears to be much less efficient than irrational design (i.e. random mutagenesis followed by in vitro selection) clearly indicates that we understand little about the relationships between primary sequence, three-dimensional structure and function. The use of evolutionary approaches and concepts will bring insights to this difficult question. The increasing availability of multigene family sequences that has resulted from genome projects has inspired the creation of novel in silico evolutionary methods to predict details of protein function in duplicated (paralogous) proteins. The underlying principle of all such approaches is to compare the evolutionary properties of homologous sequence positions in paralogs. It has been proposed that the positions that show switches in substitution rate over time--i.e., 'heterotachous sites'--are good indicators of functional divergence. However, it appears that heterotachy is a much more general process, since most variable sites of homologous proteins with no evidence of functional shift are heterotachous. Similarly, it appears that switches in substitution rate are as frequent when paralogous sequences are compared as when orthologous sequences are compared. Heterotachy, instead of being indicative of functional shift, may more generally reflect a less specific process related to the many intra- and inter-molecular interactions compatible with a range of more or less equally viable protein conformations. These interactions will lead to different constraints on the nature of the primary sequences, consistently with theories suggesting the non-independence of substitutions in proteins. However, a specific type of amino acid variation might constitute a good indicator of functional divergence: substitutions occurring at positions that are generally slowly evolving. Such substitutions at constrained sites are indeed much more frequent soon after gene duplication. The identification and analysis of these sites by complementing structural information with evolutionary data may represent a promising direction to future studies dealing with the functional characterization of an ever increasing number of multi-gene families identified by complete genome analysis.     

Isolation of two genomic sequences encoding the Mr = 14,000 subunit of rat prostatein

Complementary DNAs to rat ventral prostate poly(A) RNA were cloned into pBR322 by the "dG-dC tailing" procedure. Clones containing cDNAs to the mRNAs coding for each of the three subunits of a major secretory protein (prostatein) were identified by hybrid-arrested translation. A 457-nucleotide base pair cDNA (E45) and a portion of a 365-base pair cDNA (E85) were analyzed to determine the composite complete DNA coding sequence for the Mr = 14,000 (C3) subunit of prostatein. A sequence of 12-nucleotide bases (TTTGCTGCTATG) in the signal peptide of C3 was noted to be homologous to signal peptide nucleotide sequences reported in cDNAs coding for the other two prostatein subunits, Mr = 6,000 (C1) and 10,000 (C2). Complementary DNA coding for the C3 subunit was used as a hybridization probe to screen an EcoRI rat genomic DNA library. Two unique 12-kilobase genomic clones, each containing mRNA coding sequences within 2.5-3-kilobase fragments, were identified by restriction enzyme mapping and Southern blot analysis. Restriction enzyme sites within the coding regions of both genes were analogous to the cDNA. Differences in restriction enzyme sites in regions of intervening sequences and flanking DNA established the uniqueness of the two genes. It is suggested that both genes may be transcribed in vivo.     

Use of synthetic DNA in expression of foreign proteins

Synthetic DNAs and oligonucleotides, which can be prepared conveniently by combining chemical synthesis and enzymatic methods, have been used extensively in recombinant DNA research. Examples include total gene synthesis, probes for the isolation of specific genes from cDNA or genomic libraries, linkers containing specific restriction sites for cloning, primers for DNA and RNA sequencing, and primers for the construction of specific mutations (either deletion, insertion or point mutations) by oligonucleotide-directed site-specific mutagenesis.This article reviews recent advances in the chemical and enzymatic synthesis of oligo- and polynucleotides and the application of synthetic DNA to the expression of foreign proteins. The synthesis of genes, including structural genes and regulatory genes are reviewed. Oligonucleotide-directed site-specific mutagenesis and use of synthetic DNA to optimize foreign protein expression are also discussed.     

Plasmid pU29, a vehicle for mutagenesis of the photosynthetic puf operon in Rhodopseudomonas capsulata

Plasmid pU21, which carries the reaction center and light-harvesting genes (puf operon) of Rhodopseudomonas capsulata, has been redesigned by site-specific mutagenesis. Five restriction sites have been removed and three unique restriction sites have been introduced into this 11,589-bp pBR322 derivative. The modifications divide the puf structural genes into four regions separated by five unique and nonmutagenic restriction sites. These four fragments have been subcloned into the M13-mp series of vectors to facilitate oligonucleotide-mediated site-specific mutagenesis experiments on the photosynthetic apparatus structural genes. The inserts can then be returned from the M13 replicative form to the redesigned pU21 derivative. The modified plasmid, pU29, greatly facilitates in vitro mutagenesis experiments since previously described techniques and screening procedures are more efficient with M13 derivatives carrying smaller inserts. Additionally, tandem homologous sequences (the reaction center L and M subunits) within the puf operon are now separated on different phage vectors, eliminating problems encountered in the targeting of mutagenic oligonucleotides to only one of the two homologous sites.     

Chemical synthesis of a biologically active gene for human immune interferon-gamma. Prospect for site-specific mutagenesis and structure-function studies

A 453-base pair DNA duplex consisting of a gene coding for human interferon-gamma and initiation and termination signals plus appropriate restriction enzyme sites for plasmid insertion has been totally synthesized. The synthesis involved preparation of 66 oligodeoxynucleotides by a modified, solid phase phosphite procedure and enzymatic ligation of the oligonucleotides. The gene, when inserted into a previously constructed expression vector, was expressed in Escherichia coli, demonstrating functional activity for the synthetic gene. Several strategically located restriction cleavage sites have been introduced into the sequence. This provides a convenient system for site-specific mutagenesis for structure-function studies.     

Intraspecific evolution of a gene family coding for urinary proteins

Summary The genome of the laboratory mouse contains about 35 major urinary protein (MUP) genes, many of which are clustered on chromosome 4. We have used distance and parsimony methods to estimate phylogenetic relationships between MUP genes from nucleotide sequence and restriction maps. By analyzing coding sequences we show that the genes fall into four main groups of related sequences (groups 1–4). Comparisons of restriction maps and the nucleotide sequences of hypervariable regions that lie 50 nucleotides 5 to the cap sites show that the group 1 genes and probably also the group 2 pseudogenes fall into subgroups. The most parsimonious trees are consistent with the evolution of the array of group 1 and 2 genes by mutation accompanied by a process tending toward homogenization such as unequal crossing-over or gene conversion. The phylogenetic grouping correlates with grouping according to aspects of function. The genomes of the inbred strains BALB/c and C57BL contain different MUP gene arrays that we take to be samples from the wild population of arrays.     

Isolation and characterization of maize genes coding for zein proteins of the 21000 dalton size class

Cloning in lambda gt WES phage of EcoRI fragments from maize seedlings DNA led to the isolation of four fragments containing genes coding for 21000 dalton zein proteins. The zein genes, identified by electron microscopic analysis, do not contain intervening sequences detectable by this method. The flanking sequences were analyzed by restriction sites mapping and hybridization and showed areas of homology between each other and with sequences surrounding a previously isolated gene of the 19000 dalton zein class.     

SILMUT: a computer program for the identification of regions suitable for silent mutagenesis to introduce restriction enzyme recognition sequences.

We describe a set of IBM-compatible computer programs designed to selectively identify the potential sites for silent mutagenesis within a target DNA sequence. This program is based on a novel strategy of identifying amino acid motifs compatible with each restriction site (BioTechniques 12:382-384, 1991). The programs can be used to identify the suitability for the introduction of any 6-base nucleic acid sequences, such as restriction enzyme sites in cassette mutagenesis strategies. The Table program generates a table of multiple amino acid motifs for each restriction enzyme, obtained by translating each unique recognition sequence in all three reading frames. The Silmut program, which utilizes the features of Table, will further identify the presence of a match between any amino acid motif of each restriction enzyme and the input target sequence. Minor manipulations of the data base files will enable the individual researcher to identify the potential for introduction of any 6-base sequences by silent mutagenesis.     

A Macintosh computer program for designing DNA sequences that code for specific peptides and proteins

A computer program (PINCERS) is described for use in the design of synthetic genes and mixed-probe DNA sequences. A protein sequence is reverse translated with generation of synonymous codons at each position producing a degenerate sequence. In order to locate potential restriction enzyme sites, the degenerate sequence is searched with a library of restriction enzymes for sites that utilize any combination of synonymous codons. These sites are indicated in a map so that they may be incorporated into the synthetic gene sequence. The program allows the user to select the appropriate codon usage table for the organism of interest and then to set a threshold usage frequency below which codons are not generated. PINCERS may also be used to assist in planning the synthesis of mixed-probe DNA sequences for cross-hybridization experiments. It can identify regions of specified length with the protein sequence that have the least overall degeneracy, thereby minimizing the number of probes to be synthesized and, therefore, maximizing the concentration of a given probe sequence.     

A rapid and efficient PCR-based mutagenesis method applicable to cell physiology study

PCR-based mutagenesis is a cornerstone of molecular biology and protein engineering studies. Herein we describe a rapid and highly efficient mutagenesis method using type IIs restriction enzymes. A template gene is amplified into two separate PCR fragments using two pairs of anchor and mutagenic primers. Mutated sequences are located near the recognition site of a type IIs restriction enzyme. After digestion of two fragments with a type IIs enzyme, exposed cohesive ends that are complementary to each other are then ligated together to generate a mutated gene. We applied this method to introduce multiple site-directed mutations in EGFP and Bcl-2 family genes and observed perfect mutagenesis efficiency at the desired sites. This efficient and cost-effective mutagenesis method can be applied to a wide variety of structural and functional studies in cell physiology. Type IIs restriction enzyme; enhanced green fluorescent protein; Bcl-2     

insilico.mutagenesis: a primer selection tool designed for sequence scanning applications used in directed evolution experiments

Several primer prediction programs have been developed for a variety of applications. However none of these tools allows the prediction of a large set of primers for whole gene site-directed mutagenesis experiments using the megaprimer method. We report a novel primer prediction tool (insilico.mutagenesis), accessible at www.insilico.uni-duesseldorf.de, developed for the application to high-throughput mutagenesis used in directed evolution or structure-function dependency projects, which involve the subsequent mutagenesis of a large number of amino acid positions (e.g., in whole gene saturation or gene scanning mutagenesis experiments). Furthermore, the program is suitable for all site-directed (saturation) mutagenesis approaches, such as saturation mutagenesis of promoter sequences and other types of untranslated intergenic regions. In anticipation of downstream cloning steps, the primer design tool also includes a restriction site control feature alerting the user if unwanted restriction sites have been introduced within the mutagenesis primer. The use of our tool promises to speed up the process of site-directed mutagenesis, as it instantly allows predicting a large set of primers.     

Oligonucleotide-directed mutagenesis by microscale 'shot-gun' gene synthesis.

We describe a rapid and efficient microscale method for in vitro site-directed mutagenesis by gene synthesis. Mutants are constructed by "shot-gun ligation" of overlapping synthetic oligonucleotides yielding double stranded synthetic DNA of more than 120 nucleotides in length. The terminal oligonucleotides of the DNA segment to be synthesized are designed to create sticky ends complementary to unique restriction sites of a polylinker present in an M13 vector. The oligonucleotides are hybridized and ligated to the M13 vector without any purification of the synthetic DNA segment. After cloning, about half of the progeny from such shot-gun ligations contained the predicted sequence demonstrating the efficacy of this method for gene synthesis and its potential for the extensive mutational analysis of genes.     

Synthesis, cloning and expression of recombinant aprotinin

Synthetic DNA fragments containing the coding sequence for the serine proteinase inhibitor aprotinin, also known as bovine pancreatic trypsin inhibitor (BPTI) a Kunitz type inhibitor were fused to form a synthetic aprotinin gene by the method of Khorana and cloned into E. coli. The synthetic gene is characterized by the presence of certain restriction sites. These restriction sites are unique within the used cloning system. Therefore, a great number of modifications can be achieved easily by exchange of appropriate restriction fragments. Using this method the variant [Glu52]aprotinin was obtained starting from the aprotinin gene. Both genes were successfully expressed in E. coli as fusion proteins with beta-galactosidase using vector pUR 278. No translation products could be detected in four other expression system (pUR 108, pDR 540, pKK 223-3 and pUC 8). [Glu52]aprotinin was purified and renatured after cyanogen bromide cleavage of the fusion protein. This recombinant [Glu52]aprotinin shows exactly the same trypsin-inhibitory profile as natural aprotinin.