Intense neutral drifts yield robust and evolvable consensus proteins

Toward a systems biology perspective on enzyme evolution

Large superfamilies of enzymes derived from a common progenitor have emerged by duplication and divergence of genes encoding metabolic enzymes. Division of the functions of early generalist enzymes enhanced catalytic power and control over metabolic fluxes. Later, novel enzymes evolved from inefficient secondary activities in specialized enzymes. Enzymes operate in the context of complex metabolic and regulatory networks. The potential for evolution of a new enzyme depends upon the collection of enzymes in a microbe, the topology of the metabolic network, the environmental conditions, and the net effect of trade-offs between the original and novel activities of the enzyme.

The transition of human estrogen sulfotransferase from generalist to specialist using directed enzyme evolution

Broad specificity is believed to be a property of primordial enzymes that diverged during natural protein evolution to produce highly specific and efficient enzymes. Human estrogen sulfotransferase (SULT1E1) is a broad-specificity enzyme that detoxifies a variety of chemicals, including estrogens, by the transfer of sulfate. To study the molecular basis for the broad specificity of this enzyme and to investigate the process of SULT1E1 specialization, we have adopted a directed enzyme evolution approach. Using two iterative rounds of evolution, we generated SULT1E1 mutants with increased thermostability and narrower specificity from the broadly specific wild-type enzyme. To identify mutants with enhanced specificity, we developed an unbiased screening assay to assess sulfate transfer to three different acceptors in parallel. Such an assay enabled the isolation of SULT1E1 mutants with enhanced or wild-type activity toward an estrogen acceptor and significantly reduced activity for phenol or coumarin type of acceptors, leading to up to 3 orders of magnitude increase in specificity. We found that mutations conferring novel specificity are located in the vicinity of the active site and thus may play a direct role in reshaping the acceptor-binding site. Finally, such mutations resulted in reduced SULT1E1 thermostability, revealing a trade-off between SULT1E1 thermostability and acquisition of novel function.

The molecular basis for the broad substrate specificity of human sulfotransferase 1A1

Cytosolic sulfotransferases (SULTs) are mammalian enzymes that detoxify a wide variety of chemicals through the addition of a sulfate group. Despite extensive research, the molecular basis for the broad specificity of SULTs is still not understood. Here, structural, protein engineering and kinetic approaches were employed to obtain deep understanding of the molecular basis for the broad specificity, catalytic activity and substrate inhibition of SULT1A1. We have determined five new structures of SULT1A1 in complex with different acceptors, and utilized a directed evolution approach to generate SULT1A1 mutants with enhanced thermostability and increased catalytic activity. We found that active site plasticity enables binding of different acceptors and identified dramatic structural changes in the SULT1A1 active site leading to the binding of a second acceptor molecule in a conserved yet non-productive manner. Our combined approach highlights the dominant role of SULT1A1 structural flexibility in controlling the specificity and activity of this enzyme.

Network models of TEM β-lactamase mutations coevolving under antibiotic selection show modular structure and anticipate evolutionary trajectories

Understanding how novel functions evolve (genetic adaptation) is a critical goal of evolutionary biology. Among asexual organisms, genetic adaptation involves multiple mutations that frequently interact in a non-linear fashion (epistasis). Non-linear interactions pose a formidable challenge for the computational prediction of mutation effects. Here we use the recent evolution of β-lactamase under antibiotic selection as a model for genetic adaptation. We build a network of coevolving residues (possible functional interactions), in which nodes are mutant residue positions and links represent two positions found mutated together in the same sequence. Most often these pairs occur in the setting of more complex mutants. Focusing on extended-spectrum resistant sequences, we use network-theoretical tools to identify triple mutant trajectories of likely special significance for adaptation. We extrapolate evolutionary paths (n = 3) that increase resistance and that are longer than the units used to build the network (n = 2). These paths consist of a limited number of residue positions and are enriched for known triple mutant combinations that increase cefotaxime resistance. We find that the pairs of residues used to build the network frequently decrease resistance compared to their corresponding singlets. This is a surprising result, given that their coevolution suggests a selective advantage. Thus, β-lactamase adaptation is highly epistatic. Our method can identify triplets that increase resistance despite the underlying rugged fitness landscape and has the unique ability to make predictions by placing each mutant residue position in its functional context. Our approach requires only sequence information, sufficient genetic diversity, and discrete selective pressures. Thus, it can be used to analyze recent evolutionary events, where coevolution analysis methods that use phylogeny or statistical coupling are not possible. Improving our ability to assess evolutionary trajectories will help predict the evolution of clinically relevant genes and aid in protein design.

Understanding how new biological activities evolve on the molecular level has critical implications for biotechnology and for human health. Here we collect a database of mutations that contribute to the evolution of β-lactamase resistance to inhibitors and to new β-lactam antibiotics in bacterial pathogens, such as Escherichia coli. We compiled a database of TEM β-lactamase sequences evolved under antibiotic pressure and identified functional interactions between individual residue positions. We visualized these complex molecular interactions as a network and used network theory to derive information regarding the origin of individual mutations and their contribution to the observed resistance. Our approach should help interpret sequence databases for clinically relevant proteins undergoing high mutation rates and under selective (drug, immune) pressure, such as surface proteins of pathogens (particularly of RNA viruses such as HIV) or targets for chemotherapy in microbial pathogen or tumor cells. Notably, our approach only requires sequence data; detailed phylogenetic or tertiary structure information for the target gene is not necessary. Our analysis of how individual mutations work together to produce new biological activities should help anticipate evolution driven by a variety of clinically-relevant selections such as drug resistance, virulence, and immunity.

Recombining low homology, functionally rich regions of bacterial subtilisins by combinatorial fragment exchange

