Functional and nonfunctional mutations distinguished by random recombination of homologous genes

Flexible active-site loops fine-tune substrate specificity of hyperthermophilic metallo-oxidases

Hyperthermophilic ('superheat-loving') archaea found in high-temperature environments such as Pyrobaculum aerophilum contain multicopper oxidases (MCOs) with remarkable efficiency for oxidizing cuprous and ferrous ions. In this work, directed evolution was used to expand the substrate specificity of P. aerophilum McoP for organic substrates. Six rounds of error-prone PCR and DNA shuffling followed by high-throughput screening lead to the identification of a hit variant with a 220-fold increased efficiency (kcat/Km) than the wild-type for 2,2'-azino-bis(3-ethylbenzothiazoline-6-sulphonic acid) (ABTS) without compromising its intrinsic activity for metal ions. The analysis of the X-ray crystal structure reveals four proximal mutations close to the T1Cu active site. One of these mutations is within the 23-residues loop that occludes this site, a distinctive feature of prokaryotic MCOs. The increased flexibility of this loop results in an enlarged tunnel and one additional pocket that facilitates bulky substrate-enzyme interactions. These findings underscore the synergy between mutations that modulate the dynamics of the active-site loop enabling enhanced catalytic function. This study highlights the potential of targeting loops close to the T1Cu for engineering improvements suitable for biotechnological applications.

Characterization of tumor evolution by functional clonality and phylogenetics in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a molecularly heterogeneous solid malignancy, and its fitness may be shaped by how its tumor cells evolve. However, ability to monitor tumor cell evolution is hampered by the presence of numerous passenger mutations that do not provide any biological consequences. Here we develop a strategy to determine the tumor clonality of three independent HCC cohorts of 524 patients with diverse etiologies and race/ethnicity by utilizing somatic mutations in cancer driver genes. We identify two main types of tumor evolution, i.e., linear, and non-linear models where non-linear type could be further divided into classes, which we call shallow branching and deep branching. We find that linear evolving HCC is less aggressive than other types. GTF2IRD2B mutations are enriched in HCC with linear evolution, while TP53 mutations are the most frequent genetic alterations in HCC with non-linear models. Furthermore, we observe significant B cell enrichment in linear trees compared to non-linear trees suggesting the need for further research to uncover potential variations in immune cell types within genomically determined phylogeny types. These results hint at the possibility that tumor cells and their microenvironment may collectively influence the tumor evolution process.

Engineering an efficient and bright split Corynactis californica green fluorescent protein

Split green fluorescent protein (GFP) has been used in a panoply of cellular biology applications to study protein translocation, monitor protein solubility and aggregation, detect protein–protein interactions, enhance protein crystallization, and even map neuron contacts. Recent work shows the utility of split fluorescent proteins for large scale labeling of proteins in cells using CRISPR, but sets of efficient split fluorescent proteins that do not cross-react are needed for multiplexing experiments. We present a new monomeric split green fluorescent protein (ccGFP) engineered from a tetrameric GFP found in Corynactis californica, a bright red colonial anthozoan similar to sea anemones and scleractinian stony corals. Split ccGFP from C. californica complements up to threefold faster compared to the original Aequorea victoria split GFP and enable multiplexed labeling with existing A. victoria split YFP and CFP.

Directed evolution study unveiling key sequence factors that affect translation efficiency in Escherichia coli

Synonymous mutations in protein coding genes significantly impact translation efficiency. We synthesized a pair of genes encoding green fluorescent protein that were separated by 160 synonymous mutations to investigate key factors that affect translation efficiency. One sequence was optimized for Escherichia coli (GFP(Eco)) and the other for Bacillus subtilis (GFP(Bsu)). When the genes were expressed in E. coli, GFP(Eco) fluoresced 12-fold stronger than GFP(Bsu), confirming the suboptimal nature of the GFP(Bsu) gene. We then employed directed evolution to improve the expression of GFP(Bsu). Random mutagenesis and DNA shuffling was used to generate mutant libraries, which were screened for fluorescence. A variant showing 6-fold fluorescence enhancement was identified, which contained a single mutation (G10A) in a rare codon for Gly-4. However, the substitution generated another type of rare codon, AGA, for Arg, suggesting that the improvement was caused by a factor other than the rare codon. We next applied saturation mutagenesis to Gly-4. The darkest variant contained a GGG codon (GFP(Bsu)-G) for Gly-4. Taking the location of the mutation into account, we hypothesized that destabilization of the mRNA secondary structure around the initiation codon improved the expression. We then randomized the nucleotide triplet in 5'-untranslated region (5'UTR) of GFP(Bsu), which is complementary to the Gly-4 codon. A variant showing 6-fold fluorescence enhancement was identified, which exhibited a destabilized secondary structure. When this 5'UTR sequence was combined with GFP(Bsu)-G, 22-fold fluorescent improvement was achieved. Collectively, the stability of the mRNA secondary structure around the initiation codon predominantly affected the translation efficiency.

