A compromise required by gene sharing enables survival: Implications for evolution of new enzyme activities

Transcriptomic and proteomic ramifications of segmental amplification in Escherichia coli

Gene amplification can drive adaptation by rapidly increasing the cellular dosage of critical gene products. Segmental amplifications often encompass large genomic regions surrounding the gene(s) under selection for higher dosage. Overexpression of coamplified neighboring genes imposes a substantial metabolic burden. While compensatory mutations can decrease inappropriate overexpression of coamplified genes, it takes time for such mutations to arise. The extent to which intrinsic regulatory mechanisms modulate expression of coamplified genes in the immediate aftermath of segmental amplification is largely unknown. To address the collateral effects of segmental amplification, we evolved replicate cultures of an Escherichia coli mutant under conditions that select for higher dosage of an inefficient enzyme whose weak activity limits growth rate. Segmental amplifications encompassing the gene encoding the weak-link enzyme arose in all populations. Amplified regions ranged in size (33 to 125 kb) and copy number (2 to ≥14 copies). We performed RNA-seq and label-free proteomics to quantify expression of amplified genes present at 2, 6, and 14 copies. mRNA expression generally scales with gene copy number, but protein expression scales less well with both gene copy number and mRNA expression. We characterize the molecular mechanisms underlying discrepancies between gene copy number and expression for several cases. We also show that segmental amplifications can have system-wide consequences by indirectly altering expression of nonamplified genes. Our findings indicate that the fitness benefit derived from segmental amplification depends on the combined effects of amplicon size, gene content, and copy number as well as collateral effects on nonamplified genes.

Intramolecular epistasis correlates with divergence of specificity in promiscuous and bifunctional NSAR/OSBS enzymes

Understanding the functions and evolution of specificity-determining residues is essential for improving strategies to predict and design enzyme functions. Whether the function of an amino acid residue is retained during evolution depends on intramolecular epistasis, which occurs when the same residue contributes to different phenotypes in different genetic backgrounds. This study examines the relationship between epistasis and functional divergence by investigating a conserved specificity determinant in five homologs from the N-succinylamino acid racemase (NSAR)/o-succinylbenzoate synthase (OSBS) subfamily. NSAR activity originated as a promiscuous (non-biological) activity of an ancestral OSBS. Some extant NSAR/OSBS subfamily enzymes still have OSBS activity as a biological function and NSAR as a promiscuous activity, while some use both OSBS and NSAR activities as biological functions. Others use only NSAR activity as a biological function but can still catalyze the OSBS reaction as a promiscuous activity. Previously, we determined that the conserved residue R266 in Amycolatopsis sp. T-1-60 NSAR contributes to NSAR specificity by enabling K263 to act as a general acid/base catalyst. Here, we show that mutating R266 decreased relative specificity for NSAR activity in four of five NSAR/OSBS subfamily enzymes, as predicted. However, other phenotypes exhibited epistasis related to the pleiotropy of R266, including the proton exchange rate between the catalytic lysines and the substrate, the impact on OSBS activity, and thermostability. The strength of epistasis was associated with functional and evolutionary divergence of NSAR/OSBS enzymes. These results illustrate the benefits of comparing multiple homologs for understanding mechanisms of enzyme specificity.

Frustration can Limit the Adaptation of Promiscuous Enzymes Through Gene Duplication and Specialisation

Virtually all enzymes catalyse more than one reaction, a phenomenon known as enzyme promiscuity. It is unclear whether promiscuous enzymes are more often generalists that catalyse multiple reactions at similar rates or specialists that catalyse one reaction much more efficiently than other reactions. In addition, the factors that shape whether an enzyme evolves to be a generalist or a specialist are poorly understood. To address these questions, we follow a three-pronged approach. First, we examine the distribution of promiscuity in empirical enzymes reported in the BRENDA database. We find that the promiscuity distribution of empirical enzymes is bimodal. In other words, a large fraction of promiscuous enzymes are either generalists or specialists, with few intermediates. Second, we demonstrate that enzyme biophysics is not sufficient to explain this bimodal distribution. Third, we devise a constraint-based model of promiscuous enzymes undergoing duplication and facing selection pressures favouring subfunctionalization. The model posits the existence of constraints between the catalytic efficiencies of an enzyme for different reactions and is inspired by empirical case studies. The promiscuity distribution predicted by our constraint-based model is consistent with the empirical bimodal distribution. Our results suggest that subfunctionalization is possible and beneficial only in certain enzymes. Furthermore, the model predicts that conflicting constraints and selection pressures can cause promiscuous enzymes to enter a 'frustrated' state, in which competing interactions limit the specialisation of enzymes. We find that frustration can be both a driver and an inhibitor of enzyme evolution by duplication and subfunctionalization. In addition, our model predicts that frustration becomes more likely as enzymes catalyse more reactions, implying that natural selection may prefer catalytically simple enzymes. In sum, our results suggest that frustration may play an important role in enzyme evolution.

Cell Survival Enabled by Leakage of a Labile Metabolic Intermediate

