Digestion-after-shuffling: a method for combinatorial library construction

We have developed a convenient method for family shuffling of amino acid sequences, termed digestion-after-shuffling. After DNA shuffling of homologous genes, plasmid mixture is extracted from a library and used for several double digestions with restriction enzymes. For each double digestion, two restriction enzymes are selected, corresponding to the single restriction sites of different parental genes. After digestions, fragments with expected sizes are obtained by gel purification and religated to construct recombinant plasmids. Thus, the obtained genes should be chimeras and have at least two restriction sites originating from different parental sequences.     

Origin of new genes: evidence from experimental and computational analyses

Exon shuffling is an essential molecular mechanism for the formation of new genes. Many cases of exon shuffling have been reported in vertebrate genes. These discoveries revealed the importance of exon shuffling in the origin of new genes. However, only a few cases of exon shuffling were reported from plants and invertebrates, which gave rise to the assertion that the intron-mediated recombination mechanism originated very recently. We focused on the origin of new genes by exon shuffling and retroposition. We will first summarize our experimental work, which revealed four new genes in Drosophila, plants, and humans. These genes are 10⁶ to 10⁸ million years old. The recency of these genes allows us to directly examine the origin and evolution of genes in detail. These observations show firstly the importance of exon shuffling and retroposition in the rapid creation of new gene structures. They also show that the resultant chimerical structures appearing as mosaic proteins or as retroposed coding structures with novel regulatory systems, often confer novel functions. Furthermore, these newly created genes appear to have been governed by positive Darwinian selection throughout their history, with rapid changes of amino acid sequence and gene structure in very short periods of evolution. We further analyzed the distribution of intron phases in three non-vertebrate species, Drosophila melanogaster, Caenorhabditis elegans, and Arabidosis thaliana, as inferred from their genome sequences. As in the case of vertebrate genes, we found that intron phases in these species are unevenly distributed with an excess of phase zero introns and a significant excess of symmetric exons. Both findings are consistent with the requirements for the molecular process of exon shuffling. Thus, these non-vertebrate genomes may have also been strongly impacted by exon shuffling in general.     

Protein engineering with synthetic Escherichia coli amber suppressor genes

We have constructed synthetic genes encoding different Escherichia coli suppressor tRNAs for use in amino acid substitution studies and protein engineering. We used oligonucleotides to assemble the genes for different tRNAs with the anticodon 5' CTA 3'. The suppressor genes are expressed from a synthetic promoter derived from the promoter sequence of the E. coli lipoprotein gene. The genes have been used to suppress an amber mutation in a protein coding sequence, and the resulting altered protein has been subjected to sequence analysis to determine the nature of the amino acid inserted at the amber site. Twelve amino acids can now be added in response to the amber codon. We have employed these suppressors to study amino acid substitutions in the lac repressor.     

Correlation of gene and protein structures in the FXYD family proteins

The FXYD family proteins are auxiliary subunits of the Na,K-ATPase, expressed primarily in tissues that specialize in fluid or solute transport, or that are electrically excitable. These proteins range in size from about 60 to 160 amino acid residues, and share a core homology of 35 amino acid residues in and around a single transmembrane segment. Despite their relatively small sizes, they are all encoded by genes with six to nine small exons. We show that the helical secondary structures of three FXYD family members, FXYD1, FXYD3, and FXYD4, determined in micelles by NMR spectroscopy, reflect the structures of their corresponding genes. The coincidence of helical regions, and connecting segments, with the positions of intron-exon junctions in the genes, support the hypothesis that the FXYD proteins may have been assembled from discrete structural modules through exon shuffling.     

Shuffling of amino acid sequence: an important control in synthetic peptide studies of nucleic acid-binding domains. Binding properties of fragments of a conserved eukaryotic RNA binding motif.