Revealing biases inherent in recombination protocols

Background

The recombination of homologous genes is an effective protein engineering tool to evolve proteins. DNA shuffling by gene fragmentation and reassembly has dominated the literature since its first publication, but this fragmentation-based method is labor intensive. Recently, a fragmentation-free PCR based protocol has been published, termed recombination-dependent PCR, which is easy to perform. However, a detailed comparison of both methods is still missing.

Results

We developed different test systems to compare and reveal biases from DNA shuffling and recombination-dependent PCR (RD-PCR), a StEP-like recombination protocol. An assay based on the reactivation of β-lactamase was developed to simulate the recombination of point mutations. Both protocols performed similarly here, with slight advantages for RD-PCR. However, clear differences in the performance of the recombination protocols were observed when applied to homologous genes of varying DNA identities. Most importantly, the recombination-dependent PCR showed a less pronounced bias of the crossovers in regions with high sequence identity. We discovered that template variations, including engineered terminal truncations, have significant influence on the position of the crossovers in the recombination-dependent PCR. In comparison, DNA shuffling can produce higher crossover numbers, while the recombination-dependent PCR frequently results in one crossover. Lastly, DNA shuffling and recombination-dependent PCR both produce counter-productive variants such as parental sequences and have chimeras that are over-represented in a library, respectively. Lastly, only RD-PCR yielded chimeras in the low homology situation of GFP/mRFP (45% DNA identity level).

Conclusion

By comparing different recombination scenarios, this study expands on existing recombination knowledge and sheds new light on known biases, which should improve library-creation efforts. It could be shown that the recombination-dependent PCR is an easy to perform alternative to DNA shuffling.

Extending the diversity of cytochrome P450 enzymes by DNA family shuffling

The cytochrome P450 enzymes involved in xenobiotic metabolism are an excellent starting point for the directed evolution of novel biocatalysts due to their wide substrate specificity. A shuffled library of three highly homologous mammalian genes (for P450 2C9, P450 2C11 and P450 2C19) was constructed by applying a modified DNA family shuffling procedure. The modifications made to the traditional DNA shuffling protocols involved non-random digestion via the use of different combinations of restriction enzymes (REs) followed by isolation of fragments under 300 bp by size-selective filtration. Shuffled cytochrome P450 mutants were co-expressed in Escherichia coli with their redox partner, NADPH-cytochrome P450 reductase (NPR). We report here how non-random fragmentation may help in chimeragenesis within the areas of low sequence similarity such as substrate recognition sites (SRSs) that are generally underrepresented in recombination using the random fragmentation process. Size-selective filtration was used to limit recovery of incompletely digested fragments and consequently minimize the chances for contamination of the shuffled library with parental forms. No parental forms could be detected in the shuffled library using restriction fragment length polymorphism (RFLP) analysis, suggesting the library was free of parental contamination. Sequencing of randomly selected mutants demonstrated a high level of chimeragenesis with on average of 8.0+/-2.2 crossovers and a low level of mutagenesis with 5.2+/-2.8 spontaneous mutations per approximately 1.5 kbp of the full-length P450 sequence. The proportion of properly folded protein as indicated by the observation of characteristic Fe(II).CO vs. Fe(II) difference spectra was 15% (4/27) of analysed mutants. Screening of the shuffled library for indole oxidation revealed four clones with similar or higher levels of indigo pigment production to those of the parental P450s and two clones with elevated P450 expression. In this paper we present a method for the effective family shuffling of cytochrome P450 enzymes, applicable to the creation of mutant libraries with expanded metabolic diversity and with a significant proportion of functional clones.

Selection of novel vesicular stomatitis virus glycoprotein variants from a peptide insertion library for enhanced purification of retroviral and lentiviral vectors

The introduction of new features or functions that are not present in an original protein is a significant challenge in protein engineering. For example, modifications to vesicular stomatitis virus glycoprotein (VSV-G), which is commonly used to pseudotype retroviral and lentiviral vectors for gene delivery, have been hindered by a lack of structural knowledge of the protein. We have developed a transposon-based approach that randomly incorporates designed polypeptides throughout a protein to generate saturated insertion libraries and a subsequent high-throughput selection process in mammalian cells that enables the identification of optimal insertion sites for a novel designed functionality. This method was applied to VSV-G in order to construct a comprehensive library of mutants whose combined members have a His6 tag inserted at likely every site in the original protein sequence. Selecting the library via iterative retroviral infections of mammalian cells led to the identification of several VSV-G-His6 variants that were able to package high-titer viral vectors and could be purified by Ni-nitrilotriacetic acid affinity chromatography. Column purification of vectors reduced protein and DNA impurities more than 5,000-fold and 14,000-fold, respectively, from the viral supernatant. This substantially improved purity elicited a weaker immune response in the brain, without altering the infectivity or tropism from wild-type VSV-G-pseudotyped vectors. This work applies a powerful new tool for protein engineering to construct novel viral envelope variants that can greatly improve the safety and use of retroviral and lentiviral vectors for clinical gene therapy. Furthermore, this approach of library generation and selection can readily be extended to other challenges in protein engineering.