Many metabolites are generated in one step of a biochemical pathway and consumed in a subsequent step. Such metabolic intermediates are often reactive molecules which, if allowed to freely diffuse in the intracellular milieu, could lead to undesirable side reactions and even become toxic to the cell. Therefore, metabolic intermediates are often protected as protein-bound species and directly transferred between enzyme active sites in multi-functional enzymes, multi-enzyme complexes, and metabolons. Sequestration of reactive metabolic intermediates thus contributes to metabolic efficiency. It is not known, however, whether this evolutionary adaptation can be relaxed in response to challenges to organismal survival. Here, we report evolutionary repair experiments on Escherichia coli cells in which an enzyme crucial for the biosynthesis of proline has been deleted. The deletion makes cells unable to grow in a culture medium lacking proline. Remarkably, however, cell growth is efficiently restored by many single mutations (12 at least) in the gene of glutamine synthetase. The mutations cause the leakage to the intracellular milieu of a highly reactive phosphorylated intermediate common to the biosynthetic pathways of glutamine and proline. This intermediate is generally assumed to exist only as a protein-bound species. Nevertheless, its diffusion upon mutation-induced leakage enables a new route to proline biosynthesis. Our results support that leakage of sequestered metabolic intermediates can readily occur and contribute to organismal adaptation in some scenarios. Enhanced availability of reactive molecules may enable the generation of new biochemical pathways and the potential of mutation-induced leakage in metabolic engineering is noted.

How to Recruit a Promiscuous Enzyme to Serve a New Function

Promiscuous enzymes can be recruited to serve new functions when a genetic or environmental change makes catalysis of a novel reaction important for fitness or even survival. Subsequently, gene duplication and divergence can lead to evolution of an efficient and specialized new enzyme. Every organism likely has thousands of promiscuous enzyme activities that provide a vast reservoir of catalytic potential. However, much of this potential may not be accessible. We compiled kinetic parameters for promiscuous reactions catalyzed by 108 enzymes. The median value of k_cat/K_M is a very modest 31 M^-1 s^-1. Based upon the fluxes through metabolic pathways in E. coli, we estimate that many, if not most, promiscuous activities are too inefficient to impact fitness. However, mutations can elevate the level of an insufficient promiscuous activity by increasing enzyme expression, improving k_cat/K_M, or altering concentrations of the promiscuous and native substrates and allosteric regulators. Particularly in large bacterial populations, stochastic mutations may provide a viable pathway for recruitment of even inefficient promiscuous activities.

Characterization of a novel β-alanine biosynthetic pathway consisting of promiscuous metabolic enzymes

Bacteria adapt to utilize the nutrients available in their environment through a sophisticated metabolic system composed of highly specialized enzymes. Although these enzymes can metabolize molecules other than those for which they evolved, their efficiency toward promiscuous substrates is considered too low to be of physiological relevance. Herein, we investigated the possibility that these promiscuous enzymes are actually efficient enough at metabolizing secondary substrates to modify the phenotype of the cell. For example, in the bacterium Acinetobacter baylyi ADP1 (ADP1), panD (coding for l-aspartate decarboxylase) encodes the only protein known to catalyze the synthesis of β-alanine, an obligate intermediate in CoA synthesis. However, we show that the ADP1 ΔpanD mutant could also form this molecule through an unknown metabolic pathway arising from promiscuous enzymes and grow as efficiently as the wildtype strain. Using metabolomic analyses, we identified 1,3-diaminopropane and 3-aminopropanal as intermediates in this novel pathway. We also conducted activity screening and enzyme kinetics to elucidate candidate enzymes involved in this pathway, including 2,4-diaminobutyrate aminotransferase (Dat) and 2,4-diaminobutyrate decarboxylase (Ddc) and validated this pathway in vivo by analyzing the phenotype of mutant bacterial strains. Finally, we experimentally demonstrate that this novel metabolic route is not restricted to ADP1. We propose that the occurrence of conserved genes in hundreds of genomes across many phyla suggests that this previously undescribed pathway is widespread in prokaryotes.

Substrate Dynamics Contribute to Enzymatic Specificity in Human and Bacterial Methionine Adenosyltransferases

Protein conformational changes can facilitate the binding of noncognate substrates and underlying promiscuous activities. However, the contribution of substrate conformational dynamics to this process is comparatively poorly understood. Here, we analyze human (hMAT2A) and Escherichia coli (eMAT) methionine adenosyltransferases that have identical active sites but different substrate specificity. In the promiscuous hMAT2A, noncognate substrates bind in a stable conformation to allow catalysis. In contrast, noncognate substrates sample stable productive binding modes less frequently in eMAT owing to altered mobility in the enzyme active site. Different cellular concentrations of substrates likely drove the evolutionary divergence of substrate specificity in these orthologues. The observation of catalytic promiscuity in hMAT2A led to the detection of a new human metabolite, methyl thioguanosine, that is produced at elevated levels in a cancer cell line. This work establishes that identical active sites can result in different substrate specificity owing to the effects of substrate and enzyme dynamics.

Promiscuous activity of 3-isopropylmalate dehydrogenase produced at physiological level affords Escherichia coli growth on d-malate

Promiscuous activities of enzymes may serve as starting points for the evolution of new functions. However, most experimental examples of promiscuity affording an observable phenotype necessitate the artificial overexpression of the target enzyme. Here, we show that 3-isopropylmalate dehydrogenase (IPMDH), an enzyme involved in leucine biosynthesis, has a secondary activity on d-malate, which is sufficient for d-malate assimilation under physiological conditions where the enzyme is upregulated. In vitro, the turnover constant (k_cat ) of IPMDH for d-malate is about 30-fold lower than the k_cat for 3-isopropylmalate, yet sufficiently high to support the growth on d-malate. From an evolutionary perspective, our results highlight the possibility of phenotype emergence triggered by arbitrary changes in environmental conditions and prior to any mutational event.

