In the light of directed evolution: pathways of adaptive protein evolution

A subdomain swap strategy for reengineering nonribosomal peptides

Nonribosomal peptide synthetases (NRPSs) protect microorganisms from environmental threats by producing diverse siderophores, antibiotics, and other peptide natural products. Their modular molecular structure is also attractive from the standpoint of biosynthetic engineering. Here we evaluate a methodology for swapping module specificities of these mega-enzymes that takes advantage of flavodoxin-like subdomains involved in substrate recognition. Nine subdomains encoding diverse specificities were transplanted into the Phe-specific GrsA initiation module of gramicidin S synthetase. All chimeras could be purified as soluble protein. One construct based on a Val-specific subdomain showed sizable adenylation activity and functioned as a Val-Pro diketopiperazine synthetase upon addition of the proline-specific GrsB1 module. These results suggest that subdomain swapping could be a viable alternative to previous NRPS design approaches targeting binding pockets, domains, or entire modules. The short length of the swapped sequence stretch may facilitate straightforward exploitation of the wealth of existing NRPS modules for combinatorial biosynthesis.

Environmental changes bridge evolutionary valleys

In the basic fitness landscape metaphor for molecular evolution, evolutionary pathways are presumed to follow uphill steps of increasing fitness. How evolution can cross fitness valleys is an open question. One possibility is that environmental changes alter the fitness landscape such that low-fitness sequences reside on a hill in alternate environments. We experimentally test this hypothesis on the antibiotic resistance gene TEM-15 β-lactamase by comparing four evolutionary strategies shaped by environmental changes. The strategy that included initial steps of selecting for low antibiotic resistance (negative selection) produced superior alleles compared with the other three strategies. We comprehensively examined possible evolutionary pathways leading to one such high-fitness allele and found that an initially deleterious mutation is key to the allele's evolutionary history. This mutation is an initial gateway to an otherwise relatively inaccessible area of sequence space and participates in higher-order, positive epistasis with a number of neutral to slightly beneficial mutations. The ability of negative selection and environmental changes to provide access to novel fitness peaks has important implications for natural evolutionary mechanisms and applied directed evolution.

The Evolution of New Catalytic Mechanisms for Xenobiotic Hydrolysis in Bacterial Metalloenzymes

An increasing number of bacterial metalloenzymes have been shown to catalyse the breakdown of xenobiotics in the environment, while others exhibit a variety of promiscuous xenobiotic-degrading activities. Several different evolutionary processes have allowed these enzymes to gain or enhance xenobiotic-degrading activity. In this review, we have surveyed the range of xenobiotic-degrading metalloenzymes, and discuss the molecular and catalytic basis for the development of new activities. We also highlight how our increased understanding of the natural evolution of xenobiotic-degrading metalloenzymes can be been applied to laboratory enzyme design.

Distal substitutions drive divergent DNA specificity among paralogous transcription factors through subdivision of conformational space

Many genomes contain families of paralogs--proteins with divergent function that evolved from a common ancestral gene after a duplication event. To understand how paralogous transcription factors evolve divergent DNA specificities, we examined how the glucocorticoid receptor and its paralogs evolved to bind activating response elements [(+)GREs] and negative glucocorticoid response elements (nGREs). We show that binding to nGREs is a property of the glucocorticoid receptor (GR) DNA-binding domain (DBD) not shared by other members of the steroid receptor family. Using phylogenetic, structural, biochemical, and molecular dynamics techniques, we show that the ancestral DBD from which GR and its paralogs evolved was capable of binding both nGRE and (+)GRE sequences because of the ancestral DBD's ability to assume multiple DNA-bound conformations. Subsequent amino acid substitutions in duplicated daughter genes selectively restricted protein conformational space, causing this dual DNA-binding specificity to be selectively enhanced in the GR lineage and lost in all others. Key substitutions that determined the receptors' response element-binding specificity were far from the proteins' DNA-binding interface and interacted epistatically to change the DBD's function through DNA-induced allosteric mechanisms. These amino acid substitutions subdivided both the conformational and functional space of the ancestral DBD among the present-day receptors, allowing a paralogous family of transcription factors to control disparate transcriptional programs despite high sequence identity.

Chemomimetic biocatalysis: exploiting the synthetic potential of cofactor-dependent enzymes to create new catalysts

Despite the astonishing breadth of enzymes in nature, no enzymes are known for many of the valuable catalytic transformations discovered by chemists. Recent work in enzyme design and evolution, however, gives us good reason to think that this will change. We describe a chemomimetic biocatalysis approach that draws from small-molecule catalysis and synthetic chemistry, enzymology, and molecular evolution to discover or create enzymes with non-natural reactivities. We illustrate how cofactor-dependent enzymes can be exploited to promote reactions first established with related chemical catalysts. The cofactors can be biological, or they can be non-biological to further expand catalytic possibilities. The ability of enzymes to amplify and precisely control the reactivity of their cofactors together with the ability to optimize non-natural reactivity by directed evolution promises to yield exceptional catalysts for challenging transformations that have no biological counterparts.

Evolving new protein-protein interaction specificity through promiscuous intermediates