Improving protein solubility: the use of the Escherichia coli dihydrofolate reductase gene as a fusion reporter

We have devised a strategy for screening mutant libraries for enzyme variants with enhanced solubility. The method is based on the observation that Escherichia coli can become insensitive to the antibiotic trimethoprim (TMP) if dihydrofolate reductase (DHFR) is expressed at an appropriate level. DHFR is a very soluble protein and can be expressed at levels that exceed normally lethal concentrations of TMP. In our approach, the gene encoding an insoluble target protein is placed in a vector so that the translated protein will be fused to DHFR. The resulting fusion protein will form inclusion bodies and inactivate DHFR-the cells will be susceptible to TMP. Mutations to the target protein that make it more soluble will also make the fusion protein more soluble so that DHFR will be at least partially active-the cells will be resistant to TMP. As the solubility of the target protein increases, the cells will become more resistant to TMP. The system was tested with a putative acetyltransferase (ACE) from a strain of the marine bacterium Vibrio fischerii. The gene encoding this protein was of interest since it is part of a mobile gene cassette within an integron array of the strain in question. After multiple rounds of shuffling and selection, ACE mutants were produced that had significantly improved solubility.

Directed evolution of enzyme stability

Modern enzyme development relies to an increasing extent on strategies based on diversity generation followed by screening for variants with optimised properties. In principle, these directed evolution strategies might be used for optimising any enzyme property, which can be screened for in an economically feasible way, even if the molecular basis of that property is not known. Stability is an interesting property of enzymes because (1) it is of great industrial importance, (2) it is relatively easy to screen for, and (3) the molecular basis of stability relates closely to contemporary issues in protein science such as the protein folding problem and protein folding diseases. Thus, engineering enzyme stability is of both commercial and scientific interest. Here, we review how directed evolution has contributed to the development of stable enzymes and to new insight into the principles of protein stability. Several recent examples are described. These examples show that directed evolution is an effective strategy to obtain stable enzymes, especially when used in combination with rational or semi-rational engineering strategies. With respect to the principles of protein stability, some important lessons to learn from recent efforts in directed evolution are (1) that there are many structural ways to stabilize a protein, which are not always easy to rationalize, (2) that proteins may very well be stabilized by optimizing their surfaces, and (3) that high thermal stability may be obtained without forfeiture of catalytic performance at low temperatures.

Biased mutation-assembling: an efficient method for rapid directed evolution through simultaneous mutation accumulation

We have developed an efficient optimization technique, 'biased mutation-assembling', for improving protein properties such as thermostability. In this strategy, a mutant library is constructed using the overlap extension polymerase chain reaction technique with DNA fragments from wild-type and phenotypically advantageous mutant genes, in which the number of mutations assembled in the wild-type gene is stochastically controlled by the mixing ratio of the mutant DNA fragments to wild-type fragments. A high mixing ratio results in a mutant composition biased to favor multiple-point mutants. We applied this strategy to improve the thermostability of prolyl endopeptidase from Flavobacterium meningosepticum as a case study and found that the proportion of thermostable mutants in a library increased as the mixing ratio was increased. If the proportion of thermostable mutants increases, the screening effort needed to find them should be reduced. Indeed, we isolated a mutant with a 1200-fold longer activity half-life at 60 degrees C than that of wild-type prolyl endopeptidase after screening only 2000 mutants from a library prepared with a high mixing ratio. Our results indicate that an aggressive accumulation of advantageous mutations leads to an increase in the quality of the mutant library and a reduction in the screening effort required to find superior mutants.

Modified colorimetric assay for uricase activity and a screen for mutant Bacillus subtilis uricase genes following StEP mutagenesis

This study describes a modified colorimetric assay for uricase activity in flexible 96-well microtiter plates using the uricase/uric acid/horseradish peroxidase/4-aminoantipyrine/3,5-dichloro-2-hydroxybenzene sulfonate colorimetric reaction. The utility of this assay was demonstrated in a screen for mutant uricase enzymes derived from the uricase gene of the thermophilic bacterium Bacillus subtilis by a modified staggered extension process (StEP) mutagenesis. An Escherichia coli library of StEP-derived uricase mutant clones was screened yielding two identical active mutant uricase genes. Two motifs conserved in eukaryotic and prokaryotic uricases are highly conserved in the mutant uricase. The mutant uricase protein was found to exhibit high uricase activity (13.1 U.mg(-1)). Finally, the modified colorimetric method is much more efficient than the conventional ones and greatly reduces assay time from 4 days to less than 20 h.