The physical basis and practical consequences of biological promiscuity

Proteins interact with metabolites, nucleic acids, and other proteins to orchestrate the myriad catalytic, structural and regulatory functions that support life from the simplest microbes to the most complex multicellular organisms. These molecular interactions are often exquisitely specific, but never perfectly so. Adventitious "promiscuous" interactions are ubiquitous due to the thousands of macromolecules and small molecules crowded together in cells. Such interactions may perturb protein function at the molecular level, but as long as they do not compromise organismal fitness, they will not be removed by natural selection. Although promiscuous interactions are physiologically irrelevant, they are important because they can provide a vast reservoir of potential functions that can provide the starting point for evolution of new functions, both in nature and in the laboratory.

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.

Multifaceted Stoichiometry Control of Bacterial Operons Revealed by Deep Proteome Quantification

More than half of the protein-coding genes in bacteria are organized in polycistronic operons composed of two or more genes. It remains under debate whether the operon organization maintains the stoichiometric expression of the genes within an operon. In this study, we performed a label-free data-independent acquisition hyper reaction monitoring mass-spectrometry (HRM-MS) experiment to quantify the Escherichia coli proteome in exponential phase and quantified 93.6% of the cytosolic proteins, covering 67.9% and 56.0% of the translating polycistronic operons in BW25113 and MG1655 strains, respectively. We found that the translational regulation contributes largely to the proteome complexity: the shorter operons tend to be more tightly controlled for stoichiometry than longer operons; the operons which mainly code for complexes is more tightly controlled for stoichiometry than the operons which mainly code for metabolic pathways. The gene interval (distance between adjacent genes in one operon) may serve as a regulatory factor for stoichiometry. The catalytic efficiency might be a driving force for differential expression of enzymes encoded in one operon. These results illustrated the multifaceted nature of the operon regulation: the operon unified transcriptional level and gene-specific translational level. This multi-level regulation benefits the host by optimizing the efficiency of the productivity of metabolic pathways and maintenance of different types of protein complexes.

Experimental Evolution of Innovation and Novelty

How does novelty, a new, genetically based function, evolve? A compelling answer has been elusive because there are few model systems where both the genetic mechanisms generating novel functions and the ecological conditions that govern their origin and spread can be studied in detail. This review article considers what we have learned about the evolution of novelty from microbial selection experiments. This work reveals that the genetic routes to novelty can be more highly variable than standard models have led us to believe and underscores the importance of considering both genetics and ecology in this process.

Exploring the sequence, function, and evolutionary space of protein superfamilies using sequence similarity networks and phylogenetic reconstructions

Integrative computational methods can facilitate the discovery of new protein functions and enzymatic reactions by enabling the observation and investigation of complex sequence-structure-function and evolutionary relationships within protein superfamilies. Here, we highlight the use of sequence similarity networks (SSNs) and phylogenetic reconstructions to map the functional divergence and evolutionary history of protein superfamilies. We exemplify this approach using the nitroreductase (NTR) flavoenzyme superfamily, demonstrating that SSN investigations can provide a rapid and effective means to classify groups of proteins, expose sequence similarity relationships across the global scale of a protein superfamily, and efficiently support detailed phylogenetic analyses. Integration of such approaches with systematic experimental characterization will expand our understanding of the functional diversity of enzymes, their evolution, and their associated physiological roles.

Synonymous mutations make dramatic contributions to fitness when growth is limited by a weak-link enzyme

Synonymous mutations do not alter the specified amino acid but may alter the structure or function of an mRNA in ways that impact fitness. There are few examples in the literature, however, in which the effects of synonymous mutations on microbial growth rates have been measured, and even fewer for which the underlying mechanism is understood. We evolved four populations of a strain of Salmonella enterica in which a promiscuous enzyme has been recruited to replace an essential enzyme. A previously identified point mutation increases the enzyme’s ability to catalyze the newly needed reaction (required for arginine biosynthesis) but decreases its ability to catalyze its native reaction (required for proline biosynthesis). The poor performance of this enzyme limits growth rate on glucose. After 260 generations, we identified two synonymous mutations in the first six codons of the gene encoding the weak-link enzyme that increase growth rate by 41 and 67%. We introduced all possible synonymous mutations into the first six codons and found substantial effects on growth rate; one doubles growth rate, and another completely abolishes growth. Computational analyses suggest that these mutations affect either the stability of a stem-loop structure that sequesters the start codon or the accessibility of the region between the Shine-Dalgarno sequence and the start codon. Thus, these mutations would be predicted to affect translational efficiency and thereby indirectly affect mRNA stability because translating ribosomes protect mRNA from degradation. Experimental data support these hypotheses. We conclude that the effects of the synonymous mutations are due to a combination of effects on mRNA stability and translation efficiency that alter levels of the weak-link enzyme. These findings suggest that synonymous mutations can have profound effects on fitness under strong selection and that their importance in evolution may be under-appreciated.