Combinatorial fragment exchange was utilised to recombine key structural and functional low homology regions of bacilli subtilisins to generate new active hybrid proteases with altered substrate profiles. Up to six different regions comprising mostly of loop residues from the commercially important subtilisin Savinase were exchanged with the structurally equivalent regions of six other subtilisins. The six additional subtilisins derive from diverse origins and included thermophilic and intracellular subtilisins as well as other academically and commercially relevant subtilisins. Savinase was largely tolerant to fragment exchange; rational replacement of all six regions with 5 of 6 donating subtilisin sequences preserved activity, albeit reduced compared to Savinase. A combinatorial approach was used to generate hybrid Savinase variants in which the sequences derived from all seven subtilisins at each region were recombined to generate new region combinations. Variants with different substrate profiles and with greater apparent activity compared to Savinase and the rational fragment exchange variants were generated with the substrate profile exhibited by variants dependent on the sequence combination at each region.

Slow protein evolutionary rates are dictated by surface-core association

Why do certain proteins evolve much slower than others? We compared not only rates per protein, but also rates per position within individual proteins. For ∼90% of proteins, the distribution of positional rates exhibits three peaks: a peak of slow evolving residues, with average log(2)[normalized rate], log(2)μ, of ca. -2, corresponding primarily to core residues; a peak of fast evolving residues (log(2)μ ∼ 0.5) largely corresponding to surface residues; and a very fast peak (log(2)μ ∼ 2) associated with disordered segments. However, a unique fraction of proteins that evolve very slowly exhibit not only a negligible fast peak, but also a peak with a log(2)μ ∼ -4, rather than the standard core peak of -2. Thus, a "freeze" of a protein's surface seems to stop core evolution as well. We also observed a much higher fraction of substitutions in potentially interacting residues than expected by chance, including substitutions in pairs of contacting surface-core residues. Overall, the data suggest that accumulation of surface substitutions enables the acceptance of substitutions in core positions. The underlying reason for slow evolution might therefore be a highly constrained surface due to protein-protein interactions or the need to prevent misfolding or aggregation. If the surface is inaccessible to substitutions, so becomes the core, thus resulting in very slow overall rates.

Mechanisms for the evolution of a derived function in the ancestral glucocorticoid receptor

Understanding the genetic, structural, and biophysical mechanisms that caused protein functions to evolve is a central goal of molecular evolutionary studies. Ancestral sequence reconstruction (ASR) offers an experimental approach to these questions. Here we use ASR to shed light on the earliest functions and evolution of the glucocorticoid receptor (GR), a steroid-activated transcription factor that plays a key role in the regulation of vertebrate physiology. Prior work showed that GR and its paralog, the mineralocorticoid receptor (MR), duplicated from a common ancestor roughly 450 million years ago; the ancestral functions were largely conserved in the MR lineage, but the functions of GRs-reduced sensitivity to all hormones and increased selectivity for glucocorticoids-are derived. Although the mechanisms for the evolution of glucocorticoid specificity have been identified, how reduced sensitivity evolved has not yet been studied. Here we report on the reconstruction of the deepest ancestor in the GR lineage (AncGR1) and demonstrate that GR's reduced sensitivity evolved before the acquisition of restricted hormone specificity, shortly after the GR-MR split. Using site-directed mutagenesis, X-ray crystallography, and computational analyses of protein stability to recapitulate and determine the effects of historical mutations, we show that AncGR1's reduced ligand sensitivity evolved primarily due to three key substitutions. Two large-effect mutations weakened hydrogen bonds and van der Waals interactions within the ancestral protein, reducing its stability. The degenerative effect of these two mutations is extremely strong, but a third permissive substitution, which has no apparent effect on function in the ancestral background and is likely to have occurred first, buffered the effects of the destabilizing mutations. Taken together, our results highlight the potentially creative role of substitutions that partially degrade protein structure and function and reinforce the importance of permissive mutations in protein evolution.

Constructing de novo biosynthetic pathways for chemical synthesis inside living cells

Living organisms have evolved a vast array of catalytic functions that make them ideally suited for the production of medicinally and industrially relevant small-molecule targets. Indeed, native metabolic pathways in microbial hosts have long been exploited and optimized for the scalable production of both fine and commodity chemicals. Our increasing capacity for DNA sequencing and synthesis has revealed the molecular basis for the biosynthesis of a variety of complex and useful metabolites and allows the de novo construction of novel metabolic pathways for the production of new and exotic molecular targets in genetically tractable microbes. However, the development of commercially viable processes for these engineered pathways is currently limited by our ability to quickly identify or engineer enzymes with the correct reaction and substrate selectivity as well as the speed by which metabolic bottlenecks can be determined and corrected. Efforts to understand the relationship among sequence, structure, and function in the basic biochemical sciences can advance these goals for synthetic biology applications while also serving as an experimental platform for elucidating the in vivo specificity and function of enzymes and reconstituting complex biochemical traits for study in a living model organism. Furthermore, the continuing discovery of natural mechanisms for the regulation of metabolic pathways has revealed new principles for the design of high-flux pathways with minimized metabolic burden and has inspired the development of new tools and approaches to engineering synthetic pathways in microbial hosts for chemical production.

Applying neutral drift to the directed molecular evolution of a β-glucuronidase into a β-galactosidase: Two different evolutionary pathways lead to the same variant

Background

Directed protein evolution has been used to modify protein activity and research has been carried out to enhance the production of high quality mutant libraries. Many theoretical approaches suggest that allowing a population to undergo neutral selection may be valuable in directed evolution experiments.

Findings

