Structural and functional innovations in the real-time evolution of new (βα)8 barrel enzymes

Rational Design of Loop Dynamics for a Barrel-Shaped Enzyme by Introducing Disulfide Bonds

Loop dynamics redesign is an important strategy to manipulate protein function. Cellobiose 2-epimerase (CE) and other members of its superfamily are widely used for diverse industrial applications. The structural feature of the loops connecting barrel helices contributes greatly to the differences in their functional characteristics. Inspired by the in-silico mutation with molecular dynamics (MD) simulation analysis, we propose a strategy for identifying disulfide bond mutation candidates based on the prediction of protein flexibility and residue-residue interaction. The most beneficial mutant with the newly introduced disulfide bond would simultaneously improve both its thermostability and its reaction propensity to the targeting isomerization product. The ratio of the isomerization/epimerization catalytic rate was improved from 4:103 to 9:22. MD simulation and binding free energy calculations were applied to provide insights into molecular recognition upon mutations. The comparative analysis of enzyme/substrate binding modes indicates that the altered catalytic reaction pathway is due to less efficient binding of the native product. The key residue responsible for the observed phenotype was identified by energy decomposition and was further confirmed by the mutation experiment. The rational design of the key loop region might be a promising strategy to alter the catalytic behavior of all (α/α)6-barrel-like proteins.

Catalytic Redundancies and Conformational Plasticity Drives Selectivity and Promiscuity in Quorum Quenching Lactonases

Several enzymes from the metallo-β-lactamase-like family of lactonases (MLLs) degrade N- acyl-L-homoserine lactones (AHLs). In doing so, they play a role in a microbial communication system, quorum sensing, which contributes to pathogenicity and biofilm formation. There is currently great interest in designing quorum quenching ( QQ ) enzymes that can interfere with this communication and be used in a range of industrial and biomedical applications. However, tailoring these enzymes for specific targets requires a thorough understanding of their mechanisms and the physicochemical properties that determine their substrate specificities. We present here a detailed biochemical, computational, and structural study of the MLL GcL, which is highly proficient, thermostable, and has broad substrate specificity. Strikingly, we show that GcL does not only accept a broad range of substrates but is also capable of utilizing different reaction mechanisms that are differentially used in function of the substrate structure or the remodeling of the active site via mutations. Comparison of GcL to other lactonases such as AiiA and AaL demonstrates similar mechanistic promiscuity, suggesting this is a shared feature across lactonases in this enzyme family. Mechanistic promiscuity has previously been observed in the lactonase/paraoxonase PON1, as well as with protein tyrosine phosphatases that operate via a dual general-acid mechanism. The apparent prevalence of this phenomenon is significant from both a biochemical and an engineering perspective: in addition to optimizing for specific substrates, it is possible to optimize for specific mechanisms, opening new doors not just for the design of novel quorum quenching enzymes, but also of other mechanistically promiscuous enzymes.

Back in time to the Gly-rich prototype of the phosphate binding elementary function

Binding of nucleotides and their derivatives is one of the most ancient elementary functions dating back to the Origin of Life. We review here the works considering one of the key elements in binding of (di)nucleotide-containing ligands - phosphate binding. We start from a brief discussion of major participants, conditions, and events in prebiotic evolution that resulted in the Origin of Life. Tracing back to the basic functions, including metal and phosphate binding, and, potentially, formation of primitive protein-protein interactions, we focus here on the phosphate binding. Critically assessing works on the structural, functional, and evolutionary aspects of phosphate binding, we perform a simple computational experiment reconstructing its most ancient and generic sequence prototype. The profiles of the phosphate binding signatures have been derived in form of position-specific scoring matrices (PSSMs), their peculiarities depending on the type of the ligands have been analyzed, and evolutionary connections between them have been delineated. Then, the apparent prototype that gave rise to all relevant phosphate-binding signatures had also been reconstructed. We show that two major signatures of the phosphate binding that discriminate between the binding of dinucleotide- and nucleotide-containing ligands are GxGxxG and GxxGxG, respectively. It appears that the signature archetypal for dinucleotide-containing ligands is more generic, and it can frequently bind phosphate groups in nucleotide-containing ligands as well. The reconstructed prototype's key signature GxGGxG underlies the role of glycine residues in providing flexibility and interactions necessary for binding the phosphate groups. The prototype also contains other ancient amino acids, valine, and alanine, showing versatility towards evolutionary design and functional diversification.

