Molecular determinants of protein evolvability.

The plethora of biological functions that sustain life is rooted in the remarkable evolvability of proteins. An emerging view highlights the importance of a protein's initial state in dictating evolutionary success. A deeper comprehension of the mechanisms that govern the evolvability of these initial states can provide invaluable insights into protein evolution. In this review, we describe several molecular determinants of protein evolvability, unveiled by experimental evolution and ancestral sequence reconstruction studies. We further discuss how genetic variation and epistasis can promote or constrain functional innovation and suggest putative underlying mechanisms. By establishing a clear framework for these determinants, we provide potential indicators enabling the forecast of suitable evolutionary starting points and delineate molecular mechanisms in need of deeper exploration.

Increasing the Soluble Expression and Whole-Cell Activity of the Plastic-Degrading Enzyme MHETase through Consensus Design.

The mono(2-hydroxyethyl) terephthalate hydrolase (MHETase) from Ideonella sakaiensis carries out the second step in the enzymatic depolymerization of poly(ethylene terephthalate) (PET) plastic into the monomers terephthalic acid (TPA) and ethylene glycol (EG). Despite its potential industrial and environmental applications, poor recombinant expression of MHETase has been an obstacle to its industrial application. To overcome this barrier, we developed an assay allowing for the medium-throughput quantification of MHETase activity in cell lysates and whole-cell suspensions, which allowed us to screen a library of engineered variants. Using consensus design, we generated several improved variants that exhibit over 10-fold greater whole-cell activity than wild-type (WT) MHETase. This is revealed to be largely due to increased soluble expression, which biochemical and structural analysis indicates is due to improved protein folding.

Measuring differential fitness costs and interactions between genetic cassettes using fluorescent spectrophotometry.

In this article, we present a method for designing, executing, and analyzing data from a microbial competition experiment. We use fluorescent reporters to label different competing strains and resolve individual growth curves using a fluorescent spectrophotometer. Our comprehensive data analysis pipeline integrates multiple experiments to simultaneously infer sources of variation, extract selection coefficients, and estimate the genetic contributions to fitness for various synthetic genetic cassettes (SGCs). To demonstrate the method, we employ a synthetic biological system based on Escherichia coli. Strains carry 1 of 10 different plasmids and one of three genomically integrated fluorescent markers. All strains are co-cultured to obtain real-time measurements of optical density (total population density) and fluorescence (sub-population densities). We identify challenges in calibrating between fluorescence and density and of fluorescent proteins maturing at different rates. To resolve these issues, we compare two methods of fluorescence calibration and correct for maturation by measuring in vivo maturation times. We provide evidence of genetic interactions occurring between our SGCs and further show how to use our statistical model to test some hypotheses about microbial growth and the costs of protein expression.IMPORTANCEFluorescently labeled co-cultures are becoming increasingly popular. The approach proposed here offers a high standard for experimental design and data analysis to measure selection coefficients and growth rates in competition. Measuring competitive differences is useful in many laboratory studies, allowing for fitness cost-correction of growth rates and ecological interactions and testing hypotheses in synthetic biology. Using time-resolved growth curves, rather than endpoint measurements, for competition assays allows us to construct a detailed scientific model that can be used to ask questions about fine-grained phenomena, such as bacterial growth dynamics, as well as higher-level phenomena, such as the interactions between synthetic cassette expression.

Insertions and Deletions (Indels): A Missing Piece of the Protein Engineering Jigsaw.

Over the years, protein engineers have studied nature and borrowed its tricks to accelerate protein evolution in the test tube. While there have been considerable advances, our ability to generate new proteins in the laboratory is seemingly limited. One explanation for these shortcomings may be that insertions and deletions (indels), which frequently arise in nature, are largely overlooked during protein engineering campaigns. The profound effect of indels on protein structures, by way of drastic backbone alterations, could be perceived as "saltation" events that bring about significant phenotypic changes in a single mutational step. Should we leverage these effects to accelerate protein engineering and gain access to unexplored regions of adaptive landscapes? In this Perspective, we describe the role played by indels in the functional diversification of proteins in nature and discuss their untapped potential for protein engineering, despite their often-destabilizing nature. We hope to spark a renewed interest in indels, emphasizing that their wider study and use may prove insightful and shape the future of protein engineering by unlocking unique functional changes that substitutions alone could never achieve.

