Search

Among the enzyme lineages that undoubtedly emerged prior to the last universal common ancestor is the so-called HUP, which includes Class I aminoacyl tRNA synthetases (AARSs) as well as enzymes mediating NAD, FAD, and CoA biosynthesis. Here, we provide a detailed analysis of HUP evolution, from emergence to structural and functional diversification. The HUP is a nucleotide binding domain that uniquely catalyzes adenylation via the release of pyrophosphate. In contrast to other ancient nucleotide binding domains with the αβα sandwich architecture, such as P-loop NTPases, the HUP’s most conserved feature is not phosphate binding, but rather ribose binding by backbone interactions to the tips of β1 and/or β4. Indeed, the HUP exhibits unusual evolutionary plasticity and, while ribose binding is conserved, the location and mode of binding to the base and phosphate moieties of the nucleotide, and to the substrate(s) reacting with it, have diverged with time, foremost along the emergence of the AARSs. The HUP also beautifully demonstrates how a well- packed scaffold combined with evolvable surface elements promotes evolutionary innovation. Finally, we offer a scenario for the emergence of the HUP from a seed βαβ fragment, and suggest that despite an identical architecture, the HUP and the Rossmann represent independent emergences.

Superoxide dismutases (SODs) are critical metalloenzymes mitigating the damages of the modern oxygenated world. However, the emergence of one family of SODs, the Fe/Mn SOD, has been recurrently proposed to predate the great oxygenation event (GOE). This ancient family lacks metal binding selectivity, but displays strong catalytic selectivity. Therefore, some homologues would only be active when bound to Fe or Mn, although others, dubbed cambialistic, would function when loaded with either ion. This posed the longstanding question about the identity of the cognate metal ion of the first SODs to emerge. In this work, we utilize ancestral sequence reconstruction techniques to infer the earliest SODs. We show that the "ancestors" are active in vivo and in vitro. Further, we test their metal specificity and demonstrate that they are cambialistic in nature. Our findings shed light on how the predicted Last Common Universal Ancestor was capable of dealing with decomposition of the superoxide anion, and the early relationship between life, oxygen, and metal ion availability.

Enzymes are well known for their catalytic abilities, some even reaching "catalytic perfection" in the sense that the reaction they catalyze has reached the physical bound of the diffusion rate. However, our growing understanding of enzyme superfamilies has revealed that only some share a catalytic chemistry while others share a substrate-handle binding motif, for example, for a particular phosphate group. This suggests that some families emerged through a "substrate-handle-binding-first" mechanism ("binding-first" for brevity) instead of "chemistry-first" and we are, therefore, left to wonder what the role of non-catalytic binders might have been during enzyme evolution. In the last of their eight seminal, back-to-back articles from 1976, John Albery and Jeremy Knowles addressed the question of enzyme evolution by arguing that the simplest mode of enzyme evolution is what they defined as "uniform binding" (parallel stabilization of all enzyme-bound states to the same degree). Indeed, we show that a uniform-binding proto-catalyst can accelerate a reaction, but only when catalysis is already present, that is, when the transition state is already stabilized to some degree. Thus, we sought an alternative explanation for the cases where substrate-handle-binding preceded any involvement of a catalyst. We find that evolutionary starting points that exhibit negative catalysis can redirect the reaction's course to a preferred product without need for rate acceleration or product release; that is, if they do not stabilize, or even destabilize, the transition state corresponding to an undesired product. Such a mechanism might explain the emergence of "binding-first" enzyme families like the aldolase superfamily.

Quinone methide (QM) chemistry is widely applied including in enzyme inhibitors. Typically, enzyme-mediated bond breaking releases a phenol product that rearranges into an electrophilic QM that in turn covalently modifies protein side chains. However, the factors that govern the reactivity of QM-based inhibitors and their mode of inhibition have not been systematically explored. Foremost, enzyme inactivation might occur in cis, whereby a QM molecule inactivates the very same enzyme molecule that released it, or by trans if the released QMs diffuse away and inactivate other enzyme molecules. We examined QM-based inhibitors for enzymes exhibiting phosphoester hydrolase activity. We tested different phenolic substituents and benzylic leaving groups, thereby modulating the rates of enzymatic hydrolysis, phenolate-to-QM rearrangement, and the electrophilicity of the resulting QM. By developing assays that distinguish between cis and trans inhibition, we have identified certain combinations of leaving groups and phenyl substituents that lead to inhibition in the cis mode, while other combinations gave trans inhibition. Our results suggest that cis-acting QM-based substrates could be used as activity-based probes to identify various phospho- and phosphono-ester hydrolases, and potentially other hydrolases.

Aminoacyl‐tRNA synthetases (AARSs) charge tRNA with their cognate amino acids. Many other enzymes use amino acids as substrates, yet discrimination against noncognate amino acids that threaten the accuracy of protein translation is a hallmark of AARSs. Comparing AARSs to these other enzymes allowed us to recognize patterns in molecular recognition and strategies used by evolution for exercising selectivity. Overall, AARSs are 2–3 orders of magnitude more selective than most other amino acid utilizing enzymes. AARSs also reveal the physicochemical limits of molecular discrimination. For example, amino acids smaller by a single methyl moiety present a discrimination ceiling of ~200, while larger ones can be discriminated by up to 105‐fold. In contrast, substrates larger by a hydroxyl group challenge AARS selectivity, due to promiscuous H‐ bonding with polar active site groups. This ‘hydroxyl paradox’ is resolved by editing. Indeed, when the physicochemical discrimination limits are reached, post‐transfer editing – hydrolysis of tRNAs charged with noncognate amino acids, evolved. The editing site often selectively recognizes the edited noncognate substrate using the very same feature that the synthetic site could not efficiently discriminate against. Finally, the comparison to other enzymes also reveals that the selectivity of AARSs is an explicitly evolved trait, showing some clear examples of how selection acted not only to optimize catalytic efficiency with the target substrate, but also to abolish activity with noncognate threat substrates (‘negative selection’).

