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.

Engineering a tRNA and aminoacyl-tRNA synthetase for the site-specific incorporation of unnatural amino acids into proteins in vivo

In an effort to expand the scope of protein mutagenesis, we have completed the first steps toward a general method to allow the site-specific incorporation of unnatural amino acids into proteins in vivo. Our approach involves the generation of an “orthogonal” suppressor tRNA that is uniquely acylated in Escherichia coli by an engineered aminoacyl-tRNA synthetase with the desired unnatural amino acid. To this end, eight mutations were introduced into tRNA2Gln based on an analysis of the x-ray crystal structure of the glutaminyl-tRNA aminoacyl synthetase (GlnRS)–tRNA2Gln complex and on previous biochemical data. The resulting tRNA satisfies the minimal requirements for the delivery of an unnatural amino acid: it is not acylated by any endogenous E. coli aminoacyl-tRNA synthetase including GlnRS, and it functions efficiently in protein translation. Repeated rounds of DNA shuffling and oligonucleotide-directed mutagenesis followed by genetic selection resulted in mutant GlnRS enzymes that efficiently acylate the engineered tRNA with glutamine in vitro. The mutant GlnRS and engineered tRNA also constitute a functional synthetase–tRNA pair in vivo. The nature of the GlnRS mutations, which occur both at the protein–tRNA interface and at sites further away, is discussed.

Novel folded protein domains generated by combinatorial shuffling of polypeptide segments

It has been proposed that the architecture of protein domains has evolved by the combinatorial assembly and/or exchange of smaller polypeptide segments. To investigate this proposal, we fused DNA encoding the N-terminal half of a β-barrel domain (from cold shock protein CspA) with fragmented genomic Escherichia coli DNA and cloned the repertoire of chimeric polypeptides for display on filamentous bacteriophage. Phage displaying folded polypeptides were selected by proteolysis; in most cases the protease-resistant chimeric polypeptides comprised genomic segments in their natural reading frames. Although the genomic segments appeared to have no sequence homologies with CspA, one of the originating proteins had the same fold as CspA, but another had a different fold. Four of the chimeric proteins were expressed as soluble polypeptides; they formed monomers and exhibited cooperative unfolding. Indeed, one of the chimeric proteins contained a set of very slowly exchanging amides and proved more stable than CspA itself. These results indicate that native-like proteins can be generated directly by combinatorial segment assembly from nonhomologous proteins, with implications for theories of the evolution of new protein folds, as well as providing a means of creating novel domains and architectures in vitro.

Homologous genetic recombination as an intrinsic dynamic property of a DNA structure induced by RecA/Rad51-family proteins: A possible advantage of DNA over RNA as genomic material

Heteroduplex joints are general intermediates of homologous genetic recombination in DNA genomes. A heteroduplex joint is formed between a single-stranded region (or tail), derived from a cleaved parental double-stranded DNA, and homologous regions in another parental double-stranded DNA, in a reaction mediated by the RecA/Rad51-family of proteins. In this reaction, a RecA/Rad51-family protein first forms a filamentous complex with the single-stranded DNA, and then interacts with the double-stranded DNA in a search for homology. Studies of the three-dimensional structures of single-stranded DNA bound either to Escherichia coli RecA or Saccharomyces cerevisiae Rad51 have revealed a novel extended DNA structure. This structure contains a hydrophobic interaction between the 2′ methylene moiety of each deoxyribose and the aromatic ring of the following base, which allows bases to rotate horizontally through the interconversion of sugar puckers. This base rotation explains the mechanism of the homology search and base-pair switch between double-stranded and single-stranded DNA during the formation of heteroduplex joints. The pivotal role of the 2′ methylene-base interaction in the heteroduplex joint formation is supported by comparing the recombination of RNA genomes with that of DNA genomes. Some simple organisms with DNA genomes induce homologous recombination when they encounter conditions that are unfavorable for their survival. The extended DNA structure confers a dynamic property on the otherwise chemically and genetically stable double-stranded DNA, enabling gene segment rearrangements without disturbing the coding frame (i.e., protein-segment shuffling). These properties may give an extensive evolutionary advantage to DNA.

Exon shuffling and the origin of the mitochondrial targeting function in plant cytochrome c1 precursor.

Since most of the examples of "exon shuffling" are between vertebrate genes, the view is often expressed that exon shuffling is limited to the evolutionarily recent lineage of vertebrates. Although exon shuffling in plants has been inferred from the analysis of intron phases of plant genes [Long, M., Rosenberg, C. & Gilbert, W. (1995) Proc. Natl. Acad. Sci. USA 92, 12495-12499] and from the comparison of two functionally unknown sunflower genes [Domon, C. & Steinmetz, A. (1994) Mol. Gen. Genet. 244, 312-317], clear cases of exon shuffling in plant genes remain to be uncovered. Here, we report an example of exon shuffling in two important nucleus-encoded plant genes: cytosolic glyceraldehyde-3-phosphate dehydrogenase (cytosolic GAPDH or GapC) and cytochrome c1 precursor. The intron-exon structures of the shuffled region indicate that the shuffling event took place at the DNA sequence level. In this case, we can establish a donor-recipient relationship for the exon shuffling. Three amino terminal exons of GapC have been donated to cytochrome c1, where, in a new protein environment, they serve as a source of the mitochondrial targeting function. This finding throws light upon an old important but unsolved question in gene evolution: the origin of presequences or transit peptides that generally exist in nucleus-encoded organelle genes.

Infectivity of chimeric human T-cell leukemia virus type I molecular clones assessed by naked DNA inoculation.

Two human T-cell leukemia virus type I (HTLV-I) molecular clones, K30p and K34p were derived from HTLV-I-infected rabbit cell lines. K30p and K34p differ by 18 bp with changes in the long terminal repeats (LTRs) as well as in the gag, pol, and rex but not tax or env gene products. Cells transfected with clone K30p were infectious in vitro and injection of the K30p transfectants or naked K30p DNA into rabbits leads to chronic infection. In contrast, K34p did not mediate infection in vitro or in vivo, although the cell line from which it was derived is fully infectious and K34p transfectants produce intact virus particles. To localize differences involved in the ability of the clones to cause infection, six chimeric HTLV-I clones were constructed by shuffling corresponding fragments containing the substitutions in the LTRs, the gag/pol region and the rex region between K30p and K34p. Cells transfected with any of the six chimeras produced virus, but higher levels of virus were produced by cells transfected with those constructs containing the K30p rex region. Virus production was transient except in cells transfected with K30p or with a chimera consisting of the entire protein coding region of K30p flanked by K34p LTRs; only the transfectants showing persistent virus production mediated in vitro infection. In vivo infection in rabbits following intramuscular DNA injection was mediated by K30p as well as by a chimera of K30p containing the K34p rex gene. Comparisons revealed that virus production was greater and appeared earlier in rabbits injected with K30p. These data suggest that several defects in the K34p clone preclude infectivity and furthermore, provide systems to explore functions of HTLV-I genes.