Interacting proteins typically coevolve, and the identification of coevolving amino acids can pinpoint residues required for interaction specificity. This approach often assumes that an interface-disrupting mutation in one protein drives selection of a compensatory mutation in its partner during evolution. However, this model requires a non-functional intermediate state prior to the compensatory change. Alternatively, a mutation in one protein could first broaden its specificity, allowing changes in its partner, followed by a specificity-restricting mutation. Using bacterial toxin-antitoxin systems, we demonstrate the plausibility of this second, promiscuity-based model. By screening large libraries of interface mutants, we show that toxins and antitoxins with high specificity are frequently connected in sequence space to more promiscuous variants that can serve as intermediates during a reprogramming of interaction specificity. We propose that the abundance of promiscuous variants promotes the expansion and diversification of toxin-antitoxin systems and other paralogous protein families during evolution.

The virtue of innovation: innovation through the lenses of biological evolution

We rehearse the processes of innovation and discovery in general terms, using as our main metaphor the biological concept of an evolutionary fitness landscape. Incremental and disruptive innovations are seen, respectively, as successful searches carried out locally or more widely. They may also be understood as reflecting evolution by mutation (incremental) versus recombination (disruptive). We also bring a platonic view, focusing on virtue and memory. We use 'virtue' as a measure of efforts, including the knowledge required to come up with disruptive and incremental innovations, and 'memory' as a measure of their lifespan, i.e. how long they are remembered. Fostering innovation, in the evolutionary metaphor, means providing the wherewithal to promote novelty, good objective functions that one is trying to optimize, and means to improve one's knowledge of, and ability to navigate, the landscape one is searching. Recombination necessarily implies multi- or inter-disciplinarity. These principles are generic to all kinds of creativity, novel ideas formation and the development of new products and technologies.

Focused Directed Evolution of Aryl-Alcohol Oxidase in Saccharomyces cerevisiae by Using Chimeric Signal Peptides

Aryl-alcohol oxidase (AAO) is an extracellular flavoprotein that supplies ligninolytic peroxidases with H2O2 during natural wood decay. With a broad substrate specificity and highly stereoselective reaction mechanism, AAO is an attractive candidate for studies into organic synthesis and synthetic biology, and yet the lack of suitable heterologous expression systems has precluded its engineering by directed evolution. In this study, the native signal sequence of AAO from Pleurotus eryngii was replaced by those of the mating α-factor and the K1 killer toxin, as well as different chimeras of both prepro-leaders in order to drive secretion in Saccharomyces cerevisiae. The secretion of these AAO constructs increased in the following order: preproα-AAO > preαproK-AAO > preKproα-AAO > preproK-AAO. The chimeric preαproK-AAO was subjected to focused-directed evolution with the aid of a dual screening assay based on the Fenton reaction. Random mutagenesis and DNA recombination was concentrated on two protein segments (Met[α1]-Val109 and Phe392-Gln566), and an array of improved variants was identified, among which the FX7 mutant (harboring the H91N mutation) showed a dramatic 96-fold improvement in total activity with secretion levels of 2 mg/liter. Analysis of the N-terminal sequence of the FX7 variant confirmed the correct processing of the preαproK hybrid peptide by the KEX2 protease. FX7 showed higher stability in terms of pH and temperature, whereas the pH activity profiles and the kinetic parameters were maintained. The Asn91 lies in the flavin attachment loop motif, and it is a highly conserved residue in all members of the GMC superfamily, except for P. eryngii and P. pulmonarius AAO. The in vitro involution of the enzyme by restoring the consensus ancestor Asn91 promoted AAO expression and stability.

Viewing protein fitness landscapes through a next-gen lens

High-throughput sequencing has enabled many powerful approaches in biological research. Here, we review sequencing approaches to measure frequency changes within engineered mutational libraries subject to selection. These analyses can provide direct estimates of biochemical and fitness effects for all individual mutations across entire genes (and likely compact genomes in the near future) in genetically tractable systems such as microbes, viruses, and mammalian cells. The effects of mutations on experimental fitness can be assessed using sequencing to monitor time-dependent changes in mutant frequency during bulk competitions. The impact of mutations on biochemical functions can be determined using reporters or other means of separating variants based on individual activities (e.g., binding affinity for a partner molecule can be interrogated using surface display of libraries of mutant proteins and isolation of bound and unbound populations). The comprehensive investigation of mutant effects on both biochemical function and experimental fitness provide promising new avenues to investigate the connections between biochemistry, cell physiology, and evolution. We summarize recent findings from systematic mutational analyses; describe how they relate to a field rich in both theory and experimentation; and highlight how they may contribute to ongoing and future research into protein structure-function relationships, systems-level descriptions of cell physiology, and population-genetic inferences on the relative contributions of selection and drift.

Intermolecular epistasis shaped the function and evolution of an ancient transcription factor and its DNA binding sites

Complexes of specifically interacting molecules, such as transcription factor proteins (TFs) and the DNA response elements (REs) they recognize, control most biological processes, but little is known concerning the functional and evolutionary effects of epistatic interactions across molecular interfaces. We experimentally characterized all combinations of genotypes in the joint protein-DNA sequence space defined by an historical transition in TF-RE specificity that occurred some 500 million years ago in the DNA-binding domain of an ancient steroid hormone receptor. We found that rampant epistasis within and between the two molecules was essential to specific TF-RE recognition and to the evolution of a novel TF-RE complex with unique derived specificity. Permissive and restrictive epistatic mutations across the TF-RE interface opened and closed potential evolutionary paths accessible by the other, making the evolution of each molecule contingent on its partner's history and allowing a molecular complex with novel specificity to evolve.