Genetic Perturbation Alters Functional Substates in Alkaline Phosphatase.

Enzymes inherently exhibit molecule-to-molecule heterogeneity in their conformational and functional states, which is considered to be a key to the evolution of new functions. Single-molecule enzyme assays enable us to directly observe such multiple functional states or functional substates. Here, we quantitatively analyzed functional substates in the wild-type and 69 single-point mutants of Escherichia coli alkaline phosphatase by employing a high-throughput single-molecule assay with a femtoliter reactor array device. Interestingly, many mutant enzymes exhibited significantly heterogeneous functional substates with various types, while the wild-type enzyme showed a highly homogeneous substate. We identified a correlation between the degree of functional substates and the level of improvement in promiscuous activities. Our work provides much comprehensive evidence that the functional substates can be easily altered by mutations, and the evolution toward a new catalytic activity may involve the modulation of the functional substates.

Evolution of ß-lactamase-mediated cefiderocol resistance.

BACKGROUND: Cefiderocol is a novel siderophore ß-lactam with improved hydrolytic stability toward ß-lactamases, including carbapenemases, achieved by combining structural moieties of two clinically efficient cephalosporins, ceftazidime and cefepime. Consequently, cefiderocol represents a treatment alternative for infections caused by MDR Gram-negatives. OBJECTIVES: To study the role of cefiderocol on resistance development and on the evolution of ß-lactamases from all Ambler classes, including KPC-2, CTX-M-15, NDM-1, CMY-2 and OXA-48. METHODS: Directed evolution, using error-prone PCR followed by selective plating, was utilized to investigate how the production and the evolution of different ß-lactamases cause changes in cefiderocol susceptibility determined using microbroth dilution assays (MIC and IC50). RESULTS: We found that the expression of blaOXA-48 did not affect cefiderocol susceptibility. On the contrary, the expression of blaKPC-2, blaCMY-2, blaCTX-M-15 and blaNDM-1 substantially reduced cefiderocol susceptibility by 4-, 16-, 8- and 32-fold, respectively. Further, directed evolution on these enzymes showed that, with the acquisition of only 1-2 non-synonymous mutations, all ß-lactamases were evolvable to further cefiderocol resistance by 2- (NDM-1, CTX-M-15), 4- (CMY-2), 8- (OXA-48) and 16-fold (KPC-2). Cefiderocol resistance development was often associated with collateral susceptibility changes including increased resistance to ceftazidime and ceftazidime/avibactam as well as functional trade-offs against different ß-lactam drugs. CONCLUSIONS: The expression of contemporary ß-lactamase genes can potentially contribute to cefiderocol resistance development and the acquisition of mutations in these genes results in enzymes adapting to increasing cefiderocol concentrations. Resistance development caused clinically important cross-resistance, especially against ceftazidime and ceftazidime/avibactam.

Author Correction: Higher-order epistasis shapes the fitness landscape of a xenobiotic-degrading enzyme.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Higher-order epistasis shapes the fitness landscape of a xenobiotic-degrading enzyme.

Characterizing the adaptive landscapes that encompass the emergence of novel enzyme functions can provide molecular insights into both enzymatic and evolutionary mechanisms. Here, we combine ancestral protein reconstruction with biochemical, structural and mutational analyses to characterize the functional evolution of methyl-parathion hydrolase (MPH), an organophosphate-degrading enzyme. We identify five mutations that are necessary and sufficient for the evolution of MPH from an ancestral dihydrocoumarin hydrolase. In-depth analyses of the adaptive landscapes encompassing this evolutionary transition revealed that the mutations form a complex interaction network, defined in part by higher-order epistasis, that constrained the adaptive pathways available. By also characterizing the adaptive landscapes in terms of their functional activities towards three additional organophosphate substrates, we reveal that subtle differences in the polarity of the substrate substituents drastically alter the network of epistatic interactions. Our work suggests that the mutations function collectively to enable substrate recognition via subtle structural repositioning.