Retracing the Rapid Evolution of an Herbicide-Degrading Enzyme by Protein Engineering

The mechanisms underlying the rapid evolution of novel enzymatic activities from promiscuous side activities are poorly understood. Recently emerged enzymes catalyzing the catabolic degradation of xenobiotic substances that have been spread out into the environment during the last decades provide an exquisite opportunity to study these mechanisms. A prominent example is the herbicide atrazine (2-chloro-4-ethylamino-6-isopropylamino-1,3,5-triazine), which is degraded through a number of enzymatic reactions constituting the Atz pathway. Here, we analyzed the evolution of the hydroxyatrazine ethylaminohydrolase AtzB, a Zn(II)-dependent metalloenzyme that adopts the amidohydrolase fold and catalyzes the second step of the Atz pathway. We searched for promiscuous side activities of AtzB, which might point to the identity of its progenitor. These investigations revealed that AtzB has low promiscuous guanine deaminase activity. Furthermore, we found that the two closest AtzB homologues, which have not been functionally annotated up to now, are guanine deaminases with modest promiscuous hydroxyatrazine hydrolase activity. Based on sequence comparisons with the closest AtzB homologues, the guanine deaminase activity of AtzB could be increased by three orders of magnitude through the introduction of only four active site mutations. Interestingly, introducing the inverse four mutations into the AtzB homologues significantly enhanced their hydroxyatrazine hydrolase activity, and in one case is even equivalent to that of wild-type AtzB. Molecular dynamics simulations elucidated the structural and molecular basis for the mutation-induced activity changes. The example of AtzB highlights how novel enzymes with high catalytic proficiency can evolve from low promiscuous side activities by only few mutational events within a short period of time.

Loop dynamics and the evolution of enzyme activity

In the early 2000s, Tawfik presented his 'New View' on enzyme evolution, highlighting the role of conformational plasticity in expanding the functional diversity of limited repertoires of sequences. This view is gaining increasing traction with increasing evidence of the importance of conformational dynamics in both natural and laboratory evolution of enzymes. The past years have seen several elegant examples of harnessing conformational (particularly loop) dynamics to successfully manipulate protein function. This Review revisits flexible loops as critical participants in regulating enzyme activity. We showcase several systems of particular interest: triosephosphate isomerase barrel proteins, protein tyrosine phosphatases and β-lactamases, while briefly discussing other systems in which loop dynamics are important for selectivity and turnover. We then discuss the implications for engineering, presenting examples of successful loop manipulation in either improving catalytic efficiency, or changing selectivity completely. Overall, it is becoming clearer that mimicking nature by manipulating the conformational dynamics of key protein loops is a powerful method of tailoring enzyme activity, without needing to target active-site residues.

Identification and Mutation Analysis of Nonconserved Residues on the TIM-Barrel Surface of GH5_5 Cellulases for Catalytic Efficiency and Stability Improvement

Exploring the potential functions of nonconserved residues on the outer side of α-helices and systematically optimizing them are pivotal for their application in protein engineering. Based on the evolutionary structural conservation analysis of GH5_5 cellulases, a practical molecular improvement strategy was developed. Highly variable sites on the outer side of the α-helices of the GH5_5 cellulase from Aspergillus niger (AnCel5A) were screened, and 14 out of the 34 highly variable sites were confirmed to exert a positive effect on the activity. After the modular combination of the positive mutations, the catalytic efficiency of the mutants was further improved. By using CMC-Na as the substrate, the catalytic efficiency and specific activity of variant AnCel5A_N193A/T300P/D307P were approximately 2.0-fold that of AnCel5A (227 ± 21 versus 451 ± 43 ml/s/mg and 1,726 ± 19 versus 3,472 ± 42 U/mg, respectively). The half-life (t_1/2) of variant AnCel5A_N193A/T300P/D307P at 75°C was 2.36 times that of AnCel5A. The role of these sites was successfully validated in other GH5_5 cellulases. Computational analyses revealed that the flexibility of the loop 6-loop 7-loop 8 region was responsible for the increased catalytic performance. This work not only illustrated the important role of rapidly evolving positions on the outer side of the α-helices of GH5_5 cellulases but also revealed new insights into engineering the proteins that nature left as clues for us to find. IMPORTANCE A comprehensive understanding of the residues on the α-helices of the GH5_5 cellulases is important for catalytic efficiency and stability improvement. The main objective of this study was to use the evolutionary conservation and plasticity of the TIM-barrel fold to probe the relationship between nonconserved residues on the outer side of the α-helices and the catalytic efficiency of GH5_5 cellulases by conducting structure-guided protein engineering. By using a four-step nonconserved residue screening strategy, the functional role of nonconserved residues on the outer side of the α-helices was effectively identified, and a variant with superior performance and capability was constructed. Hence, this study proved the effectiveness of this strategy in engineering GH5_5 cellulases and provided a potential competitor for industrial applications. Furthermore, this study sheds new light on engineering TIM-barrel proteins.