Natural Variants of the KPC-2 Carbapenemase have Evolved Increased Catalytic Efficiency for Ceftazidime Hydrolysis at the Cost of Enzyme Stability

The spread of β-lactamases that hydrolyze penicillins, cephalosporins and carbapenems among Gram-negative bacteria has limited options for treating bacterial infections. Initially, Klebsiella pneumoniae carbapenemase-2 (KPC-2) emerged as a widespread carbapenem hydrolyzing β-lactamase that also hydrolyzes penicillins and cephalosporins but not cephamycins and ceftazidime. In recent years, single and double amino acid substitution variants of KPC-2 have emerged among clinical isolates that show increased resistance to ceftazidime. Because it confers multi-drug resistance, KPC β-lactamase is a threat to public health. In this study, the evolution of KPC-2 function was determined in nine clinically isolated variants by examining the effects of the substitutions on enzyme kinetic parameters, protein stability and antibiotic resistance profile. The results indicate that the amino acid substitutions associated with KPC-2 natural variants lead to increased catalytic efficiency for ceftazidime hydrolysis and a consequent increase in ceftazidime resistance. Single substitutions lead to modest increases in catalytic activity while the double mutants exhibit significantly increased ceftazidime hydrolysis and resistance levels. The P104R, V240G and H274Y substitutions in single and double mutant combinations lead to the largest increases in ceftazidime hydrolysis and resistance. Molecular modeling suggests that the P104R and H274Y mutations could facilitate ceftazidime hydrolysis through increased hydrogen bonding interactions with the substrate while the V240G substitution may enhance backbone flexibility so that larger substrates might be accommodated in the active site. Additionally, we observed a strong correlation between gain of catalytic function for ceftazidime hydrolysis and loss of enzyme stability, which is in agreement with the ‘stability-function tradeoff’ phenomenon. The high T_m of KPC-2 (66.5°C) provides an evolutionary advantage as compared to other class A enzymes such as TEM (51.5°C) and CTX-M (51°C) in that it can acquire multiple destabilizing substitutions without losing the ability to fold into a functional enzyme.

The absence of new antibiotics combined with the emergence of antibiotic-resistance enzymes like KPC-2 that can inactivate most β-lactam antibiotics has resulted in a longer duration of medical treatment, higher costs of medical care, and increased mortality. In recent years, a number of amino acid sequence variants of KPC-2 have been identified in clinical isolates worldwide suggesting continued evolution of resistance in KPC-2. In this study we have characterized nine clinically isolated variants of KPC-2 (KPC-3 to -11) that differ from the initial KPC-2 isolate by one to two amino acids. The KPC variants confer increased resistance to the antibiotic ceftazidime as compared to KPC-2. This increase in resistance is correlated with improved ability of the variant enzymes to hydrolyze the antibiotic. Additionally, the changes associated with increased ceftazidime hydrolysis also reduce the thermal stability of the enzyme, indicating the mutations that assist catalysis come with a cost on the overall stability of the enzyme. The high thermal stability of KPC-2 allows destabilizing mutations that enhance catalysis to accumulate while the enzyme retains a folded, functional structure. Thus, the high stability of KPC-2 provides an evolutionary advantage to acquire multiple mutations and retain function as compared to other β-lactamase enzymes.

The food and environmental safety of Bt crops

Bacillus thuringiensis (Bt) microbial pesticides have a 50-year history of safety in agriculture. Cry proteins are among the active insecticidal ingredients in these pesticides, and genes coding for Cry proteins have been introduced into agricultural crops using modern biotechnology. The Cry gene sequences are often modified to enable effective expression in planta and several Cry proteins have been modified to increase biological activity against the target pest(s). Additionally, the domains of different but structurally conserved Cry proteins can be combined to produce chimeric proteins with enhanced insecticidal properties. Environmental studies are performed and include invertebrates, mammals, and avian species. Mammalian studies used to support the food and feed safety assessment are also used to support the wild mammal assessment. In addition to the NTO assessment, the environmental assessment includes a comparative assessment between the Bt crop and the appropriate conventional control that is genetically similar but lacks the introduced trait to address unintended effects. Specific phenotypic, agronomic, and ecological characteristics are measured in the Bt crop and the conventional control to evaluate whether the introduction of the insect resistance has resulted in any changes that might cause ecological harm in terms of altered weed characteristics, susceptibility to pests, or adverse environmental impact. Additionally, environmental interaction data are collected in field experiments for Bt crop to evaluate potential adverse effects. Further to the agronomic and phenotypic evaluation, potential movement of transgenes from a genetically modified crop plants into wild relatives is assessed for a new pest resistance gene in a new crop. This review summarizes the evidence for safety of crops containing Cry proteins for humans, livestock, and other non-target organisms.

The food and environmental safety of Bt crops