When a new enzyme is needed, microbes often recruit a pre-existing enzyme with a promiscuous activity corresponding to the newly needed activity. Such enzymes are often the “weak-link” in metabolism because they have not evolved to efficiently catalyze the new reaction. Under these circumstances, increasing the level of the weak-link enzyme can improve fitness. We evolved a strain of S. enterica in which a weak-link enzyme–E383A ProA–serves essential functions in synthesis of proline and arginine for 260 generations and then sequenced the genomes of several evolved strains. A mutation in the promoter of the operon encoding E383A ProA increased growth rate 9-fold. More surprisingly, a mutation upstream of the start codon and two synonymous mutations within the first six codons also increased growth rate by up to 68%. Introduction of all possible synonymous mutations in the first six codons showed that some doubled growth rate, while others slowed or even prevented growth. Computational and experimental data suggest that these effects were due to enhanced translational efficiency of the weak-link enzyme. These results show that synonymous mutations, once assumed to be selectively neutral, can have strong impacts on fitness when growth rate is limited by a weak-link enzyme.

Revealing Unexplored Sequence-Function Space Using Sequence Similarity Networks

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

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.

Expanding Primary Metabolism Helps Generate the Metabolic Robustness To Facilitate Antibiotic Biosynthesis in Streptomyces

The expansion of the genetic repertoire of an organism by gene duplication or horizontal gene transfer (HGT) can aid adaptation. Streptomyces bacteria are prolific producers of bioactive specialized metabolites that have adaptive functions in nature and have found extensive utility in human medicine. While the biosynthesis of these specialized metabolites is directed by dedicated biosynthetic gene clusters, little attention has been focused on how these organisms have evolved robustness in their genomes to facilitate the metabolic plasticity required to provide chemical precursors for biosynthesis during the complex metabolic transitions from vegetative growth to specialized metabolite production and sporulation. Here, we examine genetic redundancy in actinobacteria and show that specialized metabolite-producing bacterial families exhibit gene family expansion in primary metabolism. Focusing on a gene duplication event, we show that the two pyruvate kinases in the genome of Streptomyces coelicolor arose by an ancient duplication event and that each has evolved altered enzymatic kinetics, with Pyk1 having a 20-fold-higher kcat than Pyk2 (4,703 s-1 compared to 215 s-1, respectively), and yet both are constitutively expressed. The pyruvate kinase mutants were also found to be compromised in terms of fitness compared to wild-type Streptomyces These data suggest that expanding gene families can help maintain cell functionality during metabolic perturbation such as nutrient limitation and/or specialized metabolite production.IMPORTANCE The rise of antimicrobial-resistant infections has prompted a resurgence in interest in understanding the production of specialized metabolites, such as antibiotics, by Streptomyces The presence of multiple genes encoding the same enzymatic function is an aspect of Streptomyces biology that has received little attention; however, understanding how the metabolic expansion influences these organisms can help enhance production of clinically useful molecules. Here, we show that expanding the number of pyruvate kinases enables metabolic adaptation, increases strain fitness, and represents an excellent target for metabolic engineering of industrial specialized metabolite-producing bacteria and the activation of cryptic specialized metabolites.

SAR202 Genomes from the Dark Ocean Predict Pathways for the Oxidation of Recalcitrant Dissolved Organic Matter

Deep-ocean regions beyond the reach of sunlight contain an estimated 615 Pg of dissolved organic matter (DOM), much of which persists for thousands of years. It is thought that bacteria oxidize DOM until it is too dilute or refractory to support microbial activity. We analyzed five single-amplified genomes (SAGs) from the abundant SAR202 clade of dark-ocean bacterioplankton and found they encode multiple families of paralogous enzymes involved in carbon catabolism, including several families of oxidative enzymes that we hypothesize participate in the degradation of cyclic alkanes. The five partial genomes encoded 152 flavin mononucleotide/F420-dependent monooxygenases (FMNOs), many of which are predicted to be type II Baeyer-Villiger monooxygenases (BVMOs) that catalyze oxygen insertion into semilabile alicyclic alkanes. The large number of oxidative enzymes, as well as other families of enzymes that appear to play complementary roles in catabolic pathways, suggests that SAR202 might catalyze final steps in the biological oxidation of relatively recalcitrant organic compounds to refractory compounds that persist.IMPORTANCE Carbon in the ocean is massively sequestered in a complex mixture of biologically refractory molecules that accumulate as the chemical end member of biological oxidation and diagenetic change. However, few details are known about the biochemical machinery of carbon sequestration in the deep ocean. Reconstruction of the metabolism of a deep-ocean microbial clade, SAR202, led to postulation of new biochemical pathways that may be the penultimate stages of DOM oxidation to refractory forms that persist. These pathways are tied to a proliferation of oxidative enzymes. This research illuminates dark-ocean biochemistry that is broadly consequential for reconstructing the global carbon cycle.

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

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

Understanding how enzymes evolve is a fundamental question that can help us decipher not only the mechanisms of evolution on a higher level, i.e., whole organisms, but also advances our knowledge of sequence-structure-function relationships as a guide to artificial evolution in the test tube. An important yet unexplained phenomenon occurs during the evolution of a new enzymatic function; it has been observed that new and ancestral functions often trade-off only weakly, meaning the original native activity is initially maintained at a high level despite drastic improvement of the new promiscuous activity. It has previously been proposed that weak trade-offs occur because the native activity is robust to mutations while the promiscuous activity is not. However, the present work contradicts this hypothesis, based on the detailed characterization of mutational effects on both activities in two examples of enzyme evolution. We propose an alternative explanation: the weak activity trade-off is consistent with being a by-product of strong selection for the new activity rather than an intrinsic property of the native activity.

