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

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.

Coevolution of Metabolic Pathways in Blattodea and Their Blattabacterium Endosymbionts, and Comparisons with Other Insect-Bacteria Symbioses

Many insects harbor bacterial endosymbionts that supply essential nutrients and enable their hosts to thrive on a nutritionally unbalanced diet. Comparisons of the genomes of endosymbionts and their insect hosts have revealed multiple cases of mutually-dependent metabolic pathways that require enzymes encoded in 2 genomes. Complementation of metabolic reactions at the pathway level has been described for hosts feeding on unbalanced diets, such as plant sap. However, the level of collaboration between symbionts and hosts that feed on more variable diets is largely unknown. In this study, we investigated amino acid and vitamin/cofactor biosynthetic pathways in Blattodea, which comprises cockroaches and termites, and their obligate endosymbiont Blattabacterium cuenoti (hereafter Blattabacterium). In contrast to other obligate symbiotic systems, we found no clear evidence of "collaborative pathways" for amino acid biosynthesis in the genomes of these taxa, with the exception of collaborative arginine biosynthesis in 2 taxa, Cryptocercus punctulatus and Mastotermes darwiniensis. Nevertheless, we found that several gaps specific to Blattabacterium in the folate biosynthetic pathway are likely to be complemented by their host. Comparisons with other insects revealed that, with the exception of the arginine biosynthetic pathway, collaborative pathways for essential amino acids are only observed in phloem-sap feeders. These results suggest that the host diet is an important driving factor of metabolic pathway evolution in obligate symbiotic systems. IMPORTANCE The long-term coevolution between insects and their obligate endosymbionts is accompanied by increasing levels of genome integration, sometimes to the point that metabolic pathways require enzymes encoded in two genomes, which we refer to as "collaborative pathways". To date, collaborative pathways have only been reported from sap-feeding insects. Here, we examined metabolic interactions between cockroaches, a group of detritivorous insects, and their obligate endosymbiont, Blattabacterium, and only found evidence of collaborative pathways for arginine biosynthesis. The rarity of collaborative pathways in cockroaches and Blattabacterium contrasts with their prevalence in insect hosts feeding on phloem-sap. Our results suggest that host diet is a factor affecting metabolic integration in obligate symbiotic systems.

Low protein expression enhances phenotypic evolvability by intensifying selection on folding stability

Protein abundance affects the evolution of protein genotypes, but we do not know how it affects the evolution of protein phenotypes. Here we investigate the role of protein abundance in the evolvability of green fluorescent protein (GFP) towards the novel phenotype of cyan fluorescence. We evolve GFP in E. coli through multiple cycles of mutation and selection and show that low GFP expression facilitates the evolution of cyan fluorescence. A computational model whose predictions we test experimentally helps explain why: lowly expressed proteins are under stronger selection for proper folding, which facilitates their evolvability on short evolutionary time scales. The reason is that high fluorescence can be achieved by either few proteins that fold well or by many proteins that fold less well. In other words, we observe a synergy between a protein's scarcity and its stability. Because many proteins meet the essential requirements for this scarcity-stability synergy, it may be a widespread mechanism by which low expression helps proteins evolve new phenotypes and functions.

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.

Bacterial rhomboid proteases mediate quality control of orphan membrane proteins

Although multiprotein membrane complexes play crucial roles in bacterial physiology and virulence, the mechanisms governing their quality control remain incompletely understood. In particular, it is not known how unincorporated, orphan components of protein complexes are recognised and eliminated from membranes. Rhomboids, the most widespread and largest superfamily of intramembrane proteases, are known to play key roles in eukaryotes. In contrast, the function of prokaryotic rhomboids has remained enigmatic. Here, we show that the Shigella sonnei rhomboid proteases GlpG and the newly identified Rhom7 are involved in membrane protein quality control by specifically targeting components of respiratory complexes, with the metastable transmembrane domains (TMDs) of rhomboid substrates protected when they are incorporated into a functional complex. Initial cleavage by GlpG or Rhom7 allows subsequent degradation of the orphan substrate. Given the occurrence of this strategy in an evolutionary ancient organism and the presence of rhomboids in all domains of life, it is likely that this form of quality control also mediates critical events in eukaryotes and protects cells from the damaging effects of orphan proteins.

In Vitro Selection of Peptides and Proteins-Advantages of mRNA Display