Bacillus thuringiensis (Bt) microbial pesticides have a 50-year history of safety in agriculture. Cry proteins are among the active insecticidal ingredients in these pesticides, and genes coding for Cry proteins have been introduced into agricultural crops using modern biotechnology. The Cry gene sequences are often modified to enable effective expression in planta and several Cry proteins have been modified to increase biological activity against the target pest(s). Additionally, the domains of different but structurally conserved Cry proteins can be combined to produce chimeric proteins with enhanced insecticidal properties. Environmental studies are performed and include invertebrates, mammals, and avian species. Mammalian studies used to support the food and feed safety assessment are also used to support the wild mammal assessment. In addition to the NTO assessment, the environmental assessment includes a comparative assessment between the Bt crop and the appropriate conventional control that is genetically similar but lacks the introduced trait to address unintended effects. Specific phenotypic, agronomic, and ecological characteristics are measured in the Bt crop and the conventional control to evaluate whether the introduction of the insect resistance has resulted in any changes that might cause ecological harm in terms of altered weed characteristics, susceptibility to pests, or adverse environmental impact. Additionally, environmental interaction data are collected in field experiments for Bt crop to evaluate potential adverse effects. Further to the agronomic and phenotypic evaluation, potential movement of transgenes from a genetically modified crop plants into wild relatives is assessed for a new pest resistance gene in a new crop. This review summarizes the evidence for safety of crops containing Cry proteins for humans, livestock, and other non-target organisms.

Protein folding and binding can emerge as evolutionary spandrels through structural coupling

Binding interactions between proteins and other molecules mediate numerous cellular processes, including metabolism, signaling, and gene regulation. These interactions often evolve in response to changes in the protein's chemical or physical environment (such as the addition of an antibiotic). Several recent studies have shown the importance of folding stability in constraining protein evolution. Here we investigate how structural coupling between folding and binding--the fact that most proteins can only bind their targets when folded--gives rise to an evolutionary coupling between the traits of folding stability and binding strength. Using a biophysical and evolutionary model, we show how these protein traits can emerge as evolutionary "spandrels" even if they do not confer an intrinsic fitness advantage. In particular, proteins can evolve strong binding interactions that have no functional role but merely serve to stabilize the protein if its misfolding is deleterious. Furthermore, such proteins may have divergent fates, evolving to bind or not bind their targets depending on random mutational events. These observations may explain the abundance of apparently nonfunctional interactions among proteins observed in high-throughput assays. In contrast, for proteins with both functional binding and deleterious misfolding, evolution may be highly predictable at the level of biophysical traits: adaptive paths are tightly constrained to first gain extra folding stability and then partially lose it as the new binding function is developed. These findings have important consequences for our understanding of how natural and engineered proteins evolve under selective pressure.

Engineering the ligninolytic enzyme consortium

The ligninolytic enzyme consortium is one of the most-efficient oxidative systems found in nature, playing a pivotal role during wood decay and coal formation. Typically formed by high redox-potential oxidoreductases, this array of enzymes can be used within the emerging lignocellulose biorefineries in processes that range from the production of bioenergy to that of biomaterials. To ensure that these versatile enzymes meet industry standards and needs, they have been subjected to directed evolution and hybrid approaches that surpass the limits imposed by nature. This Opinion article analyzes recent achievements in this field, including the incipient groundbreaking research into the evolution of resurrected enzymes, and the engineering of ligninolytic secretomes to create consolidated bioprocessing microbes with synthetic biology applications.

A thermostable transketolase evolved for aliphatic aldehyde acceptors

Directed evolution of a thermostable transketolase yields catalysts with significant improvement in activity, enantioselectivity and substrate scope.

The amino acid alphabet and the architecture of the protein sequence-structure map. I. Binary alphabets

The correspondence between protein sequences and structures, or sequence-structure map, relates to fundamental aspects of structural, evolutionary and synthetic biology. The specifics of the mapping, such as the fraction of accessible sequences and structures, or the sequences' ability to fold fast, are dictated by the type of interactions between the monomers that compose the sequences. The set of possible interactions between monomers is encapsulated by the potential energy function. In this study, I explore the impact of the relative forces of the potential on the architecture of the sequence-structure map. My observations rely on simple exact models of proteins and random samples of the space of potential energy functions of binary alphabets. I adopt a graph perspective and study the distribution of viable sequences and the structures they produce, as networks of sequences connected by point mutations. I observe that the relative proportion of attractive, neutral and repulsive forces defines types of potentials, that induce sequence-structure maps of vastly different architectures. I characterize the properties underlying these differences and relate them to the structure of the potential. Among these properties are the expected number and relative distribution of sequences associated to specific structures and the diversity of structures as a function of sequence divergence. I study the types of binary potentials observed in natural amino acids and show that there is a strong bias towards only some types of potentials, a bias that seems to characterize the folding code of natural proteins. I discuss implications of these observations for the architecture of the sequence-structure map of natural proteins, the construction of random libraries of peptides, and the early evolution of the natural amino acid alphabet.

Deep mutational scanning: a new style of protein science

Mutagenesis provides insight into proteins, but only recently have assays that couple genotype to phenotype been used to assess the activities of as many as 1 million mutant versions of a protein in a single experiment. This approach-'deep mutational scanning'-yields large-scale data sets that can reveal intrinsic protein properties, protein behavior within cells and the consequences of human genetic variation. Deep mutational scanning is transforming the study of proteins, but many challenges must be tackled for it to fulfill its promise.

Tracking the molecular evolution of calcium permeability in a nicotinic acetylcholine receptor