A Synonymous Mutation Upstream of the Gene Encoding a Weak-Link Enzyme Causes an Ultrasensitive Response in Growth Rate

When microbes are faced with an environmental challenge or opportunity, preexisting enzymes with promiscuous secondary activities can be recruited to provide newly important functions. Mutations that increase the efficiency of a new activity often compromise the original activity, resulting in an inefficient bifunctional enzyme. We have investigated the mechanisms by which growth of Escherichia coli can be improved when fitness is limited by such an enzyme, E383A ProA (ProA*). ProA* can serve the functions of both ProA (required for synthesis of proline) and ArgC (required for synthesis of arginine), albeit poorly. We identified four genetic changes that improve the growth rate by up to 6.2-fold. Two point mutations in the promoter of the proBA* operon increase expression of the entire operon. Massive amplification of a genomic segment around the proBA* operon also increases expression of the entire operon. Finally, a synonymous point mutation in the coding region of proB creates a new promoter for proA* This synonymous mutation increases the level of ProA* by 2-fold but increases the growth rate by 5-fold, an ultrasensitive response likely arising from competition between two substrates for the active site of the inefficient bifunctional ProA*.

IMPORTANCE

The high-impact synonymous mutation we discovered in proB is remarkable for two reasons. First, most polar effects documented in the literature are detrimental. This finding demonstrates that polar effect mutations can have strongly beneficial effects, especially when an organism is facing a difficult environmental challenge for which it is poorly adapted. Furthermore, the consequence of the synonymous mutation in proB is a 2-fold increase in the level of ProA* but a disproportionately large 5.1-fold increase in growth rate. While ultrasensitive responses are often found in signaling networks and genetic circuits, an ultrasensitive response to an adaptive mutation has not been previously reported.

Rapid bursts and slow declines: on the possible evolutionary trajectories of enzymes

The evolution of enzymes is often viewed as following a smooth and steady trajectory, from barely functional primordial catalysts to the highly active and specific enzymes that we observe today. In this review, we summarize experimental data that suggest a different reality. Modern examples, such as the emergence of enzymes that hydrolyse human-made pesticides, demonstrate that evolution can be extraordinarily rapid. Experiments to infer and resurrect ancient sequences suggest that some of the first organisms present on the Earth are likely to have possessed highly active enzymes. Reconciling these observations, we argue that rapid bursts of strong selection for increased catalytic efficiency are interspersed with much longer periods in which the catalytic power of an enzyme erodes, through neutral drift and selection for other properties such as cellular energy efficiency or regulation. Thus, many enzymes may have already passed their catalytic peaks.

Crystal Structure of the LysY·LysW Complex from Thermus thermophilus

Several bacteria and archaea utilize the amino group-carrier protein, LysW, for lysine biosynthesis, in which an isopeptide bond is formed between the C-terminal Glu of LysW and an amino group of α-aminoadipate (AAA). The resulting LysW-γ-AAA is phosphorylated by LysZ to form LysW-γ-AAA phosphate, which is subsequently reduced to LysW-γ-aminoadipic semialdehyde (LysW-γ-AASA) through a reaction catalyzed by LysY. In this study, we determined the crystal structures of LysY from Thermus thermophilus HB27 (TtLysY) complexed with TtLysW-γ-AASA and TtLysW-γ-AAA, respectively. In both structures, the globular domain of TtLysW was recognized by positively charged residues on helix α9 and the β11-α10 loop of TtLysY through conformational changes. A mutational analysis confirmed that the interactions observed between TtLysY and TtLysW are important for the function of TtLysY. The extended LysW recognition loop and conserved arginine residue were identified as signatures to discriminate LysY from ArgC, which is involved in arginine biosynthesis. Combined with the previously determined TtLysZ·TtLysW complex structure, TtLysW may simultaneously bind TtLysZ and TtLysY. These structural insights suggest the formation of a TtLysWZY ternary complex, in which the flexible C-terminal extension of TtLysW promotes the efficient transfer of the labile intermediate from the active site of TtLysZ to that of TtLysY during the sequential reactions catalyzed by TtLysZY.

Evidence that Environmental Heterogeneity Maintains a Detoxifying Enzyme Polymorphism in Drosophila melanogaster

Environmental heterogeneity is thought to be an important process maintaining genetic variation in populations [1-4]: if alternative alleles are favored in different environments, a stable polymorphism can be maintained [1, 5, 6]. This situation has been hypothesized to occur in genes encoding multi-substrate enzymes [7], in which changes that increase activity with one substrate typically decrease activity with others [8-10], but examples of polymorphisms maintained by this mechanism are rare. Here, we present evidence that a polymorphism in an enzyme gene in Drosophila melanogaster is maintained by such a trade-off. The mitochondrially localized aldehyde dehydrogenase in D. melanogaster has two important functions: detoxifying acetaldehyde derived from dietary ethanol [11] and detoxifying larger aldehydes produced as byproducts of oxidative phosphorylation [12]. A derived variant of the enzyme, Leu479Phe, is present in moderate frequencies in most temperate populations but is rare in more ethanol-averse tropical populations. Using purified recombinant protein, we show that the Leu-Phe substitution increases turnover rate of acetaldehyde but decreases turnover rate of larger aldehydes. Furthermore, using transgenic fly lines, we show that the substitution increases lifetime fitness on medium supplemented with an ecologically relevant ethanol concentration but decreases fitness on medium lacking ethanol. The strong, opposing selection pressures, coupled with documented highly variable ethanol concentrations in breeding sites of temperate populations, implicate an essential role for environmental heterogeneity in maintaining the polymorphism.