Although molecular oxygen is a relative newcomer to the biosphere, it has had a profound impact on metabolism. About 700 oxygen-dependent enzymatic reactions are known, the vast majority of which emerged only after the appearance of oxygen in the biosphere, circa 3 billion years ago. Oxygen was a major driving force for evolutionary innovation-~60% of all known oxygen-dependent enzyme families emerged as such; that is, the founding ancestor was an O

SummaryMulti-enzyme assemblies composed of metabolic enzymes catalyzing sequential reactions are being increasingly studied. Here, we report the discovery of a 1.6 megadalton multi-enzyme complex from Bacillus subtilis composed of two enzymes catalyzing opposite rather than sequential reactions (“counter-enzymes”): glutamate synthase (GltAB), and glutamate dehydrogenase (GudB), that make and break glutamate, respectively. In vivo and in vitro studies show that the primary role of complex formation is to inhibit GudB’s activity as this enzyme is constitutively expressed including in glutamate-limiting conditions. Using cryo-electron microscopy, we elucidated the structure of the complex and the basis of GudB’s inhibition. Finally, we show that this complex that exhibits unusual oscillatory progress curves is a necessity for planktonic growth in glutamate-limiting conditions, but is also essential for biofilm growth in glutamate-rich media, suggesting a regulatory role at fluctuating glutamate concentrations.

The vast majority of theoretically possible polypeptide chains do not fold, let alone confer function. Hence, protein evolution from preexisting building blocks has clear potential advantages over ab initio emergence from random sequences. In support of this view, sequence similarities between different proteins is generally indicative of common ancestry, and we collectively refer to such homologous sequences as ‘themes’. At the domain level, sequence homology is routinely detected. However, short themes which are segments, or fragments of intact domains, are particularly interesting because they may provide hints about the emergence of domains, as opposed to divergence of preexisting domains, or their mixing-and-matching to form multi-domain proteins. Here we identified 525 representative short themes, comprising 20-to-80 residues, that are unexpectedly shared between domains considered to have emerged independently. Among these ‘bridging themes’ are ones shared between the most ancient domains, e.g., Rossmann, P-loop NTPase, TIM-barrel, Flavodoxin, and Ferredoxin-like. We elaborate on several particularly interesting cases, where the bridging themes mediate ligand binding. Ligand binding may have contributed to the stability and the plasticity of these building blocks, and to their ability to invade preexisting domains or serve as starting points for completely new domains.

Oligomeric proteins are central to life. Duplication and divergence of their genes is a key evolutionary driver, also because duplications can yield very different outcomes. Given a homomeric ancestor, duplication can yield two paralogs that form two distinct homomeric complexes, or a heteromeric complex comprising both paralogs. Alternatively, one paralog remains a homomer while the other acquires a new partner. However, so far, conflicting trends have been noted with respect to which fate dominates, primarily because different methods and criteria are being used to assign the interaction status of paralogs. Here, we systematically analyzed all Saccharomyces cerevisiae and Escherichia coli oligomeric complexes that include paralogous proteins. We found that the proportions of homo-hetero duplication fates strongly depend on a variety of factors, yet that nonetheless, rigorous filtering gives a consistent picture. In E. coli about 50%, of the paralogous pairs appear to have retained the ancestral homomeric interaction, whereas in S. cerevisiae only ~10% retained a homomeric state. This difference was also observed when unique complexes were counted instead of paralogous gene pairs. We further show that this difference is accounted for by multiple cases of heteromeric yeast complexes that share common ancestry with homomeric bacterial complexes. Our analysis settles contradicting trends and conflicting previous analyses, and provides a systematic and rigorous pipeline for delineating the fate of duplicated oligomers in any organism for which protein-protein interaction data are available., Author summary About half of all proteins assemble as oligomers, either by self-interaction (homomers) or via interaction with another protein (heteromers). The latter can be unrelated, yet, quite commonly, the interacting proteins are paralogs, namely two genes that arose by gene duplication. Indeed, while a homomer is encoded by a single gene, heteromers demand two genes as a minimum. Duplication can therefore yield two discrete homomeric complexes or a single heteromer. Do paralogs tend to retain the ancestral homomeric interaction, or do they mostly diverge into heteromeric complexes? Despite several studies addressing this question, to date, we lack a systematic, rigorous approach for delineating the oligomeric fates of paralogs on an organism scale. To this end, we developed a new pipeline for analysis of molecular interaction databases that includes various filtering steps and unambiguous definitions of all possible oligomeric fates. Applying this method to Escherichia coli and Saccharomyces cerevisiae we noted that paralogous pairs tend to remain homomeric in the former while in the latter heteromeric complexes dominate. We consequently note a systematic trend of homomeric bacterial proteins diverging into heteromeric complexes in eukaryotes.

This Special Issue is composed of 10 reviews that delve into the intricacies behind enzyme promiscuity and evolution, an area that is of increasing interest in the biological research community. In particular, the reviews in this Special Issue explore enzyme promiscuity and evolution in the context of cellular metabolism, as discussed in this introductory Editorial. It is our hope that you enjoy these fascinating and informative reviews and we wish to thank the authors for their compelling contributions to The FEBS Journal. doi: 10.1111/febs.12650.

Evolutionary trajectories are deemed largely irreversible. In a newly diverged protein, reversion of mutations that led to the functional switch typically results in loss of both the new and the ancestral functions. Nonetheless, evolutionary transitions where reversions are viable have also been described. The structural and mechanistic causes of reversion compatibility versus incompatibility therefore remain unclear. We examined two laboratory evolution trajectories of mammalian paraoxonase-1, a lactonase with promiscuous organophosphate hydrolase (OPH) activity. Both trajectories began with the same active-site mutant, His115Trp, which lost the native lactonase activity and acquired higher OPH activity. A neo-functionalization trajectory amplified the promiscuous OPH activity, whereas the re-functionalization trajectory restored the native activity, thus generating a new lactonase that lacks His115. The His115 revertants of these trajectories indicated opposite trends. Revertants of the neo-functionalization trajectory lost both the evolved OPH and the original lactonase activity. Revertants of the trajectory that restored the original lactonase function were, however, fully active. Crystal structures and molecular simulations show that in the newly diverged OPH, the reverted His115 and other catalytic residues are displaced, thus causing loss of both the original and the new activity. In contrast, in the re-functionalization trajectory, reversion compatibility of the original lactonase activity derives from mechanistic versatility whereby multiple residues can fulfill the same task. This versatility enables unique sequence-reversible compositions that are inaccessible when the active site was repurposed toward a new function.