High-level expression of Aliciclobacillus acidocaldarius thioredoxin in Pichia pastoris and Bacillus subtilis

Thioredoxins are ubiquitous proteins which catalyze the reduction of disulfide bridges on target proteins and are involved in many cellular reactions. In a previous work, a thioredoxin from the thermophilic organism Aliciclobacillus acidocaldarius (Alitrx) was purified, characterized, and its gene expressed in Escherichia coli. In order to produce larger quantities of Alitrx, the protein has been expressed in the methylotrophic yeast Pichia pastoris and in the gram positive bacteria Bacillus subtilis. The growth conditions of strains showing high-level expression of Alitrx were optimized for both systems in shake-flask cultures. Active proteins were secreted in the culture media at a level of approximately 0.9 and 0.5 g/l, respectively, for P. pastoris and B. subtilis. The proteins were purified almost to homogeneity by a thermal precipitation procedure, with a 90-fold and 50-fold higher total yield with respect to that obtained with the same protein expressed in E. coli. The results indicate that either of these two systems could be utilized as a host for large-scale production of recombinant Alitrx.

Computational and experimental analysis of DNA shuffling

We describe a computational model of DNA shuffling based on the thermodynamics and kinetics of this process. The model independently tracks a representative ensemble of DNA molecules and records their states at every stage of a shuffling reaction. These data can subsequently be analyzed to yield information on any relevant metric, including reassembly efficiency, crossover number, type and distribution, and DNA sequence length distributions. The predictive ability of the model was validated by comparison to three independent sets of experimental data, and analysis of the simulation results led to several unique insights into the DNA shuffling process. We examine a tradeoff between crossover frequency and reassembly efficiency and illustrate the effects of experimental parameters on this relationship. Furthermore, we discuss conditions that promote the formation of useless "junk" DNA sequences or multimeric sequences containing multiple copies of the reassembled product. This model will therefore aid in the design of optimal shuffling reaction conditions.

Molecular basis for the subtype discrimination of the estrogen receptor-beta-selective ligand, diarylpropionitrile

Although the two subtypes of the human estrogen receptor (ER), ERalpha and ERbeta, share only 56% amino acid sequence identity in their ligand binding domain (LBD), the residues that surround the ligand are nearly identical; nevertheless, subtype-selective ligands are known. To understand the molecular basis by which diarylpropionitrile (DPN), an ERbeta-selective ligand, is able to discriminate between the two ERs, we examined its activity on ER mutants and chimeric constructs generated by DNA shuffling. The N-terminal region of the ERbeta LBD (through helix 6) appears to be fully responsible for the ERbeta selectivity of DPN. In fact, a single ERalpha point mutation (L384M) was largely sufficient to switch the DPN response of this ER to that of the ERbeta type, but residues in helix 3 are also important in achieving the full ERbeta selectivity of DPN. Using molecular modeling, we found an energetically favorable fit for the S-DPN enantiomer in ERbeta, in which the proximal phenol mimics the A ring of estradiol, and the nitrile engages in stabilizing interactions with residues in the ligand-binding pocket of ERbeta. Our findings highlight that a limited number of critical interactions of DPN with the ERbeta ligand-binding pocket underlie its ER subtype-selective character.

Directed evolution to increase camptothecin sensitivity of human DNA topoisomerase I

BACKGROUND

Human DNA topoisomerase I (top1) relaxes DNA supercoiling during basic nuclear processes. The enzyme is the main target of antitumor agents, such as camptothecins (CPT), that transform top1 into a DNA-damaging agent.

RESULTS

By directed evolution of a C-terminal portion, we selected human top1 mutants that were 22-28-fold more CPT-sensitive than wild-type top1 in Saccharomyces cerevisiae cells. The evolved enzymes showed unique mutation patterns and were more processive in plasmid relaxation assays. A top1 mutant had only two amino acid changes in the linker domain, one of which may change a linker/core domain contact surface. The mutant stimulated DNA cleavage to higher levels than the wild-type enzyme and was more sensitive to CPT in a cleavage assay. Moreover, the mutant was more CPT-sensitive than wild-type top1 in a repair-deficient yeast strain.

CONCLUSIONS

Mutations in the linker domain can affect DNA binding and CPT sensitivity of human top1. Such drug-hypersensitive topoisomerases may be useful in developing DNA cutters with high cell lethality and in new drug discovery programs.

Applied molecular evolution of O6-benzylguanine-resistant DNA alkyltransferases in human hematopoietic cells