We have used synthetic peptides to study a conserved RNA binding motif in yeast poly(A)-binding protein. Two peptides, 45 and 44 amino acids in length, corresponding to amino and carboxyl halves of a 90-amino acid RNA-binding domain in the protein were synthesized. While the amino-terminal peptide had no significant affinity for nucleic acids, the carboxyl-terminal peptide-bound nucleic acids with similar characteristics to that for the entire 577 residue yeast poly(A)-binding protein. In 100 mM NaCl, the latter peptide retained over 50% of the intrinsic binding free energy of the protein, as well as, similar RNA versus DNA binding specificity. However, shuffling of the sequence of this 44 residue peptide had surprisingly little effect on its nucleic acid binding properties suggesting the overriding importance of amino acid composition as opposed to primary sequence. Deletion studies on the 44 residue peptide with the "correct" sequence succeeded in identifying amino acids important for conferring RNA specificity and for increasing our understanding of the molecular basis for nucleic acid binding by synthetic peptides. The shuffled peptide study, however, clearly indicates that considerable caution must be exercised before extrapolating results of structure/function studies on synthetic peptide analogues to the parent protein.     

Correcting errors in synthetic DNA through consensus shuffling

Although efficient methods exist to assemble synthetic oligonucleotides into genes and genomes, these suffer from the presence of 1–3 random errors/kb of DNA. Here, we introduce a new method termed consensus shuffling and demonstrate its use to significantly reduce random errors in synthetic DNA. In this method, errors are revealed as mismatches by re-hybridization of the population. The DNA is fragmented, and mismatched fragments are removed upon binding to an immobilized mismatch binding protein (MutS). PCR assembly of the remaining fragments yields a new population of full-length sequences enriched for the consensus sequence of the input population. We show that two iterations of consensus shuffling improved a population of synthetic green fluorescent protein (GFPuv) clones from ~60 to >90% fluorescent, and decreased errors 3.5- to 4.3-fold to final values of ~1 error per 3500 bp. In addition, two iterations of consensus shuffling corrected a population of GFPuv clones where all members were non-functional, to a population where 82% of clones were fluorescent. Consensus shuffling should facilitate the rapid and accurate synthesis of long DNA sequences.     

Sequence shuffle controls morphological consequences in a self-assembling tetrapeptide.

Peptide and protein self-assembly is a well-studied phenomenon in chemistry and biology, where nanoscopic building blocks exhibit rapid self-association to reveal supramolecular aggregates of defined structural features. These superstructures are stabilized by hydrophobic interactions, hydrogen bonding and a host of other noncovalent interactions. Thus, amino acid side chains in the primary structure hold importance in dictating secondary structures and preference for particular conformational signatures in peptide aggregates. This report describes contrasting nanoscale morphologies in antamanide-derived synthetic tetrapeptide mutants, which are composed by shuffling only two amino acids: phenylalanine and proline. Remarkable differences in ultrastructures in primary sequence-shuffled tetrapeptides suggest dissimilar aggregational pathways due to context-dependent location of proline and phenylalanine residues with respect to one another.     

Directed evolution of alpha-aspartyl dipeptidase from Salmonella typhimurium.

Model-free approaches (error-prone PCR to introduce random mutations, DNA shuffling to combine positive mutations, and screening of the resultant mutant libraries) have been used to enhance the catalytic activity and thermostability of alpha-aspartyl dipeptidase from Salmonella typhimurium, which is uniquely able to hydrolyze Asp-X dipeptides (where X is any amino acid) and one tripeptide (Asp-Gly-Gly). Under double selective pressures of activity and thermostability, through two rounds of error-prone PCR and three sequential generations of DNA shuffling, coupled with screening, a mutant pepEM3074 with approximately 47-fold increased enzyme activity compared with its wild-type parent was obtained. Moreover, the stability of pepEM3074 is increased significantly. Three amino acid substitutions (Asn89His, Gln153Glu, and Leu205Arg), two of them are near the active site and substrate binding pocket, were identified by sequencing the genes encoding this evolved enzyme. The mechanism of the enhancement of activity and stability was analyzed in this paper.     

