Strategies and applications of in vitro mutagenesis

Random Mutagenesis as a Promising Tool for Microalgal Strain Improvement towards Industrial Production

Microalgae have become a promising novel and sustainable feedstock for meeting the rising demand for food and feed. However, microalgae-based products are currently hindered by high production costs. One major reason for this is that commonly cultivated wildtype strains do not possess the robustness and productivity required for successful industrial production. Several strain improvement technologies have been developed towards creating more stress tolerant and productive strains. While classical methods of forward genetics have been extensively used to determine gene function of randomly generated mutants, reverse genetics has been explored to generate specific mutations and target phenotypes. Site-directed mutagenesis can be accomplished by employing different gene editing tools, which enable the generation of tailor-made genotypes. Nevertheless, strategies promoting the selection of randomly generated mutants avoid the introduction of foreign genetic material. In this paper, we review different microalgal strain improvement approaches and their applications, with a primary focus on random mutagenesis. Current challenges hampering strain improvement, selection, and commercialization will be discussed. The combination of these approaches with high-throughput technologies, such as fluorescence-activated cell sorting, as tools to select the most promising mutants, will also be discussed.

DNA gap repair in Escherichia coli for multiplex site-directed mutagenesis

Site-directed mutagenesis allows the generation of novel DNA sequences that can be used for a variety of important applications such as the functional analysis of genetic variants. To overcome the limitations of existing site-directed mutagenesis approaches, we explored in vivo DNA gap repair. We found that site-specific mutations in plasmid DNA can be generated in Escherichia coli using mutant single-stranded oligonucleotides to target PCR-derived linear double-stranded plasmid DNA. We called this method DeGeRing, and we characterized its advantages, including non-biased multiplex mutagenesis, over existing site-directed mutagenesis methods such as recombineering (recombination-mediated genetic engineering), single DNA break repair (SDBR, introduced by W. Mandecki), and QuikChange (Agilent Technologies, La Jolla, CA). We determined the efficiency of DeGeRing to induce site-directed mutations with and without a phenotype in three K-12 E coli strains using multiple single-stranded oligonucleotides containing homological and heterological parts of various sizes. Virtual lack of background made the isolation of mutants with frequencies up to 10^-6 unnecessary. Our data show that endogenous DNA gap repair in E coli supports efficient multiplex site-directed mutagenesis. DeGeRing might facilitate the generation of mutant DNA sequences for protein engineering and the functional analysis of genetic variants in reverse genetics.

An Efficient Approach for Two Distal Point Site-Directed Mutagenesis from Randomly Ligated PCR Products

Site-directed mutagenesis is one of the most important tools in molecular biology. The majority of the mutagenesis methods have been developed to mutate one region of target DNA in each cycle of mutagenesis, while in some cases there is a need to mutate several distal points. We used a new method to simultaneously mutate two distal points in the target DNA. Different regions of the target DNA were amplified in three separate PCR reactions. The PCR products were back-to-back and together they made the complete length of the template DNA. Mutations were introduced to PCR products by middle mutagenic primers. PCR products were mixed and ligated with random blunt ligation, and then the desired mutated DNA fragments were selected in two steps by flanking restriction enzyme digestion and size selection. Selected fragments were amplified in another PCR reaction using flanking primers and finally cloned into the plasmid vector. This mutagenesis process is simple, there is no need to use modified primers and long or difficult PCR reactions.

Random Scanning Mutagenesis

In vitro oligonucleotide and polymerase chain reaction (PCR)-based mutagenesis is generally used for altering the nucleotide sequence of genes to study their functional importance and the products they encode. A thorough approach to this problem is to systematically change each successive amino acid residue in the protein to alanine (i.e., alanine-scanning mutagenesis) or to a limited number of alternative amino acids. Although these strategies can provide useful information, it is sometimes desirable to test a broader spectrum of amino acid changes at the targeted positions. "Random scanning mutagenesis" was developed to examine the functional importance of individual amino acid residues in the conserved structural motif of human immunodeficiency virus (HIV) reverse transcriptase, and this protocol is adapted from that method. This strategy is an oligonucleotide-based method for generating all 19 possible replacements at individual amino acid sites within a protein.

Finding the Needle in the Haystack-the Use of Microfluidic Droplet Technology to Identify Vitamin-Secreting Lactic Acid Bacteria

Efficient screening technologies aim to reduce both the time and the cost required for identifying rare mutants possessing a phenotype of interest in a mutagenized population. In this study, we combined a mild mutagenesis strategy with high-throughput screening based on microfluidic droplet technology to identify Lactococcus lactis variants secreting vitamin B₂ (riboflavin). Initially, we used a roseoflavin-resistant mutant of L. lactis strain MG1363, JC017, which secreted low levels of riboflavin. By using fluorescence-activated droplet sorting, several mutants that secreted riboflavin more efficiently than JC017 were readily isolated from the mutagenesis library. The screening was highly efficient, and candidates with as few as 1.6 mutations per million base pairs (Mbp) were isolated. The genetic characterization revealed that riboflavin production was triggered by mutations inhibiting purine biosynthesis, which is surprising since the purine nucleotide GTP is a riboflavin precursor. Purine starvation in the mutants induced overexpression of the riboflavin biosynthesis cluster ribABGH. When the purine starvation was relieved by purine supplementation in the growth medium, the outcome was an immediate downregulation of the riboflavin biosynthesis cluster and a reduction in riboflavin production. Finally, by applying the new isolates in milk fermentation, the riboflavin content of milk (0.99 mg/liter) was improved to 2.81 mg/liter, compared with 0.66 mg/liter and 1.51 mg/liter by using the wild-type strain and the original roseoflavin-resistant mutant JC017, respectively. The results obtained demonstrate how powerful classical mutagenesis can be when combined with droplet-based microfluidic screening technology for obtaining microorganisms with useful attributes.