Here we report on an investigation into the value of neutral selection in a classical model system for directed evolution, the conversion of the E. coli β-glucuronidase to a β-galactosidase activity. We find that neutral selection, i.e. selection for retaining glucuronidase activity, can efficiently identify the majority of sites of mutation that have been identified as beneficial for galactosidase activity in previous experiments. Each variant demonstrating increased galactosidase activity identified by our neutral drift experiments contained a mutation at one of four sites, T509, S557, N566 or W529. All of these sites have previously been identified using direct selection for beta galactosidase activity.

Conclusions

Our results are consistent with others that show that a neutral selection approach can be effective in selecting improved variants. However, we interpret our results to show that neutral selection is, in this case, not a more efficient approach than conventional directed evolution approaches. However, the neutral approach is likely to be beneficial when the resulting library can be screened for a range of related activities. More detailed statistical studies to resolve the apparent differences between this system and others are likely to be a fruitful avenue for future research.

Protein design in metabolic engineering and synthetic biology

Starting from experimental data on sequence, structure or biochemical properties of enzymes, protein design seeks to construct enzymes with desired activity, stability, specificity and selectivity. Two strategies are widely used to investigate sequence-structure-function relationships: statistical methods to analyse protein families or mutant libraries, and molecular modelling methods to study proteins and their interaction with ligands or substrates. On the basis of these methods, protein design has been successfully applied to fine-tune bottleneck enzymes in metabolic engineering and to design enzymes with new substrate spectra and new functions. However, constructing efficient metabolic pathways by integrating individual enzymes into a complex system is challenging. The field of synthetic biology is still in its infancy, but promising results have demonstrated the feasibility and usefulness of the concept.

Utilizing natural diversity to evolve protein function: applications towards thermostability

Protein evolution relies on designing a library of sequences that capture meaningful functional diversity in a limited number of protein variants. Several approaches take advantage of the sequence space already explored through natural selection by incorporating sequence diversity available from modern genomes (and their ancestors) when designing these libraries. The success of these approaches is, partly, owing to the fact that modern sequence diversity has already been subjected to evolutionary selective forces and thus the diversity has already been deemed 'fit to survive'. Five of these approaches will be discussed in this review to highlight how protein engineers can use evolutionary sequence history/diversity of homologous proteins in unique ways to design protein libraries.

The contribution of coevolving residues to the stability of KDO8P synthase

Background

The evolutionary tree of 3-deoxy-D-manno-octulosonate 8-phosphate (KDO8P) synthase (KDO8PS), a bacterial enzyme that catalyzes a key step in the biosynthesis of bacterial endotoxin, is evenly divided between metal and non-metal forms, both having similar structures, but diverging in various degrees in amino acid sequence. Mutagenesis, crystallographic and computational studies have established that only a few residues determine whether or not KDO8PS requires a metal for function. The remaining divergence in the amino acid sequence of KDO8PSs is apparently unrelated to the underlying catalytic mechanism.

Methodology/Principal Findings

The multiple alignment of all known KDO8PS sequences reveals that several residue pairs coevolved, an indication of their possible linkage to a structural constraint. In this study we investigated by computational means the contribution of coevolving residues to the stability of KDO8PS. We found that about 1/4 of all strongly coevolving pairs probably originated from cycles of mutation (decreasing stability) and suppression (restoring it), while the remaining pairs are best explained by a succession of neutral or nearly neutral covarions.

Conclusions/Significance

Both sequence conservation and coevolution are involved in the preservation of the core structure of KDO8PS, but the contribution of coevolving residues is, in proportion, smaller. This is because small stability gains or losses associated with selection of certain residues in some regions of the stability landscape of KDO8PS are easily offset by a large number of possible changes in other regions. While this effect increases the tolerance of KDO8PS to deleterious mutations, it also decreases the probability that specific pairs of residues could have a strong contribution to the thermodynamic stability of the protein.

Initial mutations direct alternative pathways of protein evolution

Whether evolution is erratic due to random historical details, or is repeatedly directed along similar paths by certain constraints, remains unclear. Epistasis (i.e. non-additive interaction between mutations that affect fitness) is a mechanism that can contribute to both scenarios. Epistasis can constrain the type and order of selected mutations, but it can also make adaptive trajectories contingent upon the first random substitution. This effect is particularly strong under sign epistasis, when the sign of the fitness effects of a mutation depends on its genetic background. In the current study, we examine how epistatic interactions between mutations determine alternative evolutionary pathways, using in vitro evolution of the antibiotic resistance enzyme TEM-1 β-lactamase. First, we describe the diversity of adaptive pathways among replicate lines during evolution for resistance to a novel antibiotic (cefotaxime). Consistent with the prediction of epistatic constraints, most lines increased resistance by acquiring three mutations in a fixed order. However, a few lines deviated from this pattern. Next, to test whether negative interactions between alternative initial substitutions drive this divergence, alleles containing initial substitutions from the deviating lines were evolved under identical conditions. Indeed, these alternative initial substitutions consistently led to lower adaptive peaks, involving more and other substitutions than those observed in the common pathway. We found that a combination of decreased enzymatic activity and lower folding cooperativity underlies negative sign epistasis in the clash between key mutations in the common and deviating lines (Gly238Ser and Arg164Ser, respectively). Our results demonstrate that epistasis contributes to contingency in protein evolution by amplifying the selective consequences of random mutations.