An evolutionary biochemist's perspective on promiscuity

Evolutionary biochemists define enzyme promiscuity as the ability to catalyze secondary reactions that are physiologically irrelevant, either because they are too inefficient to affect fitness or because the enzyme never encounters the substrate. Promiscuous activities are common because evolution of a perfectly specific active site is both difficult and unnecessary; natural selection ceases when the performance of a protein is 'good enough' that it no longer affects fitness. Although promiscuous functions are accidental and physiologically irrelevant, they are of great importance because they provide opportunities for the evolution of new functions in nature and in the laboratory, as well as targets for therapeutic drugs and tools for a wide range of technological applications.

Differential effects of a mutation on the normal and promiscuous activities of orthologs: implications for natural and directed evolution

Neutral drift occurring over millions or billions of years results in substantial sequence divergence among enzymes that catalyze the same reaction. Although natural selection maintains the primary activity of orthologous enzymes, there is, by definition, no selective pressure to maintain physiologically irrelevant promiscuous activities. Thus, the levels and the evolvabilities of promiscuous activities may vary among orthologous enzymes. Consistent with this expectation, we have found that the levels of a promiscuous activity in nine gamma-glutamyl phosphate reductase (ProA) orthologs vary by about 50-fold. Remarkably, a single amino acid change from Glu to Ala near the active site appeared to be critical for improvement of the promiscuous activity in every ortholog. The effects of this change varied dramatically. The improvement in the promiscuous activity varied from 50- to 770-fold, and, importantly, was not correlated with the initial level of the promiscuous activity. The decrease in the original activity varied from 190- to 2,100-fold. These results suggest that evolution of a novel enzyme may be possible in some microbes, but not in others. Further, these results underscore the importance of using multiple orthologs as starting points for directed evolution of novel enzyme activities.

Escherichia coli D-malate dehydrogenase, a generalist enzyme active in the leucine biosynthesis pathway

The enzymes of the β-decarboxylating dehydrogenase superfamily catalyze the oxidative decarboxylation of D-malate-based substrates with various specificities. Here, we show that, in addition to its natural function affording bacterial growth on D-malate as a carbon source, the D-malate dehydrogenase of Escherichia coli (EcDmlA) naturally expressed from its chromosomal gene is capable of complementing leucine auxotrophy in a leuB(-) strain lacking the paralogous isopropylmalate dehydrogenase enzyme. To our knowledge, this is the first example of an enzyme that contributes with a physiologically relevant level of activity to two distinct pathways of the core metabolism while expressed from its chromosomal locus. EcDmlA features relatively high catalytic activity on at least three different substrates (L(+)-tartrate, D-malate, and 3-isopropylmalate). Because of these properties both in vivo and in vitro, EcDmlA may be defined as a generalist enzyme. Phylogenetic analysis highlights an ancient origin of DmlA, indicating that the enzyme has maintained its generalist character throughout evolution. We discuss the implication of these findings for protein evolution.

Improved cyclopropanation activity of histidine-ligated cytochrome P450 enables the enantioselective formal synthesis of levomilnacipran

Engineering enzymes capable of modes of activation unprecedented in nature will increase the range of industrially important molecules that can be synthesized through biocatalysis. However, low activity for a new function is often a limitation in adopting enzymes for preparative-scale synthesis, reaction with demanding substrates, or when a natural substrate is also present. By mutating the proximal ligand and other key active-site residues of the cytochrome P450 enzyme from Bacillus megaterium (P450-BM3), a highly active His-ligated variant of P450-BM3 that can be employed for the enantioselective synthesis of the levomilnacipran core was engineered. This enzyme, BM3-Hstar, catalyzes the cyclopropanation of N,N-diethyl-2-phenylacrylamide with an estimated initial rate of over 1000 turnovers per minute and can be used under aerobic conditions. Cyclopropanation activity is highly dependent on the electronic properties of the P450 proximal ligand, which can be used to tune this non-natural enzyme activity.

Substrate ambiguous enzymes within the Escherichia coli proteome offer different evolutionary solutions to the same problem