The food industry prefers to use classical approaches, e.g., random mutagenesis followed by screening, to improve microorganisms used in food production, as the use of recombinant DNA technologies is still not widely accepted. Although modern automated screening platforms are widely accessible, screening remains as a bottleneck in strain development, especially when a mild mutagenesis approach is applied to reduce the chance of accumulating unintended mutations, which may cause unwanted phenotypic changes. Here, we incorporate a droplet-based high-throughput screening method into the strain development process and readily capture L. lactis variants with more efficient vitamin secretion from low-error-rate mutagenesis libraries. This study shows that useful mutants showing strong phenotypes but without extensive mutations can be identified with efficient screening technologies. It is therefore possible to avoid accumulating detrimental mutations while enriching beneficial ones through iterative mutagenesis screening. Due to the low mutation rates, the genetic determinants are also readily identified.

Rapid generation of drug-resistance alleles at endogenous loci using CRISPR-Cas9 indel mutagenesis

Genetic alterations conferring resistance to the effects of chemical inhibitors are valuable tools for validating on-target effects in cells. Unfortunately, for many therapeutic targets such alleles are not available. To address this issue, we evaluated whether CRISPR-Cas9-mediated insertion/deletion (indel) mutagenesis can produce drug-resistance alleles at endogenous loci. This method takes advantage of the heterogeneous in-frame alleles produced following Cas9-mediated DNA cleavage, which we show can generate rare alleles that confer resistance to the growth-arrest caused by chemical inhibitors. We used this approach to identify novel resistance alleles of two lysine methyltransferases, DOT1L and EZH2, which are each essential for the growth of MLL-fusion leukemia cells. We biochemically characterized the DOT1L mutation, showing that it is significantly more active than the wild-type enzyme. These findings validate the on-target anti-leukemia activities of existing DOT1L and EZH2 inhibitors and reveal a simple method for deriving drug-resistance alleles for novel targets, which may have utility during early stages of drug development.

The power of multiplexed functional analysis of genetic variants

New technologies have recently enabled saturation mutagenesis and functional analysis of nearly all possible variants of regulatory elements or proteins of interest in single experiments. Here we discuss the past, present, and future of such multiplexed (functional) assays for variant effects (MAVEs). MAVEs provide detailed insight into sequence-function relationships, and they may prove critical for the prospective clinical interpretation of genetic variants.

Massively parallel single-amino-acid mutagenesis

Random mutagenesis methods only partially cover the mutational space, and are constrained by DNA synthesis length limitations. Here, we demonstrate PALS, a single-volume, site-directed mutagenesis approach using microarray-programmed oligonucleotides. We created libraries including nearly every missense mutation as singleton events for the yeast transcription factor Gal4 (99.9% coverage) and human tumor suppressor p53 (93.5%). PALS-based comprehensive missense mutational scans may aid structure-function studies, protein engineering, and the interpretation of variants identified by clinical sequencing.

Saturation editing of genomic regions by multiplex homology-directed repair

Saturation mutagenesis^1,2 – coupled to an appropriate biological assay – represents a fundamental means of achieving a high-resolution understanding of regulatory³ and protein-coding⁴ nucleic acid sequences of interest. However, mutagenized sequences introduced in trans on episomes or via random or “safe-harbor” integration fail to capture the native context of the endogenous chromosomal locus⁵. This shortcoming markedly limits the interpretability of the resulting measurements of mutational impact. Here, we couple CRISPR/Cas9 RNA-guided cleavage⁶ with multiplex homology-directed repair (HDR) using a complex library of donor templates to demonstrate saturation editing of genomic regions. In exon 18 of BRCA1, we replace a six base-pair (bp) genomic region with all possible hexamers, or the full exon with all possible single nucleotide variants (SNVs), and measure strong effects on transcript abundance attributable to nonsense-mediated decay and exonic splicing elements. We similarly perform saturation genome editing of a well-conserved coding region of an essential gene, DBR1, and measure relative effects on growth that correlate with functional impact. Measurement of the functional consequences of large numbers of mutations with saturation genome editing will potentially facilitate high-resolution functional dissection of both cis-regulatory elements and trans-acting factors, as well as the interpretation of variants of uncertain significance observed in clinical sequencing.

Empirical complexities in the genetic foundations of lethal mutagenesis

From population genetics theory, elevating the mutation rate of a large population should progressively reduce average fitness. If the fitness decline is large enough, the population will go extinct in a process known as lethal mutagenesis. Lethal mutagenesis has been endorsed in the virology literature as a promising approach to viral treatment, and several in vitro studies have forced viral extinction with high doses of mutagenic drugs. Yet only one empirical study has tested the genetic models underlying lethal mutagenesis, and the theory failed on even a qualitative level. Here we provide a new level of analysis of lethal mutagenesis by developing and evaluating models specifically tailored to empirical systems that may be used to test the theory. We first quantify a bias in the estimation of a critical parameter and consider whether that bias underlies the previously observed lack of concordance between theory and experiment. We then consider a seemingly ideal protocol that avoids this bias-mutagenesis of virions-but find that it is hampered by other problems. Finally, results that reveal difficulties in the mere interpretation of mutations assayed from double-strand genomes are derived. Our analyses expose unanticipated complexities in testing the theory. Nevertheless, the previous failure of the theory to predict experimental outcomes appears to reside in evolutionary mechanisms neglected by the theory (e.g., beneficial mutations) rather than from a mismatch between the empirical setup and model assumptions. This interpretation raises the specter that naive attempts at lethal mutagenesis may augment adaptation rather than retard it.

Massively parallel functional dissection of mammalian enhancers in vivo

The functional consequences of genetic variation in mammalian regulatory elements are poorly understood. We report the in vivo dissection of three mammalian enhancers at single-nucleotide resolution through a massively parallel reporter assay. For each enhancer, we synthesized a library of >100,000 mutant haplotypes with 2-3% divergence from the wild-type sequence. Each haplotype was linked to a unique sequence tag embedded within a transcriptional cassette. We introduced each enhancer library into mouse liver and measured the relative activities of individual haplotypes en masse by sequencing the transcribed tags. Linear regression analysis yielded highly reproducible estimates of the effect of every possible single-nucleotide change on enhancer activity. The functional consequence of most mutations was modest, with ∼22% affecting activity by >1.2-fold and ∼3% by >2-fold. Several, but not all, positions with higher effects showed evidence for purifying selection, or co-localized with known liver-associated transcription factor binding sites, demonstrating the value of empirical high-resolution functional analysis.

