Multicopy suppression underpins metabolic evolvability

Connectivity between catalytic landscapes of the metallo-β-lactamase superfamily

The expansion of functions in an enzyme superfamily is thought to occur through recruitment of latent promiscuous functions within existing enzymes. Thus, the promiscuous activities of enzymes represent connections between different catalytic landscapes and provide an additional layer of evolutionary connectivity between functional families alongside their sequence and structural relationships. Functional connectivity has been observed between individual functional families; however, little is known about how catalytic landscapes are connected throughout a highly diverged superfamily. Here, we describe a superfamily-wide analysis of evolutionary and functional connectivity in the metallo-β-lactamase (MBL) superfamily. We investigated evolutionary connections between functional families and related evolutionary to functional connectivity; 24 enzymes from 15 distinct functional families were challenged against 10 catalytically distinct reactions. We revealed that enzymes of this superfamily are generally promiscuous, as each enzyme catalyzes on average 1.5 reactions in addition to its native one. Catalytic landscapes in the MBL superfamily overlap substantially; each reaction is connected on average to 3.7 other reactions whereas some connections appear to be unrelated to recent evolutionary events and occur between chemically distinct reactions. These findings support the idea that the highly distinct reactions in the MBL superfamily could have evolved from a common ancestor traversing a continuous network via promiscuous enzymes. Several functional connections (e.g., the lactonase/phosphotriesterase and phosphonatase/phosphodiesterase/arylsulfatase reactions) are also observed in structurally and evolutionary distinct superfamilies, suggesting that these catalytic landscapes are substantially connected. Our results show that new enzymatic functions could evolve rapidly from the current diversity of enzymes and range of promiscuous activities.

k-OptForce: integrating kinetics with flux balance analysis for strain design

Computational strain design protocols aim at the system-wide identification of intervention strategies for the enhanced production of biochemicals in microorganisms. Existing approaches relying solely on stoichiometry and rudimentary constraint-based regulation overlook the effects of metabolite concentrations and substrate-level enzyme regulation while identifying metabolic interventions. In this paper, we introduce k-OptForce, which integrates the available kinetic descriptions of metabolic steps with stoichiometric models to sharpen the prediction of intervention strategies for improving the bio-production of a chemical of interest. It enables identification of a minimal set of interventions comprised of both enzymatic parameter changes (for reactions with available kinetics) and reaction flux changes (for reactions with only stoichiometric information). Application of k-OptForce to the overproduction of L-serine in E. coli and triacetic acid lactone (TAL) in S. cerevisiae revealed that the identified interventions tend to cause less dramatic rearrangements of the flux distribution so as not to violate concentration bounds. In some cases the incorporation of kinetic information leads to the need for additional interventions as kinetic expressions render stoichiometry-only derived interventions infeasible by violating concentration bounds, whereas in other cases the kinetic expressions impart flux changes that favor the overproduction of the target product thereby requiring fewer direct interventions. A sensitivity analysis on metabolite concentrations shows that the required number of interventions can be significantly affected by changing the imposed bounds on metabolite concentrations. Furthermore, k-OptForce was capable of finding non-intuitive interventions aiming at alleviating the substrate-level inhibition of key enzymes in order to enhance the flux towards the product of interest, which cannot be captured by stoichiometry-alone analysis. This study paves the way for the integrated analysis of kinetic and stoichiometric models and enables elucidating system-wide metabolic interventions while capturing regulatory and kinetic effects.

Enzyme recruitment and its role in metabolic expansion

Although more than 10(9) years have passed since the existence of the last universal common ancestor, proteins have yet to reach the limits of divergence. As a result, metabolic complexity is ever expanding. Identifying and understanding the mechanisms that drive and limit the divergence of protein sequence space impact not only evolutionary biologists investigating molecular evolution but also synthetic biologists seeking to design useful catalysts and engineer novel metabolic pathways. Investigations over the past 50 years indicate that the recruitment of enzymes for new functions is a key event in the acquisition of new metabolic capacity. In this review, we outline the genetic mechanisms that enable recruitment and summarize the present state of knowledge regarding the functional characteristics of extant catalysts that facilitate recruitment. We also highlight recent examples of enzyme recruitment, both from the historical record provided by phylogenetics and from enzyme evolution experiments. We conclude with a look to the future, which promises fruitful consequences from the convergence of molecular evolutionary theory, laboratory-directed evolution, and synthetic biology.

Copy number change: evolving views on gene amplification