Applied molecular evolution is a rapidly developing technology that can be used to create and identify novel enzymes that nature has not selected. An important application of this technology is the creation of highly drug-resistant enzymes for cancer gene therapy. Seventeen O(6)-alkylguanine-DNA alkyltransferase (AGT) mutants highly resistant to O(6)-benzylguanine (BG) were identified previously by screening 8 million variants, using genetic complementation in Escherichia coli. To examine the potential of these mutants for use in humans, the sublibrary of AGT clones was introduced to human hematopoietic cells and stringently selected for resistance to killing by the combination of BG and 1,3-bis(2-chloroethyl)-1-nitrosourea. This competitive analysis between the mutants in human cells revealed three AGT mutants that conferred remarkable resistance to the combination of BG and 1,3-bis(2-chloroethyl)-1-nitrosourea. Of these, one was recovered significantly more frequently than the others. Upon further analysis, this mutant displayed a level of BG resistance in human hematopoietic cells greater than that of any previously reported mutant.

Predicting crossover generation in DNA shuffling

We introduce a quantitative framework for assessing the generation of crossovers in DNA shuffling experiments. The approach uses free energy calculations and complete sequence information to model the annealing process. Statistics obtained for the annealing events then are combined with a reassembly algorithm to infer crossover allocation in the reassembled sequences. The fraction of reassembled sequences containing zero, one, two, or more crossovers and the probability that a given nucleotide position in a reassembled sequence is the site of a crossover event are estimated. Comparisons of the predictions against experimental data for five example systems demonstrate good agreement despite the fact that no adjustable parameters are used. An in silico case study of a set of 12 subtilases examines the effect of fragmentation length, annealing temperature, sequence identity and number of shuffled sequences on the number, type, and distribution of crossovers. A computational verification of crossover aggregation in regions of near-perfect sequence identity and the presence of synergistic reassembly in family DNA shuffling is obtained.

Protein engineering of subtilisin

The serine protease subtilisin is an important industrial enzyme as well as a model for understanding the enormous rate enhancements affected by enzymes. For these reasons along with the timely cloning of the gene, ease of expression and purification and availability of atomic resolution structures, subtilisin became a model system for protein engineering studies in the 1980s. Fifteen years later, mutations in well over 50% of the 275 amino acids of subtilisin have been reported in the scientific literature. Most subtilisin engineering has involved catalytic amino acids, substrate binding regions and stabilizing mutations. Stability has been the property of subtilisin which has been most amenable to enhancement, yet perhaps least understood. This review will give a brief overview of the subtilisin engineering field, critically review what has been learned about subtilisin stability from protein engineering experiments and conclude with some speculation about the prospects for future subtilisin engineering.

Approaches for deciphering the structural basis of low temperature enzyme activity

An increasing number of enzymes active at low temperature are being studied to help determine the structural features important for cold-activity. This review examines the diversity of prokaryotic cold-active enzymes and the features proposed to account for low temperature activity. We then consider the difficulty of identifying the key structural features needed for cold-activity and the need to compare enzymes having different temperature optima from phylogenetically related organisms to determine features responsible for low temperature activity. In addition to studying naturally occurring enzymes, directed evolution experiments are discussed as methods for examining the proposed mechanisms influencing the thermal dependence of activity.

Subtilisins of Bacillus spp. hydrolyze keratin and allow growth on feathers

Keratinase is a serine protease produced by Bacillus licheniformis PWD-1 that effectively degrades keratin and confers the ability to grow on feathers to a protease-deficient B. subtilis strain. Studies presented herein demonstrate that B. licheniformis Carlsberg strain NCIMB 6816, which produces the well-characterized serine protease subtilisin Carlsberg, also degrades and grows on feathers. The PWD-1 and Carlsberg strains showed a similar time-course of enzyme production, and the purified serine proteases have similar enzymatic properties on insoluble azokeratin and soluble FITC-casein. Kinetic analysis of both enzymes demonstrated that they have high specificity for aromatic and hydrophobic amino acids in the P1 substrate position, although keratinase discriminates more than subtilisin Carlsberg against charged residues at this site. Nucleotide sequence analysis of the serine protease genes from B. licheniformis strains PWD-1, Carlsberg NCIMB 6816, ATCC 12759, and NCIMB 10689 showed that the kerA-encoded protease of PWD-1 differs from the others only by having V222, rather than A222, near the active site serine S220. Further, high-level expression of subE-encoded subtilisin from B. subtilis (78% similar to subtilisin Carlsberg) also confers growth on feathers on a protease-deficient B. subtilis strain. While strain PWD-1 and the kerA protease efficiently degrade keratin, keratin hydrolysis and growth on feathers is a property that can be conferred by appropriate expression of the major subtilisins, including the industrially produced enzymes.

Directed evolution of an enantioselective lipase

BACKGROUND

The biocatalytic production of enantiopure compounds is of steadily increasing importance to the chemical and biotechnological industry. In most cases, however, it is impossible to identify an enzyme that possesses the desired enantioselectivity. Therefore, there is a strong need to create by molecular biological methods novel enzymes which display high enantioselectivity.