Evolutionary repurposing of a sulfatase: A new Michaelis complex leads to efficient transition state charge offset.

The recruitment and evolutionary optimization of promiscuous enzymes is key to the rapid adaptation of organisms to changing environments. Our understanding of the precise mechanisms underlying enzyme repurposing is, however, limited: What are the active-site features that enable the molecular recognition of multiple substrates with contrasting catalytic requirements? To gain insights into the molecular determinants of adaptation in promiscuous enzymes, we performed the laboratory evolution of an arylsulfatase to improve its initially weak phenylphosphonate hydrolase activity. The evolutionary trajectory led to a 100,000-fold enhancement of phenylphosphonate hydrolysis, while the native sulfate and promiscuous phosphate mono- and diester hydrolyses were only marginally affected (≤50-fold). Structural, kinetic, and in silico characterizations of the evolutionary intermediates revealed that two key mutations, T50A and M72V, locally reshaped the active site, improving access to the catalytic machinery for the phosphonate. Measured transition state (TS) charge changes along the trajectory suggest the creation of a new Michaelis complex (E•S, enzyme-substrate), with enhanced leaving group stabilization in the TS for the promiscuous phosphonate (ßleavinggroup from -1.08 to -0.42). Rather than altering the catalytic machinery, evolutionary repurposing was achieved by fine-tuning the molecular recognition of the phosphonate in the Michaelis complex, and by extension, also in the TS. This molecular scenario constitutes a mechanistic alternative to adaptation solely based on enzyme flexibility and conformational selection. Instead, rapid functional transitions between distinct chemical reactions rely on the high reactivity of permissive active-site architectures that allow multiple substrate binding modes.

Evolution of cyclohexadienyl dehydratase from an ancestral solute-binding protein.

The emergence of enzymes through the neofunctionalization of noncatalytic proteins is ultimately responsible for the extraordinary range of biological catalysts observed in nature. Although the evolution of some enzymes from binding proteins can be inferred by homology, we have a limited understanding of the nature of the biochemical and biophysical adaptations along these evolutionary trajectories and the sequence in which they occurred. Here we reconstructed and characterized evolutionary intermediate states linking an ancestral solute-binding protein to the extant enzyme cyclohexadienyl dehydratase. We show how the intrinsic reactivity of a desolvated general acid was harnessed by a series of mutations radiating from the active site, which optimized enzyme-substrate complementarity and transition-state stabilization and minimized sampling of noncatalytic conformations. Our work reveals the molecular evolutionary processes that underlie the emergence of enzymes de novo, which are notably mirrored by recent examples of computational enzyme design and directed evolution.

Evolutionary and molecular foundations of multiple contemporary functions of the nitroreductase superfamily.

Insight regarding how diverse enzymatic functions and reactions have evolved from ancestral scaffolds is fundamental to understanding chemical and evolutionary biology, and for the exploitation of enzymes for biotechnology. We undertook an extensive computational analysis using a unique and comprehensive combination of tools that include large-scale phylogenetic reconstruction to determine the sequence, structural, and functional relationships of the functionally diverse flavin mononucleotide-dependent nitroreductase (NTR) superfamily (>24,000 sequences from all domains of life, 54 structures, and >10 enzymatic functions). Our results suggest an evolutionary model in which contemporary subgroups of the superfamily have diverged in a radial manner from a minimal flavin-binding scaffold. We identified the structural design principle for this divergence: Insertions at key positions in the minimal scaffold that, combined with the fixation of key residues, have led to functional specialization. These results will aid future efforts to delineate the emergence of functional diversity in enzyme superfamilies, provide clues for functional inference for superfamily members of unknown function, and facilitate rational redesign of the NTR scaffold.