Nicotinic acetylcholine receptors are a family of ligand-gated nonselective cationic channels that participate in fundamental physiological processes at both the central and the peripheral nervous system. The extent of calcium entry through ligand-gated ion channels defines their distinct functions. The α9α10 nicotinic cholinergic receptor, expressed in cochlear hair cells, is a peculiar member of the family as it shows differences in the extent of calcium permeability across species. In particular, mammalian α9α10 receptors are among the ligand-gated ion channels which exhibit the highest calcium selectivity. This acquired differential property provides the unique opportunity of studying how protein function was shaped along evolutionary history, by tracking its evolutionary record and experimentally defining the amino acid changes involved. We have applied a molecular evolution approach of ancestral sequence reconstruction, together with molecular dynamics simulations and an evolutionary-based mutagenesis strategy, in order to trace the molecular events that yielded a high calcium permeable nicotinic α9α10 mammalian receptor. Only three specific amino acid substitutions in the α9 subunit were directly involved. These are located at the extracellular vestibule and at the exit of the channel pore and not at the transmembrane region 2 of the protein as previously thought. Moreover, we show that these three critical substitutions only increase calcium permeability in the context of the mammalian but not the avian receptor, stressing the relevance of overall protein structure on defining functional properties. These results highlight the importance of tracking evolutionarily acquired changes in protein sequence underlying fundamental functional properties of ligand-gated ion channels.

The robustness and innovability of protein folds

Assignment of protein folds to functions indicates that >60% of folds carry out one or two enzymatic functions, while few folds, for example, the TIM-barrel and Rossmann folds, exhibit hundreds. Are there structural features that make a fold amenable to functional innovation (innovability)? Do these features relate to robustness--the ability to readily accumulate sequence changes? We discuss several hypotheses regarding the relationship between the architecture of a protein and its evolutionary potential. We describe how, in a seemingly paradoxical manner, opposite properties, such as high stability and rigidity versus conformational plasticity and structural order versus disorder, promote robustness and/or innovability. We hypothesize that polarity--differentiation and low connectivity between a protein's scaffold and its active-site--is a key prerequisite for innovability.

Inferring fitness landscapes by regression produces biased estimates of epistasis

The genotype-fitness map plays a fundamental role in shaping the dynamics of evolution. However, it is difficult to directly measure a fitness landscape in practice, because the number of possible genotypes is astronomical. One approach is to sample as many genotypes as possible, measure their fitnesses, and fit a statistical model of the landscape that includes additive and pairwise interactive effects between loci. Here, we elucidate the pitfalls of using such regressions by studying artificial but mathematically convenient fitness landscapes. We identify two sources of bias inherent in these regression procedures, each of which tends to underestimate high fitnesses and overestimate low fitnesses. We characterize these biases for random sampling of genotypes as well as samples drawn from a population under selection in the Wright-Fisher model of evolutionary dynamics. We show that common measures of epistasis, such as the number of monotonically increasing paths between ancestral and derived genotypes, the prevalence of sign epistasis, and the number of local fitness maxima, are distorted in the inferred landscape. As a result, the inferred landscape will provide systematically biased predictions for the dynamics of adaptation. We identify the same biases in a computational RNA-folding landscape as well as regulatory sequence binding data treated with the same fitting procedure. Finally, we present a method to ameliorate these biases in some cases.

The role of substrate binding pocket residues phenylalanine 176 and phenylalanine 196 on Pseudomonas sp. OX1 toluene o-xylene monooxygenase activity and regiospecificity

Saturation mutagenesis was used to generate eleven substitutions of toluene-o-xylene monooxygenase (ToMO) at alpha subunit (TouA) positions F176 and F196 among which nine were novel: F176H, F176N, F176S, F176T, F196A, F196L, F196T, F196Y, F196H, F196I, and F196V. By testing the substrates phenol, toluene, and naphthalene, these positions were found to influence ToMO oxidation activity and regiospecificity. Specifically, TouA variant F176H was identified that had 4.7-, 4.3-, and 1.8-fold faster hydroxylation activity towards phenol, toluene, and naphthalene, respectively, compared to native ToMO. The F176H variant also produced the novel product hydroquinone (61%) from phenol, made twofold more 2-naphthol from naphthalene (34% vs. 16% by the wild-type ToMO), and had the regiospecificity of toluene changed from 51% to 73% p-cresol. The TouA F176N variant had the most para-hydroxylation capability, forming p-cresol (92%) from toluene and hydroquinone (82%) from phenol as the major product, whereas native ToMO formed 30% o-cresol, 19% m-cresol, and 51% of p-cresol from toluene and 100% catechol from phenol. For naphthalene oxidation, TouA variant F176S exhibited the largest shift in the product distribution by producing threefold more 2-naphthol. Among the other F196 variants, F196L produced catechol from phenol two times faster than the wild-type enzyme. The TouA F196I variant produced twofold less o-cresol and 19% more p-cresol from toluene, and the TouA F196A variant produced 62% more 2-naphthol from naphthalene compared to wild-type ToMO. Both of these positions have never been studied through the saturation mutagenesis and some of the best substitutions uncovered here have never been predicted and characterized for aromatics hydroxylation.

Mutagenic Organized Recombination Process by Homologous IN vivo Grouping (MORPHING) for directed enzyme evolution