A long-term goal of evolutionary biology is to understand the factors that govern the outcome of evolution. Epistasis (i.e. the situation in which the fitness effect of a mutation depends on its genetic background) is one such factor. Epistasis not only affects the dynamics of evolution, it may also direct its outcome by affecting the type and order of selected mutations. This effect is particularly strong under sign epistasis, which occurs when the sign of a mutation's fitness effect depends on its genetic background. Here, we demonstrate how epistasis causes divergence of mutational pathways of an antibiotic resistance enzyme, TEM-1 β-lactamase. First, we use in vitro mutagenesis followed by selection for cefotaxime resistance to demonstrate that alternative mutational pathways towards highly resistant variants exist in addition to the main pathway that was previously described. Next, to test whether negative interactions between alternative initial substitutions govern this diversification, we start identical evolution experiments with alleles containing initial substitutions from the deviating lines. These alleles consistently evolve to lower adaptive peaks and acquire different mutations than those in the main pathway. Our results demonstrate that sign epistasis between alternative initial substitutions may force evolution to follow different mutational pathways.

Natural evolution of TEM-1 β-lactamase: experimental reconstruction and clinical relevance

TEM-1 β-lactamase is one of the most well-known antibiotic resistance determinants around. It confers resistance to penicillins and early cephalosporins and has shown an astonishing functional plasticity in response to the introduction of novel drugs derived from these antibiotics. Since its discovery in the 1960s, over 170 variants of TEM-1 - with different amino acid sequences and often resistance phenotypes - have been isolated in hospitals and clinics worldwide. Next to this well-documented 'natural' evolution, the in vitro evolution of TEM-1 has been the focus of attention of many experimental studies. In this review, we compare the natural and laboratory evolution of TEM-1 in order to address the question to what extent the evolution of antibiotic resistance can be repeated, and hence might have been predicted, under laboratory conditions. We also use the comparison to gain an insight into the adaptive relevance of hitherto uncharacterized substitutions present in clinical isolates and to predict substitutions not yet observed in nature. Based on new structural insights, we review what is known about substitutions in TEM-1 that contribute to the extension of its resistance phenotype. Finally, we address the clinical relevance of TEM alleles during the past decade, which has been dominated by the emergence of another β-lactamase, CTX-M.

Switching from an esterase to a hydroxynitrile lyase mechanism requires only two amino acid substitutions

The alpha/beta hydrolase superfamily contains mainly esterases, which catalyze hydrolysis, but also includes hydroxynitrile lyases, which catalyze addition of cyanide to aldehydes, a carbon-carbon bond formation. Here, we convert a plant esterase, SABP2, into a hydroxynitrile lyase using just two amino acid substitutions. Variant SABP2-G12T-M239K lost the ability to catalyze ester hydrolysis (<0.9 mU/mg) and gained the ability to catalyze the release of cyanide from mandelonitrile (20 mU/mg, k(cat)/K(M) = 70 min(-1)M(-1)). This variant also catalyzed the reverse reaction, formation of mandelonitrile with low enantioselectivity: 20% ee (S), E = 1.5. The specificity constant for the lysis of mandelontrile is 13,000-fold faster than the uncatalyzed reaction and only 1300-fold less efficient (k(cat/)K(M)) than hydroxynitrile lyase from rubber tree.

Status of protein engineering for biocatalysts: how to design an industrially useful biocatalyst

Recent advances in the development of both experimental and computational protein engineering tools have enabled a number of further successes in the development of biocatalysts ready for large-scale applications. Key tools are first, the targeting of libraries, leading to far smaller but more useful libraries than in the past, second, the combination of structural, mechanistic, and sequence-based knowledge often based on prior successful cases, and third, the advent of structurally based algorithms allowing the design of novel functions. Based on these tools, a number of improved biocatalysts for pharmaceutical applications have been presented, such as an (R)-transaminase for the synthesis of active pharmaceutical ingredients (APIs) of sitagliptin (Januvia®) and ketoreductases, glucose dehydrogenases, and haloalkane dehalogenases for the API synthesis toward atorvastatin (Lipitor®) and montelukast (Singulair®).

Experimental evolution of adenylate kinase reveals contrasting strategies toward protein thermostability

Success in evolution depends critically upon the ability of organisms to adapt, a property that is also true for the proteins that contribute to the fitness of an organism. Successful protein evolution is enhanced by mutational pathways that generate a wide range of physicochemical mechanisms to adaptation. In an earlier study, we used a weak-link method to favor changes to an essential but maladapted protein, adenylate kinase (AK), within a microbial population. Six AK mutants (a single mutant followed by five double mutants) had success within the population, revealing a diverse range of adaptive strategies that included changes in nonpolar packing, protein folding dynamics, and formation of new hydrogen bonds and electrostatic networks. The first mutation, AK(BSUB) Q199R, was essential in defining the structural context that facilitated subsequent mutations as revealed by a considerable mutational epistasis and, in one case, a very strong dependence upon the order of mutations. Namely, whereas the single mutation AK(BSUB) G213E decreases protein stability by >25 degrees C, the same mutation in the background of AK(BSUB) Q199R increases stability by 3.4 degrees C, demonstrating that the order of mutations can play a critical role in favoring particular molecular pathways to adaptation. In turn, protein folding kinetics shows that four of the five AK(BSUB) double mutants utilize a strategy in which an increase in the folding rate accompanied by a decrease in the unfolding rate results in additional stability. However, one mutant exhibited a dramatic increase in the folding relative to a modest increase in the unfolding rate, suggesting a different adaptive strategy for thermostability. In all cases, an increase in the folding rates for the double mutants appears to be the preferred mechanism in conferring additional stability and may be an important aspect of protein evolution. The range of overlapping as well as contrasting strategies for success illustrates both the power and subtlety of adaptation at even the smallest unit of change, a single amino acid.

Pervasive cryptic epistasis in molecular evolution