The rapid pace of genomic sequence analysis is increasing the awareness of intrinsically dynamic genetic landscapes. Gene duplication and amplification (GDA) contribute to adaptation and evolution by allowing DNA regions to expand and contract in an accordion-like fashion. This process affects diverse aspects of bacterial infection, including antibiotic resistance and host-pathogen interactions. In this review, microbial GDA is discussed, primarily using recent bacterial examples that demonstrate medical and evolutionary consequences. Interplay between GDA and horizontal gene transfer further impact evolutionary trajectories. Complementing the discovery of gene duplication in clinical and environmental settings, experimental evolution provides a powerful method to document genetic change over time. New methods for GDA detection highlight both its importance and its potential application for genetic engineering, synthetic biology and biotechnology.

Addition of Escherichia coli K-12 growth observation and gene essentiality data to the EcoCyc database

The sets of compounds that can support growth of an organism are defined by the presence of transporters and metabolic pathways that convert nutrient sources into cellular components and energy for growth. A collection of known nutrient sources can therefore serve both as an impetus for investigating new metabolic pathways and transporters and as a reference for computational modeling of known metabolic pathways. To establish such a collection for Escherichia coli K-12, we have integrated data on the growth or nongrowth of E. coli K-12 obtained from published observations using a variety of individual media and from high-throughput phenotype microarrays into the EcoCyc database. The assembled collection revealed a substantial number of discrepancies between the high-throughput data sets, which we investigated where possible using low-throughput growth assays on soft agar and in liquid culture. We also integrated six data sets describing 16,119 observations of the growth of single-gene knockout mutants of E. coli K-12 into EcoCyc, which are relevant to antimicrobial drug design, provide clues regarding the roles of genes of unknown function, and are useful for validating metabolic models. To make this information easily accessible to EcoCyc users, we developed software for capturing, querying, and visualizing cellular growth assays and gene essentiality data.

Copy-number changes in evolution: rates, fitness effects and adaptive significance

Gene copy-number differences due to gene duplications and deletions are rampant in natural populations and play a crucial role in the evolution of genome complexity. Per-locus analyses of gene duplication rates in the pre-genomic era revealed that gene duplication rates are much higher than the per nucleotide substitution rate. Analyses of gene duplication and deletion rates in mutation accumulation lines of model organisms have revealed that these high rates of copy-number mutations occur at a genome-wide scale. Furthermore, comparisons of the spontaneous duplication and deletion rates to copy-number polymorphism data and bioinformatic-based estimates of duplication rates from sequenced genomes suggest that the vast majority of gene duplications are detrimental and removed by natural selection. The rate at which new gene copies appear in populations greatly influences their evolutionary dynamics and standing gene copy-number variation in populations. The opportunity for mutations that result in the maintenance of duplicate copies, either through neofunctionalization or subfunctionalization, also depends on the equilibrium frequency of additional gene copies in the population, and hence on the spontaneous gene duplication (and loss) rate. The duplication rate may therefore have profound effects on the role of adaptation in the evolution of duplicated genes as well as important consequences for the evolutionary potential of organisms. We further discuss the broad ramifications of this standing gene copy-number variation on fitness and adaptive potential from a population-genetic and genome-wide perspective.

Limited influence of oxygen on the evolution of chemical diversity in metabolic networks

Oxygen is thought to promote species and biomolecule diversity. Previous studies have suggested that oxygen expands metabolic networks by acquiring metabolites with different chemical properties (higher hydrophobicity, for example). However, such conclusions are typically based on biased evaluation, and are therefore non-conclusive. Thus, we re-investigated the effect of oxygen on metabolic evolution using a phylogenetic comparative method and metadata analysis to reduce the bias as much as possible. Notably, we found no difference in metabolic network expansion between aerobes and anaerobes when evaluating phylogenetic relationships. Furthermore, we showed that previous studies have overestimated or underestimated the degrees of differences in the chemical properties (e.g., hydrophobicity) between oxic and anoxic metabolites in metabolic networks of unicellular organisms; however, such overestimation was not observed when considering the metabolic networks of multicellular organisms. These findings indicate that the contribution of oxygen to increased chemical diversity in metabolic networks is lower than previously thought; rather, phylogenetic signals and cell-cell communication result in increased chemical diversity. However, this conclusion does not contradict the effect of oxygen on metabolic evolution; instead, it provides a deeper understanding of how oxygen contributes to metabolic evolution despite several limitations in data analysis methods.

Quantitative assessment of thermodynamic constraints on the solution space of genome-scale metabolic models