Approaches that depend on directed evolution require reliable methods to generate DNA diversity so that mutant libraries can focus on specific target regions. We took advantage of the high frequency of homologous DNA recombination in Saccharomyces cerevisiae to develop a strategy for domain mutagenesis aimed at introducing and in vivo recombining random mutations in defined segments of DNA. Mutagenic Organized Recombination Process by Homologous IN vivo Grouping (MORPHING) is a one-pot random mutagenic method for short protein regions that harnesses the in vivo recombination apparatus of yeast. Using this approach, libraries can be prepared with different mutational loads in DNA segments of less than 30 amino acids so that they can be assembled into the remaining unaltered DNA regions in vivo with high fidelity. As a proof of concept, we present two eukaryotic-ligninolytic enzyme case studies: i) the enhancement of the oxidative stability of a H₂O₂-sensitive versatile peroxidase by independent evolution of three distinct protein segments (Leu28-Gly57, Leu149-Ala174 and Ile199-Leu268); and ii) the heterologous functional expression of an unspecific peroxygenase by exclusive evolution of its native 43-residue signal sequence.

Dynamic New World: Refining Our View of Protein Structure, Function and Evolution

Proteins are crucial to the functioning of all lifeforms. Traditional understanding posits that a single protein occupies a single structure ("fold"), which performs a single function. This view is radically challenged with the recognition that high structural dynamism-the capacity to be extra "floppy"-is more prevalent in functional proteins than previously assumed. As reviewed here, this dynamic take on proteins affects our understanding of protein "structure", function, and evolution, and even gives us a glimpse into protein origination. Specifically, this review will discuss historical developments concerning protein structure, and important new relationships between dynamism and aspects of protein sequence, structure, binding modes, binding promiscuity, evolvability, and origination. Along the way, suggestions will be provided for how key parts of textbook definitions-that so far have excluded membership to intrinsically disordered proteins (IDPs)-could be modified to accommodate our more dynamic understanding of proteins.

New generation of biocatalysts for organic synthesis

The use of enzymes as catalysts for the preparation of novel compounds has received steadily increasing attention over the past few years. High demands are placed on the identification of new biocatalysts for organic synthesis. The catalysis of more ambitious reactions reflects the high expectations of this field of research. Enzymes play an increasingly important role as biocatalysts in the synthesis of key intermediates for the pharmaceutical and chemical industry, and new enzymatic technologies and processes have been established. Enzymes are an important part of the spectrum of catalysts available for synthetic chemistry. The advantages and applications of the most recent and attractive biocatalysts--reductases, transaminases, ammonia lyases, epoxide hydrolases, and dehalogenases--will be discussed herein and exemplified by the syntheses of interesting compounds.

Metalloenzyme design and engineering through strategic modifications of native protein scaffolds

Metalloenzymes are among the major targets of protein design and engineering efforts aimed at attaining novel and efficient catalysis for biochemical transformation and biomedical applications, due to the diversity of functions imparted by the metallo-cofactors along with the versatility of the protein environment. Naturally evolved protein scaffolds can often serve as robust foundations for sustaining artificial active sites constructed by rational design, directed evolution, or a combination of the two strategies. Accumulated knowledge of structure-function relationship and advancement of tools such as computational algorithms and unnatural amino acids incorporation all contribute to the design of better metalloenzymes with catalytic properties approaching the needs of practical applications.

Computationally designed libraries for rapid enzyme stabilization

The ability to engineer enzymes and other proteins to any desired stability would have wide-ranging applications. Here, we demonstrate that computational design of a library with chemically diverse stabilizing mutations allows the engineering of drastically stabilized and fully functional variants of the mesostable enzyme limonene epoxide hydrolase. First, point mutations were selected if they significantly improved the predicted free energy of protein folding. Disulfide bonds were designed using sampling of backbone conformational space, which tripled the number of experimentally stabilizing disulfide bridges. Next, orthogonal in silico screening steps were used to remove chemically unreasonable mutations and mutations that are predicted to increase protein flexibility. The resulting library of 64 variants was experimentally screened, which revealed 21 (pairs of) stabilizing mutations located both in relatively rigid and in flexible areas of the enzyme. Finally, combining 10-12 of these confirmed mutations resulted in multi-site mutants with an increase in apparent melting temperature from 50 to 85°C, enhanced catalytic activity, preserved regioselectivity and a >250-fold longer half-life. The developed Framework for Rapid Enzyme Stabilization by Computational libraries (FRESCO) requires far less screening than conventional directed evolution.

Stability-activity tradeoffs constrain the adaptive evolution of RubisCO