The functional effects of most amino acid replacements accumulated during molecular evolution are unknown, because most are not observed naturally and the possible combinations are too numerous. We created 168 single mutations in wild-type Escherichia coli isopropymalate dehydrogenase (IMDH) that match the differences found in wild-type Pseudomonas aeruginosa IMDH. 104 mutant enzymes performed similarly to E. coli wild-type IMDH, one was functionally enhanced, and 63 were functionally compromised. The transition from E. coli IMDH, or an ancestral form, to the functional wild-type P. aeruginosa IMDH requires extensive epistasis to ameliorate the combined effects of the deleterious mutations. This result stands in marked contrast with a basic assumption of molecular phylogenetics, that sites in sequences evolve independently of each other. Residues that affect function are scattered haphazardly throughout the IMDH structure. We screened for compensatory mutations at three sites, all of which lie near the active site and all of which are among the least active mutants. No compensatory mutations were found at two sites indicating that a single site may engage in compound epistatic interactions. One complete and three partial compensatory mutations of the third site are remote and lie in a different domain. This demonstrates that epistatic interactions can occur between distant (>20Å) sites. Phylogenetic analysis shows that incompatible mutations were fixed in different lineages.

Evolutionary fates within a microbial population highlight an essential role for protein folding during natural selection

Systems biology can offer a great deal of insight into evolution by quantitatively linking complex properties such as protein structure, folding, and function to the fitness of an organism. Although the link between diseases such as Alzheimer's and misfolding is well appreciated, directly showing the importance of protein folding to success in evolution has been more difficult. We show here that predicting success during adaptation can depend critically on enzyme kinetic and folding models. We used a 'weak link' method to favor mutations to an essential, but maladapted, adenylate kinase gene within a microbial population that resulted in the identification of five mutants that arose nearly simultaneously and competed for success. Physicochemical characterization of these mutants showed that, although steady-state enzyme activity is important, success within the population is critically dependent on resistance to denaturation and aggregation. A fitness function based on in vitro measurements of enzyme activity, reversible and irreversible unfolding, and the physiological context reproduces in vivo evolutionary fates in the population linking organismal adaptation to its physical basis.

Multiple global suppressors of protein stability defects facilitate the evolution of extended-spectrum TEM β-lactamases

The introduction of extended-spectrum cephalosporins and β-lactamase inhibitors has driven the evolution of extended-spectrum β-lactamases (ESBLs) that possess the ability to hydrolyze these drugs. The evolved TEM ESBLs from clinical isolates of bacteria often contain substitutions that occur in the active site and alter the catalytic properties of the enzyme to provide an increased hydrolysis of extended-spectrum cephalosporins or an increased resistance to inhibitors. These active-site substitutions often result in a cost in the form of reduced enzyme stability. The evolution of TEM ESBLs is facilitated by mutations that act as global suppressors of protein stability defects in that they allow the enzyme to absorb multiple amino acid changes despite incremental losses in stability associated with the substitutions. The best-studied example is the M182T substitution, which corrects protein stability defects and is commonly found in TEM ESBLs or inhibitor-resistant variants from clinical isolates. In this study, a genetic selection for second-site mutations that could partially restore function to a severely destabilized primary mutant enabled the identification of A184V, T265M, R275Q, and N276D, which are known to occur in TEM ESBLs from clinical isolates, as suppressors of TEM-1 protein stability defects. Further characterization demonstrated that these substitutions increased the thermal stability of TEM-1 and were able to correct the stability defects of two different sets of destabilizing mutations. The acquisition of compensatory global suppressors of stability costs associated with active-site mutations may be a common mechanism for the evolution of novel protein function.

Robustness and evolvability

Why isn't random variation always deleterious? Are there factors that sometimes make adaptation easier? Biological systems are extraordinarily robust to perturbation by mutations, recombination and the environment. It has been proposed that this robustness might make them more evolvable. Robustness to mutation allows genetic variation to accumulate in a cryptic state. Switching mechanisms known as evolutionary capacitors mean that the amount of heritable phenotypic variation available can be correlated to the degree of stress and hence to the novelty of the environment and remaining potential for adaptation. There have been two somewhat separate literatures relating robustness to evolvability. One has focused on molecular phenotypes and new mutations, the other on morphology and cryptic genetic variation. Here, we review both literatures, and show that the true distinction is whether recombination rates are high or low. In both cases, the evidence supports the claim that robustness promotes evolvability.

Mutational effects and the evolution of new protein functions

The divergence of new genes and proteins occurs through mutations that modulate protein function. However, mutations are pleiotropic and can have different effects on organismal fitness depending on the environment, as well as opposite effects on protein function and dosage. We review the pleiotropic effects of mutations. We discuss how they affect the evolution of gene and protein function, and how these complex mutational effects dictate the likelihood and mechanism of gene duplication and divergence. We propose several factors that can affect the divergence of new protein functions, including mutational trade-offs and hidden, or apparently neutral, variation.

Enzyme promiscuity: a mechanistic and evolutionary perspective

Many, if not most, enzymes can promiscuously catalyze reactions, or act on substrates, other than those for which they evolved. Here, we discuss the structural, mechanistic, and evolutionary implications of this manifestation of infidelity of molecular recognition. We define promiscuity and related phenomena and also address their generality and physiological implications. We discuss the mechanistic enzymology of promiscuity--how enzymes, which generally exert exquisite specificity, catalyze other, and sometimes barely related, reactions. Finally, we address the hypothesis that promiscuous enzymatic activities serve as evolutionary starting points and highlight the unique evolutionary features of promiscuous enzyme functions.

Identification and characterization of beta-lactamase inhibitor protein-II (BLIP-II) interactions with beta-lactamases using phage display