mRNA display is a robust in vitro selection technique that allows the selection of peptides and proteins with desired functions from libraries of trillions of variants. mRNA display relies upon a covalent linkage between a protein and its encoding mRNA molecule; the power of the technique stems from the stability of this link, and the large degree of control over experimental conditions afforded to the researcher. This article describes the major advantages that make mRNA display the method of choice among comparable in vivo and in vitro methods, including cell-surface display, phage display, and ribosomal display. We also describe innovative techniques that harness mRNA display for directed evolution, protein engineering, and drug discovery.

Library Selection with a Randomized Repertoire of (βα)8-Barrel Enzymes Results in Unexpected Induction of Gene Expression

The potential of the frequently encountered (βα)₈-barrel fold to acquire new functions was tested by an approach combining random mutagenesis and selection in vivo. For this purpose, the genes encoding 52 different phosphate-binding (βα)₈-barrel proteins were subjected to error-prone PCR and cloned into an expression plasmid. The resulting mixed repertoire was used to transform different auxotrophic Escherichia coli strains, each lacking an enzyme with a phosphate-containing substrate. After plating of the different transformants on minimal medium, growth was observed only for two strains, lacking either the gene for the serine phosphatase SerB or the phosphoserine aminotransferase SerC. The same mutants of the E. coli genes nanE (encoding a putative N-acetylmannosamine-6-phosphate 2-epimerase) and pdxJ (encoding the pyridoxine 5'-phosphate synthase) were responsible for rescuing both ΔserB and ΔserC. Unexpectedly, the complementing NanE and PdxJ variants did not catalyze the SerB or SerC reactions in vitro. Instead, RT-qPCR, RNAseq, and transcriptome analysis showed that they rescue the deletions by enlisting the help of endogenous E. coli enzymes HisB and HisC through exclusive up-regulation of histidine operon transcription. While the promiscuous SerB activity of HisB is well-established, our data indicate that HisC is promiscuous for the SerC reaction, as well. The successful rescue of ΔserB and ΔserC through point mutations and recruitment of additional amino acids in NanE and PdxJ provides another example for the adaptability of the (βα)₈-barrel fold.

Cryptic genetic variation shapes the adaptive evolutionary potential of enzymes

Genetic variation among orthologous proteins can cause cryptic phenotypic properties that only manifest in changing environments. Such variation may impact the evolvability of proteins, but the underlying molecular basis remains unclear. Here, we performed comparative directed evolution of four orthologous metallo-β-lactamases toward a new function and found that different starting genotypes evolved to distinct evolutionary outcomes. Despite a low initial fitness, one ortholog reached a significantly higher fitness plateau than its counterparts, via increasing catalytic activity. By contrast, the ortholog with the highest initial activity evolved to a less-optimal and phenotypically distinct outcome through changes in expression, oligomerization and activity. We show how cryptic molecular properties and conformational variation of active site residues in the initial genotypes cause epistasis, that could lead to distinct evolutionary outcomes. Our work highlights the importance of understanding the molecular details that connect genetic variation to protein function to improve the prediction of protein evolution.

Reframing gene essentiality in terms of adaptive flexibility

BACKGROUND

Essentiality assays are important tools commonly utilized for the discovery of gene functions. Growth/no growth screens of single gene knockout strain collections are also often utilized to test the predictive power of genome-scale models. False positive predictions occur when computational analysis predicts a gene to be non-essential, however experimental screens deem the gene to be essential. One explanation for this inconsistency is that the model contains the wrong information, possibly an incorrectly annotated alternative pathway or isozyme reaction. Inconsistencies could also be attributed to experimental limitations, such as growth tests with arbitrary time cut-offs. The focus of this study was to resolve such inconsistencies to better understand isozyme activities and gene essentiality.

RESULTS

In this study, we explored the definition of conditional essentiality from a phenotypic and genomic perspective. Gene-deletion strains associated with false positive predictions of gene essentiality on defined minimal medium for Escherichia coli were targeted for extended growth tests followed by population sequencing and transcriptome analysis. Of the twenty false positive strains available and confirmed from the Keio single gene knock-out collection, 11 strains were shown to grow with longer incubation periods making these actual true positives. These strains grew reproducibly with a diverse range of growth phenotypes. The lag phase observed for these strains ranged from less than one day to more than 7 days. It was found that 9 out of 11 of the false positive strains that grew acquired mutations in at least one replicate experiment and the types of mutations ranged from SNPs and small indels associated with regulatory or metabolic elements to large regions of genome duplication. Comparison of the detected adaptive mutations, modeling predictions of alternate pathways and isozymes, and transcriptome analysis of KO strains suggested agreement for the observed growth phenotype for 6 out of the 9 cases where mutations were observed.

CONCLUSIONS

Longer-term growth experiments followed by whole genome sequencing and transcriptome analysis can provide a better understanding of conditional gene essentiality and mechanisms of adaptation to such perturbations. Compensatory mutations are largely reproducible mechanisms and are in agreement with genome-scale modeling predictions to loss of function gene deletion events.

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.

Determinants and Prediction of Esterase Substrate Promiscuity Patterns

Esterases receive special attention because of their wide distribution in biological systems and environments and their importance for physiology and chemical synthesis. The prediction of esterases' substrate promiscuity level from sequence data and the molecular reasons why certain such enzymes are more promiscuous than others remain to be elucidated. This limits the surveillance of the sequence space for esterases potentially leading to new versatile biocatalysts and new insights into their role in cellular function. Here, we performed an extensive analysis of the substrate spectra of 145 phylogenetically and environmentally diverse microbial esterases, when tested with 96 diverse esters. We determined the primary factors shaping their substrate range by analyzing substrate range patterns in combination with structural analysis and protein-ligand simulations. We found a structural parameter that helps rank (classify) the promiscuity level of esterases from sequence data at 94% accuracy. This parameter, the active site effective volume, exemplifies the topology of the catalytic environment by measuring the active site cavity volume corrected by the relative solvent accessible surface area (SASA) of the catalytic triad. Sequences encoding esterases with active site effective volumes (cavity volume/SASA) above a threshold show greater substrate spectra, which can be further extended in combination with phylogenetic data. This measure provides also a valuable tool for interrogating substrates capable of being converted. This measure, found to be transferred to phosphatases of the haloalkanoic acid dehalogenase superfamily and possibly other enzymatic systems, represents a powerful tool for low-cost bioprospecting for esterases with broad substrate ranges, in large scale sequence data sets.

Enzyme engineering: reaching the maximal catalytic efficiency peak

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

Refined experimental annotation reveals conserved corrinoid autotrophy in chloroform-respiring Dehalobacter isolates

Two novel chlorinated alkane-respiring Dehalobacter restrictus strains CF and DCA were isolated from the same enrichment culture, ACT-3, and characterized. The closed genomes of these highly similar sister strains were previously assembled from metagenomic sequence data and annotated. The isolation of the strains enabled experimental verification of predicted annotations, particularly focusing on irregularities or predicted gaps in central metabolic pathways and cofactor biosynthesis. Similar to D. restrictus strain PER-K23, strains CF and DCA require arginine, histidine and threonine for growth, although the corresponding biosynthesis pathways are predicted to be functional. Using strain CF to experimentally verify annotations, we determined that the predicted defective serine biosynthesis pathway can be rescued with a promiscuous serine hydroxymethyltransferase. Strain CF grew without added thiamine although the thiamine biosynthesis pathway is predicted to be absent; intracellular thiamine diphosphate, the cofactor of carboxylases in central metabolism, was not detected in cell extracts. Thus, strain CF may use amino acids to replenish central metabolites, portending entangled metabolite exchanges in ACT-3. Consistent with annotation, strain CF possesses a functional corrinoid biosynthesis pathway, demonstrated by increasing corrinoid content during growth and guided cobalamin biosynthesis in corrinoid-free medium. Chloroform toxicity to corrinoid-producing methanogens and acetogens may drive the conservation of corrinoid autotrophy in Dehalobacter strains. Heme detection in strain CF cell extracts suggests the 'archaeal' heme biosynthesis pathway also functions in anaerobic Firmicutes. This study reinforces the importance of incorporating enzyme promiscuity and cofactor availability in genome-scale functional predictions and identifies essential nutrient interdependencies in anaerobic dechlorinating microbial communities.

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.

A protein constructed de novo enables cell growth by altering gene regulation

Recent advances in protein design rely on rational and computational approaches to create novel sequences that fold and function. In contrast, natural systems selected functional proteins without any design a priori. In an attempt to mimic nature, we used large libraries of novel sequences and selected for functional proteins that rescue Escherichia coli cells in which a conditionally essential gene has been deleted. In this way, the de novo protein SynSerB3 was selected as a rescuer of cells in which serB, which encodes phosphoserine phosphatase, an enzyme essential for serine biosynthesis, was deleted. However, SynSerB3 does not rescue the deleted activity by catalyzing hydrolysis of phosphoserine. Instead, SynSerB3 up-regulates hisB, a gene encoding histidinol phosphate phosphatase. This endogenous E. coli phosphatase has promiscuous activity that, when overexpressed, compensates for the deletion of phosphoserine phosphatase. Thus, the de novo protein SynSerB3 rescues the deletion of serB by altering the natural regulation of the His operon.

Evolution of new functions de novo and from preexisting genes

How the enormous structural and functional diversity of new genes and proteins was generated (estimated to be 10(10)-10(12) different proteins in all organisms on earth [Choi I-G, Kim S-H. 2006. Evolution of protein structural classes and protein sequence families. Proc Natl Acad Sci 103: 14056-14061] is a central biological question that has a long and rich history. Extensive work during the last 80 years have shown that new genes that play important roles in lineage-specific phenotypes and adaptation can originate through a multitude of different mechanisms, including duplication, lateral gene transfer, gene fusion/fission, and de novo origination. In this review, we focus on two main processes as generators of new functions: evolution of new genes by duplication and divergence of pre-existing genes and de novo gene origination in which a whole protein-coding gene evolves from a noncoding sequence.

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.

A substrate ambiguous enzyme facilitates genome reduction in an intracellular symbiont

BACKGROUND

Genome evolution in intracellular microbial symbionts is characterized by gene loss, generating some of the smallest and most gene-poor genomes known. As a result of gene loss these genomes commonly contain metabolic pathways that are fragmented relative to their free-living relatives. The evolutionary retention of fragmented metabolic pathways in the gene-poor genomes of endosymbionts suggests that they are functional. However, it is not always clear how they maintain functionality. To date, the fragmented metabolic pathways of endosymbionts have been shown to maintain functionality through complementation by host genes, complementation by genes of another endosymbiont and complementation by genes in host genomes that have been horizontally acquired from a microbial source that is not the endosymbiont. Here, we demonstrate a fourth mechanism.

RESULTS

We investigate the evolutionary retention of a fragmented pathway for the essential nutrient pantothenate (vitamin B5) in the pea aphid, Acyrthosiphon pisum endosymbiosis with Buchnera aphidicola. Using quantitative analysis of gene expression we present evidence for complementation of the Buchnera pantothenate biosynthesis pathway by host genes. Further, using complementation assays in an Escherichia coli mutant we demonstrate functional replacement of a pantothenate biosynthesis enzyme, 2-dehydropantoate 2-reductase (E.C. 1.1.1.169), by an endosymbiont gene, ilvC, encoding a substrate ambiguous enzyme.

CONCLUSIONS

Earlier studies have speculated that missing enzyme steps in fragmented endosymbiont metabolic pathways are completed by adaptable endosymbiont enzymes from other pathways. Here, we experimentally demonstrate completion of a fragmented endosymbiont vitamin biosynthesis pathway by recruitment of a substrate ambiguous enzyme from another pathway. In addition, this work extends host/symbiont metabolic collaboration in the aphid/Buchnera symbiosis from amino acid metabolism to include vitamin biosynthesis.

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.

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.

The logic of metabolism and its fuzzy consequences

Intermediary metabolism molecules are orchestrated into logical pathways stemming from history (L-amino acids, D-sugars) and dynamic constraints (hydrolysis of pyrophosphate or amide groups is the driving force of anabolism). Beside essential metabolites, numerous variants derive from programmed or accidental changes. Broken down, variants enter standard pathways, producing further variants. Macromolecule modification alters enzyme reactions specificity. Metabolism conform thermodynamic laws, precluding strict accuracy. Hence, for each regular pathway, a wealth of variants inputs and produces metabolites that are similar to but not the exact replicas of core metabolites. As corollary, a shadow, paralogous metabolism, is associated to standard metabolism. We focus on a logic of paralogous metabolism based on diversion of the core metabolic mimics into pathways where they are modified to minimize their input in the core pathways where they create havoc. We propose that a significant proportion of paralogues of well-characterized enzymes have evolved as the natural way to cope with paralogous metabolites. A second type of denouement uses a process where protecting/deprotecting unwanted metabolites - conceptually similar to the procedure used in the laboratory of an organic chemist - is used to enter a completely new catabolic pathway.