Constraint-based methods provide powerful computational techniques to allow understanding and prediction of cellular behavior. These methods rely on physiochemical constraints to eliminate infeasible behaviors from the space of available behaviors. One such constraint is thermodynamic feasibility, the requirement that intracellular flux distributions obey the laws of thermodynamics. The past decade has seen several constraint-based methods that interpret this constraint in different ways, including those that are limited to small networks, rely on predefined reaction directions, and/or neglect the relationship between reaction free energies and metabolite concentrations. In this work, we utilize one such approach, thermodynamics-based metabolic flux analysis (TMFA), to make genome-scale, quantitative predictions about metabolite concentrations and reaction free energies in the absence of prior knowledge of reaction directions, while accounting for uncertainties in thermodynamic estimates. We applied TMFA to a genome-scale network reconstruction of Escherichia coli and examined the effect of thermodynamic constraints on the flux space. We also assessed the predictive performance of TMFA against gene essentiality and quantitative metabolomics data, under both aerobic and anaerobic, and optimal and suboptimal growth conditions. Based on these results, we propose that TMFA is a useful tool for validating phenotypes and generating hypotheses, and that additional types of data and constraints can improve predictions of metabolite concentrations.

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.

Does habitat variability really promote metabolic network modularity?

The hypothesis that variability in natural habitats promotes modular organization is widely accepted for cellular networks. However, results of some data analyses and theoretical studies have begun to cast doubt on the impact of habitat variability on modularity in metabolic networks. Therefore, we re-evaluated this hypothesis using statistical data analysis and current metabolic information. We were unable to conclude that an increase in modularity was the result of habitat variability. Although horizontal gene transfer was also considered because it may contribute for survival in a variety of environments, closely related to habitat variability, and is known to be positively correlated with network modularity, such a positive correlation was not concluded in the latest version of metabolic networks. Furthermore, we demonstrated that the previously observed increase in network modularity due to habitat variability and horizontal gene transfer was probably due to a lack of available data on metabolic reactions. Instead, we determined that modularity in metabolic networks is dependent on species growth conditions. These results may not entirely discount the impact of habitat variability and horizontal gene transfer. Rather, they highlight the need for a more suitable definition of habitat variability and a more careful examination of relationships of the network modularity with horizontal gene transfer, habitats, and environments.

FOCAL: an experimental design tool for systematizing metabolic discoveries and model development

Current computational tools can generate and improve genome-scale models based on existing data; however, for many organisms, the data needed to test and refine such models are not available. To facilitate model development, we created the forced coupling algorithm, FOCAL, to identify genetic and environmental conditions such that a reaction becomes essential for an experimentally measurable phenotype. This reaction's conditional essentiality can then be tested experimentally to evaluate whether network connections occur or to create strains with desirable phenotypes. FOCAL allows network connections to be queried, which improves our understanding of metabolism and accuracy of developed models.

EcoCyc: fusing model organism databases with systems biology