Protein-protein interactions are critical to cellular processes yet the ability to predict and rationally design interactions is limited because of incomplete knowledge of the principles governing these interactions. The beta-lactamase inhibitory protein (BLIP)/beta-lactamase interaction has become a model system to investigate protein-protein interactions and has been the focus of several structural, thermodynamic and binding specificity studies. BLIP-II also inhibits beta-lactamase but has no sequence homology with BLIP. The structure of BLIP-II in complex with TEM-1 beta-lactamase revealed that BLIP-II has a completely different structure than BLIP but it interacts with the same protruding loop-helix region of TEM-1 as does BLIP. The significance of the individual interacting residues in molecular recognition by BLIP-II is currently unknown. Therefore, a phage display vector was developed with the purpose of expressing BLIP-II onto the surface of the M13 filamentous bacteriophage. The BLIP-II displayed phage bound to TEM-1 with picomolar affinity indicating that BLIP-II is properly folded while on the surface of the phage. The phage system, as well as enzyme inhibition assays with purified proteins, revealed that BLIP-II is a more potent inhibitor than BLIP for several class A beta-lactamases with K(i) values in the low picomolar range.

USER friendly DNA recombination (USERec): a simple and flexible near homology-independent method for gene library construction

USER friendly DNA recombination (USERec) is introduced as a near homology-independent method that allows the simultaneous recombination of an unprecedented number of 10 DNA fragments (approximately 40-400 bp) within a day. The large number of fragments and their ease of preparation enables the creation of libraries of much larger genetic diversity (potentially approximately 10(10)-10(11) sequences) than current alternative methods based on DNA truncation (ITCHY, SCRATCHY and SHIPREC) or type IIb restriction enzymes (SISDC). At the same time, the frequency of frameshifts in the recombined library is low (90% of the recombined sequences are in frame). Compared to overlap extension PCR, USERec also requires much reduced crossover sequence constraints (only a 5'-AN(4-8)T-3' motif) and fewer experimental steps. Based on its simplicity and flexibility, and the accessibility of large and high quality recombined DNA libraries, USERec is established as a convenient alternative for the combinatorial assembly of gene fragments (e.g. exon or domain shuffling) and for a number of applications in gene library construction, such as loop grafting and multi-site-directed or random mutagenesis.

Ultrahigh-throughput screening in drop-based microfluidics for directed evolution

The explosive growth in our knowledge of genomes, proteomes, and metabolomes is driving ever-increasing fundamental understanding of the biochemistry of life, enabling qualitatively new studies of complex biological systems and their evolution. This knowledge also drives modern biotechnologies, such as molecular engineering and synthetic biology, which have enormous potential to address urgent problems, including developing potent new drugs and providing environmentally friendly energy. Many of these studies, however, are ultimately limited by their need for even-higher-throughput measurements of biochemical reactions. We present a general ultrahigh-throughput screening platform using drop-based microfluidics that overcomes these limitations and revolutionizes both the scale and speed of screening. We use aqueous drops dispersed in oil as picoliter-volume reaction vessels and screen them at rates of thousands per second. To demonstrate its power, we apply the system to directed evolution, identifying new mutants of the enzyme horseradish peroxidase exhibiting catalytic rates more than 10 times faster than their parent, which is already a very efficient enzyme. We exploit the ultrahigh throughput to use an initial purifying selection that removes inactive mutants; we identify approximately 100 variants comparable in activity to the parent from an initial population of approximately 10(7). After a second generation of mutagenesis and high-stringency screening, we identify several significantly improved mutants, some approaching diffusion-limited efficiency. In total, we screen approximately 10(8) individual enzyme reactions in only 10 h, using < 150 microL of total reagent volume; compared to state-of-the-art robotic screening systems, we perform the entire assay with a 1,000-fold increase in speed and a 1-million-fold reduction in cost.

Evolutionary trajectories of beta-lactamase CTX-M-1 cluster enzymes: predicting antibiotic resistance

Extended-spectrum beta-lactamases (ESBL) constitute a key antibiotic-resistance mechanism affecting Gram-negative bacteria, and also an excellent model for studying evolution in real time. A shift in the epidemiology of ESBLs is being observed, which is characterized by the explosive diversification and increase in frequency of the CTX-M-type beta-lactamases in different settings. This provides a unique opportunity for studying a protein evolutionary radiation by the sequential acquisition of specific mutations enhancing protein efficiency and fitness concomitantly. The existence of driver antibiotic molecules favoring protein divergence has been investigated by combining evolutionary analyses and experimental site-specific mutagenesis. Phylogenetic reconstruction with all the CTX-M variants described so far provided a hypothetical evolutionary scenario showing at least three diversification events. CTX-M-3 was likely the enzyme at the origin of the diversification in the CTX-M-1 cluster, which was coincident with positive selection acting on several amino acid positions. Sixty-three CTX-M-3 derivatives containing all combinations of mutations under positively selected positions were constructed, and their phenotypic efficiency was evaluated. The CTX-M-3 diversification process can only be explained in a complex selective landscape with at least two antibiotics (cefotaxime and ceftazidime), indicating the need to invoke mixtures of selective drivers in order to understand the final evolutionary outcome. Under this hypothesis, we found congruent results between the in silico and in vitro analyses of evolutionary trajectories. Three pathways driving the diversification of CTX-M-3 towards the most complex and efficient variants were identified. Whereas the P167S pathway has limited possibilities of further diversification, the D240G route shows a robust diversification network. In the third route, drift may have played a role in the early stages of CTX-M-3 evolution. Antimicrobial agents should not be considered only as selectors for efficient mechanisms of resistance but also as diversifying agents of the evolutionary trajectories. Different trajectories were identified using a combination of phylogenetic reconstructions and directed mutagenesis analyses, indicating that such an approach might be useful to fulfill the desirable goal of predicting evolutionary trajectories in antimicrobial resistance.

Exploring protein fitness landscapes by directed evolution

