Laboratory-Directed Protein Evolution

Systematic approaches to directed evolution of proteins have been documented since the 1970s. The ability to recruit new protein functions arises from the considerable substrate ambiguity of many proteins. The substrate ambiguity of a protein can be interpreted as the evolutionary potential that allows a protein to acquire new specificities through mutation or to regain function via mutations that differ from the original protein sequence. All organisms have evolutionarily exploited this substrate ambiguity. When exploited in a laboratory under controlled mutagenesis and selection, it enables a protein to “evolve” in desired directions. One of the most effective strategies in directed protein evolution is to gradually accumulate mutations, either sequentially or by recombination, while applying selective pressure. This is typically achieved by the generation of libraries of mutants followed by efficient screening of these libraries for targeted functions and subsequent repetition of the process using improved mutants from the previous screening. Here we review some of the successful strategies in creating protein diversity and the more recent progress in directed protein evolution in a wide range of scientific disciplines and its impacts in chemical, pharmaceutical, and agricultural sciences.     

Natural Selection and Drift as Individual-Level Causes of Evolution

In this paper I critically evaluate Reisman and Forber’s (Philos Sci 72(5):1113–1123, 2005) arguments that drift and natural selection are population-level causes of evolution based on what they call the manipulation condition. Although I agree that this condition is an important step for identifying causes for evolutionary change, it is insufficient. Following Woodward, I argue that the invariance of a relationship is another crucial parameter to take into consideration for causal explanations. Starting from Reisman and Forber’s example on drift and after having briefly presented the criterion of invariance, I show that once both the manipulation condition and the criterion of invariance are taken into account, drift, in this example, should better be understood as an individual-level rather than a population-level cause. Later, I concede that it is legitimate to interpret natural selection and drift as population-level causes when they rely on genuinely indeterministic events and some cases of frequency-dependent selection.     

Adaptive Evolution of p53 Thermodynamic Stability

The thermodynamic stability of a protein plays an important role during evolution and adaptation in order to maintain a folded and active conformation. p53 is a tumour suppressor involved in the regulation of numerous genes. Human p53 has an unusually low thermodynamic stability and is frequently inactivated by oncogenic missense mutations. Here, we examined the thermodynamic and kinetic stability of p53 DNA binding domains from selected invertebrate and vertebrate species by differential scanning calorimetry and equilibrium urea denaturation. There is a correlation in the apparent melting temperature of p53 with the body temperature of homeotherm vertebrates. We found that p53 from these organisms has a half-life for spontaneous unfolding at organismal body temperature of 10-20 min. We also found that p53 from invertebrates has higher stability, bearing more resemblance towards p63 and p73 from humans. Using structure-guided mutagenesis on the human p53 scaffold, we demonstrated that the amino acid changes on the protein surface and in the protein interior lead to the elevated stability of p53 orthologs. We propose a model in which the p53 DNA binding domain has been shaped by the complex interplay of different selective pressures and underwent adaptive evolution leading to pronounced effects on its stability. p53 from vertebrates has evolved to have a low thermodynamic stability and similarly short spontaneous half-life at organismal body temperature, which is related to function.     

Thermodynamic Framework of the Interaction between Protein and Solvent Drives Protein Folding

Abstract

We use internal coordinate molecular mechanics calculations to study the impact of abasic sites on the conformation and the mechanics of the DNA double helix. Abasic sites, which are common mutagenic lesions, are shown to locally modify both the groove geometry and the curvature of DNA in a sequence dependent manner. By controlled twisting and bending, it is also shown that these lesions modify the deformability of the duplex, generally increasing its flexibility, but again to an extent which depends on the nature of the abasic site and on the surrounding base sequence. Both the conformational and mechanical influence of this type of DNA damage may be significant for recognition and repair mechanisms.     

Genetic Constraints on Protein Evolution

ABSTRACT

Evolution requires the generation and optimization of new traits (“adaptation”) and involves the selection of mutations that improve cellular function. These mutations were assumed to arise by selection of neutral mutations present at all times in the population. Here we review recent evidence that indicates that deleterious mutations are more frequent in the population than previously recognized and that these mutations play a significant role in protein evolution through continuous positive selection. Positively selected mutations include adaptive mutations, i.e. mutations that directly affect enzymatic function, and compensatory mutations, which suppress the pleiotropic effects of adaptive mutations. Compensatory mutations are by far the most frequent of the two and would allow potentially adaptive but deleterious mutations to persist long enough in the population to be positively selected during episodes of adaptation. Compensatory mutations are, by definition, context-dependent and thus constrain the paths available for evolution. This provides a mechanistic basis for the examples of highly constrained evolutionary landscapes and parallel evolution reported in natural and experimental populations. The present review article describes these recent advances in the field of protein evolution and discusses their implications for understanding the genetic basis of disease and for protein engineering in vitro.     

Evolution of Protein Binding Modes in Homooligomers