Needles in stacks of needles: finding disease-causal variants in a wealth of genomic data

Genome and exome sequencing yield extensive catalogues of human genetic variation. However, pinpointing the few phenotypically causal variants among the many variants present in human genomes remains a major challenge, particularly for rare and complex traits wherein genetic information alone is often insufficient. Here, we review approaches to estimate the deleteriousness of single nucleotide variants (SNVs), which can be used to prioritize disease-causal variants. We describe recent advances in comparative and functional genomics that enable systematic annotation of both coding and non-coding variants. Application and optimization of these methods will be essential to find the genetic answers that sequencing promises to hide in plain sight.

PCRless library mutagenesis via oligonucleotide recombination in yeast

The directed evolution of biomolecules with new functions is largely performed in vitro, with PCR mutagenesis followed by high-throughput assays for desired activities. As synthetic biology creates impetus for generating biomolecules that function in living cells, new technologies are needed for performing mutagenesis and selection for directed evolution in vivo. Homologous recombination, routinely exploited for targeted gene alteration, is an attractive tool for in vivo library mutagenesis, yet surprisingly is not routinely used for this purpose. Here, we report the design and characterization of a yeast-based system for library mutagenesis of protein loops via oligonucleotide recombination. In this system, a linear vector is co-transformed with single-stranded mutagenic oligonucleotides. Using repair of nonsense codons engineered in three different active-site loops in the selectable marker TRP1 as a model system, we first optimized the recombination efficiency. Single-loop recombination was highly efficient, averaging 5%, or 4.0×10(5) recombinants. Multiple loops could be simultaneously mutagenized, although the efficiencies dropped to 0.2%, or 6.0×10(3) recombinants, for two loops and 0.01% efficiency, or 1.5×10(2) recombinants, for three loops. Finally, the utility of this system for directed evolution was tested explicitly by selecting functional variants from a mock library of 1:10(6) wild-type:nonsense codons. Sequencing showed that oligonucleotide recombination readily covered this large library, mutating not only the target codon but also encoded silent mutations on either side of the library cassette. Together these results establish oligonucleotide recombination as a simple and powerful library mutagenesis technique and advance efforts to engineer the cell for fully in vivo directed evolution.

Random-scanning mutagenesis

Oligonucleotide-mediated mutagenesis is a useful tool for engineering nucleotide changes at defined positions in a DNA sequence. Oligonucleotide-based approaches are commonly used to introduce missense mutations at individual codons in a gene or gene segment, thereby revealing the functional importance of specific amino acid residues in a protein. For mutagenesis studies involving tracts of polypeptide sequence, investigators typically change each successive residue to alanine or to a limited number of alternative amino acids. Although these strategies can provide useful information, it is sometimes desirable to test a broader spectrum of amino acid changes at the targeted positions. This article describes a facile, oligonucleotide-based method for generating all 19 possible replacements at individual amino acid sites within a protein. This technique is known as "random-scanning mutagenesis" and is illustrated herein using examples from our studies of a conserved polymerase motif in HIV-1 reverse transcriptase.

Prediction of domain interactive motif pairs

Protein domain-domain interaction pairs supply functional information about the interacting proteins; and finding interaction motif pairs in protein-protein interaction database can deeply disclose the essence of the protein interaction. Up to now, there is little research work on prediction of interaction motif pairs within domain-domain interaction pairs. In this paper, we propose a new method to predict domain interaction motif pairs. We start from collecting contact segment pairs in the PDB protein complexes, and then use the contact segment pairs as seeds to iteratively cluster the protein-protein interaction database with the help of functional domains, finally we generalize the similar segment pair clusters to produce motif pairs. Using our method, we find 528 motif pairs.

Toward Assessing the Position-Dependent Contributions of Backbone Hydrogen Bonding to β-Sheet Folding Thermodynamics Employing Amide-to-Ester Perturbations

An amide-to-ester backbone substitution in a protein is accomplished by replacing an alpha-amino acid residue with the corresponding alpha-hydroxy acid, preserving stereochemistry, and conformation of the backbone and the structure of the side chain. This substitution replaces the amide NH (a hydrogen bond donor) with an ester O (which is not a hydrogen bond donor) and the amide carbonyl (a strong hydrogen bond acceptor) with an ester carbonyl (a weaker hydrogen bond acceptor), thus perturbing folding energetics. Amide-to-ester perturbations were used to evaluate the thermodynamic contribution of each hydrogen bond in the PIN WW domain, a three-stranded beta-sheet protein. Our results reveal that removing a hydrogen bond donor destabilizes the native state more than weakening a hydrogen bond acceptor and that the degree of destabilization is strongly dependent on the location of the amide bond replaced. Hydrogen bonds near turns or at the ends of beta-strands are less influential than hydrogen bonds that are protected within a hydrophobic core. Beta-sheet destabilization caused by an amide-to-ester substitution cannot be directly related to hydrogen bond strength because of differences in the solvation and electrostatic interactions of amides and esters. We propose corrections for these differences to obtain approximate hydrogen bond strengths from destabilization energies. These corrections, however, do not alter the trends noted above, indicating that the destabilization energy of an amide-to-ester mutation is a good first-order approximation of the free energy of formation of a backbone amide hydrogen bond.

Discovery of stable and significant binding motif pairs from PDB complexes and protein interaction datasets

MOTIVATION

Discovery of binding sites is important in the study of protein-protein interactions. In this paper, we introduce stable and significant motif pairs to model protein-binding sites. The stability is the pattern's resistance to some transformation. The significance is the unexpected frequency of occurrence of the pattern in a sequence dataset comprising known interacting protein pairs. Discovery of stable motif pairs is an iterative process, undergoing a chain of changing but converging patterns. Determining the starting point for such a chain is an interesting problem. We use a protein complex dataset extracted from the Protein Data Bank to help in identifying those starting points, so that the computational complexity of the problem is much released.

RESULTS

We found 913 stable motif pairs, of which 765 are significant. We evaluated these motif pairs using comprehensive comparison results against random patterns. Wet-experimentally discovered motifs reported in the literature were also used to confirm the effectiveness of our method.