RESULTS

A bacterial lipase from Pseudomonas aeruginosa (PAL) was evolved to catalyze with high enantioselectivity the hydrolysis of the chiral model substrate 2-methyldecanoic acid p-nitrophenyl ester. Successive rounds of random mutagenesis by ep-PCR and saturation mutagenesis resulted in an increase in enantioselectivity from E=1.1 for the wild-type enzyme to E=25.8 for the best variant which carried five amino acid substitutions. The recently solved three-dimensional structure of PAL allowed us to analyze the structural consequences of these substitutions.

CONCLUSIONS

A highly enantioselective lipase was created by increasing the flexibility of distinct loops of the enzyme. Our results demonstrate that enantioselective enzymes can be created by directed evolution, thereby opening up a large area of novel applications in biotechnology.

Modeling DNA mutation and recombination for directed evolution experiments

Directed evolution experiments rely on the cyclical application of mutagenesis, screening and amplification in a test tube. They have led to the creation of novel proteins for a wide range of applications. However, directed evolution currently requires an uncertain, typically large, number of labor intensive and expensive experimental cycles before proteins with improved function are identified. This paper introduces predictive models for quantifying the outcome of the experiments aiding in the setup of directed evolution for maximizing the chances of obtaining DNA sequences encoding enzymes with improved activities. Two methods of DNA manipulation are analysed: error-prone PCR and DNA recombination. Error-prone PCR is a DNA replication process that intentionally introduces copying errors by imposing mutagenic reaction conditions. The proposed model calculates the probability of producing a specific nucleotide sequence after a number of PCR cycles. DNA recombination methods rely on the mixing and concatenation of genetic material from a number of parent sequences. This paper focuses on modeling a specific DNA recombination protocol, DNA shuffling. Three aspects of the DNA shuffling procedure are modeled: the fragment size distribution after random fragmentation by DNase I, the assembly of DNA fragments, and the probability of assembling specific sequences or combinations of mutations. Results obtained with the proposed models compare favorably with experimental data.

Evolution of evolvability

Genomic sequence data provide evidence for a common origin of life and for its evolution by genetic variation via mutation and recombination. This paper discusses the fundamental dialectic paradigm of evolution--stability versus variability--at the crossroads of molecular genetics, population genetics, ecology, and the emerging science of experimental evolution. Experimental evolution of molecules, viruses, and bacteria can be used not only to test some basic evolutionary hypotheses but also to create new organisms for applications in biotechnology, agriculture, and medicine.

Modeling DNA shuffling

In vitro evolution is a new, important laboratory method to evolve molecules with desired properties. It has been used in a variety of biological studies and drug development. In this paper, we study one important mutagenesis method used in in vitro evolution experiments called DNA shuffling. We construct a mathematical model for DNA shuffling and study the properties of molecules after DNA shuffling experiments based on this model. The model for DNA shuffling consists of two parts. First we apply the Lander-Waterman model for physical mapping by fingerprinting random clones to model the distribution of regions that can be reassembled through DNA shuffling. Then we present a model for recombination between different DNA species with different mutations. We compare our theoretical results with experimental data. Finally we propose novel applications of the theoretical results to the optimal design of DNA shuffling experiments and to physical mapping using DNA shuffling.

In vitro evolution of thermostable p53 variants

The tumor suppressor p53 is conformationally unstable at physiological temperature. Even the activated p53delta30 variant, which lacks the self-inhibiting carboxy terminal domain, has a half-life of only 8 min at 37 degrees C in vitro. We have developed a genetic approach to identify p53 variants that stabilize the active conformation. The human p53delta30 gene was randomly mutated, and the resulting library was expressed in Escherichia coli under conditions that apparently denatured the parental protein. Stable p53 variants were identified based on their ability to specifically bind a p53 consensus site. The initial thermostable variants were randomly recombined by DNA shuffling, and substitutions that were functionally additive or synergistic were identified in a second more stringent round of screening. The DNA binding activity of N239Y/N268D/E336V p53delta30 variant has a half-life of 100 min at 37 degrees C, 12 times longer than that of the parental protein. The thermostable variants should be more amenable to crystallographic studies and more effective in gene therapies than the wild-type protein.

Evolution of differential substrate specificities in Mu class glutathione transferases probed by DNA shuffling