Uniform binding and negative catalysis at the origin of enzymes

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.

Enabling Role of Ligand-Driven Conformational Changes in Enzyme Evolution

Many enzymes that show a large specificity in binding the enzymatic transition state with a higher affinity than the substrate utilize substrate binding energy to drive protein conformational changes to form caged substrate complexes. These protein cages provide strong stabilization of enzymatic transition states. Using part of the substrate binding energy to drive the protein conformational change avoids a similar strong stabilization of the Michaelis complex and irreversible ligand binding. A seminal step in the development of modern enzyme catalysts was the evolution of enzymes that couple substrate binding to a conformational change. These include enzymes that function in glycolysis (triosephosphate isomerase), the biosynthesis of lipids (glycerol phosphate dehydrogenase), the hexose monophosphate shunt (6-phosphogluconate dehydrogenase), and the mevalonate pathway (isopentenyl diphosphate isomerase), catalyze the final step in the biosynthesis of pyrimidine nucleotides (orotidine monophosphate decarboxylase), and regulate the cellular levels of adenine nucleotides (adenylate kinase). The evolution of enzymes that undergo ligand-driven conformational changes to form active protein-substrate cages is proposed to proceed by selection of variants, in which the selected side chain substitutions destabilize a second protein conformer that shows compensating enhanced binding interactions with the substrate. The advantages inherent to enzymes that incorporate a conformational change into the catalytic cycle provide a strong driving force for the evolution of flexible protein folds such as the TIM barrel. The appearance of these folds represented a watershed event in enzyme evolution that enabled the rapid propagation of enzyme activities within enzyme superfamilies.

Complex Loop Dynamics Underpin Activity, Specificity, and Evolvability in the (βα)8 Barrel Enzymes of Histidine and Tryptophan Biosynthesis

Enzymes are conformationally dynamic, and their dynamical properties play an important role in regulating their specificity and evolvability. In this context, substantial attention has been paid to the role of ligand-gated conformational changes in enzyme catalysis; however, such studies have focused on tremendously proficient enzymes such as triosephosphate isomerase and orotidine 5'-monophosphate decarboxylase, where the rapid (μs timescale) motion of a single loop dominates the transition between catalytically inactive and active conformations. In contrast, the (βα)₈-barrels of tryptophan and histidine biosynthesis, such as the specialist isomerase enzymes HisA and TrpF, and the bifunctional isomerase PriA, are decorated by multiple long loops that undergo conformational transitions on the ms (or slower) timescale. Studying the interdependent motions of multiple slow loops, and their role in catalysis, poses a significant computational challenge. This work combines conventional and enhanced molecular dynamics simulations with empirical valence bond simulations to provide rich details of the conformational behavior of the catalytic loops in HisA, PriA, and TrpF, and the role of their plasticity in facilitating bifunctionality in PriA and evolved HisA variants. In addition, we demonstrate that, similar to other enzymes activated by ligand-gated conformational changes, loops 3 and 4 of HisA and PriA act as gripper loops, facilitating the isomerization of the large bulky substrate ProFAR, albeit now on much slower timescales. This hints at convergent evolution on these different (βα)₈-barrel scaffolds. Finally, our work reemphasizes the potential of engineering loop dynamics as a tool to artificially manipulate the catalytic repertoire of TIM-barrel proteins.