SUPPLEMENTARY INFORMATION

http://sdmc.i2r.a-star.edu.sg/BindingMotifPairs.

Structural double-mutant cycle analysis of a hydrogen bond network in ketosteroid isomerase from Pseudomonas putida biotype B

KSI (ketosteroid isomerase) catalyses an allylic isomerization reaction at a diffusion-controlled rate. A hydrogen bond network, Asp(99).Water(504).Tyr(14).Tyr(55).Tyr(30), connects two critical catalytic residues, Tyr(14) and Asp(99), with Tyr(30), Tyr(55) and a water molecule in the highly apolar active site of the Pseudomonas putida KSI. In order to characterize the interactions among these amino acids in the hydrogen bond network of KSI, double-mutant cycle analysis was performed, and the crystal structure of each mutant protein within the cycle was determined respectively to interpret the coupling energy. The DeltaDeltaG(o) values of Y14F/D99L (Tyr(14)-->Phe/Asp(99)-->Leu) KSI, 25.5 kJ/mol for catalysis and 28.9 kJ/mol for stability, were smaller than the sums (i.e. 29.7 kJ/mol for catalysis and 34.3 kJ/mol for stability) for single mutant KSIs respectively, indicating that the effect of the Y14F/D99L mutation was partially additive for both catalysis and stability. The partially additive effect of the Y14F/D99L mutation suggests that Tyr(14) and Asp(99) should interact positively for the stabilization of the transition state during the catalysis. The crystal structure of Y14F/D99L KSI indicated that the Y14F/D99L mutation increased the hydrophobic interaction while disrupting the hydrogen bond network. The DeltaDeltaG(o) values of both Y30F/D99L and Y55F/D99L KSIs for the catalysis and stability were larger than the sum of single mutants, suggesting that either Tyr(30) and Asp(99) or Tyr(55) and Asp(99) should interact negatively for the catalysis and stability. These synergistic effects of both Y30F/D99L and Y55F/D99L mutations resulted from the disruption of the hydrogen bond network. The synergistic effect of the Y55F/D99L mutation was larger than that of the Y30F/D99L mutation, since the former mutation impaired the proper positioning of a critical catalytic residue, Tyr(14), involved in the catalysis of KSI. The present study can provide insight into interpreting the coupling energy measured by double-mutant cycle analysis based on the crystal structures of the wild-type and mutant proteins.

Site-directed mutagenesis and generation of chimeric viruses by homologous recombination in yeast to facilitate analysis of plant-virus interactions

A yeast homologous recombination system was used to generate mutants and chimeras in the genome of Potato leafroll virus (PLRV). A yeast-bacteria shuttle vector was developed that allows mutants and chimeras generated in yeast to be transformed into Escherichia coli for confirmation of the mutations and transformed into Agrobacterium tumefaciens to facilitate agroinfection of plants by the mutant PLRV genomes. The advantages of the system include the high frequency of recovered mutants generated by yeast homologous recombination, the ability to generate over 20 mutants and chimeras using only two restriction endonuclease sites, the ability to introduce multiple additional sequences using three and four DNA fragments, and the mobilization of the same plasmid from yeast to E. coli, A. tumefaciens, and plants. The wild-type PLRV genome showed no loss of virulence after sequential propagation in yeast, E. coli, and A. tumefaciens. Moreover, many PLRV clones with mutations generated in the capsid protein and readthrough domain of the capsid protein replicated and moved throughout plants. This approach will facilitate the analysis of plant-virus interactions of in vivo-generated mutants for many plant viruses, especially those not transmissible mechanically to plants.

Site-specific genomic (SSG) and random domain-localized (RDL) mutagenesis in yeast

Background

A valuable weapon in the arsenal available to yeast geneticists is the ability to introduce specific mutations into yeast genome. In particular, methods have been developed to introduce deletions into the yeast genome using PCR fragments. These methods are highly efficient because they do not require cloning in plasmids.

Results

We have modified the existing method for introducing deletions in the yeast (S. cerevisiae) genome using PCR fragments in order to target point mutations to this genome. We describe two PCR-based methods for directing point mutations into the yeast genome such that the final product contains no other disruptions. In the first method, site-specific genomic (SSG) mutagenesis, a specific point mutation is targeted into the genome. In the second method, random domain-localized (RDL) mutagenesis, a mutation is introduced at random within a specific domain of a gene. Both methods require two sequential transformations, the first transformation integrates the URA3 marker into the targeted locus, and the second transformation replaces URA3 with a PCR fragment containing one or a few mutations. This PCR fragment is synthesized using a primer containing a mutation (SSG mutagenesis) or is synthesized by error-prone PCR (RDL mutagenesis). In SSG mutagenesis, mutations that are proximal to the URA3 site are incorporated at higher frequencies than distal mutations, however mutations can be introduced efficiently at distances of at least 500 bp from the URA3 insertion. In RDL mutagenesis, to ensure that incorporation of mutations occurs at approximately equal frequencies throughout the targeted region, this region is deleted at the same time URA3 is integrated.

Conclusion

SSG and RDL mutagenesis allow point mutations to be easily and efficiently incorporated into the yeast genome without disrupting the native locus.

Construction of block-shuffled libraries of DNA for evolutionary protein engineering: Y-ligation-based block shuffling

Evolutionary protein engineering is now proceeding to a new stage in which novel technologies, besides the conventional point mutations, to generate a library of proteins, are required. In this context, a novel method for shuffling and rearranging DNA blocks (leading to protein libraries) is reported. A cycle of processes for producing combinatorial diversity was devised and designated Y-ligation-based block shuffling (YLBS). Methodological refinement was made by applying it to the shuffling of module-sized and amino acid-sized blocks. Running three cycles of YLBS with module-sized GFP blocks resulted in a high diversity of an eight-block shuffled library. Partial shuffling of the central four blocks of GFP was performed to obtain in-effect shuffled protein, resulting in an intact arrangement. Shuffling of amino acid monomer-sized blocks by YLBS was also performed and a diversity of more than 10(10) shuffled molecules was attained. The deletion problems encountered during these experiments were shown to be solved by additional measures which tame type IIS restriction enzymes. The frequency of appearance of each block was skewed but was within a permissible range. Therefore, YLBS is the first general method for generating a huge diversity of shuffled proteins, recombining domains, exons and modules with ease.