A well-known case of evolutionary adaptation is that of ribulose-1,5-bisphosphate carboxylase (RubisCO), the enzyme responsible for fixation of CO2 during photosynthesis. Although the majority of plants use the ancestral C3 photosynthetic pathway, many flowering plants have evolved a derived pathway named C4 photosynthesis. The latter concentrates CO2, and C4 RubisCOs consequently have lower specificity for, and faster turnover of, CO2. The C4 forms result from convergent evolution in multiple clades, with substitutions at a small number of sites under positive selection. To understand the physical constraints on these evolutionary changes, we reconstructed in silico ancestral sequences and 3D structures of RubisCO from a large group of related C3 and C4 species. We were able to precisely track their past evolutionary trajectories, identify mutations on each branch of the phylogeny, and evaluate their stability effect. We show that RubisCO evolution has been constrained by stability-activity tradeoffs similar in character to those previously identified in laboratory-based experiments. The C4 properties require a subset of several ancestral destabilizing mutations, which from their location in the structure are inferred to mainly be involved in enhancing conformational flexibility of the open-closed transition in the catalytic cycle. These mutations are near, but not in, the active site or at intersubunit interfaces. The C3 to C4 transition is preceded by a sustained period in which stability of the enzyme is increased, creating the capacity to accept the functionally necessary destabilizing mutations, and is immediately followed by compensatory mutations that restore global stability.

Tracing ancestral specificity of lectins: ancestral sequence reconstruction method as a new approach in protein engineering

Protein evolution is a process of molecular design leading to the diversity of functional proteins found in nature. Recent advances in bioinformatics and structural biology, in addition to recombinant protein expression techniques, enable us to analyze more directly the molecular evolution of proteins by a new method using ancestral sequence reconstruction (ASR), the so-called experimental molecular archaeology. ASR has been used to reveal molecular properties and structures correlating with changing geology, ecology, and physiology, and to identify the structure elements important to changing physiological functions to fill substantial gaps in the processes of protein evolution. In this chapter, we describe ASR as a new method of protein engineering studies, and their application to analyzing lectins, of which evolutionary processes and structural features contributing to molecular stability, specificity, and unique functions have been elucidated. Experimental molecular archeology using ASR and crystal structures of full-length ancestral proteins is useful in understanding the evolutionary process of the functional and structural diversified lectins by tracing ancestral specificities.

Random mutagenesis by error-prone pol plasmid replication in Escherichia coli

Directed evolution is an approach that mimics natural evolution in the laboratory with the goal of modifying existing enzymatic activities or of generating new ones. The identification of mutants with desired properties involves the generation of genetic diversity coupled with a functional selection or screen. Genetic diversity can be generated using PCR or using in vivo methods such as chemical mutagenesis or error-prone replication of the desired sequence in a mutator strain. In vivo mutagenesis methods facilitate iterative selection because they do not require cloning, but generally produce a low mutation density with mutations not restricted to specific genes or areas within a gene. For this reason, this approach is typically used to generate new biochemical properties when large numbers of mutants can be screened or selected. Here we describe protocols for an advanced in vivo mutagenesis method that is based on error-prone replication of a ColE1 plasmid bearing the gene of interest. Compared to other in vivo mutagenesis methods, this plasmid-targeted approach allows increased mutation loads and facilitates iterative selection approaches. We also describe the mutation spectrum for this mutagenesis methodology in detail, and, using cycle 3 GFP as a target for mutagenesis, we illustrate the phenotypic diversity that can be generated using our method. In sum, error-prone Pol I replication is a mutagenesis method that is ideally suited for the evolution of new biochemical activities when a functional selection is available.

Alanine scan of core positions in ubiquitin reveals links between dynamics, stability, and function

Mutations at solvent-inaccessible core positions in proteins can impact function through many biophysical mechanisms including alterations to thermodynamic stability and protein dynamics. As these properties of proteins are difficult to investigate, the impacts of core mutations on protein function are poorly understood for most systems. Here, we determined the effects of alanine mutations at all 15 core positions in ubiquitin on function in yeast. The majority (13 of 15) of alanine substitutions supported yeast growth as the sole ubiquitin. Both the two null mutants (I30A and L43A) were less stable to temperature-induced unfolding in vitro than wild type (WT) but were well folded at physiological temperatures. Heteronuclear NMR studies indicated that the L43A mutation reduces temperature stability while retaining a ground-state structure similar to WT. This structure enables L43A to bind to common ubiquitin receptors in vitro. Many of the core alanine ubiquitin mutants, including one of the null variants (I30A), exhibited an increased accumulation of high-molecular-weight species, suggesting that these mutants caused a defect in the processing of ubiquitin-substrate conjugates. In contrast, L43A exhibited a unique accumulation pattern with reduced levels of high-molecular-weight species and undetectable levels of free ubiquitin. When conjugation to other proteins was blocked, L43A ubiquitin accumulated as free ubiquitin in yeast. Based on these findings, we speculate that ubiquitin's stability to unfolding may be required for efficient recycling during proteasome-mediated substrate degradation.

Multiplex iterative plasmid engineering for combinatorial optimization of metabolic pathways and diversification of protein coding sequences

Engineering complex biological systems typically requires combinatorial optimization to achieve the desired functionality. Here, we present Multiplex Iterative Plasmid Engineering (MIPE), which is a highly efficient and customized method for combinatorial diversification of plasmid sequences. MIPE exploits ssDNA mediated λ Red recombineering for the introduction of mutations, allowing it to target several sites simultaneously and generate libraries of up to 10(7) sequences in one reaction. We also describe "restriction digestion mediated co-selection (RD CoS)", which enables MIPE to produce enhanced recombineering efficiencies with greatly simplified coselection procedures. To demonstrate this approach, we applied MIPE to fine-tune gene expression level in the 5-gene riboflavin biosynthetic pathway and successfully isolated a clone with 2.67-fold improved production in less than a week. We further demonstrated the ability of MIPE for highly multiplexed diversification of protein coding sequence by simultaneously targeting 23 codons scattered along the 750 bp sequence. We anticipate this method to benefit the optimization of diverse biological systems in synthetic biology and metabolic engineering.