Conservation and diversity in allosteric fingerprints of proteins for evolutionary-inspired engineering and design

Hand-in-hand work of physics and evolution delivered protein universe with diversity of forms, sizes, and functions. Pervasiveness and advantageous traits of allostery made it an important component of the protein function regulation, calling for thorough investigation of its structural determinants and evolution. Learning directly from nature, we explored here allosteric communication in several major folds and repeat proteins, including α/β and β-barrels, β-propellers, Ig-like fold, ankyrin and α/β leucine-rich repeat proteins, which provide structural platforms for many different enzymatic and signalling functions. We obtained a picture of conserved allosteric communication characteristic in different fold types, modifications of the structure-driven signalling patterns via sequence-determined divergence to specific functions, as well as emergence and potential diversification of allosteric regulation in multi-domain proteins and oligomeric assemblies. Our observations will be instrumental in facilitating the engineering and de novo design of proteins with allosterically regulated functions, including development of therapeutic biologics. In particular, results described here may guide the identification of the optimal structural platforms (e.g. fold type, size, and oligomerization states) and the types of diversifications/perturbations, such as mutations, effector binding, and order-disorder transition. The tunable allosteric linkage across distant regions can be used as a pivotal component in the design/engineering of modular biological systems beyond the traditional scaffolding function.

Simulations reveal the key role of Arg15 in the promiscuous activity in the HisA enzyme

The HisA enzyme catalyzes the first step of histidine biosynthesis via the Amadori rearrangement of the substrate ProFAR. Since it possesses the most conserved and ancient TIM-barrel fold, it provides an ideal framework for bioengineering of a new function from ancestral enzymes. In the present study, first, the catalytic mechanism of HisA biosynthesis was elucidated using hybrid Quantum Mechanical/Molecular Mechanical calculations, and thereafter, key residues contributing towards the promiscuity for TrpF activity were revealed using several MD simulations of a wild type enzyme and its variant with the native (ProFAR) and promiscuous (PRA) substrates. Our study reveals that the two loops (βα)₁ and (βα)₅ on the catalytic site of the HisA enzyme have incredible adaptability for the native and promiscuous substrates. The conformational interplay between these two loops is substrate driven and precise bioengineering targeting these loops is key to the emergence of new functions. Furthermore, the study reveals a key role of the Arg 15 residue which is close to the catalytic center of the enzyme in the bifunctionality of the HisA enzyme by increasing the loop flexibility. Therefore, our study provides crucial information for future bioengineering work to use the HisA enzyme as a scaffold for new enzymatic activity.

Evolution, folding, and design of TIM barrels and related proteins

Proteins are chief actors in life that perform a myriad of exquisite functions. This diversity has been enabled through the evolution and diversification of protein folds. Analysis of sequences and structures strongly suggest that numerous protein pieces have been reused as building blocks and propagated to many modern folds. This information can be traced to understand how the protein world has diversified. In this review, we discuss the latest advances in the analysis of protein evolutionary units, and we use as a model system one of the most abundant and versatile topologies, the TIM-barrel fold, to highlight the existing common principles that interconnect protein evolution, structure, folding, function, and design.

Improving the activity and thermostability of GH2 β-glucuronidases via domain reassembly

Glycoside hydrolase family 2 (GH2) enzymes are generally composed of three domains: TIM-barrel domain (TIM), immunoglobulin-like β-sandwich domain (ISD), and sugar-binding domain (SBD). The combination of these three domains yields multiple structural combinations with different properties. Theoretically, the drawbacks of a given GH2 fold may be circumvented by efficiently reassembling the three domains. However, very few successful cases have been reported. In this study, we used six GH2 β-glucuronidases (GUSs) from bacteria, fungi, or humans as model enzymes and constructed a series of mutants by reassembling the domains from different GUSs. The mutants PGUS-At, GUS-PAA, and GUS-PAP, with reassembled domains from fungal GUSs, showed improved expression levels, activity, and thermostability, respectively. Specifically, compared to the parental enzyme, the mutant PGUS-At displayed 3.8 times higher expression, the mutant GUS-PAA displayed 1.0 time higher catalytic efficiency (k_cat /K_m ), and the mutant GUS-PAP displayed 7.5 times higher thermostability at 65°C. Furthermore, two-hybrid mutants, GUS-AEA and GUS-PEP, were constructed with the ISD from a bacterial GUS and SBD and TIM domain from fungal GUSs. GUS-AEA and GUS-PEP showed 30.4% and 23.0% higher thermostability than GUS-PAP, respectively. Finally, molecular dynamics simulations were conducted to uncover the molecular reasons for the increased thermostability of the mutant.