A library of variant enzymes was created by combined shuffling of the DNA encoding the human Mu class glutathione transferases GST M1-1 and GST M2-2. The parental GSTs are 84 % sequence identical at the protein level, but their specific activities with the substrates aminochrome and 2-cyano-1,3-dimethyl-1-nitrosoguanidine (cyanoDMNG) differ by more than 100-fold. Aminochrome is of particular interest as an oxidation product of dopamine and of possible significance in the etiology of Parkinson's disease, and cyanoDMNG is a model for genotoxic and potentially carcinogenic nitroso compounds. GST M2-2 has at least two orders of magnitude higher catalytic activity with both of the substrates than any of the other known GSTs, including GST M1-1. The DNA library of variant Mu class GST sequences contained "mosaic" structures composed of alternating segments of both parental sequences. All clones contained the 5'-end of a GST M1-1 clone optimized for high-level expression in Escherichia coli. The remainder of the sequences derived from segments of GST M2-2 and GST M1-1 DNA. All of the clones analyzed contained between two and seven distinct DNA segments. In addition, each clone contained an average of approximately one point mutation. None of the library clones analyzed was identical with either of the two parental structures. Variant GST sequences were expressed in E. coli, and their enzymatic activities with aminochrome, cyanoDMNG, and 1-chloro-2,4-dinitrobenzene (CDNB) were determined in bacterial lysates. Such screening of more than 70 clones demonstrated a continuous range of activities covering at least two orders of magnitude for each of the substrates. For a given clone, the activities with aminochrome and cyanoDMNG, in spite of their different chemistries, were clearly correlated, whereas no strong correlation was found with CDNB. This functional correlation suggests a common structural basis for the enzymatic mechanisms for conjugation of aminochrome and denitrosation of cyanoDMNG. From an evolutionary perspective, the results show that recombination of segments from homologous proteins gives rise to a large proportion of functionally competent proteins with a range of activities. The data support the proposal that natural evolution of protein functions may involve recombination of DNA segments followed by selection for advantageous functional properties of the resulting proteins. Clearly, the same approach can be utilized in the engineering of proteins displaying novel functions by in vitro evolution.

Improving enzymes for cancer gene therapy

New techniques now make it feasible to tailor enzymes for cancer gene therapy. Novel enzymes with desired properties can be created and selected from vast libraries of mutants containing random substitutions within catalytic domains. In this review, we first consider genes for the ablation of tumors, namely, genes that have been mutated (or potentially can be mutated) to afford enhanced activation of prodrugs and increased sensitization of tumors to specific chemotherapeutic agents. We then consider genes that have been mutated to provide better protection of normal host tissues, such as bone marrow, against the toxicity of specific chemotherapeutic agents. Expression of the mutant enzyme could render sensitive tissues, such as bone marrow, more resistant to specific cytotoxic agents.

Directed evolution converts subtilisin E into a functional equivalent of thermitase

We used directed evolution to convert Bacillus subtilis subtilisin E into an enzyme functionally equivalent to its thermophilic homolog thermitase from Thermoactinomyces vulgaris. Five generations of random mutagenesis, recombination and screening created subtilisin E 5-3H5, whose half-life at 83 degrees C (3.5 min) and temperature optimum for activity (Topt, 76 degrees C) are identical with those of thermitase. The Topt of the evolved enzyme is 17 degrees C higher and its half-life at 65 degrees C is >200 times that of wild-type subtilisin E. In addition, 5-3H5 is more active towards the hydrolysis of succinyl-Ala-Ala-Pro-Phe-p-nitroanilide than wild-type at all temperatures from 10 to 90 degrees C. Thermitase differs from subtilisin E at 157 amino acid positions. However, only eight amino acid substitutions were sufficient to convert subtilisin E into an enzyme equally thermostable. The eight substitutions, which include known stabilizing mutations (N218S, N76D) and also several not previously reported, are distributed over the surface of the enzyme. Only two (N218S, N181D) are found in thermitase. Directed evolution provides a powerful tool to unveil mechanisms of thermal adaptation and is an effective and efficient approach to increasing thermostability without compromising enzyme activity.

Delineation of an evolutionary salvage pathway by compensatory mutations of a defective lysozyme

Model-free approaches (random mutagenesis, DNA shuffling) in combination with more "rational," three-dimensional information-guided randomization have been used for directed evolution of lysozyme activity in a defective T4 lysozyme mutant. A specialized lysozyme cloning vector phage, derived from phage lambda, depends upon T4 lysozyme function for its ability to form plaques. The substitution W138P in T4 lysozyme totally abolishes its plaque-forming ability. Compensating mutations in W138P T4 lysozyme after sequential random mutagenesis of the whole gene as well as after targeted randomization of residues in the vicinity of Trp138 were selected. In a second stage, these mutations were randomly recombined by the recombinatorial PCR method of DNA shuffling. Shuffled and selected W138P T4 lysozyme variants provide the hybrid lambda phage with sufficient lysozyme activity to produce normal-size plaques, even at elevated temperature (42 degrees C). The individual mutations with the highest compensatory information for W138P repair are the substitutions A146F and A146M, selected after targeted randomization of three residues in the neighborhood of Trp138 by combinatorial mutagenesis. The best evolved W138P T4 lysozymes, however, accumulated mutations originating from both randomly mutagenized as well as target-randomized variants.