Directed evolution circumvents our profound ignorance of how a protein's sequence encodes its function by using iterative rounds of random mutation and artificial selection to discover new and useful proteins. Proteins can be tuned to adapt to new functions or environments by simple adaptive walks involving small numbers of mutations. Directed evolution studies have shown how rapidly some proteins can evolve under strong selection pressures and, because the entire 'fossil record' of evolutionary intermediates is available for detailed study, they have provided new insight into the relationship between sequence and function. Directed evolution has also shown how mutations that are functionally neutral can set the stage for further adaptation.

Changing the DNA recognition specificity of the EcoDam DNA-(adenine-N6)-methyltransferase by directed evolution

EcoDam is an adenine-N6 DNA methyltransferase that methylates the GATC sites in the Escherichia coli genome. We have changed the target specificity of EcoDam from GATC to GATT by directed evolution, combining different random mutagenesis methods with restriction protection at GATT sites for selection and screening. By co-evolution of an enzyme library and a substrate library, we identified GATT as the best non-GATC site and discover a double mutation, R124S/P134S, as the first step to increase enzyme activity at GATT sites. After four generations of mutagenesis and selection, we obtained enzyme variants with new specificity for GATT. While the wild-type EcoDam shows no detectable activity at GATT sites in E. coli cells, some variants prefer methylation at GATT over GATC sites by about 10-fold in cells. In vitro DNA methylation kinetics carried out under single-turnover conditions using a hemimethylated GATC and a GATT oligonucleotide substrate confirmed that the evolved proteins prefer methylation of GATT sites to a similar degree. They show up to 1600-fold change in specificity in vitro and methylate the new GATT target site with 20% of the rate of GATC methylation by the wild-type enzyme, indicating good activity. We conclude that the new methyltransferases are fully functional in vivo and in vitro but show a new target-site specificity.

Following evolutionary paths to protein-protein interactions with high affinity and selectivity

How do intricate multi-residue features such as protein-protein interfaces evolve? To address this question, we evolved a new colicin-immunity binding interaction. We started with Im9, which inhibits its cognate DNase ColE9 at 10(-14) M affinity, and evolved it toward ColE7, which it inhibits promiscuously (Kd > 10(-8) M). Iterative rounds of random mutagenesis and selection toward higher affinity for ColE7, and selectivity (against ColE9 inhibition), led to an approximately 10(5)-fold increase in affinity and a 10(8)-fold increase in selectivity. Analysis of intermediates along the evolved variants revealed that changes in the binding configuration of the Im protein uncovered a latent set of interactions, thus providing the key to the rapid divergence of new Im7 variants. Overall, protein-protein interfaces seem to share the evolvability features of enzymes, that is, the exploitation of promiscuous interactions and alternative binding configurations via 'generalist' intermediates, and the key role of compensatory stabilizing mutations in facilitating the divergence of new functions.

Directed evolution of proteins in vitro using compartmentalization in emulsions

This unit describes a protocol for the directed evolution of proteins utilizing in vitro compartmentalization. This method uses a large number of independent in vitro transcription and translation (IVTT) reactions in water droplets suspended in an oil emulsion to enable selection of proteins that bind a target molecule. Protein variants that bind the target also bind to and allow recovery of the genes that encoded them. This protocol serves as a basis for carrying out selections in emulsions, and can potentially be modified to select for other functionalities, including catalysis. This selection method is advantageous compared to alternative selection protocols due to the ability to screen through very large-size libraries and the ability to express and screen or select for functions that would otherwise be toxic or inaccessible to in vivo selections and screens.

Directed evolution of serum paraoxonase PON3 by family shuffling and ancestor/consensus mutagenesis, and its biochemical characterization

Serum paraoxonases (PONs) are calcium-dependent lactonases with anti-atherogenic and detoxification functions. Here we describe the directed evolution and characterization of recombinant variants of serum paraoxonase PON3 that express in an active and soluble manner in Escherichia coli. These variants were obtained by combining family shuffling and phylogeny-based mutagenesis: the limited diversity of accessible, cloned PON3 genes was complemented by spiking the shuffling reaction with ancestor/consensus mutations, mutations to residues that comprise the consensus or appear in the predicted ancestors of the PON family. We screened the resulting libraries for PON3's lactonase activity while ensuring that the selected variants retained the substrate specificity of wild-type mammalian PON3s. The availability of highly stable, recombinant PON3 that is free of all other serum components enabled us to explore unknown biochemical features of PON3, including its binding to HDL particles, the effect of HDL on PON3's stability and enzymatic activity, and ex vivo tests of its anti-atherogenic properties. Overall, it appears that PON3 possesses properties very similar to those of PON1: the enzyme's lactonase activity is selectively stimulated by binding to apoAI-HDL, with a concomitant increase in its stability. PON3 also exhibits potentially anti-atherogenic functions, although at levels lower than those of PON1.

Chaperonin overexpression promotes genetic variation and enzyme evolution

Most protein mutations, and mutations that alter protein functions in particular, undermine stability and are therefore deleterious. Chaperones, or heat-shock proteins, are often implicated in buffering mutations, and could thus facilitate the acquisition of neutral genetic diversity and the rate of adaptation. We examined the ability of the Escherichia coli GroEL/GroES chaperonins to buffer destabilizing and adaptive mutations. Here we show that mutational drifts performed in vitro with four different enzymes indicated that GroEL/GroES overexpression doubled the number of accumulating mutations, and promoted the folding of enzyme variants carrying mutations in the protein core and/or mutations with higher destabilizing effects (destabilization energies of >3.5 kcal mol(-)(1), on average, versus approximately 1 kcal mol(-)(1) in the absence of GroEL/GroES). The divergence of modified enzymatic specificity occurred much faster under GroEL/GroES overexpression, in terms of the number of adapted variants (>or=2-fold) and their improved specificity and activity (>or=10-fold). These results indicate that protein stability is a major constraint in protein evolution, and buffering mechanisms such as chaperonins are key in alleviating this constraint.