Many enzymes exhibit some catalytic promiscuity or substrate ambiguity. These weak activities do not affect the fitness of the organism under ordinary circumstances, but can serve as potential evolutionary precursors of new catalytic functions. We wondered whether different proteins with the same substrate ambiguous activity evolve differently under identical selection conditions. Patrick et al. (Patrick WM, Quandt EM, Swartzlander DB, Matsumura I. 2007. Multicopy suppression underpins metabolic evolvability. Mol Biol Evol. 24:2716-2722.) previously showed that three multicopy suppressors, gph, hisB, and ytjC, rescue ΔserB Escherichia coli cells from starvation on minimal media. We directed the evolution of variants of Gph, histidinol phosphatase (HisB), and YtjC that complemented ΔserB more efficiently, and characterized the effects of the amino acid changes, alone and in combination, upon the evolved phosphoserine phosphatase (PSP) activity. Gph and HisB are members of the HAD superfamily of hydrolases, but they adapted through different, kinetically distinguishable, biochemical mechanisms. All of the selected mutations, except N102T in YtjC, proved to be beneficial in isolation. They exhibited a pattern of antagonistic epistasis, as their effects in combination upon the kinetic parameters of the three proteins in reactions with phosphoserine were nonmultiplicative. The N102T mutation exhibited sign epistasis, as it was deleterious in isolation but beneficial in the context of other mutations. We also showed that the D57N mutation in the chromosomal copy of hisB is sufficient to suppress the ΔserB deletion. These results in combination show that proteomes can offer multiple mechanistic solutions to a molecular recognition problem.

Probing the mutational interplay between primary and promiscuous protein functions: a computational-experimental approach

Protein promiscuity is of considerable interest due its role in adaptive metabolic plasticity, its fundamental connection with molecular evolution and also because of its biotechnological applications. Current views on the relation between primary and promiscuous protein activities stem largely from laboratory evolution experiments aimed at increasing promiscuous activity levels. Here, on the other hand, we attempt to assess the main features of the simultaneous modulation of the primary and promiscuous functions during the course of natural evolution. The computational/experimental approach we propose for this task involves the following steps: a function-targeted, statistical coupling analysis of evolutionary data is used to determine a set of positions likely linked to the recruitment of a promiscuous activity for a new function; a combinatorial library of mutations on this set of positions is prepared and screened for both, the primary and the promiscuous activities; a partial-least-squares reconstruction of the full combinatorial space is carried out; finally, an approximation to the Pareto set of variants with optimal primary/promiscuous activities is derived. Application of the approach to the emergence of folding catalysis in thioredoxin scaffolds reveals an unanticipated scenario: diverse patterns of primary/promiscuous activity modulation are possible, including a moderate (but likely significant in a biological context) simultaneous enhancement of both activities. We show that this scenario can be most simply explained on the basis of the conformational diversity hypothesis, although alternative interpretations cannot be ruled out. Overall, the results reported may help clarify the mechanisms of the evolution of new functions. From a different viewpoint, the partial-least-squares-reconstruction/Pareto-set-prediction approach we have introduced provides the computational basis for an efficient directed-evolution protocol aimed at the simultaneous enhancement of several protein features and should therefore open new possibilities in the engineering of multi-functional enzymes.

Interpretations of evolutionary processes at the molecular level have been determined to a significant extent by the concept of “trade-off”, the idea that improving a given feature of a protein molecule by mutation will likely bring about deterioration in other features. For instance, if a protein is able to carry out two different molecular tasks based on the same functional site (competing tasks), optimization for one task could be naively expected to impair its performance for the other task. In this work, we report a computational/experimental approach to assess the potential patterns of modulation of two competing molecular tasks in the course of natural evolution. Contrary to the naïve expectation, we find that diverse modulation patterns are possible, including the simultaneous optimization of the two tasks. We show, however, that this simultaneous optimization is not in conflict with the trade-offs expected for two competing tasks: using the language of the theory of economic efficiency, trade-offs are realized in the Pareto set of optimal variants for the two tasks, while most protein variants do not belong to such Pareto set. That is, most protein variants are not Pareto-efficient and can potentially be improved in terms of several features.

A trait-based approach for modelling microbial litter decomposition

Trait-based models are an emerging tool in ecology with the potential to link community dynamics, environmental responses and ecosystem processes. These models represent complex communities by defining taxa with trait combinations derived from prior distributions that may be constrained by trade-offs. Herein I develop a model that links microbial community composition with physiological and enzymatic traits to predict litter decomposition rates. This approach allows for trade-offs among traits that represent alternative microbial strategies for resource acquisition. The model predicts that optimal strategies depend on the level of enzyme production in the whole community, which determines resource availability and decomposition rates. There is also evidence for facilitation and competition among microbial taxa that co-occur on decomposing litter. These interactions vary with community investment in extracellular enzyme production and the magnitude of trade-offs affecting enzyme biochemical traits. The model accounted for 69% of the variation in decomposition rates of 15 Hawaiian litter types and up to 26% of the variation in enzyme activities. By explicitly representing diversity, trait-based models can predict ecosystem processes based on functional trait distributions in a community. The model developed herein illustrates that traits influencing microbial enzyme production are some of the key controls on litter decomposition rates.

Enzyme functional evolution through improved catalysis of ancestrally nonpreferred substrates

In this study, we investigated the role for ancestral functional variation that may be selected upon to generate protein functional shifts using ancestral protein resurrection, statistical tests for positive selection, forward and reverse evolutionary genetics, and enzyme functional assays. Data are presented for three instances of protein functional change in the salicylic acid/benzoic acid/theobromine (SABATH) lineage of plant secondary metabolite-producing enzymes. In each case, we demonstrate that ancestral nonpreferred activities were improved upon in a daughter enzyme after gene duplication, and that these functional shifts were likely coincident with positive selection. Both forward and reverse mutagenesis studies validate the impact of one or a few sites toward increasing activity with ancestrally nonpreferred substrates. In one case, we document the occurrence of an evolutionary reversal of an active site residue that reversed enzyme properties. Furthermore, these studies show that functionally important amino acid replacements result in substrate discrimination as reflected in evolutionary changes in the specificity constant (k(cat)/K(M)) for competing substrates, even though adaptive substitutions may affect K(M) and k(cat) separately. In total, these results indicate that nonpreferred, or even latent, ancestral protein activities may be coopted at later times to become the primary or preferred protein activities.