The evolution of protein interactions cannot be deciphered without a detailed analysis of interaction interfaces and binding modes. We performed a large-scale study of protein homooligomers in terms of their symmetry, interface sizes, and conservation of binding modes. We also focused specifically on the evolution of protein binding modes from nine families of homooligomers and mapped 60 different binding modes and oligomerization states onto the phylogenetic trees of these families. We observed a significant tendency for the same binding modes to be clustered together and conserved within clades on phylogenetic trees; this trend is especially pronounced for close homologs with 70% sequence identity or higher. Some binding modes are conserved among very distant homologs, pointing to their ancient evolutionary origin, while others are very specific for a certain phylogenetic group. Moreover, we found that the most ancient binding modes have a tendency to involve symmetrical (isologous) homodimer binding arrangements with larger interfaces, while recently evolved binding modes more often exhibit asymmetrical arrangements and smaller interfaces.     

Episodic Evolution of Protein Hormones in Mammals

Pituitary growth hormone (GH) and prolactin have been shown previously to display a pattern of evolution in which episodes of rapid change are imposed on a low underlying basal rate (near-stasis). This study was designed to explore whether a similar pattern is seen in the evolution of other protein hormones in mammals. Seven protein hormones were examined (with the common α-subunit of the glycoprotein hormones providing an additional polypeptide for analysis)—those for which sequences from at least four eutherian orders are available with a suitable non-eutherian outgroup. Six of these (GH, prolactin, insulin, parathyroid hormone, glycoprotein hormone α-subunit, and luteinizing hormone β-subunit) showed markedly variable evolutionary rates in each case with a pattern of a slow basal rate and bursts of rapid change, the precise positions of the bursts varying from protein to protein. Two protein hormones (follicle-stimulating hormone β-subunit and thyroid-stimulating hormone β-subunit) showed no significant rate variation. Based on the sequences currently available, and pooling data from all eight proteins, the phase of slow basal change occupied about 85% of the sampled evolutionary time, but most evolutionary change (about 62% of the substitutions accepted) occurred during the episodes of rapid change. It is concluded that, in mammals at least, a pattern of prolonged periods of near-stasis with occasional episodes of rapid change provides a better model of evolutionary change for protein hormones than the one of constant evolutionary rates that is commonly favored. The mechanisms underlying this episodic evolution are not yet clear, and it may be that they vary from one group to another; in some cases, positive selection appears to underlie bursts of rapid change. Where gene duplication is associated with a period of accelerated evolution this often occurs at the end rather than the beginning of the episode. To what extent the type of pattern seen for protein hormones can be extended to other proteins remains to be established. Received: 10 October 2000 / Accepted: 18 December 2000     

Thermodynamic Characterization of the Unfolding of the Prion Protein

The prion protein appears to be unusually susceptible to conformational change, and unlike nearly all other proteins, it can easily be made to convert to alternative misfolded conformations. To understand the basis of this structural plasticity, a detailed thermodynamic characterization of two variants of the mouse prion protein (moPrP), the full-length moPrP (23–231) and the structured C-terminal domain, moPrP (121–231), has been carried out. All thermodynamic parameters governing unfolding, including the changes in enthalpy, entropy, free energy, and heat capacity, were found to be identical for the two protein variants. The N-terminal domain remains unstructured and does not interact with the C-terminal domain in the full-length protein at pH 4. Moreover, the enthalpy and entropy of unfolding of moPrP (121–231) are similar in magnitude to values reported for other proteins of similar size. However, the protein has an unusually high native-state heat capacity, and consequently, the change in heat capacity upon unfolding is much lower than that expected for a protein of similar size. It appears, therefore, that the native state of the prion protein undergoes substantial fluctuations in enthalpy and hence, in structure.     

System Relaxation and Thermodynamic Integration

Abstract

The thermodynamic integration (TI) method for calculating free energy differences has three inherent problems: statistics, numerical integration, and relaxation. In this paper the latter is analyzed for the nonequilibrium TI method introduced by Postma which combines molecular dynamics simulation and TI in a very effective way. A nontrivial extrapolation technique is presented to remove the relaxation error and to calculate the underlying relaxation time. It is shown that the optimal choice of grid points, derived in a previous paper, for minimizing statistical errors not only removes integration errors, but also minimizes relaxation errors. The methods are applied in a calculation of the free energy of cavity formation in water.     

Protein Design and Evolution for Biocatalysis

Cover illustration: Soon designing proteins on demand? This is the vision expressed by the designer pencils pointing at the logo of the meeting held in September 2006 in Greifswald, Germany. This special issue was edited by Prof. Uwe Bornscheuer from Greifswald, who selected papers from keynote speakers at this meeting. With special thanks to Prof. Romas Kazlauskas, University of Minnesota, for the design of the conference logo. Pencils © FOTOLIA.     

蛋白质定向进化技术的研究进展

蛋白质定向进化技术是蛋白质分子改造的一个重要策略.重点介绍了易错PCR、DNA改组等对编码蛋白质的基因进行随机突变和重组的技术,以及构建突变体库和高通量筛选的方法,并探讨了定向进化技术在蛋白质工程中的应用及前景.     