The vast majority of theoretically possible polypeptide chains do not fold, let alone confer function. Hence, protein evolution from preexisting building blocks has clear potential advantages over ab initio emergence from random sequences. In support of this view, sequence similarities between different proteins is generally indicative of common ancestry, and we collectively refer to such homologous sequences as "themes." At the domain level, sequence homology is routinely detected. However, short themes which are segments, or fragments of intact domains, are particularly interesting because they may provide hints about the emergence of domains, as opposed to divergence of preexisting domains, or their mixing-and-matching to form multi-domain proteins. Here we identified 525 representative short themes, comprising 20-80 residues that are unexpectedly shared between domains considered to have emerged independently. Among these "bridging themes" are ones shared between the most ancient domains, for example, Rossmann, P-loop NTPase, TIM-barrel, flavodoxin, and ferredoxin-like. We elaborate on several particularly interesting cases, where the bridging themes mediate ligand binding. Ligand binding may have contributed to the stability and the plasticity of these building blocks, and to their ability to invade preexisting domains or serve as starting points for completely new domains.

Under stress, metabolism is changing: specific up- or down-regulation of proteins and metabolites occurs as well as side effects. Distinguishing specific stress-signaling metabolites (alarmones) from side products (damage metabolites) is not trivial. One example is diadenosine tetraphosphate (Ap4A) - a side product of aminoacyl-tRNA synthetases found in all domains of life. The earliest observations suggested that Ap4A serves as an alarmone for heat stress in Escherichia coli. However, despite 50 years of research, the signaling mechanisms associated with Ap4A remain unknown. We defined a set of criteria for distinguishing alarmones from damage metabolites to systematically classify Ap4A. In a nutshell, no indications for a signaling cascade that is triggered by Ap4A were found; rather, we found that Ap4A is efficiently removed in a constitutive, nonregulated manner. Several fold perturbations in Ap4A concentrations have no effect, yet accumulation at very high levels is toxic due to disturbance of zinc homeostasis, and also because Ap4A's structural overlap with ATP can result in spurious binding and inactivation of ATP-binding proteins. Overall, Ap4A met all criteria for a damage metabolite. While we do not exclude any role in signaling, our results indicate that the damage metabolite option should be considered as the null hypothesis when examining Ap4A and other metabolites whose levels change upon stress.

The linkage between regulatory elements of transcription, such as promoters, and their protein products is central to gene function. Promoter-protein coevolution is therefore expected, but rarely observed, and the manner by which these two regulatory levels are linked remains largely unknown. We study glutamate dehydrogenase-a hub of carbon and nitrogen metabolism. In Bacillus subtilis, two paralogues exist: GudB is constitutively transcribed whereas RocG is tightly regulated. In their active, oligomeric states, both enzymes show similar enzymatic rates. However, swaps of enzymes and promoters cause severe fitness losses, thus indicating promoter-enzyme coevolution. Characterization of the proteins shows that, compared to RocG, GudB's enzymatic activity is highly dependent on glutamate and pH Promoter-enzyme swaps therefore result in excessive glutamate degradation when expressing a constitutive enzyme under a constitutive promoter, or insufficient activity when both the enzyme and its promoter are tightly regulated. Coevolution of transcriptional and enzymatic regulation therefore underlies paralogue-specific spatio-temporal control, especially under diverse growth conditions.

Improving an enzyme's initially low catalytic efficiency with a new target substrate by an order of magnitude or two may require only a few rounds of mutagenesis and screening or selection. However, subsequent rounds of optimization tend to yield decreasing degrees of improvement (diminishing returns) eventually leading to an optimization plateau. We aimed to optimize the catalytic efficiency of bacterial phosphotriesterase (PTE) toward V-type nerve agents. Previously, we improved the catalytic efficiency of wild-type PTE toward the nerve agent VX by 500-fold, to a catalytic efficiency (kcat/KM) of 5 × 106 M-1 min-1. However, effective in vivo detoxification demands an enzyme with a catalytic efficiency of >107 M-1 min-1. Here, following eight additional rounds of directed evolution and the computational design of a stabilized variant, we evolved PTE variants that detoxify VX with a kcat/KM ≥ 5 × 107 M-1 min-1 and Russian VX (RVX) with a kcat/KM ≥ 107 M-1 min-1. These final 10-fold improvements were the most time consuming and laborious, as most libraries yielded either minor or no improvements. Stabilizing the evolving enzyme, and avoiding tradeoffs in activity with different substrates, enabled us to obtain further improvements beyond the optimization plateau and evolve PTE variants that were overall improved by >5000-fold with VX and by >17 000-fold with RVX. The resulting variants also hydrolyze G-type nerve agents with high efficiency (GA, GB at kcat/KM > 5 × 107 M-1 min-1) and can thus serve as candidates for broad-spectrum nerve-agent prophylaxis and post-exposure therapy using low enzyme doses.