Linker-mediated recombinational subcloning of large DNA fragments using yeast

The homologous recombination pathway in yeast is an ideal tool for the sequence-specific assembly of plasmids. Complementary 80-nucleotide oligonucleotides that overlap a vector and a target fragment were found to serve as efficient recombination linkers for fragment subcloning. Using electroporation, single-stranded 80-mers were adequate for routine plasmid construction. A cycloheximide-based counterselection was introduced to increase the specificity of cloning by homologous recombination relative to nonspecific vector background. Reconstruction experiments suggest this counterselection increased cloning specificity by 100-fold. Cycloheximide counterselection was used in conjunction with 80-bp linkers to subclone targeted regions from bacterial artificial chromosomes. This technology may find broad application in the final stages of completing the Human Genome Sequencing Project and in applications of BAC clones to the functional analysis of complex genomes.

Alanine-scanning mutagenesis of Aspergillus gamma-tubulin yields diverse and novel phenotypes

We have created 41 clustered charged-to-alanine scanning mutations of the mipA, gamma-tubulin, gene of Aspergillus nidulans and have created strains carrying these mutations by two-step gene replacement and by a new procedure, heterokaryon gene replacement. Most mutant alleles confer a wild-type phenotype, but others are lethal or conditionally lethal. The conditionally lethal alleles exhibit a variety of phenotypes under restrictive conditions. Most have robust but highly abnormal mitotic spindles and some have abnormal cytoplasmic microtubule arrays. Two alleles appear to have reduced amounts of gamma-tubulin at the spindle pole bodies and nucleation of spindle microtubule assembly may be partially inhibited. One allele inhibits germ tube formation. The cold sensitivity of two alleles is strongly suppressed by the antimicrotubule agents benomyl and nocodazole and a third allele is essentially dependent on these compounds for growth. Together our data indicate that gamma-tubulin probably carries out functions essential to mitosis and organization of cytoplasmic microtubules in addition to its well-documented role in microtubule nucleation. We have also placed our mutations on a model of the structure of gamma-tubulin and these data give a good initial indication of the functionally important regions of the molecule.

Orthogonal combinatorial mutagenesis: a codon-level combinatorial mutagenesis method useful for low multiplicity and amino acid-scanning protocols