Evolution of new enzymes by gene duplication and divergence

Thousands of new metabolic and regulatory enzymes have evolved by gene duplication and divergence since the dawn of life. New enzyme activities often originate from promiscuous secondary activities that have become important for fitness due to a change in the environment or a mutation. Mutations that make a promiscuous activity physiologically relevant can occur in the gene encoding the promiscuous enzyme itself, but can also occur elsewhere, resulting in increased expression of the enzyme or decreased competition between the native and novel substrates for the active site. If a newly useful activity is inefficient, gene duplication/amplification will set the stage for divergence of a new enzyme. Even a few mutations can increase the efficiency of a new activity by orders of magnitude. As efficiency increases, amplified gene arrays will shrink to provide two alleles, one encoding the original enzyme and one encoding the new enzyme. Ultimately, genomic rearrangements eliminate co-amplified genes and move newly evolved paralogs to a distant region of the genome.

Structure and kinetics of indole-3-glycerol phosphate synthase from Pseudomonas aeruginosa: Decarboxylation is not essential for indole formation

In tryptophan biosynthesis, the reaction catalyzed by the enzyme indole-3-glycerol phosphate synthase (IGPS) starts with a condensation step in which the substrate's carboxylated phenyl group makes a nucleophilic attack to form the pyrrole ring of the indole, followed by a decarboxylation that restores the aromaticity of the phenyl. IGPS from Pseudomonas aeruginosa has the highest turnover number of all characterized IGPS enzymes, providing an excellent model system to test the necessity of the decarboxylation step. Since the 1960s, this step has been considered to be mechanistically essential based on studies of the IGPS–phosphoribosylanthranilate isomerase fusion protein from Escherichia coli. Here, we present the crystal structure of P. aeruginosa IGPS in complex with reduced CdRP, a nonreactive substrate analog, and using a sensitive discontinuous assay, we demonstrate weak promiscuous activity on the decarboxylated substrate 1-(phenylamino)-1-deoxyribulose-5-phosphate, with an ∼1000× lower rate of IGP formation than from the native substrate. We also show that E. coli IGPS, at an even lower rate, can produce IGP from decarboxylated substrate. Our structure of P. aeruginosa IGPS has eight molecules in the asymmetric unit, of which seven contain ligand and one displays a previously unobserved conformation closer to the reactive state. One of the few nonconserved active-site residues, Phe²⁰¹ in P. aeruginosa IGPS, is by mutagenesis demonstrated to be important for the higher turnover of this enzyme on both substrates. Our results demonstrate that despite IGPS's classification as a carboxy-lyase (i.e. decarboxylase), decarboxylation is not a completely essential step in its catalysis.

Mutational Pathways and Trade-Offs Between HisA and TrpF Functions: Implications for Evolution via Gene Duplication and Divergence

When a new activity evolves by changes in a pre-existing enzyme this is likely to reduce the original activity, generating a functional trade-off. The properties of this trade-off will affect the continued evolution of both functions. If the trade-off is strong, gene duplication and subsequent divergence would be favored whereas if the trade-off is weak a bi-functional enzyme could evolve that performs both functions. We previously showed that when a bi-functional HisA enzyme was evolved under selection for both HisA and TrpF functions, evolution mainly proceeded via duplication-divergence and specialization, implying that the trade-off is strong between these two functions. Here, we examined this hypothesis by identifying the mutational pathways (i.e., the mutational landscape) in the Salmonella enterica HisA enzyme that conferred a TrpF-like activity, and examining the trade-offs between the original and new activity. For the HisA enzyme there are many different paths toward the new TrpF function, each with its own unique trade-off. A total of 16 single mutations resulted in HisA enzyme variants that acquired TrpF activity and only three of them maintained HisA activity. Twelve mutants were evolved further toward increased TrpF activity and during evolution toward improved TrpF activity the original HisA activity was completely lost in all lineages. We propose that, aside from various relevant ecological factors, two main genetic factors influence whether evolution of a new function proceeds via duplication – divergence (specialization) or by evolution of a generalist: (i) the relative mutation supply of the two pathways and (ii) the shape of the trade-off curve between the native and new function.