Nerve agents are organophosphates that potently inhibit acetylcholinesterase and their enzymatic detoxification has been a long-standing goal. Nerve agents vary widely in size, charge, hydrophobicity, and the cleavable ester bond. A single enzyme is therefore unlikely to efficiently hydrolyze all agents. Here, we describe a mixture of three previously developed variants of the bacterial phosphotriesterase (Bd-PTE) that are highly stable and nearly sequence identical. This mixture enables effective detoxification of a broad spectrum of known threat agents – GA (tabun), GB (sarin), GD (soman), GF (cyclosarin), VX, and Russian-VX. The potential for dimer dissociation and exchange that could inactivateBd-PTE has minimal impact, and the three enzyme variants are as active in a mixture as they are individually. To our knowledge, this engineered enzyme ‘cocktail’ comprises the first solution for enzymatic detoxification of the entire range of threat nerve agents.

Aminoacyl-tRNA synthetases (aaRSs), the enzymes responsible for coupling tRNAs to their cognate amino acids, minimize translational errors by intrinsic hydrolytic editing. Here, we compared norvaline (Nva), a linear amino acid not coded for protein synthesis, to the proteinogenic, branched valine (Val) in their propensity to mistranslate isoleucine (Ile) in proteins. We show that in the synthetic site of isoleucyl- tRNA synthetase (IleRS), Nva and Val are activated and transferred to tRNA at similar rates. The efficiency of the synthetic site in pre-transfer editing of Nva and Val also appears to be similar. Post-transfer editing was, however, more rapid with Nva and consequently IleRS misaminoacylates Nva- tRNA(Ile) at slower rate than Val-tRNA(Ile). Accordingly, an Escherichia coli strain lacking IleRS post-transfer editing misincorporated Nva and Val in the proteome to a similar extent and at the same Ile positions. However, Nva mistranslation inflicted higher toxicity than Val, in agreement with IleRS editing being optimized for hydrolysis of Nva-tRNA(Ile). Furthermore, we found that the evolutionary- related IleRS, leucyl- and valyl-tRNA synthetases (I/L/VRSs), all efficiently hydrolyze Nva-tRNAs even when editing of Nva seems redundant. We thus hypothesize that editing of Nva-tRNAs had already existed in the last common ancestor of I/L/VRSs, and that the editing domain of I/L/VRSs had primarily evolved to prevent infiltration of Nva into modern proteins.

Ancestral reconstruction provides instrumental insights regarding the biochemical and biophysical characteristics of past proteins. A striking observation relates to the remarkably high thermostability of reconstructed ancestors. The latter has been linked to high environmental temperatures in the Precambrian era, the era relating to most reconstructed proteins. We found that inferred ancestors of the serum paraoxonase (PON) enzyme family, including the mammalian ancestor, exhibit dramatically increased thermostabilities compared with the extant, human enzyme (up to 30 °C higher melting temperature). However, the environmental temperature at the time of emergence of mammals is presumed to be similar to the present one. Additionally, the mammalian PON ancestor has superior folding properties (kinetic stability)-unlike the extant mammalian PONs, it expresses in E. coli in a soluble and functional form, and at a high yield. We discuss two potential origins of this unexpectedly high stability. First, ancestral stability may be overestimated by a "consensus effect," whereby replacing amino acids that are rare in contemporary sequences with the amino acid most common in the family increases protein stability. Comparison to other reconstructed ancestors indicates that the consensus effect may bias some but not all reconstructions. Second, we note that high stability may relate to factors other than high environmental temperature such as oxidative stress or high radiation levels. Foremost, intrinsic factors such as high rates of genetic mutations and/or of transcriptional and translational errors, and less efficient protein quality control systems, may underlie the high kinetic and thermodynamic stability of past proteins.

The last decade has seen a growing number of experiments aimed at systematically mapping the effects of mutations in different proteins, and of attempting to correlate their biophysical and biochemical effects with organismal fitness. While insightful, systematic laboratory measurements of fitness effects present challenges and difficulties. Here, we discuss the limitations associated with such measurements, and in particular the challenge of correlating the effects of mutations at the single protein level ("protein fitness") with their effects on organismal fitness. A variety of experimental setups are used, with some measuring the direct effects on protein function and others monitoring the growth rate of a model organism carrying the protein mutants. The manners by which fitness effects are calculated and presented also vary, and the conclusions, including the derived distributions of fitness effects of mutations, vary accordingly. The comparison of the effects of mutations in the laboratory to the natural protein diversity, namely to amino acid changes that have fixed in the course of millions of years of evolution, is also debatable. The results of laboratory experiments may, therefore, be less relevant to understanding long-term inter-species variations yet insightful with regard to short-term polymorphism, for example, in the study of the effects of human SNPs.