Combinatorial protein design by in vitro recombination

DNA recombination is a powerful engine for the creation of new phenotypes. Recently, methods for in vitro DNA recombination (DNA shuffling) have been developed and applied to the evolution of novel molecules in the laboratory. An exciting new development is the shuffling of homologous genes to create diversity for directed evolution.

Molecular evolution by staggered extension process (StEP) in vitro recombination

We have developed a simple and efficient method for in vitro mutagenesis and recombination of polynucleotide sequences. The staggered extension process (StEP) consists of priming the template sequence(s) followed by repeated cycles of denaturation and extremely abbreviated annealing/polymerase-catalyzed extension. In each cycle the growing fragments anneal to different templates based on sequence complementarity and extend further. This is repeated until full-length sequences form. Due to template switching, most of the polynucleotides contain sequence information from different parental sequences. The method is demonstrated by the recombination of two genes encoding thermostable subtilisins carrying two phenotypic markers separated by 113 base pairs and eight other point mutation markers. To demonstrate its utility for directed evolution, we have used StEP to recombine a set of five thermostabilized subtilisin E variants identified during a single round of error-prone PCR mutagenesis and screening. Screening the StEP-recombined library yielded an enzyme whose half-life at 65 degrees C is 50 times that of wild-type subtilisin E.

Engineering of a cold-adapted protease by sequential random mutagenesis and a screening system

A cold-adapted protease subtilisin was successfully isolated by evolutionary engineering based on sequential in vitro random mutagenesis and an improved method of screening (H. Kano, S. Taguchi, and H. Momose, Appl. Microbiol. Biotechnol. 47:46-51, 1997). The mutant subtilisin, termed m-63, exhibited a catalytic efficiency (expressed as the kcat/Km value) 100% higher than that of the wild type at 10 degrees C when N-succinyl-L-Ala-L-Ala-L-Pro-L-Phe-p-nitroanilide was used as a synthetic substrate. This cold adaptation was achieved with three mutations, Val to Ile at position 72 (V72I), Ala to Thr at position 92 (A92T), and Gly to Asp at position 131 (G131D), and it was found that an increase in substrate affinity (i.e., a decreased Km value) was mostly responsible for the increased activity. Analysis of kinetic parameters revealed that the V72I mutation contributed negatively to the activity but that the other two mutations, A92T and G131D, overcame the negative contribution to confer the 100% increase in activity. Besides suppression of the activity-negative mutation (V72I) by A92T and G131D, suppression of structural stability was observed in measurements of activity retention at 60 degrees C and circular dichroism spectra at 10 degrees C.

Random-priming in vitro recombination: an effective tool for directed evolution

A simple and efficient method for in vitro mutagenesis and recombination of polynucleotide sequences is reported. The method involves priming template polynucleotide(s) with random-sequence primers and extending to generate a pool of short DNA fragments which contain a controllable level of point mutations. The fragments are reassembled during cycles of denaturation, annealing and further enzyme-catalyzed DNA polymerization to produce a library of full-length sequences. Screening or selecting the expressed gene products leads to new variants with improved functions, as demonstrated by the recombination of genes encoding different thermostable subtilisins in order to obtain enzymes more stable than either parent.

Strategies for the in vitro evolution of protein function: enzyme evolution by random recombination of improved sequences

Sets of genes improved by directed evolution can be recombined in vitro to produce further improvements in protein function. Recombination is particularly useful when improved sequences are available; costs of generating such sequences, however, must be weighed against the costs of further evolution by sequential random mutagenesis. Four genes encoding para-nitrobenzyl (pNB) esterase variants exhibiting enhanced activity were recombined in two cycles of high-fidelity DNA shuffling and screening. Genes encoding enzymes exhibiting further improvements in activity were analyzed in order to elucidate evolutionary processes at the DNA level and begin to provide an experimental basis for choosing in vitro evolution strategies and setting key parameters for recombination. DNA sequencing of improved variants from the two rounds of DNA shuffling confirmed important features of the recombination process: rapid fixation and accumulation of beneficial mutations from multiple parent sequences as well as removal of silent and deleterious mutations. The five to sixfold further enhancement of total activity towards the para-nitrophenyl (pNP) ester of loracarbef was obtained through recombination of mutations from several parent sequences as well as new point mutations. Computer simulations of recombination and screening illustrate the trade-offs between recombining fewer parent sequences (in order to reduce screening requirements) and lowering the potential for further evolution. Search strategies which may substantially reduce screening requirements in certain situations are described.

Combinatorial protein design: strategies for screening protein libraries

Powerful strategies for screening protein libraries further strengthen the arguments for applying 'irrational' approaches to understanding and designing new proteins. Developments during the past year include the application of functional complementation and automation to reduce screening loads, as well as the use of computerized data acquisition to characterize whole protein libraries rather than just selected individuals.