Functional Trade-Offs in Promiscuous Enzymes Cannot Be Explained by Intrinsic Mutational Robustness of the Native Activity.

The extent to which an emerging new function trades off with the original function is a key characteristic of the dynamics of enzyme evolution. Various cases of laboratory evolution have unveiled a characteristic trend; a large increase in a new, promiscuous activity is often accompanied by only a mild reduction of the native, original activity. A model that associates weak trade-offs with "evolvability" was put forward, which proposed that enzymes possess mutational robustness in the native activity and plasticity in promiscuous activities. This would enable the acquisition of a new function without compromising the original one, reducing the benefit of early gene duplication and therefore the selection pressure thereon. Yet, to date, no experimental study has examined this hypothesis directly. Here, we investigate the causes of weak trade-offs by systematically characterizing adaptive mutations that occurred in two cases of evolutionary transitions in enzyme function: (1) from phosphotriesterase to arylesterase, and (2) from atrazine chlorohydrolase to melamine deaminase. Mutational analyses in various genetic backgrounds revealed that, in contrast to the prevailing model, the native activity is less robust to mutations than the promiscuous activity. For example, in phosphotriesterase, the deleterious effect of individual mutations on the native phosphotriesterase activity is much larger than their positive effect on the promiscuous arylesterase activity. Our observations suggest a revision of the established model: weak trade-offs are not caused by an intrinsic robustness of the native activity and plasticity of the promiscuous activity. We propose that upon strong adaptive pressure for the new activity without selection against the original one, selected mutations will lead to the largest possible increases in the new function, but whether and to what extent they decrease the old function is irrelevant, creating a bias towards initially weak trade-offs and the emergence of generalist enzymes.

Revealing Unexplored Sequence-Function Space Using Sequence Similarity Networks.

The rapidly expanding number of protein sequences found in public databases can improve our understanding of how protein functions evolve. However, our current knowledge of protein function likely represents a small fraction of the diverse repertoire that exists in nature. Integrative computational methods can facilitate the discovery of new protein functions and enzymatic reactions through the observation and investigation of the complex sequence-structure-function relationships within protein superfamilies. Here, we highlight the use of sequence similarity networks (SSNs) to identify previously unexplored sequence and function space. We exemplify this approach using the nitroreductase (NTR) superfamily. We demonstrate that SSN investigations can provide a rapid and effective means to classify groups of proteins, therefore exposing experimentally unexplored sequences that may exhibit novel functionality. Integration of such approaches with systematic experimental characterization will expand our understanding of the functional diversity of enzymes and their associated physiological roles.

The role of protein dynamics in the evolution of new enzyme function.

Enzymes must be ordered to allow the stabilization of transition states by their active sites, yet dynamic enough to adopt alternative conformations suited to other steps in their catalytic cycles. The biophysical principles that determine how specific protein dynamics evolve and how remote mutations affect catalytic activity are poorly understood. Here we examine a 'molecular fossil record' that was recently obtained during the laboratory evolution of a phosphotriesterase from Pseudomonas diminuta to an arylesterase. Analysis of the structures and dynamics of nine protein variants along this trajectory, and three rationally designed variants, reveals cycles of structural destabilization and repair, evolutionary pressure to 'freeze out' unproductive motions and sampling of distinct conformations with specific catalytic properties in bi-functional intermediates. This work establishes that changes to the conformational landscapes of proteins are an essential aspect of molecular evolution and that change in function can be achieved through enrichment of preexisting conformational sub-states.

Conformational Tinkering Drives Evolution of a Promiscuous Activity through Indirect Mutational Effects.