Toward a systems biology perspective on enzyme evolution

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

The evolution of cefotaximase activity in the TEM β-lactamase

The development of a molecular-level understanding of drug resistance through β-lactamase is critical not only in designing newer-generation antibacterial agents but also in providing insight into the evolutionary mechanisms of enzymes in general. In the present study, we have evaluated the effect of four drug resistance mutations (A42G, E104K, G238S, and M182T) on the cefotaximase activity of the TEM-1 β-lactamase. Using computational methods, including docking and molecular mechanics calculations, we have been able to correctly identify the relative order of catalytic activities associated with these four single point mutants. Further analyses suggest that the changes in catalytic efficiency for mutant enzymes are correlated to structural changes within the binding site. Based on the energetic and structural analyses of the wild-type and mutant enzymes, structural rearrangement is suggested as a mechanism of evolution of drug resistance through TEM β-lactamase. The present study not only provides molecular-level insight into the effect of four drug resistance mutations on the structure and function of the TEM β-lactamase but also establishes a foundation for a future molecular-level analysis of complete evolutionary trajectory for this class of enzymes.

Molecular mechanisms of epistasis within and between genes

'Disease-causing' mutations do not cause disease in all individuals. One possible important reason for this is that the outcome of a mutation can depend upon other genetic variants in a genome. These epistatic interactions between mutations occur both within and between molecules, and studies in model organisms show that they are extremely prevalent. However, epistatic interactions are still poorly understood at the molecular level, and consequently difficult to predict de novo. Here I provide an overview of our current understanding of the molecular mechanisms that can cause epistasis, and areas where more research is needed. A more complete understanding of epistasis will be vital for making accurate predictions about the phenotypes of individuals.

Mutational effects and the evolution of new protein functions

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

Enzyme promiscuity: a mechanistic and evolutionary perspective

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

Mechanism and plasticity of isochorismate pyruvate lyase: a computational study

The isochorismate pyruvate lyase (IPL) from Pseudomonas aeruginosa, designated as PchB, catalyzes the transformation of isochorismate into pyruvate and salicylate, but it also catalyzes the rearrangement of chorismate into prephenate, suggesting that both reactions may proceed by a pericyclic mechanism. In this work, molecular dynamics simulations employing hybrid quantum mechanics/molecular mechanics methods have been carried out to get a detailed knowledge of the reaction mechanism of PchB. The results provide a theoretical rate constant enhancement by comparison with the reaction in solution, in agreement with the experimental data, and confirm the pericyclic nature of the reaction mechanism. The robustness of this promiscuous enzyme has been checked by considering the impact of Ala37Ile mutation, previously proposed by us to improve the secondary chorismate mutase (CM) activity. The effect of this mutation, which was shown to increase the rate constant for the CM activity by a factor of 10(3), also increases the IPL catalytic efficiency, although only by a factor of 6.

Comparative study of human mitochondrial proteome reveals extensive protein subcellular relocalization after gene duplications

Background

Gene and genome duplication is the principle creative force in evolution. Recently, protein subcellular relocalization, or neolocalization was proposed as one of the mechanisms responsible for the retention of duplicated genes. This hypothesis received support from the analysis of yeast genomes, but has not been tested thoroughly on animal genomes. In order to evaluate the importance of subcellular relocalizations for retention of duplicated genes in animal genomes, we systematically analyzed nuclear encoded mitochondrial proteins in the human genome by reconstructing phylogenies of mitochondrial multigene families.

Results

The 456 human mitochondrial proteins selected for this study were clustered into 305 gene families including 92 multigene families. Among the multigene families, 59 (64%) consisted of both mitochondrial and cytosolic (non-mitochondrial) proteins (mt-cy families) while the remaining 33 (36%) were composed of mitochondrial proteins (mt-mt families). Phylogenetic analyses of mt-cy families revealed three different scenarios of their neolocalization following gene duplication: 1) relocalization from mitochondria to cytosol, 2) from cytosol to mitochondria and 3) multiple subcellular relocalizations. The neolocalizations were most commonly enabled by the gain or loss of N-terminal mitochondrial targeting signals. The majority of detected subcellular relocalization events occurred early in animal evolution, preceding the evolution of tetrapods. Mt-mt protein families showed a somewhat different pattern, where gene duplication occurred more evenly in time. However, for both types of protein families, most duplication events appear to roughly coincide with two rounds of genome duplications early in vertebrate evolution. Finally, we evaluated the effects of inaccurate and incomplete annotation of mitochondrial proteins and found that our conclusion of the importance of subcellular relocalization after gene duplication on the genomic scale was robust to potential gene misannotation.

Conclusion

Our results suggest that protein subcellular relocalization is an important mechanism for the retention and gain of function of duplicated genes in animal genome evolution.

Chaperonin overexpression promotes genetic variation and enzyme evolution

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