We describe here a method to generate combinatorial libraries of oligonucleotides mutated at the codon-level, with control of the mutagenesis rate so as to create predictable binomial distributions of mutants. The method allows enrichment of the libraries with single, double or larger multiplicity of amino acid replacements by appropriate choice of the mutagenesis rate, depending on the concentration of synthetic precursors. The method makes use of two sets of deoxynucleoside-phosphoramidites bearing orthogonal protecting groups [4,4'-dimethoxytrityl (DMT) and 9-fluorenylmethoxycarbonyl (Fmoc)] in the 5' hydroxyl. These phosphoramidites are divergently combined during automated synthesis in such a way that wild-type codons are assembled with commercial DMT-deoxynucleoside-methyl-phosphoramidites while mutant codons are assembled with Fmoc-deoxynucleoside-methyl-phosphoramidites in an NNG/C fashion in a single synthesis column. This method is easily automated and suitable for low mutagenesis rates and large windows, such as those required for directed evolution and alanine scanning. Through the assembly of three oligonucleotide libraries at different mutagenesis rates, followed by cloning at the polylinker region of plasmid pUC18 and sequencing of 129 clones, we concluded that the method performs essentially as intended.

Efficient gene targeted random mutagenesis in genetically stable Escherichia coli strains

We describe a method to generate in vivo collections of mutants orders of magnitude larger than previously possible. The method favors accumulation of mutations in the target gene, rather than in the host chromosome. This is achieved by propagating the target gene on a plasmid, in Escherichia coli cells, within the region preferentially replicated by DNA polymerase I (Pol I), which replicates only a minor fraction of the chromosome. Mutagenesis is enhanced by a conjunction of a Pol I variant that has a low replication fidelity and the absence of the mutHLS system that corrects replication errors. The method was tested with two reporter genes, encoding lactose repressor or lipase. The proportion of mutants in the collection was estimated to reach 1% after one cycle of growth and 10% upon prolonged cell cultivation, resulting in collections of 10(12)-10(13) mutants per liter of cell culture. The extended cultivation did not affect growth properties of the cells. We suggest that our method is well suited for generating protein variants too rare to be present in the collections established by methods used previously and for isolating the genes that encode such variants by submitting the cells of the collections to appropriate selection protocols.

The essential role of Saccharomyces cerevisiae CDC6 nucleotide-binding site in cell growth, DNA synthesis, and Orc1 association

Saccharomyces cerevisiae Cdc6 is a protein required for the initiation of DNA replication. The biochemical function of the protein is unknown, but the primary sequence contains motifs characteristic of nucleotide-binding sites. To study the requirement of the nucleotide-binding site for the essential function of Cdc6, we have changed the conserved Lys114 at the nucleotide-binding site to five other amino acid residues. We have used these mutants to investigate in vivo roles of the conserved lysine in the growth rate of transformant cells and the complementation of cdc6 temperature-sensitive mutant cells. Our results suggest that replacement of Lys with Glu (K114E) and Pro (K114P) leads to loss-of-function in supporting cell growth, replacement of the Lys with Gln (K114Q) or Leu (K114L) yields partially functional proteins, and replacement with Arg yields a phenotype equivalent to wild-type, a silent mutation. To investigate what leads to the growth defects derived from the mutations at the nucleotide-binding site, we evaluated its gene functions in DNA replication by the assays of the plasmid stability and chromosomal DNA synthesis. Indeed, the K114P and K114E mutants showed the complete retraction of DNA synthesis. In order to test its effect on the G1/S transition of the cell cycle, we have carried out the temporal and spatial studies of yeast replication complex. To do this, yeast chromatin fractions from synchronized culture were prepared to detect the Mcm5 loading onto the chromatin in the presence of the wild-type Cdc6 or mutant cdc6(K114E) proteins. We found that cdc6(K114E) is defective in the association with chromatin and in the loading of Mcm5 onto chromatin origins. To further investigate the molecular mechanism of nucleotide-binding function, we have demonstrated that the Cdc6 protein associates with Orc1 in vitro and in vivo. Intriguingly, the interaction between Orc1 and Cdc6 is disrupted when the cdc6(K114E) protein is used. Our results suggest that a proper molecular interaction between Orc1 and Cdc6 depends on the functional ATP-binding of Cdc6, which may be a prerequisite step to assemble the operational replicative complex at the G1/S transition.

Contributions of the TATA box sequence to rate-limiting steps in transcription initiation by RNA polymerase II

We have examined the role of the TATA box in determining transcription initiation frequency in vitro by studying a collection of promoters containing different TATA sequences in the context of the adenovirus major late promoter. In addition to measuring transcription rates, we have determined how the sequence changes affected the association and dissociation kinetics and the affinity of TBP binding. We observed that transcription from promoters containing the highest affinity TATA boxes is limited by the rate with which TBP associates with the promoter. In contrast, transcription from promoters containing lower affinity TATA boxes appears to be limited both by how much TBP is bound and by the relatively low occupancy of the conformation that can undergo subsequent steps in preinitiation complex assembly. The implications of these results in understanding the mechanism of transcription enhancement by transcriptional activators is discussed.

Scoring functions for computational algorithms applicable to the design of spiked oligonucleotides

Protein engineering by inserting stretches of random DNA sequences into target genes in combination with adequate screening or selection methods is a versatile technique to elucidate and improve protein functions. Established compounds for generating semi-random DNA sequences are spiked oligonucleotides which are synthesised by interspersing wild type (wt) nucleotides of the target sequence with certain amounts of other nucleotides. Directed spiking strategies reduce the complexity of a library to a manageable format compared with completely random libraries. Computational algorithms render feasible the calculation of appropriate nucleotide mixtures to encode specified amino acid subpopulations. The crucial element in the ranking of spiked codons generated during an iterative algorithm is the scoring function. In this report three scoring functions are analysed: the sum-of-square-differences function s, a modified cubic function c, and a scoring function m derived from maximum likelihood considerations. The impact of these scoring functions on calculated amino acid distributions is demonstrated by an example of mutagenising a domain surrounding the active site serine of subtilisin-like proteases. At default weight settings of one for each amino acid, the new scoring function m is superior to functions s and c in finding matches to a given amino acid population.

Approaches to DNA mutagenesis: an overview

In the last several years, the use of double-stranded DNA templates together with thermostable-polymerase PCR has essentially replaced the use single-stranded DNA templates using the thermolabile polymerase for in vitro mutagenesis. Numerous PCR methods are now available, such as overlap-extension PCR, megaprimer PCR, and inverse PCR. All these PCR methods are reliable, effective, and convenient, although they are more prone to high rates of spontaneous error in mutant DNAs than are methods using thermolabile polymerases. Some improvements, such as the introduction of methylated templates, have been employed to minimize PCR errors. On the other hand, because of the introduction of many selection measures (e.g., restoration of antibiotic resistance, restoration of replication origin and unique site elimination), both double-stranded and single-stranded DNAs can now be used as templates for mutagenesis using thermolabile polymerase methods. For PCR methods, selection measures such as nested PCR has developed. All these selection measures have greatly improved the efficiency of mutagenesis by removing wild-type templates prior to transformation. Many efficient methods are available for both SDM and REM. Mutations can be introduce in vitro or in vivo, either by mutagenic primers or by erroneous DNA synthesis. Thus, choices largely depend on the experimental needs and resources of the investigator.

Expression of bacterial Rubisco genes in Escherichia coli

The structural genes for three forms of Rubisco have been isolated from bacteria and introduced into various plasmids. Apart from details of the sequences which have been obtained from these constructs, they are now being exploited for mutagenesis to determine the identity and specific function of the individual amino acid residues that compose the active site. These methods have been applied to a plasmid that contains the structural gene for the simplest form of Rubisco from Rhodospirillum rubrum to obtain mutant enzymes with altered activity. The construct pRR2119 is also expressed to very high levels in Escherichia coli and enough recombinant protein of both the wild-type and m utant enzymes can be obtained for detailed physico-chemical studies. Other vectors have now been constructed, containing the genes of prokaryotic Rubisco that assemble into an active form I enzyme. The levels of expression are acceptable and the product is similar to the authentic enzyme. These constructs are now being used for mutagenesis in vitro to attempt to alter the relative rates of the oxygenase and carboxylase activities.

Intracellular location of the Saccharomyces cerevisiae CDC6 gene product

The CDC6 gene product from Saccharomyces cerevisiae is required for transition from late G1 to S phase of the cell cycle. We have investigated the subcellular localization of the CDC6 protein in yeast to explore where Cdc6p exerts its gene function (s). Using affinity-purified sera we localized Cdc6p to the cytoplasm and the nuclear matrix by both subcellular fractionation and indirect immunofluorescence microscopy. The nuclear localization was confirmed to be in the nuclear scaffold by the low-salt extraction method. The Cdc6p cannot be detected in the mitochondrial or plasma membrane fractions. Using indirect immunofluorescence, we found that a subpopulation of Cdc6p migrated into the nucleus after G1/S transition and diminished after M phase, suggesting its temporal role in nuclear DNA replication. The predicted Cdc6p polypeptide contains a conserved nuclear localization, 27PLKRKKL33, similar to that of the SV40 large T antigen and other nuclear proteins. To test whether this peptide segment plays a role in mediating nuclear transport, we have carried out site-directed mutagenesis to alter the conserved 29Lys to Thr and Arg. The wild-type nuclear localization signal of Cdc6p was found to mediate the LacZ reporter gene fused to CDC6 efficiently to the nucleus, but not the mutated versions of the nuclear localization motif. The results suggested that 29Lys is important in mediating nuclear localization, the 29Thr and 29Arg mutant versions of the CDC6 gene were also unable to complement the cdc6 temperature-sensitive mutant. However, when these mutants were expressed from a multicopy plasmid, the mutated genes could complement the mutation. Similar results were obtained in the cdc6-disrupted cells. Taken together, we suggest that (i) Cdc6p is predominantly located in the cytoplasm, (ii) the nuclear entry of Cdc6p is cell cycle dependent, and (iii) nuclear entry of Cdc6p is mediated by its nuclear localization signal. The presence of Cdc6p in both the nucleus and the cytoplasm suggests a model that Cdc6p exerts its gene function in DNA replication and mitotic restraint in the cell cycle.

Codon cassette mutagenesis: a general method to insert or replace individual codons by using universal mutagenic cassettes

We describe codon cassette mutagenesis, a simple method of mutagenesis that uses universal mutagenic cassettes to deposit single codons at specific sites in double-stranded DNA. A target molecule is first constructed that contains a blunt, double-strand break at the site targeted for mutagenesis. A double-stranded mutagenic codon cassette is then inserted at the target site. Each mutagenic codon cassette contains a three base pair direct terminal repeat and two head-to-head recognition sequences for the restriction endonuclease Sapl, an enzyme that cleaves outside of its recognition sequence. The intermediate molecule containing the mutagenic cassette is then digested with Sapl, thereby removing most of the mutagenic cassette, leaving only a three base cohesive overhang that is ligated to generate the final insertion or substitution mutation. A general method for constructing blunt-end target molecules suitable for this approach is also described. Because the mutagenic cassette is excised during this procedure and alters the target only by introducing the desired mutation, the same cassette can be used to introduce a particular codon at all target sites. Each cassette can deposit two different codons, depending on the orientation in which it is inserted into the target molecule. Therefore, a series of eleven cassettes is sufficient to insert all possible amino acids at any constructed target site. Thus codon cassettes are 'off-the-shelf' reagents, and this methodology should be a particularly useful and inexpensive approach for subjecting multiple different positions in a protein sequence to saturation mutagenesis.

RAMHA: a PC-based Monte-Carlo simulation of random saturation mutagenesis

Random mutagenesis is a powerful tool in protein structure-function analyses. One approach to random mutagenesis is the de novo synthesis of polypeptide-encoding oligodeoxy-nucleotides using doped nucleoside phosphoramidites. A Turbo PASCAL program, RAMHA, is described for modeling such mutagenesis. Upon entering the target sequence and the desired level of nucleotide contamination, RAMHA performs a Monte Carlo simulation of the mutagenesis, compiling statistics on the similarity of resultant mutant polypeptides to the wild-type sequence, the frequency of premature open-reading frame terminations, and other relevant outcomes. Simulated mutagenesis of two DNA targets has led to the development of two different strategies to avoid the random introduction of stop codons within mutagenized gene segments.

A random mutagenesis procedure: application to the POL3 gene of Saccharomyces cerevisiae

To obtain a broad spectrum of mutations in the POL3 gene, we have developed an efficient random mutagenesis procedure. Partially extended, primed, single-stranded templates were used for forced misincorporation of non-complementary nucleotides, extended to completion and ligated. Linear fragments of the resulting amplified mutagenized library were then used to transform a Saccharomyces cerevisiae strain by the marker replacement technique. This procedure has proven to be very efficient when applied to the C-terminal moiety of POL3, yielding 24 temperature-sensitive mutants and six extragenic revertants.

Random mutagenesis of the gene for bacteriophage T7 RNA polymerase

Random mutagenesis of the gene for bacteriophage T7 RNA polymerase was used to identify functionally essential amino acid residues of the enzyme. A two-plasmid system was developed that permits the straightforward isolation of T7 RNA polymerase mutants that had lost almost all catalytic activity. It was shown that substitutions of Thr and Ala for Pro at the position 563, Ser for Tyr571, Pro for Thr636, Asp for Tyr639 and of Cys for Phe646 resulted in inactivation of the enzyme. It is noteworthy that all these mutations are limited to two short regions that are highly conservative in sequences of monomeric RNA polymerases.

Random mutagenesis of a plasmid-borne glycosidase gene and phenotypic selection of mutants in Escherichia coli

The bglA gene from Bacillus polymyxa encodes a beta-glucosidase able to hydrolyze p-nitrophenyl-beta-D-glucopyranoside (PNPG), and 5-bromo-4-chloro-3-indolyl-beta-D-galactopyranoside (X-gal), chromogenic substrates for beta-glucosidases and beta-galactosidases respectively. A plasmid carrying the blgA gene inserted in vector pUC18 was mutagenized in vitro with hydroxylamine, and subsequently used to transform E. coli selecting for the ampicillin resistance conferred by the cloning vector. The transformants were screened on petri dishes for mutations causing inability to hydrolyze either one of the two substrates, and for mutations increasing resistance of the enzyme to thermal inactivation. The isolation of several mutants with such characteristics suggests that the simple procedure used here can be applied to generate modifications of enzymatic properties that fit specific industrial requirements.

Adenovirus VARNA1 gene B block promoter element sequences required for transcription and for interaction with transcription factors

We constructed mutants with a deletion of either half of the 18 base-pair B block palindrome in the VARNA1 gene, mutants with different intra-palindromic spacings, a complete set of mutants with single base substitutions, and mutants with double and triple base substitutions in the palindrome. The transcription efficiencies of these mutants were determined in human KB cell-free cytoplasmic S100 extracts. The relative competing strength of each mutant, as determined by a sequential competition experiment, was used to assess each mutant's ability to sequester factors into formation of a stable preinitiation complex. The ability of each mutant to assemble transcriptionally active preinitiation complexes was also determined by direct transcription of the isolated complexes. Finally, the ability of each mutant to interact with the transcription factor(s) TFIIIC and form a distinct gel-resolved complex was also determined. From the results of the above assays, we concluded that the two seemingly identical halves of the palindrome did not contribute equally to transcription, or to assembly of the functional preinitiation complex, nor to interaction with TFIIIC. The anterior half (B1) of the B block palindrome, which is proximal to the A block promoter element, played a stronger role in transcription and in assembly of the functional preinitiation complex than the posterior half (B2) of the palindrome. Consistent with this observation, the point mutations in four base-pairs, GTTC, from +60 to +63 in the anterior half of the B block palindrome, has the most severe effect on transcription. In contrast, we showed that the central sequence and the posterior half (B2) played a stronger role than the anterior half (B1) of the B block palindrome in the interaction of the promoter with TFIIIC. This was corroborated by the observation that base substitutions in the central four base-pair sequence of the palindrome, TCGA, from +62 to +65, had the most severe effect on interaction with TFIIIC, and that mutations in most of the sequences in the posterior half of the B block palindrome had more drastic effects than mutations in the anterior half of the palindrome in this interaction. Furthermore, the spacing between the two halves of the B block palindrome had a drastic effect on the overall transcription efficiency and the interaction of the promoter with TFIIIC, suggesting that the interaction between the two halves of the B block palindrome is not only essential, but also synergistic for the interaction with TFIIIC as well as the assembly of a transcriptionally active preinitiation complex and efficient transcription.

Heterologous protein production in yeast

The exploitation of recombinant DNA technology to engineer expression systems for heterologous proteins represented a major task within the field of biotechnology during the last decade. Yeasts attracted the attention of molecular biologists because of properties most favourable for their use as hosts in heterologous protein production. Yeasts follow the general eukaryotic posttranslational modification pattern of expressed polypeptides, exhibit the ability to secrete heterologous proteins and benefit from an established fermentation technology. Aside from the baker's yeast Saccharomyces cerevisiae, an increasing number of alternative non-Saccharomyces yeast species are used as expression systems in basic research and for an industrial application. In the following review a selection from the different yeast systems is described and compared.

Random PCR mutagenesis screening of secreted proteins by direct expression in mammalian cells

We have developed a general method for screening randomly mutagenized expression libraries in mammalian cells by using fluorescence-activated cell sorting (FACS). The cDNA sequence of a secreted protein is randomly mutagenized by PCR under conditions of reduced Taq polymerase fidelity. The mutated DNA is inserted into an expression vector encoding the membrane glycophospholipid anchor sequence of decay-accelerating factor (DAF) fused to the C terminus of the secreted protein. This results in expression of the protein on the cell surface in transiently transfected mammalian cells, which can then be screened by FACS. This method was used to isolate mutants in the kringle 1 (K1) domain of tissue plasminogen activator (t-PA) that would no longer be recognized by a specific monoclonal antibody (mAb387) that inhibits binding of t-PA to its clearance receptor. DNA sequence analysis of the mutants and localization of the mutated residues on a three-dimensional model of the K1 domain identified three key discontinuous amino acid residues that are essential for mAb387 binding. Mutants with changes in any of these three residues were found to have reduced binding to the t-PA receptor on human hepatoma HepG2 cells but to retain full clot lysis activity.

Quantitative interpretations of double mutations of enzymes

The quantitative effect of a second mutation on a mutant enzyme may be antagonistic, absent, partially additive, additive, or synergistic with respect to the first mutation. Depending on which kinetic or thermodynamic parameter of an enzyme is measured, the same two mutations can interact differently in the double mutant. Additive effects of two mutations on an equilibrium constant, such as the dissociation constant of the enzyme-substrate complex (KS), occur when noninteracting residues which facilitate the same step (substrate binding) are mutated. Partially additive effects result from the cooperative interaction with the substrate of the two residues mutated, and synergistic effects result from the anticooperative interaction with the substrate of the two residues mutated. An alternative explanation for synergy is extensive unfolding of the enzyme. Antagonistic effects on an equilibrium constant such as KS result from opposing structural effects of the two mutations on substrate binding. No additional effect of the second mutation in the double mutant represents a limiting case of either partial additivity or antagonism [corrected]. The interactions of the effects of two mutations on a rate constant such as kcat have the same explanations as those given above for equilibrium constants since the binding of a rate-limiting transition state is occurring. However, due to kinetic complexity, the following exceptions and additions exist. Additive effects of two mutations on kcat occur when noninteracting residues which facilitate the same step are mutated, provided this step is rate limiting. If the affected step is not rate limiting then synergistic effects of the two mutations are observed as each mutation causes the step to become progressively more rate limiting. Additive effects on kcat also occur when the two mutations affect consecutive steps, provided one of them is rate limiting. Partially additive effects on kcat also occur when noninteracting residues facilitating consecutive, non-rate-limiting steps are mutated. These concepts, when applied to published data on double mutants of delta 5-3-ketosteroid isomerase, staphylococcal nuclease, tyrosyl-tRNA synthetase, glutathione reductase, and subtilisin, provide deeper insights into the independent, cooperative, anticooperative, or antagonistic interactions of amino acid residues in the binding of substrates, activators, and inhibitors and in promoting catalysis.

A general strategy for random insertion and substitution mutagenesis: substoichiometric coupling of trinucleotide phosphoramidites

Results from a number of recent studies suggest that amino acid insertion mutations may provide an important alternative to substitution mutations for modifying protein structures and functional activities. To facilitate the use of single-amino acid insertions, we have developed a general strategy for inducing random, in-phase codon insertions across a defined segment of a cloned gene. In brief, a mixture of blocked and protected trinucleotide phosphoramidites is coupled at substoichiometric levels after every third monomer coupling on a conventional solid-state synthesizer. From the heterogeneous mixture of oligonucleotide sequences thus generated, those oligonucleotides that have acquired a single additional codon are purified by urea/PAGE. By using equimolar amounts of GCT and GGT trinucleotides in the oligonucleotide synthesis plus standard oligonucleotide-directed mutagenesis techniques, we have induced as many as 13 different single alanine and glycine insertion mutations into the gene for staphylococcal nuclease in one experiment. On replacement of the 5'-dimethoxytrityl blocking group on the trinucleotide phosphoramidite with an acid-stable blocking group, such as levulinate or fluoren-9-ylmethoxycarbonyl (Fmoc), this same strategy of substoichiometric couplings at codon boundaries should permit the synthesis of complex pools of oligonucleotides for the introduction, with constant efficiency, of every type of amino acid substitution at each codon across a gene segment.