Autosomal Similarity Revealed by Eukaryotic Genomic Comparison

To describe eukaryotic autosomes quantitatively and determine differences between them in terms of amino acid sequences of genes, functional classification of proteins, and complete DNA sequences, we applied two theoretical methods, the Proteome-vector method and the function of degree of disagreement (FDOD) method, that are based on function and sequence similarity respectively, to autosomes from nine eukaryotes. No matter what aspect of the autosome is considered, the autosomal differences within each organism were less than that between species. Our results show that eukaryotic autosomes resemble each other within a species while those from different organisms differ. We propose a hypothesis (named intra-species autosomal random shuffling) as an explanation for our results and suggest that lateral gene transfer (LGT) did not occur frequently during the evolution of eukarya.     

Construction of a chimeric aminopeptidase by a combination of gene shuffling and mutagenesis

Three chimeric genes were constructed by gene shuffling of aminopeptidases from Aeromonas caviae and Vibrio proteolyticus. Although expressed chimeric enzymes formed inclusion bodies in Escherichia coli, the introduction of two amino acid mutations into the chimeric genes by site-saturated mutagenesis and a random mutation on error-prone PCR resulted in solubilization of the chimeric enzyme. In addition, active chimeric enzyme showed a different thermostability and thermoactivity to parental enzymes.     

Heterologous overexpression of a mutant termite cellulase gene in Escherichia coli by DNA shuffling of four orthologous parental cDNAs

Among cellulase genes, those of animals are known for their difficulty in overexpression. We constructed a chimeric library by family shuffling of endo-beta-1,4-glucanase genes from four different termite species (Reticulitermes speratus, Nasutitermes takasagoensis, Coptotermes formosanus, and Coptotermes acinaciformis) sharing 78.5-96% homology in amino acid sequence. The constructed library was screened by Congo red plate assay combined with 96-well micro-enzyme assay, and clones showing enhanced CMCase activities were obtained. The mutated genes were overexpressed in Escherichia coli intracellularly as an active form. The endo-beta-1,4-glucanase (CMCase) activity in soluble fractions of E. coli harboring the mutant genes was 20-30 fold higher than that of wild-type genes. The mutant enzyme showed high activity against CMC and properties similar to those of the native enzymes.     

Identification of amino acids important for substrate specificity in sucrose transporters using gene shuffling

Plant sucrose transporters (SUTs) are H(+)-coupled uptake transporters. Type I and II (SUTs) are phylogenetically related but have different substrate specificities. Type I SUTs transport sucrose, maltose, and a wide range of natural and synthetic α- and β-glucosides. Type II SUTs are more selective for sucrose and maltose. Here, we investigated the structural basis for this difference in substrate specificity. We used a novel gene shuffling method called synthetic template shuffling to introduce 62 differentially conserved amino acid residues from type I SUTs into OsSUT1, a type II SUT from rice. The OsSUT1 variants were tested for their ability to transport the fluorescent coumarin β-glucoside esculin when expressed in yeast. Fluorescent yeast cells were selected using fluorescence-activated cell sorting (FACS). Substitution of five amino acids present in type I SUTs in OsSUT1 was found to be sufficient to confer esculin uptake activity. The changes clustered in two areas of the OsSUT1 protein: in the first loop and the top of TMS2 (T80L and A86K) and in TMS5 (S220A, S221A, and T224Y). The substrate specificity of this OsSUT1 variant was almost identical to that of type I SUTs. Corresponding changes in the sugarcane type II transporter ShSUT1 also changed substrate specificity, indicating that these residues contribute to substrate specificity in type II SUTs in general.     

Three amino acid changes contribute markedly to the thermostability of β-glucosidase BglC from Thermobifida fusca