Toxicological evaluation of proteins introduced into food crops

This manuscript focuses on the toxicological evaluation of proteins introduced into GM crops to impart desired traits. In many cases, introduced proteins can be shown to have a history of safe use. Where modifications have been made to proteins, experience has shown that it is highly unlikely that modification of amino acid sequences can make a non-toxic protein toxic. Moreover, if the modified protein still retains its biological function, and this function is found in related proteins that have a history of safe use (HOSU) in food, and the exposure level is similar to functionally related proteins, then the modified protein could also be considered to be "as-safe-as" those that have a HOSU. Within nature, there can be considerable evolutionary changes in the amino acid sequence of proteins within the same family, yet these proteins share the same biological function. In general, food crops such as maize, soy, rice, canola etc. are subjected to a variety of processing conditions to generate different food products. Processing conditions such as cooking, modification of pH conditions, and mechanical shearing can often denature proteins in these crops resulting in a loss of functional activity. These same processing conditions can also markedly lower human dietary exposure to (functionally active) proteins. Safety testing of an introduced protein could be indicated if its biological function was not adequately characterized and/or it was shown to be structurally/functionally related to proteins that are known to be toxic to mammals.

Systems approaches for synthetic biology: a pathway toward mammalian design

We review methods of understanding cellular interactions through computation in order to guide the synthetic design of mammalian cells for translational applications, such as regenerative medicine and cancer therapies. In doing so, we argue that the challenges of engineering mammalian cells provide a prime opportunity to leverage advances in computational systems biology. We support this claim systematically, by addressing each of the principal challenges to existing synthetic bioengineering approaches—stochasticity, complexity, and scale—with specific methods and paradigms in systems biology. Moreover, we characterize a key set of diverse computational techniques, including agent-based modeling, Bayesian network analysis, graph theory, and Gillespie simulations, with specific utility toward synthetic biology. Lastly, we examine the mammalian applications of synthetic biology for medicine and health, and how computational systems biology can aid in the continued development of these applications.

Polishing the craft of genetic diversity creation in directed evolution

Genetic diversity creation is a core technology in directed evolution where a high quality mutant library is crucial to its success. Owing to its importance, the technology in genetic diversity creation has seen rapid development over the years and its application has diversified into other fields of scientific research. The advances in molecular cloning and mutagenesis since 2008 were reviewed. Specifically, new cloning techniques were classified based on their principles of complementary overhangs, homologous sequences, overlapping PCR and megaprimers and the advantages, drawbacks and performances of these methods were highlighted. New mutagenesis methods developed for random mutagenesis, focused mutagenesis and DNA recombination were surveyed. The technical requirements of these methods and the mutational spectra were compared and discussed with references to commonly used techniques. The trends of mutant library preparation were summarised. Challenges in genetic diversity creation were discussed with emphases on creating "smart" libraries, controlling the mutagenesis spectrum and specific challenges in each group of mutagenesis methods. An outline of the wider applications of genetic diversity creation includes genome engineering, viral evolution, metagenomics and a study of protein functions. The review ends with an outlook for genetic diversity creation and the prospective developments that can have future impact in this field.

Mutants with higher stability and specific activity from a single thermosensitive variant of T7 RNA polymerase

A single strategy to select RNA polymerase from bacteriophage T7 (T7 RNAP) mutants in Escherichia coli with enhanced thermostability or enzymatic activity is described. T7 RNAP has the ability to specifically transcribe genes under control of T7 phage promoter. By using random mutagenesis of the T7 RNAP gene in combination with an appropriate screening at 25 and 42°C, we have generated and selected E.coli clones with temperature-sensitive phenotype in the presence of chloramphenicol. The resistance to chloramphenicol used to select these clones results from expression control of the chloramphenicol acetyl transferase gene by the T7 promoter. In a second phase, and using the thermosensitive T7 RNAP variants as template, a new round of random mutagenesis was performed. Combined to an appropriate screening strategy, 11 mutations (second-site T7 RNAP revertants) that restore the initial resistance to chloramphenicol at 42°C were identified. Nine of these mutations increase the thermal resistance of the wild-type T7 RNA. They include the five mutations previously described using different approaches and four novel mutations. One improves T7 RNA catalytic activity and one has no positive effect on the natural enzyme but increases the activity of some combined mutants. Additive effects of mutations amount to an increase of as much as 10°C in T1/2 compared with the wild-type enzyme and up to a 2-fold activity enhancement.

Path-based approach to random walks on networks characterizes how proteins evolve new functions

We develop a path-based approach to continuous-time random walks on networks with arbitrarily weighted edges. We describe an efficient numerical algorithm for calculating statistical properties of the stochastic path ensemble. After demonstrating our approach on two reaction rate problems, we present a biophysical model that describes how proteins evolve new functions while maintaining thermodynamic stability. We use our methodology to characterize dynamics of evolutionary adaptation, reproducing several key features observed in directed evolution experiments. We find that proteins generally fall into two qualitatively different regimes of adaptation depending on their binding and folding energetics.