EcoCyc (http://EcoCyc.org) is a model organism database built on the genome sequence of Escherichia coli K-12 MG1655. Expert manual curation of the functions of individual E. coli gene products in EcoCyc has been based on information found in the experimental literature for E. coli K-12-derived strains. Updates to EcoCyc content continue to improve the comprehensive picture of E. coli biology. The utility of EcoCyc is enhanced by new tools available on the EcoCyc web site, and the development of EcoCyc as a teaching tool is increasing the impact of the knowledge collected in EcoCyc.

Metabolic network modularity arising from simple growth processes

Metabolic networks consist of linked functional components, or modules. The mechanism underlying metabolic network modularity is of great interest not only to researchers of basic science but also to those in fields of engineering. Previous studies have suggested a theoretical model, which proposes that a change in the evolutionary goal (system-specific purpose) increases network modularity, and this hypothesis was supported by statistical data analysis. Nevertheless, further investigation has uncovered additional possibilities that might explain the origin of network modularity. In this work we propose an evolving network model without tuning parameters to describe metabolic networks. We demonstrate, quantitatively, that metabolic network modularity can arise from simple growth processes, independent of the change in the evolutionary goal. Our model is applicable to a wide range of organisms and appears to suggest that metabolic network modularity can be more simply determined than previously thought. Nonetheless, our proposition does not serve to contradict the previous model; it strives to provide an insight from a different angle in the ongoing efforts to understand metabolic evolution, with the hope of eventually achieving the synthetic engineering of metabolic networks.

Current understanding of the formation and adaptation of metabolic systems based on network theory

Formation and adaptation of metabolic networks has been a long-standing question in biology. With recent developments in biotechnology and bioinformatics, the understanding of metabolism is progressively becoming clearer from a network perspective. This review introduces the comprehensive metabolic world that has been revealed by a wide range of data analyses and theoretical studies; in particular, it illustrates the role of evolutionary events, such as gene duplication and horizontal gene transfer, and environmental factors, such as nutrient availability and growth conditions, in evolution of the metabolic network. Furthermore, the mathematical models for the formation and adaptation of metabolic networks have also been described, according to the current understanding from a perspective of metabolic networks. These recent findings are helpful in not only understanding the formation of metabolic networks and their adaptation, but also metabolic engineering.

Evolution of a new function in an esterase: simple amino acid substitutions enable the activity present in the larger paralog, BioH

Gene duplication and divergence are essential processes for the evolution of new activities. Divergence may be gradual, involving simple amino acid residue substitutions, or drastic, such that larger structural elements are inserted, deleted or rearranged. Vast protein sequence comparisons, supported by some experimental evidence, argue that large structural modifications have been necessary for certain catalytic activities to evolve. However, it is not clear whether these activities could not have been attained by gradual changes. Interestingly, catalytic promiscuity could play a fundamental evolutionary role: a preexistent secondary activity could be increased by simple amino acid residue substitutions that do not affect the enzyme's primary activity. The promiscuous profile of the enzyme may be modified gradually by genetic drift, making a pool of potentially useful activities that can be selected before duplication. In this work, we used random mutagenesis and in vivo selection to evolve the Pseudomonas aeruginosa PAO1 carboxylesterase PA3859, a small protein, to attain the function of BioH, a much larger paralog involved in biotin biosynthesis. BioH was chosen as a target activity because it provides a highly sensitive selection for evolved enzymatic activities by auxotrophy complementation. After only two cycles of directed evolution, mutants with the ability to efficiently complement biotin auxotrophy were selected. The in vivo and in vitro characterization showed that the activity of one of our mutant proteins was similar to that of the wild-type BioH enzyme. Our results demonstrate that it is possible to evolve enzymatic activities present in larger proteins by discrete amino acid substitutions.

Out of fuzzy chemistry: from prebiotic chemistry to metabolic networks

The origin of life on Earth was a chemical affair. So how did primitive biochemical systems originate from geochemical and cosmochemical processes on the young planet? Contemporary research into the origins of life subscribes to the Darwinian principle of material causes operating in an evolutionary context, as advocated by A. I. Oparin and J. B. S. Haldane in the 1920s. In its simplest form (e.g., a bacterial cell) extant biological complexity relies on the functional integration of metabolic networks and replicative genomes inside a lipid boundary. Different research programmes have explored the prebiotic plausibility of each of these autocatalytic subsystems and combinations thereof: self-maintained networks of small molecules, template chemistry, and self-reproductive vesicles. This tutorial review focuses on the debates surrounding the origin of metabolism and offers a brief overview of current studies on the evolution of metabolic networks. I suggest that a leitmotif in the origin and evolution of metabolism is the role played by catalysers' substrate ambiguity and multifunctionality.

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.

Purine biosynthesis in archaea: variations on a theme

Background

The ability to perform de novo biosynthesis of purines is present in organisms in all three domains of life, reflecting the essentiality of these molecules to life. Although the pathway is quite similar in eukaryotes and bacteria, the archaeal pathway is more variable. A careful manual curation of genes in this pathway demonstrates the value of manual curation in archaea, even in pathways that have been well-studied in other domains.

Results

We searched the Integrated Microbial Genome system (IMG) for the 17 distinct genes involved in the 11 steps of de novo purine biosynthesis in 65 sequenced archaea, finding 738 predicted proteins with sequence similarity to known purine biosynthesis enzymes. Each sequence was manually inspected for the presence of active site residues and other residues known or suspected to be required for function.

Many apparently purine-biosynthesizing archaea lack evidence for a single enzyme, either glycinamide ribonucleotide formyltransferase or inosine monophosphate cyclohydrolase, suggesting that there are at least two more gene variants in the purine biosynthetic pathway to discover. Variations in domain arrangement of formylglycinamidine ribonucleotide synthetase and substantial problems in aminoimidazole carboxamide ribonucleotide formyltransferase and inosine monophosphate cyclohydrolase assignments were also identified.

Manual curation revealed some overly specific annotations in the IMG gene product name, with predicted proteins without essential active site residues assigned product names implying enzymatic activity (21 proteins, 2.8% of proteins inspected) or Enzyme Commission (E. C.) numbers (57 proteins, 7.7%). There were also 57 proteins (7.7%) assigned overly generic names and 78 proteins (10.6%) without E.C. numbers as part of the assigned name when a specific enzyme name and E. C. number were well-justified.

Conclusions

The patchy distribution of purine biosynthetic genes in archaea is consistent with a pathway that has been shaped by horizontal gene transfer, duplication, and gene loss. Our results indicate that manual curation can improve upon automated annotation for a small number of automatically-annotated proteins and can reveal a need to identify further pathway components even in well-studied pathways.

Reviewers

This article was reviewed by Dr. Céline Brochier-Armanet, Dr Kira S Makarova (nominated by Dr. Eugene Koonin), and Dr. Michael Galperin.

Three serendipitous pathways in E. coli can bypass a block in pyridoxal-5'-phosphate synthesis

Overexpression of seven different genes restores growth of a ΔpdxB strain of E. coli, which cannot make pyridoxal phosphate (PLP), on M9/glucose.

None of the enzymes encoded by these genes has a promiscuous 4-phosphoerythronate dehydrogenase activity that can replace the activity of PdxB.

Overexpression of these genes restores PLP synthesis by three different serendipitous pathways that feed into the normal PLP synthesis pathway downstream of the blocked step.

Reactions in one of these pathways are catalyzed by low-level activities of enzymes of unknown function and a promiscuous activity of an enzyme that normally has a role in another pathway; one reaction appears to be non-enzymatic.

Most metabolic enzymes are prodigious catalysts that have evolved to accelerate chemical reactions with high efficiency and specificity. However, many enzymes have inefficient promiscuous activities, as well, as a result of the assemblage of highly reactive catalytic residues and cofactors in active sites. Although promiscuous activities are generally orders of magnitude less efficient than well-evolved activities (O'Brien and Herschlag, 1998, 2001; Wang et al, 2003; Taylor Ringia et al, 2004), they often enhance reaction rates by orders of magnitude relative to those of uncatalyzed reactions (O'Brien and Herschlag, 1998, 2001). Thus, promiscuous activities provide a reservoir of novel catalytic activities that can be recruited to serve new functions.

The evolutionary potential of promiscuous enzymes extends beyond the recruitment of single enzymes to serve new functions. Microbes contain hundreds of enzymes—E. coli contains about 1700 (Freilich et al, 2005)—raising the possibility that promiscuous enzymes can be patched together to generate ‘serendipitous' pathways that are not part of normal metabolism. We distinguish serendipitous pathways from latent or cryptic pathways, which are bona fide pathways involving dedicated enzymes that are produced only under particular environmental circumstances. In contrast, serendipitous pathways are patched together from enzymes that normally serve other functions and are not regulated in a coordinated manner in response to the need to synthesize or degrade a metabolite.

In this study, we describe the discovery of three serendipitous pathways that allow synthesis of pyridoxal phosphate (PLP) in a strain of E. coli that lacks 4-phosphoerythronate dehydrogenase (PdxB) when one of the seven different genes is overexpressed. These genes were identified in a multicopy suppression experiment in which a library of E. coli genes (from the ASKA collection) was introduced into a ΔpdxB strain of E. coli that is unable to synthesize PLP. Surprisingly, none of the enzymes encoded by these genes has a promiscuous 4-phosphoerythronate (4PE) dehydrogenase activity that can substitute for the missing PdxB. Rather, overproduction of these enzymes appears to facilitate at least three serendipitous pathways that draw material from other metabolic pathways and feed into the normal PLP synthesis pathway downstream of the blocked step (Figure 1).

We have characterized one of these pathways in detail (Figure 3). The first step, dephosphorylation of 3-phosphohydroxypyruvate, is catalyzed by YeaB, a predicted NUDIX hydrolase of unknown function. Although catalysis is inefficient (k_cat=5.7×10⁻⁵ s⁻¹ and k_cat/K_M>0.028 M⁻¹ s⁻¹), the enzymatic rate is 4×10⁷-fold faster than the rate of the uncatalyzed reaction, and is sufficient to support PLP synthesis when YeaB is overproduced. The second step in the pathway is decarboxylation of 3-hydroxypyruvate (3HP). Although we found two enzymes (1-deoxyxylulose-5-phosphate synthase and the catalytic domain of α-ketoglutarate dehydrogenase) that catalyze this reaction with low but respectable activity in vitro, their involvement in pathway 1 was ruled out by genetic methods. Surprisingly, the non-enzymatic rate of decarboxylation of 3HP appears to be sufficient to support PLP synthesis. The third step in the pathway, condensation of glycolaldehyde and glycine to form 4-hydroxy-L-threonine, is catalyzed by LtaE, a low-specificity threonine aldolase whose physiological role is not known. The final step, phosphorylation of 4-hydroxy-L-threonine, is catalyzed by homoserine kinase (ThrB), which is required for synthesis of threonine. The promiscuous phosphorylation of 4-hydroxy-L-threonine is 80-fold slower than the physiological phosphorylation of homoserine. The involvement of LtaE and ThrB in pathway 1 was confirmed by genetic experiments showing that overexpression of yeaB no longer restores growth of ΔpdxB strains lacking either ltaE or thrB.

Although pathway 1 is inefficient, it provides the ΔpdxB strain with the ability to grow under conditions in which survival is otherwise impossible. In general, serendipitous assembly of an inefficient pathway from promiscuous activities of available enzymes will be important whenever the pathway provides increased fitness. This might occur when a critical metabolite is no longer available from the environment, and survival depends on assembly of a new biosynthetic pathway. A second circumstance in which an inefficient serendipitous pathway might improve fitness is the appearance of a novel compound in the environment that can be exploited as a source of carbon, nitrogen or phosphorous. Finally, chemotherapeutic agents that block metabolic pathways in bacteria or cancer cells could provide selective pressure for assembly of serendipitous pathways that allow synthesis of the end product of the blocked pathway and thus a previously unappreciated source of drug resistance. In all of these cases, even an inefficient pathway can provide a selective advantage over other cells in a particular environmental niche, allowing survival and subsequent mutations that elevate the efficiency of the pathway.

Our work is consistent with the hypothesis that the recognized metabolic network of E. coli is underlain by a denser network of reactions due to promiscuous enzymes that use and generate recognized metabolites, but also unusual metabolites that normally have no physiological role. The findings reported here highlight the abundance of cryptic capabilities in the E. coli proteome that can be drawn on to generate novel pathways. Such pathways could provide a starting place for assembly of more efficient pathways, both in nature and in the hands of metabolic engineers.

Bacterial genomes encode hundreds to thousands of enzymes, most of which are specialized for particular functions. However, most enzymes have inefficient promiscuous activities, as well, that generally serve no purpose. Promiscuous reactions can be patched together to form multistep metabolic pathways. Mutations that increase expression or activity of enzymes in such serendipitous pathways can elevate flux through the pathway to a physiologically significant level. In this study, we describe the discovery of three serendipitous pathways that allow synthesis of pyridoxal-5′-phosphate (PLP) in a strain of E. coli that lacks 4-phosphoerythronate (4PE) dehydrogenase (PdxB) when one of seven different genes is overexpressed. We have characterized one of these pathways in detail. This pathway diverts material from serine biosynthesis and generates an intermediate in the normal PLP synthesis pathway downstream of the block caused by lack of PdxB. Steps in the pathway are catalyzed by a protein of unknown function, a broad-specificity enzyme whose physiological role is unknown, and a promiscuous activity of an enzyme that normally serves another function. One step in the pathway may be non-enzymatic.

Experimental evolution, loss-of-function mutations, and "the first rule of adaptive evolution"

Adaptive evolution can cause a species to gain, lose, or modify a function; therefore, it is of basic interest to determine whether any of these modes dominates the evolutionary process under particular circumstances. Because mutation occurs at the molecular level, it is necessary to examine the molecular changes produced by the underlying mutation in order to assess whether a given adaptation is best considered as a gain, loss, or modification of function. Although that was once impossible, the advance of molecular biology in the past half century has made it feasible. In this paper, I review molecular changes underlying some adaptations, with a particular emphasis on evolutionary experiments with microbes conducted over the past four decades. I show that by far the most common adaptive changes seen in those examples are due to the loss or modification of a pre-existing molecular function, and I discuss the possible reasons for the prominence of such mutations.

Artificial gene amplification reveals an abundance of promiscuous resistance determinants in Escherichia coli

Duplicated genes provide an important raw material for adaptive evolution. However, the relationship between gene duplication and the emergence of new biochemical functions is complicated, and it has been difficult to quantify the likelihood of evolving novelty in any systematic manner. Here, we describe a comprehensive search for artificially amplified genes that are able to impart new phenotypes on Escherichia coli, provided their expression is up-regulated. We used a high-throughput, library-on-library strategy to screen for resistance to antibiotics and toxins. Cells containing a complete E. coli ORF library were exposed to 237 toxin-containing environments. From 86 of these environments, we identified a total of 115 cases where overexpressed ORFs imparted improved growth. Of the overexpressed ORFs that we tested, most conferred small but reproducible increases in minimum inhibitory concentration (≤16-fold) for their corresponding antibiotics. In many cases, proteins were acting promiscuously to impart resistance. In the absence of toxins, most strains bore no fitness cost associated with ORF overexpression. Our results show that even the genome of a nonpathogenic bacterium harbors a substantial reservoir of resistance genes, which can be readily accessed through overexpression mutations. During the growth of a population under selection, these mutations are most likely to be gene amplifications. Therefore, our work provides validation and biochemical insight into the innovation, amplification, and divergence model of gene evolution under continuous selection [Bergthorsson U, Andersson DI, Roth JR (2007) Proc Natl Acad Sci USA 104:17004-17009], and also illustrates the high frequency at which novel traits can evolve in bacterial populations.

De novo designed proteins from a library of artificial sequences function in Escherichia coli and enable cell growth

A central challenge of synthetic biology is to enable the growth of living systems using parts that are not derived from nature, but designed and synthesized in the laboratory. As an initial step toward achieving this goal, we probed the ability of a collection of >10⁶de novo designed proteins to provide biological functions necessary to sustain cell growth. Our collection of proteins was drawn from a combinatorial library of 102-residue sequences, designed by binary patterning of polar and nonpolar residues to fold into stable 4-helix bundles. We probed the capacity of proteins from this library to function in vivo by testing their abilities to rescue 27 different knockout strains of Escherichia coli, each deleted for a conditionally essential gene. Four different strains – ΔserB, ΔgltA, ΔilvA, and Δfes – were rescued by specific sequences from our library. Further experiments demonstrated that a strain simultaneously deleted for all four genes was rescued by co-expression of four novel sequences. Thus, cells deleted for ∼0.1% of the E. coli genome (and ∼1% of the genes required for growth under nutrient-poor conditions) can be sustained by sequences designed de novo.

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.

Identification of a suppressor gene for the arginine-auxotrophic argJ mutation in Corynebacterium glutamicum

We recently proposed a metabolic engineering strategy for L-ornithine production based on the hypothesis that an increased intracellular supply of N-acetylglutamate may further enhance L-ornithine production in a well-defined recombinant strain of Corynebacterium glutamicum. In this work, an argJ-deficient arginine auxotrophic mutant of C. glutamicum is suppressed by a different locus of C. glutamicum ATCC13032. Overexpression of the NCgl1469 open reading frame (ORF), exhibiting N-acetylglutamate synthase (NAGS) activity, was able to complement the C. glutamicum arginine-auxotrophic argJ strain and showed increased NAGS activity from 0.03 to 0.17 units mg(-1) protein. Additionally, overexpression of the NCgl1469 ORF resulted in a 39% increase in excreted L-ornithine. These results indicate that the intracellular supply of N-acetylglutamate is a rate-limiting step during L-ornithine production in C. glutamicum.

Gene amplification and adaptive evolution in bacteria

Gene duplication-amplification (GDA) processes are highly relevant biologically because they generate extensive and reversible genetic variation on which adaptive evolution can act. Whenever cellular growth is restricted, escape from these growth restrictions often occurs by GDA events that resolve the selective problem. In addition, GDA may facilitate subsequent genetic change by allowing a population to grow and increase in number, thereby increasing the probability for subsequent adaptive mutations to occur in the amplified genes or in unrelated genes. Mathematical modeling of the effect of GDA on the rate of adaptive evolution shows that GDA will facilitate adaptation, especially when the supply of mutations in the population is rate-limiting. GDA can form via several mechanisms, both RecA-dependent and RecA-independent, including rolling-circle amplification and nonequal crossing over between sister chromatids. Due to the high intrinsic instability and fitness costs associated with GDAs, they are generally transient in nature, and consequently their evolutionary and medical importance is often underestimated.

Bacterial gene amplification: implications for the evolution of antibiotic resistance

Recent data suggest that, in response to the presence of antibiotics, gene duplication and amplification (GDA) constitutes an important adaptive mechanism in bacteria. For example, resistance to sulphonamide, trimethoprim and beta-lactams can be conferred by increased gene dosage through GDA of antibiotic hydrolytic enzymes, target enzymes or efflux pumps. Furthermore, most types of antibiotic resistance mechanism are deleterious in the absence of antibiotics, and these fitness costs can be ameliorated by increased gene dosage of limiting functions. In this Review, we highlight the dynamic properties of gene amplifications and describe how they can facilitate adaptive evolution in response to toxic drugs.

Genesis, effects and fates of repeats in prokaryotic genomes

DNA repeats are causes and consequences of genome plasticity. Repeats are created by intrachromosomal recombination or horizontal transfer. They are targeted by recombination processes leading to amplifications, deletions and rearrangements of genetic material. The identification and analysis of repeats in nearly 700 genomes of bacteria and archaea is facilitated by the existence of sequence data and adequate bioinformatic tools. These have revealed the immense diversity of repeats in genomes, from those created by selfish elements to the ones used for protection against selfish elements, from those arising from transient gene amplifications to the ones leading to stable duplications. Experimental works have shown that some repeats do not carry any adaptive value, while others allow functional diversification and increased expression. All repeats carry some potential to disorganize and destabilize genomes. Because recombination and selection for repeats vary between genomes, the number and types of repeats are also quite diverse and in line with ecological variables, such as host-dependent associations or population sizes, and with genetic variables, such as the recombination machinery. From an evolutionary point of view, repeats represent both opportunities and problems. We describe how repeats are created and how they can be found in genomes. We then focus on the functional and genomic consequences of repeats that dictate their fate.

Mapping catalytic promiscuity in the alkaline phosphatase superfamily

"Promiscuous" enzymes possess activities in addition to their native ones. Promiscuous activities could be remnants from an evolutionary ancestor that has been adapted to fulfil a new function following gene duplication. Alternatively, the observation of promiscuity could indicate that an enzyme has the potential to evolve into a new catalyst. Thus, the observation of promiscuity defines functional relationships in enzyme superfamilies. Crosswise promiscuity can provide an additional layer of connectivity between members of a - usually structurally defined - superfamily to establish a system for tracking the emergence and interconversion of enzymatic function. The systematic analysis of measured promiscuous rates may serve as a basis for drawing up phylogenetic relationships based on the potential for catalysis and may be useful for active use in directed evolution, suggesting evolutionary "short cuts". We review recent observations of catalytic promiscuity in members of the alkaline phosphatase (AP) superfamily that exhibit reciprocal relationships of crosswise promiscuity with rate accelerations (kcat/KM)/k2 between 106 and 1018. Specifically, we focus on the mechanistic features that appear to form the basis of catalytic promiscuity in this superfamily.

A study in molecular contingency: glutamine phosphoribosylpyrophosphate amidotransferase is a promiscuous and evolvable phosphoribosylanthranilate isomerase

The prevalence of paralogous enzymes implies that novel catalytic functions can evolve on preexisting protein scaffolds. The weak secondary activities of proteins, which reflect catalytic promiscuity and substrate ambiguity, are plausible starting points for this evolutionary process. In this study, we observed the emergence of a new enzyme from the ASKA (A Complete Set of E. coli K-12 ORF Archive) collection of Escherichia coli open reading frames. The overexpression of (His)(6)-tagged glutamine phosphoribosylpyrophosphate amidotransferase (PurF) unexpectedly rescued a Delta trpF E. coli strain from starvation on minimal media. The wild-type PurF and TrpF enzymes are unrelated in sequence, tertiary structure and catalytic mechanism. The promiscuous phosphoribosylanthranilate isomerase activity of the ASKA PurF variant apparently stems from a preexisting affinity for phosphoribosylated substrates. The relative fitness of the (His)(6)-PurF/Delta trpF strain was improved 4.8-fold to nearly wild-type levels by random mutagenesis of purF and genetic selection. The evolved and ancestral PurF proteins were purified and reacted with phosphoribosylanthranilate in vitro. The best evolvant (k(cat)/K(M)=0.3 s(-1) M(-1)) was approximately 25-fold more efficient than its ancestor but >10(7)-fold less efficient than the wild-type phosphoribosylanthranilate isomerase. These observations demonstrate in quantitative terms that the weak secondary activities of promiscuous enzymes can dramatically improve the fitness of contemporary organisms.

Sensitive genome-wide screen for low secondary enzymatic activities: the YjbQ family shows thiamin phosphate synthase activity

Contemporary enzymes are highly efficient and selective catalysts. However, due to the intrinsically very reactive nature of active sites, gratuitous secondary reactions are practically unavoidable. Consequently, even the smallest cell, with its limited enzymatic repertoire, has the potential to carry out numerous additional, very likely inefficient, secondary reactions. If selectively advantageous, secondary reactions could be the basis for the evolution of new fully functional enzymes. Here, we investigated if Escherichia coli has cryptic enzymatic activities related to thiamin biosynthesis. We selected this pathway because this vitamin is essential, but the cell's requirements are very small. Therefore, enzymes with very low activity could complement the auxotrophy of strains deleted of some bona fide thiamin biosynthetic genes. By overexpressing the E. coli's protein repertoire, we selected yjbQ, a gene that complemented a strain deleted of the thiamin phosphate synthase (TPS)-coding gene thiE. In vitro studies confirmed TPS activity, and by directed evolution experiments, this activity was enhanced. Structurally oriented mutagenesis allowed us to identify the putative active site. Remote orthologs of YjbQ from Thermotoga, Sulfolobus, and Pyrococcus were cloned and also showed thiamin auxotrophy complementation, indicating that the cryptic TPS activity is a property of this protein family. Interestingly, the thiE- and yjbQ-coded TPSs are analog enzymes with no structural similarity, reflecting distinct evolutionary origin. These results support the hypothesis that the enzymatic repertoire of a cell such as E. coli has the potential to perform vast amounts of alternative reactions, which could be exploited to evolve novel or more efficient catalysts.