A mechanistic view of enzyme evolution

New enzyme functions often evolve through the recruitment and optimization of latent promiscuous activities. How do mutations alter the molecular architecture of enzymes to enhance their activities? Can we infer general mechanisms that are common to most enzymes, or does each enzyme require a unique optimization process? The ability to predict the location and type of mutations necessary to enhance an enzyme's activity is critical to protein engineering and rational design. In this review, via the detailed examination of recent studies that have shed new light on the molecular changes underlying the optimization of enzyme function, we provide a mechanistic perspective of enzyme evolution. We first present a global survey of the prevalence of activity-enhancing mutations and their distribution within protein structures. We then delve into the molecular solutions that mediate functional optimization, specifically highlighting several common mechanisms that have been observed across multiple examples. As distinct protein sequences encounter different evolutionary bottlenecks, different mechanisms are likely to emerge along evolutionary trajectories toward improved function. Identifying the specific mechanism(s) that need to be improved upon, and tailoring our engineering efforts to each sequence, may considerably improve our chances to succeed in generating highly efficient catalysts in the future.

Collateral fitness effects of mutations

The distribution of fitness effects of mutation plays a central role in constraining protein evolution. The underlying mechanisms by which mutations lead to fitness effects are typically attributed to changes in protein specific activity or abundance. Here, we reveal the importance of a mutation's collateral fitness effects, which we define as effects that do not derive from changes in the protein's ability to perform its physiological function. We comprehensively measured the collateral fitness effects of missense mutations in the Escherichia coli TEM-1 β-lactamase antibiotic resistance gene using growth competition experiments in the absence of antibiotic. At least 42% of missense mutations in TEM-1 were deleterious, indicating that for some proteins collateral fitness effects occur as frequently as effects on protein activity and abundance. Deleterious mutations caused improper posttranslational processing, incorrect disulfide-bond formation, protein aggregation, changes in gene expression, and pleiotropic effects on cell phenotype. Deleterious collateral fitness effects occurred more frequently in TEM-1 than deleterious effects on antibiotic resistance in environments with low concentrations of the antibiotic. The surprising prevalence of deleterious collateral fitness effects suggests they may play a role in constraining protein evolution, particularly for highly expressed proteins, for proteins under intermittent selection for their physiological function, and for proteins whose contribution to fitness is buffered against deleterious effects on protein activity and protein abundance.

How enzyme promiscuity and horizontal gene transfer contribute to metabolic innovation

Promiscuity is the coincidental ability of an enzyme to catalyze its native reaction and additional reactions that are not biological functions in the same active site. Promiscuity plays a central role in enzyme evolution and is thus a useful property for protein and metabolic engineering. This review examines enzyme evolution holistically, beginning with evaluating biochemical support for four enzyme evolution models. As expected, there is strong biochemical support for the subfunctionalization and innovation-amplification-divergence models, in which promiscuity is a central feature. In many cases, however, enzyme evolution is more complex than the models indicate, suggesting much is yet to be learned about selective pressures on enzyme function. A complete understanding of enzyme evolution must also explain the ability of metabolic networks to integrate new enzyme activities. Hidden within metabolic networks are underground metabolic pathways constructed from promiscuous activities. We discuss efforts to determine the diversity and pervasiveness of underground metabolism. Remarkably, several studies have discovered that some metabolic defects can be repaired via multiple underground routes. In prokaryotes, metabolic innovation is driven by connecting enzymes acquired by horizontal gene transfer (HGT) into the metabolic network. Thus, we end the review by discussing how the combination of promiscuity and HGT contribute to evolution of metabolism in prokaryotes. Future studies investigating the contribution of promiscuity to enzyme and metabolic evolution will need to integrate deeper probes into the influence of evolution on protein biophysics, enzymology, and metabolism with more complex and realistic evolutionary models. ENZYMES: lactate dehydrogenase (EC 1.1.1.27), malate dehydrogenase (EC 1.1.1.37), OSBS (EC 4.2.1.113), HisA (EC 5.3.1.16), TrpF, PriA (EC 5.3.1.24), R-mandelonitrile lyase (EC 4.1.2.10), Maleylacetate reductase (EC 1.3.1.32).