Substantial improvements in enzyme activity demand multiple mutations at spatially proximal positions in the active site. Such mutations, however, often exhibit unpredictable epistatic (non-additive) effects on activity. Here we describe FuncLib, an automated method for designing multipoint mutations at enzyme active sites using phylogenetic analysis and Rosetta design calculations. We applied FuncLib to two unrelated enzymes, a phosphotriesterase and an acetyl-CoA synthetase. All designs were active, and most showed activity profiles that significantly differed from the wild-type and from one another. Several dozen designs with only 3-6 active-site mutations exhibited 10- to 4,000-fold higher efficiencies with a range of alternative substrates, including hydrolysis of the toxic organophosphate nerve agents soman and cyclosarin and synthesis of butyryl-CoA. FuncLib is implemented as a web server (http://FuncLib.weizmann.ac.il); it circumvents iterative, high-throughput experimental screens and opens the way to designing highly efficient and diverse catalytic repertoires.

The emergence of catalysis in a noncatalytic protein scaffold is a rare, unexplored event. Chalcone isomerase (CHI), a key enzyme in plant flavonoid biosynthesis, is presumed to have evolved from a nonenzymatic ancestor related to the widely distributed fatty-acid binding proteins (FAPs) and a plant protein family with no isomerase activity (CHILs). Ancestral inference supported the evolution of CHI from a protein lacking isomerase activity. Further, we identified four alternative founder mutations, i.e., mutations that individually instated activity, including a mutation that is not phylogenetically traceable. Despite strong epistasis in other cases of protein evolution, CHI's laboratory reconstructed mutational trajectory shows weak epistasis. Thus, enantioselective CHI activity could readily emerge despite a catalytically inactive starting point. Accordingly, X-ray crystallography, NMR, and molecular dynamics simulations reveal reshaping of the active site toward a productive substrate-binding mode and repositioning of the catalytic arginine that was inherited from the ancestral fatty-acid binding proteins.

Ancestral reconstruction is a powerful tool for studying protein evolution as well as for protein design and engineering. However, in many positions alternative predictions with relatively high marginal probabilities exist, and thus the prediction comprises an ensemble of near-ancestor sequences that relate to the historical ancestor. The ancestral phenotype should therefore be explored for the entire ensemble, rather than for the sequence comprising the most probable amino acid at all positions [the most probable ancestor (mpa)]. To this end, we constructed libraries that sample ensembles of near-ancestor sequences. Specifically, we identified positions where alternatively predicted amino acids are likely to affect the ancestor's structure and/or function. Using the serum paraoxonases (PONs) enzyme family as a test case, we constructed libraries that combinatorially sample these alternatives. We next characterized these libraries, reflecting the vertebrate and mammalian PON ancestors. We found that the mpa of vertebrate PONs represented only one out of many different enzymatic phenotypes displayed by its ensemble. The mammalian ancestral library, however, exhibited a homogeneous phenotype that was well represented by the mpa. Our library design strategy that samples near-ancestor ensembles at potentially critical positions therefore provides a systematic way of examining the robustness of inferred ancestral phenotypes.

Despite the abundance of membrane-associated enzymes, the mechanism by which membrane binding stabilizes these enzymes and stimulates their catalysis remains largely unknown. Serum paraoxonase-1 (PON1) is a lipophilic lactonase whose stability and enzymatic activity are dramatically stimulated when associated with high-density lipoprotein (HDL) particles. Our mutational and structural analyses, combined with empirical valence bond simulations, reveal a network of hydrogen bonds that connect HDL binding residues with Asn168--a key catalytic residue residing >15A from the HDL contacting interface. This network ensures precise alignment of N168, which, in turn, ligates PON1's catalytic calcium and aligns the lactone substrate for catalysis. HDL binding restrains the overall motion of the active site and particularly of N168, thus reducing the catalytic activation energy barrier. We demonstrate herein that disturbance of this network, even at its most far-reaching periphery, undermines PON1's activity. Membrane binding thus immobilizes long-range interactions via second- and third-shell residues that reduce the active site's floppiness and pre-organize the catalytic residues. Although this network is critical for efficient catalysis, as demonstrated here, unraveling these long-rage interaction networks is challenging, let alone their implementation in artificial enzyme design.

To carry out their activities, biological macromolecules balance different physical traits, such as stability, interaction affinity, and selectivity. How such often opposing traits are encoded in a macromolecular system is critical to our understanding of evolutionary processes and ability to design new molecules with desired functions. We present a framework for constraining design simulations to balance different physical characteristics. Each trait is represented by the equilibrium fractional occupancy of the desired state relative to its alternatives, ranging from none to full occupancy, and the different traits are combined using Boolean operators to effect a “fuzzy”-logic language for encoding any combination of traits. In another paper, we presented a new combinatorial backbone design algorithm AbDesign where the fuzzy-logic framework was used to optimize protein backbones and sequences for both stability and binding affinity in antibody-design simulation. We now extend this framework and find that fuzzy-logic design simulations reproduce sequence and structure design principles seen in nature to underlie exquisite specificity on the one hand and multispecificity on the other hand. The fuzzy-logic language is broadly applicable and could help define the space of tolerated and beneficial mutations in natural biomolecular systems and design artificial molecules that encode complex characteristics., Graphical abstract, Highlights • Complex combinations of physical constraints shape the evolution of biomolecules. • An analytical system based on equilibrium thermodynamics and fuzzy logic is proposed to embody these constraints. • Computational design simulations subject to fuzzy-logic design recapitulate elements observed in natural proteins. • This design language can be used to analyze natural biomolecular systems and design new ones exhibiting complex traits.

The practical need for highly efficient enzymes presents new challenges in enzyme engineering, in particular, the need to improve catalytic turnover (kcat) or efficiency (kcat/KM) by several orders of magnitude. However, optimizing catalysis demands navigation through complex and rugged fitness landscapes, with optimization trajectories often leading to strong diminishing returns and dead-ends. When no further improvements are observed in library screens or selections, it remains unclear whether the maximal catalytic efficiency of the enzyme (the catalytic 'fitness peak') has been reached; or perhaps, an alternative combination of mutations exists that could yield additional improvements. Here, we discuss fundamental aspects of the process of catalytic optimization, and offer practical solutions with respect to overcoming optimization plateaus.

Organismal adaptation to extreme temperatures yields enzymes with distinct configurational stabilities, including thermophilic and psychrophilic enzymes, which are adapted to high and low temperatures, respectively. These enzymes are widely assumed to also have unique rate–temperature dependencies. Thermophilic enzymes, for example, are considered optimal at high temperatures and effectively inactive at low temperatures due to excess rigidity. Surveying published data, we find that thermophilic, mesophilic, and psychrophilic enzymes exhibit indistinguishable rate–temperature dependencies. Furthermore, given the nonenzymatic rate–temperature dependency, all enzymes, regardless of their operation temperatures, become >10-fold less powerful catalysts per 25°C temperature increase. Among other factors, this loss of rate acceleration may be ascribed to thermally induced vibrations compromising the active-site catalytic configuration, suggesting that many enzymes are in fact insufficiently rigid.

The last decade has seen a growing number of experiments aimed at systematically mapping the effects of mutations in different proteins, and of attempting to correlate their biophysical and biochemical effects with organismal fitness. While insightful, systematic laboratory measurements of fitness effects present challenges and difficulties. Here, we discuss the limitations associated with such measurements, and in particular the challenge of correlating the effects of mutations at the single protein level ("protein fitness") with their effects on organismal fitness. A variety of experimental setups are used, with some measuring the direct effects on protein function and others monitoring the growth rate of a model organism carrying the protein mutants. The manners by which fitness effects are calculated and presented also vary, and the conclusions, including the derived distributions of fitness effects of mutations, vary accordingly. The comparison of the effects of mutations in the laboratory to the natural protein diversity, namely to amino acid changes that have fixed in the course of millions of years of evolution, is also debatable. The results of laboratory experiments may, therefore, be less relevant to understanding long-term inter-species variations yet insightful with regard to short-term polymorphism, for example, in the study of the effects of human SNPs.

Serum paraoxonases (PONs) are detoxifying lactonases that were first identified in mammals. Three mammalian families are known, PON1, 2, and 3 that reside primarily in the liver. They catalyze essentially the same reaction, lactone hydrolysis, but differ in their substrate specificity. Although some members are highly specific, others have a broad specificity profile. The evolutionary origins and substrate specificities of PONs therefore remain poorly understood. Here, we report a newly identified family of bacterial PONs, and the reconstruction of the ancestor of the three families of mammalian PONs. Both the mammalian ancestor and the characterized bacterial PONX_OCCAL were found to efficiently hydrolyze N-acyl homoserine lactones that mediate quorum sensing in many bacteria, including pathogenic ones. The mammalian PONs may therefore relate to a newly identified family of bacterial, PON-like “quorum-quenching” lactonases. The appearance of PONs in metazoa is likely to relate to innate immunity rather than detoxification. Unlike the bacterial PON, the mammalian ancestor also hydrolyzes, with low efficiency, lactones other than homoserine lactones, thus preceding the detoxifying functions that diverged later in two of the three mammalian families. The bifunctionality of the mammalian ancestor and the trade-off between the quorum-quenching and detoxifying lactonase activities explain the broad and overlapping specificities of some mammalian PONs versus the singular specificity of others.

Protein evolvability includes two elements—robustness (or neutrality, mutations having no effect) and innovability (mutations readily inducing new functions). How are these two conflicting demands bridged? Does the ability to bridge them relate to the observation that certain folds, such as TIM barrels, accommodate numerous functions, whereas other folds support only one? Here, we hypothesize that the key to innovability is polarity—an active site composed of flexible, loosely packed loops alongside a well-separated, highly ordered scaffold. We show that highly stabilized variants of TEM-1 β-lactamase exhibit selective rigidification of the enzyme's scaffold while the active-site loops maintained their conformational plasticity. Polarity therefore results in stabilizing, compensatory mutations not trading off, but instead promoting the acquisition of new activities. Indeed, computational analysis indicates that in folds that accommodate only one function throughout evolution, for example, dihydrofolate reductase, ≥ 60% of the active-site residues belong to the scaffold. In contrast, folds associated with multiple functions such as the TIM barrel show high scaffold–active-site polarity (~ 20% of the active site comprises scaffold residues) and > 2-fold higher rates of sequence divergence at active-site positions. Our work suggests structural measures of fold polarity that appear to be correlated with innovability, thereby providing new insights regarding protein evolution, design, and engineering.

The field of directed evolution has progressed to the point where it is feasible to engineer enzymes for unnatural substrates and reactions with catalytic efficiencies and regio-specificity or stereo-specificity that rival those of natural enzymes. Here, we describe the conceptual and methodological advances that have enabled this progress. We address methodologies based on small libraries enriched with improved variants and carrying compensatory stabilizing mutations. Such libraries can be combined with low-throughput screens that provide high accuracy and directly target the desired substrate and reaction conditions, and thereby provide highly improved variants.

The origins of enzyme specificity are well established. However, the molecular details underlying the ability of a single active site to promiscuously bind different substrates and catalyze different reactions remain largely unknown. To better understand the molecular basis of enzyme promiscuity, we studied the mammalian serum paraoxonase 1 (PON1) whose native substrates are lipophilic lactones. We describe the crystal structures of PON1 at a catalytically relevant pH and of its complex with a lactone analogue. The various PON1 structures and the analysis of active-site mutants guided the generation of docking models of the various substrates and their reaction intermediates. The models suggest that promiscuity is driven by coincidental overlaps between the reactive intermediate for the native lactonase reaction and the ground and/or intermediate states of the promiscuous reactions. This overlap is also enabled by different active-site conformations: the lactonase activity utilizes one active-site conformation whereas the promiscuous phosphotriesterase activity utilizes another. The hydrolysis of phosphotriesters, and of the aromatic lactone dihydrocoumarin, is also driven by an alternative catalytic mode that uses only a subset of the active-site residues utilized for lactone hydrolysis. Indeed, PON1's active site shows a remarkable level of networking and versatility whereby multiple residues share the same task and individual active-site residues perform multiple tasks (e.g., binding the catalytic calcium and activating the hydrolytic water). Overall, the coexistence of multiple conformations and alternative catalytic modes within the same active site underlines PON1's promiscuity and evolutionary potential.

We discuss the basic features of divergent versus convergent evolution and of the common scenario of parallel evolution. The example of quorum-quenching lactonases is subsequently described. Three different quorum-quenching lactonase families are known, and they belong to three different superfamilies. Their key active-site architectures have converged and are strikingly similar. Curiously, a promiscuous organophosphate hydrolase activity is observed in all three families. We describe the structural and mechanistic features that underline this converged promiscuity and how this promiscuity drove the parallel divergence of organophosphate hydrolases within these lactonase families by either natural or laboratory evolution.

Although de novo computational enzyme design has been shown to be feasible, the field is still in its infancy: the kinetic parameters of designed enzymes are still orders of magnitude lower than those of naturally occurring ones. Nonetheless, designed enzymes can be improved by directed evolution, as recently exemplified for the designed Kemp eliminase KE07. Random mutagenesis and screening resulted in variants with > 200-fold higher catalytic efficiency and provided insights about features missing in the designed enzyme. Here we describe the optimization of KE70, another designed Kemp eliminase. Amino acid substitutions predicted to improve catalysis in design calculations involving extensive backbone sampling were individually tested. Those proven beneficial were combinatorially incorporated into the originally designed KE70 along with random mutations, and the resulting libraries were screened for improved eliminase activity. Nine rounds of mutation and selection resulted in > 400-fold improvement in the catalytic efficiency of the original KE70 design, reflected in both higher k cat values and lower K m values, with the best variants exhibiting k cat / K m values of > 5 × 10 4 s − 1 M − 1 . The optimized KE70 variants were characterized structurally and biochemically, providing insights into the origins of the improvements in catalysis. Three primary contributions were identified: first, the reshaping of the active-site cavity to achieve tighter substrate binding; second, the fine-tuning of electrostatics around the catalytic His-Asp dyad; and, third, the stabilization of the active-site dyad in a conformation optimal for catalysis.

Organophosphate nerve agents are extremely lethal compounds. Rapid in vivo organophosphate clearance requires bioscavenging enzymes with catalytic efficiencies of >10(7) (M(-1) min(-1)). Although serum paraoxonase (PON1) is a leading candidate for such a treatment, it hydrolyzes the toxic S(p) isomers of G-agents with very slow rates. We improved PON1's catalytic efficiency by combining random and targeted mutagenesis with high-throughput screening using fluorogenic analogs in emulsion compartments. We thereby enhanced PON1's activity toward the coumarin analog of S(p)-cyclosarin by ∼10(5)-fold. We also developed a direct screen for protection of acetylcholinesterase from inactivation by nerve agents and used it to isolate variants that degrade the toxic isomer of the coumarin analog and cyclosarin itself with k(cat)/K(M) ∼ 10(7) M(-1) min(-1). We then demonstrated the in vivo prophylactic activity of an evolved variant. These evolved variants and the newly developed screens provide the basis for engineering PON1 for prophylaxis against other G-type agents.

Biological messiness relates to infidelity, heterogeneity, stochastic noise and variation--both genetic and phenotypic--at all levels, from single proteins to organisms. Messiness comes from the complexity and evolutionary history of biological systems and from the high cost of accuracy. For better or for worse, messiness is inherent to biology. It also provides the raw material for physiological and evolutionary adaptations to new challenges.

In the version of this article originally published, the number for the equal contributions footnote was missing for Miriam Kaltenbach and Jason R. Burke in the author list. The error has been corrected in the PDF and print versions of this article.

Summary Upon heterologous overexpression, many proteins misfold or aggregate, thus resulting in low functional yields. Human acetylcholinesterase (hAChE), an enzyme mediating synaptic transmission, is a typical case of a human protein that necessitates mammalian systems to obtain functional expression. We developed a computational strategy and designed an AChE variant bearing 51 mutations that improved core packing, surface polarity, and backbone rigidity. This variant expressed at ∼2,000-fold higher levels in E. coli compared to wild-type hAChE and exhibited 20°C higher thermostability with no change in enzymatic properties or in the active-site configuration as determined by crystallography. To demonstrate broad utility, we similarly designed four other human and bacterial proteins. Testing at most three designs per protein, we obtained enhanced stability and/or higher yields of soluble and active protein in E. coli. Our algorithm requires only a 3D structure and several dozen sequences of naturally occurring homologs, and is available at http://pross.weizmann.ac.il., Graphical Abstract, Highlights • A new computational method is used to stabilize five recalcitrant proteins • Designed variants show higher expression and stability with unmodified function • A designed human acetylcholinesterase variant expresses solubly in bacteria • The method is fully automated and implemented on a webserver, Heterologous expression of proteins and their mutants often results in misfolding and aggregation. Goldenzweig et al. (2016) developed an automated algorithm for protein stabilization requiring minimal experimental testing; for instance, the five tested variants of human acetylcholinesterase showed ≥100-fold higher soluble bacterial expression and higher melting temperatures than wild-type.

Understanding enzyme catalysis through the analysis of natural enzymes is a daunting challenge—their active sites are complex and combine numerous interactions and catalytic forces that are finely coordinated. Study of more rudimentary (wo)man-made enzymes provides a unique opportunity for better understanding of enzymatic catalysis. KE07, a computationally designed Kemp eliminase that employs a glutamate side chain as the catalytic base for the critical proton abstraction step and an apolar binding site to guide substrate binding, was optimized by seven rounds of random mutagenesis and selection, resulting in a > 200-fold increase in catalytic efficiency. Here, we describe the directed evolution process in detail and the biophysical and crystallographic studies of the designed KE07 and its evolved variants. The optimization of KE07's activity to give a k cat / K M value of ∼ 2600 s − 1 M − 1 and an ∼ 10 6 -fold rate acceleration ( k cat / k uncat ) involved the incorporation of up to eight mutations. These mutations led to a marked decrease in the overall thermodynamic stability of the evolved KE07s and in the configurational stability of their active sites. We identified two primary contributions of the mutations to KE07's improved activity: (i) the introduction of new salt bridges to correct a mistake in the original design that placed a lysine for leaving-group protonation without consideration of its “quenching” interactions with the catalytic glutamate, and (ii) the tuning of the environment, the p K a of the catalytic base, and its interactions with the substrate through the evolution of a network of hydrogen bonds consisting of several charged residues surrounding the active site.

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.

Proteins are renowned for their specificity of function. There is, however, accumulating evidence that many proteins, from enzymes to antibodies, are functionally promiscuous. Promiscuity is of considerable physiological importance. In the immune system, cross-reactive or multispecific antibodies are implicated in autoimmune and allergy conditions. In most cases, however, the mechanism behind promiscuity and the relationship between specific and promiscuous activities are unknown. Are the two contradictory? Or can a protein exhibit several unrelated activities each of which is highly specific? To address these questions, we studied a multispecific IgE antibody (SPE7) elicited against a 2,4-dinitrophenyl hapten (DNP). SPE7 is able to distinguish between closely related derivatives such as NP (nitrophenol) and DNP, yet it can also bind a number of unrelated ligands. We find that, like DNP, the cross-reactants are themselves bound specifically-close derivatives of these cross-reactants show very low or no binding to SPE7. It has been suggested that cross-reactivity is simply due to "hydrophobic stickiness", nonspecific interactions between hydrophobic ligands and binding sites. However, partitioning experiments reveal that affinity for SPE7 is unrelated to ligand hydrophobicity. These data, combined with crystal structures of SPE7 in complex with four different ligands, demonstrate that each cross-reactant is bound specifically, forming different hydrogen bonds dependant upon its particular chemistry and the availability of complementary antibody residues. SPE7 is highly homologous to the germline antinitrophenol (NP) antibody B1-8. By comparing the sequences and binding patterns of SPE7 and B1-8, we address the relationship between affinity maturation, specificity, and cross-reactivity.

The directed evolution of proteins has benefited greatly from site-specific methods of diversification such as saturation mutagenesis. These techniques target diversity to a number of chosen positions that are usually non-contiguous in the protein's primary structure. However, the number of targeted positions can be large, thus leading to impractically large library size, wherein almost all library variants are inactive and the likelihood of selecting desirable properties is extremely small. We describe a versatile combinatorial method for the partial diversification of large sets of residues. Our library oligonucleotides comprise randomized codons that are flanked by wild-type sequences. Adding these oligonucleotides to an assembly PCR of wild-type gene fragments incorporates the randomized cassettes, at their target sites, into the reassembled gene. Varying the oligonucleotides concentration resulted in library variants that carry a different average number of mutated positions that comprise a random subset of the entire set of diversified codons. This method, dubbed Incorporating Synthetic Oligos via Gene Reassembly (ISOR), was used to create libraries of a cytosine-C5 methyltransferase wherein 45 individual positions were randomized. One library, containing an average of 5.6 mutated residues per gene, was selected, and mutants with wild-type-like activities isolated. We also created libraries of serum paraoxonase PON1 harboring insertions and deletions (indels) in various areas surrounding the active site. Screening these libraries yielded a range of mutants with altered substrate specificities and indicated that certain regions of this enzyme have a surprisingly high tolerance to indels.

Small libraries for directed evolution can be obtained by neutral drifts that maintain the protein's original function, yielding highly polymorphic, stable and evolvable variants. We describe methods for preparing such libraries, using serum paraoxonase (PON1). An optimized GFP variant fused to PON1 reported levels of soluble, functional enzyme, enabling selection by flow cytometry and identification of enzyme variants exhibiting improved specific and total activities toward several substrates, including toxic organophosphates.

When generating novel tailor-made proteins, protein engineers routinely apply the principles of 'Darwinian' evolution. However, laboratory evolution of proteins also has the potential to test evolutionary theories and reproduce evolutionary scenarios, thus reconstructing putative protein intermediates and providing a glimpse of 'protein fossils'. This commentary describes research at the interface of applied and fundamental molecular evolution, and provides a personal view of how synergy between fundamental and applied experiments indicates novel and more efficient ways of generating new proteins in the laboratory.

This article describes a set of standard control experiments for the authentication of new protein variants isolated through library selection and site-directed mutagenesis. These controls are specifically designed to rule out artifacts derived from 'double transformants' -- i.e. cells transformed with, or infected by, two different plasmids simultaneously. These seem to have been the source of past artifacts and, as demonstrated here, are far more common than generally recognized. By following standard protocols for cloning, plasmid isolation, subcloning, in combination with functional assays, the presence of such artifacts can be ruled out. This protocol needs to be applied for any new variant isolated from heterogeneous gene repertoires, and in particular for variants isolated by selection for either enzymatic activity, or binding.

The modular nature of protein folds suggests that present day proteins evolved via duplication and recombination of smaller functional elements. However, the reconstruction of these putative evolutionary pathways after many millions of years of evolutionary drift has thus far proven difficult, with all attempts to date failing to produce a functional protein. Tachylecin-2 is a monomeric 236 amino acid, five-bladed beta-propeller with five sugar-binding sites. This protein was isolated from a horseshoe crab that emerged ca 500 million years ago. The modular, yet ancient, nature of Tachylectin-2 makes it an excellent model for exploring the evolution of proteins from smaller subunits. To this end, we generated genetically diverse libraries by incremental truncation of the Tachylectin-2 gene and screened them for functional lectins. A number of approximately 100 amino acid residue segments were isolated with the ability to assemble into active homo-pentamers. The topology of most of these segments follows a "hidden" module that differs from the modules observed in wild-type Tachylectin-2, yet their biophysical properties and sugar binding activities resemble the wild-type's. Since the pentamer's molecular mass is twofold higher than the wild-type (approximately 500 amino acid residues), the structure of these oligomeric forms is likely to also differ. Our laboratory evolution experiments highlight the versatility and modularity of the beta-propeller fold, while substantiating the hypothesis that proteins with high internal symmetry, such as beta-propellers, evolved from short, functional gene segments that, at later stages, duplicated, fused, and rearranged, to yield the folds we recognise today.

Derived from the yeast whole-genome duplication, Saccharomyces cerevisiae GAL1 and GAL3 encode the catabolic enzyme galactokinase (Gal1) and its transcriptional coinducer (Gal3), whereas the ancestral, preduplicated GAL1 gene performed both functions. Previous studies indicated that divergence was primarily driven by changes in upstream promoter elements, and changes in GAL3's coding region are assumed to be the result of drift. We show that replacement of GAL3's open-reading-frame with GAL1's results in an extended lag phase upon switching to growth on galactose with up to 2.5-fold differences in the initial cell masses. Accordingly, the binding affinity of Gal3 to Gal80 was found to be greater than 10-folds higher than that of Gal1, with both a higher association rate (ka) and lower dissociation (kd) rate. Thus, while changes in the noncoding, regulatory regions were the initial driving force for GAL3's subfunctionalization as a coinducer, adaptive changes in the protein sequence seem to have followed.