How remote mutations can lead to changes in enzyme function at a molecular level is a central question in evolutionary biochemistry and biophysics. Here, we combine laboratory evolution with biochemical, structural, genetic, and computational analysis to dissect the molecular basis for the functional optimization of phosphotriesterase activity in a bacterial lactonase (AiiA) from the metallo-ß-lactamase (MBL) superfamily. We show that a 1000-fold increase in phosphotriesterase activity is caused by a more favorable catalytic binding position of the paraoxon substrate in the evolved enzyme that resulted from conformational tinkering of the active site through peripheral mutations. A nonmutated active site residue, Phe68, was displaced by ∼3 Å through the indirect effects of two second-shell trajectory mutations, allowing molecular interactions between the residue and paraoxon. Comparative mutational scanning, i.e., examining the effects of alanine mutagenesis on different genetic backgrounds, revealed significant changes in the functional roles of Phe68 and other nonmutated active site residues caused by the indirect effects of trajectory mutations. Our work provides a quantitative measurement of the impact of second-shell mutations on the catalytic contributions of nonmutated residues and unveils the underlying intramolecular network of strong epistatic mutational relationships between active site residues and more remote residues. Defining these long-range conformational and functional epistatic relationships has allowed us to better understand the subtle, but cumulatively significant, role of second- and third-shell mutations in evolution.

Noise-mean relationship in mutated promoters.

Gene expression depends on the frequency of transcription events (burst frequency) and on the number of mRNA molecules made per event (burst size). Both processes are encoded in promoter sequence, yet their dependence on mutations is poorly understood. Theory suggests that burst size and frequency can be distinguished by monitoring the stochastic variation (noise) in gene expression: Increasing burst size will increase mean expression without changing noise, while increasing burst frequency will increase mean expression and decrease noise. To reveal principles by which promoter sequence regulates burst size and frequency, we randomly mutated 22 yeast promoters chosen to span a range of expression and noise levels, generating libraries of hundreds of sequence variants. In each library, mean expression (m) and noise (coefficient of variation, η) varied together, defining a scaling curve: η(2) = b/m + η(ext)(2). This relation is expected if sequence mutations modulate burst frequency primarily. The estimated burst size (b) differed between promoters, being higher in promoter containing a TATA box and lacking a nucleosome-free region. The rare variants that significantly decreased b were explained by mutations in TATA, or by an insertion of an out-of-frame translation start site. The decrease in burst size due to mutations in TATA was promoter-dependent, but independent of other mutations. These TATA box mutations also modulated the responsiveness of gene expression to changing conditions. Our results suggest that burst size is a promoter-specific property that is relatively robust to sequence mutations but is strongly dependent on the interaction between the TATA box and promoter nucleosomes.

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.

Dynamics and constraints of enzyme evolution.

The wealth of distinct enzymatic functions found in nature is impressive and the on-going evolutionary divergence of enzymatic functions continues to generate new and efficient catalysts, which can be seen through the recent emergence of enzymes able to degrade xenobiotics. However, recreating such processes in the laboratory has been met with only moderate success. What are the factors that lead to suboptimal research outputs? In this review, we discuss constraints on enzyme evolution, which can restrict evolutionary trajectories and lead to evolutionary dead-ends. We highlight recent studies that have used experimental evolution to mimic different aspects of enzymatic adaptation under simple, controlled settings to shed light on evolutionary dynamics and constraints. A better understanding of these constraints will lead to the development of more efficient strategies for directed evolution and enzyme engineering.

Do viral proteins possess unique biophysical features?

Natural selection shapes the sequence, structure and biophysical properties of proteins to fit their environment. We hypothesize that highly thermostable proteins and viral proteins represent two opposing adaptation strategies. Thermostable proteins are highly compact and possess well-packed hydrophobic cores and intensely charged surfaces. By contrast, viral proteins, and RNA viral proteins in particular, display a high occurrence of disordered segments and loosely packed cores. These features might endow viral proteins with increased structural flexibility and effective ways to interact with the components of the host. They could also be related to high adaptability levels and mutation rates observed in viruses, thus, representing a unique strategy for buffering the deleterious effects of mutations, such that those that have little (interactions), have little to lose.