Thermostability of β-glucosidase was enhanced by family shuffling, site saturation mutagenesis, and site-directed mutagenesis. Family shuffling was carried out based on β-glucosidase BglC from Thermobifida fusca and β-glucosidase BglB from Paebibacillus polymxyxa with the help of synthetic primers. High-throughput screening revealed mutants with higher thermostability than both parental enzymes. The most thermostable mutant VM2 containing three key amino acid changes in L444Y/G447S/A433V had a 144-fold increase in half-life of inactivation as compared to the parental enzyme BglC. The mutant VM2 showed 28% and 94% increase in k(cat) towards p-nitrophenyl-β-D-glucopyranoside (pNPG) and cellobiose, respectively. The mutant with enhanced stability would facilitate the recycle of β-glucosidase in the bioconversion of cellulosic biomass.     

Growth factor engineering by degenerate homoduplex gene family recombination

There is great interest in engineering human growth factors as potential therapeutic agonists and antagonists. We approached this goal with a synthetic DNA recombination method. We aligned a pool of "top-strand" oligonucleotides incorporating polymorphisms from mammalian genes encoding epidermal growth factor (EGF) using multiple polymorphic "scaffold" oligonucleotides. Top strands were then linked by gap filling and ligation. This approach avoided heteroduplex annealing in the linkage of highly degenerate oligonucleotides and thus achieved completely random recombination. Cloned genes from a human-mouse chimeric library captured every possible permutation of the parental polymorphisms, creating an apparently complete recombined gene-family library, which has not been previously described. This library yielded a chimeric protein whose agonist activity was enhanced 123-fold. A second library from five mammalian EGF homologs possessed the highest reported recombination density (1 crossover per 12.4 bp). The five-homolog library yielded the strongest-binding hEGF variant yet reported. In addition, it contained strongly binding EGF variants with antagonist properties. Our less biased approach to DNA shuffling should be useful for the engineering of a wide variety of proteins.     

Characterization of the gene for human plasminogen, a key proenzyme in the fibrinolytic system

The organization and structure of the gene coding for plasminogen has been determined by a combination of in vitro amplification of leukocyte DNA from normal individuals and isolation of unique clones from three different human genomic libraries. These clones were characterized by restriction mapping, Southern blotting, and DNA sequencing. The gene for human plasminogen spanned about 52.5 kilobases of DNA and consisted of 19 exons separated by 18 introns. DNA sequence analysis revealed that the five kringle structures in plasminogen were coded by two exons. The nucleotides in the introns at the intron-exon boundaries were GT-AG analogous to those found in other eukaryotic genes. Three polyadenylation sites for plasminogen mRNA were also identified. When the amino acid sequences deduced from the genomic DNA and cDNAs of plasminogen were compared with that of the plasma protein determined by amino acid sequence analysis, an apparent amino acid polymorphism was observed in several positions of the polypeptide chain. Nucleotide sequence analysis of the amplified genomic DNAs and genomic clones also revealed that the plasminogen gene was very closely related to several other proteins, including apolipoprotein(a). This protein may have evolved via duplication and exon shuffling of the plasminogen gene. The presence of another plasminogen-related gene(s) in the human genomic library was also observed.     

Reduction of wobble-position GC bases in Corynebacteria genes and enhancement of PCR and heterologous expression

Corynebacteria codon usage exhibits an overall GC content of 67%, and a wobble-position GC content of 88%. Escherichia coli, on the other hand has an overall GC content of 51%, and a wobble-position GC content of 55%. The high GC content of Corynebacteria genes results in an unfavorable codon preference for heterologous expression, and can present difficulties for polymerase-based manipulations due to secondary-structure effects. Since these characteristics are due primarily to base composition at the wobble-position, synthetic genes can, in principle, be designed to eliminate these problems and retain the wild-type amino acid sequence. Such genes would obviate the need for special additives or bases during in vitro polymerase-based manipulation and mutant host strains containing uncommon tRNA's for heterologous expression. We have evaluated synthetic genes with reduced wobble-position G/C content using two variants of the enzyme 2,5-diketo-D-gluconic acid reductase (2,5-DKGR A and B) from Corynebacterium. The wild-type genes are refractory to polymerase-based manipulations and exhibit poor heterologous expression in enteric bacteria. The results indicate that a subset of codons for five amino acids (alanine, arginine, glutamate, glycine and valine) contribute the greatest contribution to reduction in G/C content at the wobble-position. Furthermore, changes in codons for two amino acids (leucine and proline) enhance bias for expression in enteric bacteria without affecting the overall G/C content. The synthetic genes are readily amplified using polymerase-based methodologies, and exhibit high levels of heterologous expression in E. coli.     