Protein Design and Evolution for Biocatalysis

Cover illustration: Protein Design and Evolution for Biocatalysis. This special issue of Biotechnology Journal contains selected contributions from scientists participating to the ESF-EMBO Symposium which took place in October 2008 in San Feliu, Spain. Guest Editor is the chair and organizer of the meeting, Jiri Damborsky from Brno (Czech Republic). He highlights a variety of topics brought up in the meeting, ranging from new methods of rational design, directed evolution, metagenomics and single-molecule techniques, to construction of useful enzymes for industrial applications. Uwe Bornscheuer (Greifswald, Germany) authored a meeting report. Image colored pencils, © PhotoDisc, Inc.; Protein logo © ESF.     

Molecular Evolution of TEPP Protein Genes in Metazoans

TEPP is a gene expressed in human reproductive organs such as testis, prostate, and placenta. Here, identification and molecular evolutionary analysis of TEPP proteins in various metazoan animals including deuterostomes (chordates, hemichordates, and echinoderms), lophotrochozoans (mollusks and annelids), and cnidarians (sea anemone and coral) are reported. A multiple sequence alignment revealed two highly conserved regions in TEPP proteins that had no similarity to any other known domains or proteins. Genomic sequence analysis showed frequent shifting of the splice sites of intron 1 in mammalian TEPP genes. A comparison of the intron positions in the coding region showed that the exon/intron structure of the TEPP gene was established in an early metazoan ancestor and that independent loss of a single intron occurred in echinoderms and in vertebrates. The urochordate tunicate TEPP genes are intronless, possibly due to replacement of the original gene by a retrogene. No homolog was detected in birds, insects, nematodes, and teleost fishes despite the extensive sequence data of these species, implying that the TEPP gene might be lost in these lineages.     

Preservation of Protein Dynamics in Dihydrofolate Reductase Evolution

The hydride transfer reaction catalyzed by dihydrofolate reductase (DHFR) is a model for examining how protein dynamics contribute to enzymatic function. The relationship between functional motions and enzyme evolution has attracted significant attention. Recent studies on N23PP Escherichia coli DHFR (ecDHFR) mutant, designed to resemble parts of the human enzyme, indicated a reduced single turnover rate. NMR relaxation dispersion experiments with that enzyme showed rigidification of millisecond Met-20 loop motions (Bhabha, G., Lee, J., Ekiert, D. C., Gam, J., Wilson, I. A., Dyson, H. J., Benkovic, S. J., and Wright, P. E. (2011) Science 332, 234–238). A more recent study of this mutant, however, indicated that fast motions along the reaction coordinate are actually more dispersed than for wild-type ecDHFR (WT). Furthermore, a double mutant (N23PP/G51PEKN) that better mimics the human enzyme seems to restore both the single turnover rates and narrow distribution of fast dynamics (Liu, C. T., Hanoian, P., French, T. H., Hammes-Schiffer, S., and Benkovic, S. J. (2013) Proc. Natl. Acad. Sci. U.S.A. 110, 10159–11064). Here, we measured intrinsic kinetic isotope effects for both N23PP and N23PP/G51PEKN double mutant DHFRs over a temperature range. The findings indicate that although the C-H→C transfer and dynamics along the reaction coordinate are impaired in the altered N23PP mutant, both seem to be restored in the N23PP/G51PEKN double mutant. This indicates that the evolution of G51PEKN, although remote from the Met-20 loop, alleviated the loop rigidification that would have been caused by N23PP, enabling WT-like H-tunneling. The correlation between the calculated dynamics, the nature of C-H→C transfer, and a phylogenetic analysis of DHFR sequences are consistent with evolutionary preservation of the protein dynamics to enable H-tunneling from well reorganized active sites.     

Protein Evolution in the Context of Drosophila Development

The tempo at which a protein evolves depends not only on the rate at which mutations arise but also on the selective effects that those mutations have at the organismal level. It is intuitive that proteins functioning during different stages of development may be predisposed to having mutations of different selective effects. For example, it has been hypothesized that changes to proteins expressed during early development should have larger phenotypic consequences because later stages depend on them. Conversely, changes to proteins expressed much later in development should have smaller consequences at the organismal level. Here we assess whether proteins expressed at different times during Drosophila development vary systematically in their rates of evolution. We find that proteins expressed early in development and particularly during mid–late embryonic development evolve unusually slowly. In addition, proteins expressed in adult males show an elevated evolutionary rate. These two trends are independent of each other and cannot be explained by peculiar rates of mutation or levels of codon bias. Moreover, the observed patterns appear to hold across several functional classes of genes, although the exact developmental time of the slowest protein evolution differs among each class. We discuss our results in connection with data on the evolution of development.This revised article was published online in June 2005 with the final version of the content.Jerel C. Davis and Onn Brandman contributed equally to this work. Reviewing Editor: ManYuan Long