Selection for chaperone-like mediated genetic robustness at low mutation rate: impact of drift, epistasis and complexity

Genetic robustness is defined as the constancy of a phenotype in the face of deleterious mutations. Overexpression of chaperones, to assist the folding of proteins carrying deleterious mutations, is so far one of the most accepted molecular mechanisms enhancing genetic robustness. Most theories on the evolution of robustness have focused on the implications of high mutation rate. Here we show that genetic drift, which is modulated by population size, organism complexity, and epistasis, can be a sufficient force to select for chaperone-mediated genetic robustness. Using an exact analytical solution, we also show that selection for costly genetic robustness leads to a paradox: the decrease of population fitness on long timescales and the long-term dependency on robustness mechanisms. We suggest that selection for genetic robustness could be universal and not restricted to high mutation rate organisms such as RNA viruses. The evolution of the endosymbiont Buchnera illustrates this selection mechanism and its paradox: the increased dependency on chaperones mediating genetic robustness. Our model explains why most chaperones might have become essential even in optimal growth conditions.

Potential role of phenotypic mutations in the evolution of protein expression and stability

Phenotypic mutations (errors occurring during protein synthesis) are orders of magnitude more frequent than genetic mutations. Consequently, the sequences of individual protein molecules transcribed and translated from the same gene can differ. To test the effects of such mutations, we established a bacterial system in which an antibiotic resistance gene (TEM-1 beta-lactamase) was transcribed by either a high-fidelity RNA polymerase or its error-prone mutant. This setup enabled the analysis of individual mRNA transcripts that were synthesized under normal or error-prone conditions. We found that an increase of approximately 20-fold in the frequency of transcription errors promoted the evolution of higher TEM-1 expression levels and of more stable enzyme variants. The stabilized variants exhibited a distinct advantage under error-prone transcription, although under normal transcription they conferred resistance similar to wild-type TEM-1. They did so, primarily, by increasing TEM-1's tolerance to destabilizing deleterious mutations that arise from transcriptional errors. The stabilized TEM-1 variants also showed increased tolerance to genetic mutations. Thus, although phenotypic mutations are not individually subjected to inheritance and natural selection, as are genetic mutations, they collectively exert a direct and immediate effect on protein fitness. They may therefore play a role in shaping protein traits such as expression levels, stability, and tolerance to genetic mutations.

Applications of flow cytometry in environmental microbiology and biotechnology

Flow cytometry (FCM) is a technique for counting, examining and sorting microscopic particles suspended in a stream of fluid. It uses the principles of light scattering, light excitation and the emission from fluorescent molecules to generate specific multiparameter data from particles and cells. The cells are hydrodynamically focussed in a sheath solution before being intercepted by a focused light source provided by a laser. FCM has been used primarily in medical applications but is being used increasingly for the examination of individual cells from environmental samples. It has found uses in the isolation of both culturable and hitherto non-culturable bacteria present infrequently in environmental samples using appropriate growth conditions. FCM lends itself to high-throughput applications in directed evolution for the analysis of single cells or cell populations carrying mutant genes. It is also suitable for encapsulation studies where individual bacteria are compartmentalised with substrate in water-in-oil-in-water emulsions or with individual genes in transcriptional/translational mixtures for the production of mutant enzymes. The sensitivity of the technique has allowed the examination of gene optimisation by a procedure known as random or neutral drift where screening and selection is based on the retention of some predetermined level of activity through multiple rounds of mutagenesis.

Directed enzyme evolution: climbing fitness peaks one amino acid at a time

Directed evolution can generate a remarkable range of new enzyme properties. Alternate substrate specificities and reaction selectivities are readily accessible in enzymes from families that are naturally functionally diverse. Activities on new substrates can be obtained by improving variants with broadened specificities or by step-wise evolution through a sequence of more and more challenging substrates. Evolution of highly specific enzymes has been demonstrated, even with positive selection alone. It is apparent that many solutions exist for any given problem, and there are often many paths that lead uphill, one step at a time.

Neutral evolution of proteins: The superfunnel in sequence space and its relation to mutational robustness

Following Kimura's neutral theory of molecular evolution [M. Kimura, The Neutral Theory of Molecular Evolution (Cambridge University Press, Cambridge, 1983) (reprinted in 1986)], it has become common to assume that the vast majority of viable mutations of a gene confer little or no functional advantage. Yet, in silico models of protein evolution have shown that mutational robustness of sequences could be selected for, even in the context of neutral evolution. The evolution of a biological population can be seen as a diffusion on the network of viable sequences. This network is called a "neutral network." Depending on the mutation rate mu and the population size N, the biological population can evolve purely randomly (muN<<1) or it can evolve in such a way as to select for sequences of higher mutational robustness (muN>>1). The stringency of the selection depends not only on the product muN but also on the exact topology of the neutral network, the special arrangement of which was named "superfunnel." Even though the relation between mutation rate, population size, and selection was thoroughly investigated, a study of the salient topological features of the superfunnel that could affect the strength of the selection was wanting. This question is addressed in this study. We use two different models of proteins: on lattice and off lattice. We compare neutral networks computed using these models to random networks. From this, we identify two important factors of the topology that determine the stringency of the selection for mutationally robust sequences. First, the presence of highly connected nodes ("hubs") in the network increases the selection for mutationally robust sequences. Second, the stringency of the selection increases when the correlation between a sequence's mutational robustness and its neighbors' increases. The latter finding relates a global characteristic of the neutral network to a local one, which is attainable through experiments or molecular modeling.