Mutations that improve efficiency of a weak-link enzyme are rare compared to adaptive mutations elsewhere in the genome

New enzymes often evolve by gene amplification and divergence. Previous experimental studies have followed the evolutionary trajectory of an amplified gene, but have not considered mutations elsewhere in the genome when fitness is limited by an evolving gene. We have evolved a strain of Escherichia coli in which a secondary promiscuous activity has been recruited to serve an essential function. The gene encoding the ‘weak-link’ enzyme amplified in all eight populations, but mutations improving the newly needed activity occurred in only one. Most adaptive mutations occurred elsewhere in the genome. Some mutations increase expression of the enzyme upstream of the weak-link enzyme, pushing material through the dysfunctional metabolic pathway. Others enhance production of a co-substrate for a downstream enzyme, thereby pulling material through the pathway. Most of these latter mutations are detrimental in wild-type E. coli, and thus would require reversion or compensation once a sufficient new activity has evolved.

In Vivo Titration of Folate Pathway Enzymes

How enzymes behave in cells is likely different from how they behave in the test tube. Previous in vitro studies find that osmolytes interact weakly with folate. Removal of the osmolyte from the solvation shell of folate is more difficult than removal of water, which weakens binding of folate to its enzyme partners. To examine if this phenomenon occurs in vivo, osmotic stress titrations were performed with Escherichia coli Two strategies were employed: resistance to an antibacterial drug and complementation of a knockout strain by the appropriate gene cloned into a plasmid that allows tight control of expression levels as well as labeling by a degradation tag. The abilities of the knockout and complemented strains to grow under osmotic stress were compared. Typically, the knockout strain could grow to high osmolalities on supplemented medium, while the complemented strain stopped growing at lower osmolalities on minimal medium. This pattern was observed for an R67 dihydrofolate reductase clone rescuing a ΔfolA strain, for a methylenetetrahydrofolate reductase clone rescuing a ΔmetF strain, and for a serine hydroxymethyltransferase clone rescuing a ΔglyA strain. Additionally, an R67 dihydrofolate reductase clone allowed E. coli DH5α to grow in the presence of trimethoprim until an osmolality of ∼0.81 is reached, while cells in a control titration lacking antibiotic could grow to 1.90 osmol.IMPORTANCEE. coli can survive in drought and flooding conditions and can tolerate large changes in osmolality. However, the cell processes that limit bacterial growth under high osmotic stress conditions are not known. In this study, the dose of four different enzymes in E. coli was decreased by using deletion strains complemented by the gene carried in a tunable plasmid. Under conditions of limiting enzyme concentration (lower than that achieved by chromosomal gene expression), cell growth can be blocked by osmotic stress conditions that are normally tolerated. These observations indicate that E. coli has evolved to deal with variations in its osmotic environment and that normal protein levels are sufficient to buffer the cell from environmental changes. Additional factors involved in the osmotic pressure response may include altered protein concentration/activity levels, weak solute interactions with ligands which can make it more difficult for proteins to bind their substrates/inhibitors/cofactors in vivo, and/or viscosity effects.

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.

Metabolite–Enzyme Coevolution: From Single Enzymes to Metabolic Pathways and Networks

How individual enzymes evolved is relatively well understood. However, individual enzymes rarely confer a physiological advantage on their own. Judging by its current state, the emergence of metabolism seemingly demanded the simultaneous emergence of many enzymes. Indeed, how multicomponent interlocked systems, like metabolic pathways, evolved is largely an open question. This complexity can be unlocked if we assume that survival of the fittest applies not only to genes and enzymes but also to the metabolites they produce. This review develops our current knowledge of enzyme evolution into a wider hypothesis of pathway and network evolution. We describe the current models for pathway evolution and offer an integrative metabolite–enzyme coevolution hypothesis. Our hypothesis addresses the origins of new metabolites and of new enzymes and the order of their recruitment. We aim to not only survey established knowledge but also present open questions and potential ways of addressing them.