Directed evolution of thymidine kinase for AZT phosphorylation using DNA family shuffling

The thymidine kinase (TK) genes from herpes simplex virus (HSV) types 1 and 2 were recombined in vitro with a technique called DNA family shuffling. A high-throughput robotic screen identified chimeras with an enhanced ability to phosphorylate zidovudine (AZT). Improved clones were combined, reshuffled, and screened on increasingly lower concentrations of AZT. After four rounds of shuffling and screening, two clones were isolated that sensitize Escherichia coli to 32-fold less AZT compared with HSV-1 TK and 16,000-fold less than HSV-2 TK. Both clones are hybrids derived from several crossover events between the two parental genes and carry several additional amino acid substitutions not found in either parent, including active site mutations. Kinetic measurements show that the chimeric enzymes had acquired reduced K(M) for AZT as well as decreased specificity for thymidine. In agreement with the kinetic data, molecular modeling suggests that the active sites of both evolved enzymes better accommodate the azido group of AZT at the expense of thymidine. Despite the overall similarity of the two chimeric enzymes, each contains key contributions from different parents in positions influencing substrate affinity. Such mutants could be useful for anti-HIV gene therapy, and similar directed-evolution approaches could improve other enzyme-prodrug combinations.     

Improved assembly of multimeric genes for the biosynthetic production of protein polymers

We report a general method for the construction of highly repetitive synthetic genes and their use in the biosynthetic production of artificial protein polymers. Through the application of improved recombinant DNA techniques and high-throughput screening methods, we have developed a facile approach to rapid gene assembly and cloning which is widely applicable in the biosynthesis of novel protein polymers. Using this technique, synthetic genes encoding tandem repeats of the beta-sheet forming amino acid sequence AEAEAKAK were constructed and subsequently cloned into a bacterial expression host for inducible protein production. A 17-kDa fusion protein, poly-EAK9, was isolated from Escherichia coli and purified to homogeneity by immobilized metal affinity chromatography. The amino acid sequence and molecular weight were confirmed by amino acid analysis, N-terminal sequencing, and MALDI-TOF mass spectrometry. Circular dichroism studies on the artificial protein poly-EAK9 demonstrate the formation of a beta-sheet structure in aqueous solution.     

分段DNA shuffling: 一种大分子海藻糖合酶有效的定向进化方法

[目的]红色亚栖热菌(Meiothermus ruber)海藻糖合酶(Trehalose synthase,M-TreS)将麦芽糖转化生成海藻糖只需一步反应,且具有很好的热稳定性及pH耐受性,是潜在的工业生产海藻糖的酶源.为了提高该酶的性能,有必要对其进行定向进化.[方法]M-TreS基因(M-treS)大小为2 889bp.该蛋白质分子本身具有很大的进化空间,但是却不宜进行全长基因Shuffling.分段DNA shuffling是为大分子蛋白质(基因≥2 000 bp)的进化而设计的一种方法.该方法分为三步:(1)用两对引物分别扩增目的基因的上游片段和下游片段;(2)上下游片段各自进行Shuffling; (3)利用重叠延伸PCR连接上下游突变群,建立完整基因的突变文库.[结果]结合易错PCR,通过该方法经一轮进化获得一株酶活力是野生型1.6倍、催化效率是野生型2倍的突变株.序列分析表明,该突变株共有6个位点发生了氨基酸的替代,其中一个来自易错突变,2个来自同源重组,3个为随机突变.[结论]分段DNA shuffling是进化大分子蛋白质的有效方法.