Toward a systems biology perspective on enzyme evolution

Mimicomes: Mimicking Multienzyme System by Artificial Design

Enzymes are widely distributed in organelles of cells, which are capable of carrying out specific catalytic reactions. In general, several enzymes collaborate to facilitate complex reactions and engage in vital biochemical processes within cells, which are also called cascade systems. The cascade systems are highly efficient, and their dysfunction is associated with a multitude of endogenous diseases. The advent of nanotechnology makes it possible to mimic these cascade systems in nature and realize partial functions of natural biological processes both in vitro and in vivo. To emphasize the significance of artificial cascade systems, mimicomes is first proposed, a new concept that refers to the artificial cascade catalytic systems. Typically, mimicomes are able to mimic specific natural biochemical catalytic processes or facilitate the overall catalytic efficiency of cascade systems. Subsequently, the evolution and development of different types of mimicomes in recent decades are elucidated exhaustedly, from the natural enzyme-based mimicomes (immobilized enzyme and vesicle mimicomes) to the nanozyme-based mimicomes and enzyme-nanozyme hybrid mimicomes. In conclusion, the remaining challenges in the design of multifunctional mimicomes and their potential applications are summarized, offering insights into their future prospects.

High-throughput selection of (new) enzymes: phage display-mediated isolation of alkyl halide hydrolases from a library of active-site mutated epoxide hydrolases

Epoxide hydrolase StEH1, from potato, is similar in overall structural fold and catalytic mechanism to haloalkane dehalogenase DhlA from Xanthobacter autotrophicus. StEH1 displays low (promiscuous) hydrolytic activity with (2-chloro)- and (2-bromo)ethanebenzene producing 2-phenylethanol. To investigate possibilities to amplify these very low dehalogenase activities, StEH1 was subjected to targeted randomized mutagenesis at five active-site amino acid residues and the resulting protein library was challenged for reactivity towards a bait chloride substrate. Enzymes catalyzing the first half-reaction of a hydrolytic cycle were isolated following monovalent phage display of the mutated proteins. Several StEH1 derived enzymes were identified with enhanced dehalogenase activities.

Specialization restricts the evolutionary paths available to yeast sugar transporters

Functional innovation at the protein level is a key source of evolutionary novelties. The constraints on functional innovations are likely to be highly specific in different proteins, which are shaped by their unique histories and the extent of global epistasis that arises from their structures and biochemistries. These contextual nuances in the sequence-function relationship have implications both for a basic understanding of the evolutionary process and for engineering proteins with desirable properties. Here, we have investigated the molecular basis of novel function in a model member of an ancient, conserved, and biotechnologically relevant protein family. These Major Facilitator Superfamily sugar porters are a functionally diverse group of proteins that are thought to be highly plastic and evolvable. By dissecting a recent evolutionary innovation in an α-glucoside transporter from the yeast Saccharomyces eubayanus, we show that the ability to transport a novel substrate requires high-order interactions between many protein regions and numerous specific residues proximal to the transport channel. To reconcile the functional diversity of this family with the constrained evolution of this model protein, we generated new, state-of-the-art genome annotations for 332 Saccharomycotina yeast species spanning approximately 400 million years of evolution. By integrating phylogenetic and phenotypic analyses across these species, we show that the model yeast α-glucoside transporters likely evolved from a multifunctional ancestor and became subfunctionalized. The accumulation of additive and epistatic substitutions likely entrenched this subfunction, which made the simultaneous acquisition of multiple interacting substitutions the only reasonably accessible path to novelty.

The mechanistic insights into different aspects of promiscuity in metalloenzymes

Enzymes are nature's ultimate machinery to catalyze complex reactions. Though enzymes are evolved to catalyze specific reactions, they also show significant promiscuity in reactions and substrate selection. Metalloenzymes contain a metal ion or metal cofactor in their active site, which is crucial in their catalytic activity. Depending on the metal and its coordination environment, the metal ion or cofactor may function as a Lewis acid or base and a redox center and thus can catalyze a plethora of natural reactions. In fact, the versatility in the oxidation state of the metal ions provides metalloenzymes with a high level of catalytic adaptability and promiscuity. In this chapter, we discuss different aspects of promiscuity in metalloenzymes by using several recent experimental and theoretical works as case studies. We start our discussion by introducing the concept of promiscuity and then we delve into the mechanistic insight into promiscuity at the molecular level.

Expanded Substrate Specificity in D-Amino Acid Transaminases: A Case Study of Transaminase from Blastococcus saxobsidens

Enzymes with expanded substrate specificity are good starting points for the design of biocatalysts for target reactions. However, the structural basis of the expanded substrate specificity is still elusive, especially in the superfamily of pyridoxal-5'-phosphate-dependent transaminases, which are characterized by a conserved organization of both the active site and functional dimer. Here, we analyze the structure-function relationships in a non-canonical D-amino acid transaminase from Blastococcus saxobsidens, which is active towards D-amino acids and primary (R)-amines. A detailed study of the enzyme includes a kinetic analysis of its substrate scope and a structural analysis of the holoenzyme and its complex with phenylhydrazine-a reversible inhibitor and analogue of (R)-1-phenylethylamine-a benchmark substrate of (R)-selective amine transaminases. We suggest that the features of the active site of transaminase from B. saxobsidens, such as the flexibility of the R34 and R96 residues, the lack of bulky residues in the β-turn at the entrance to the active site, and the short O-pocket loop, facilitate the binding of substrates with and without α-carboxylate groups. The proposed structural determinants of the expanded substrate specificity can be used for the design of transaminases for the stereoselective amination of keto compounds.

Recruitment of a Middling Promiscuous Enzyme Drives Adaptive Metabolic Evolution in Escherichia coli

A key step in metabolic pathway evolution is the recruitment of promiscuous enzymes to perform new functions. Despite the recognition that promiscuity is widespread in biology, factors dictating the preferential recruitment of one promiscuous enzyme over other candidates are unknown. Escherichia coli contains four sugar kinases that are candidates for recruitment when the native glucokinase machinery is deleted-allokinase (AlsK), manno(fructo)kinase (Mak), N-acetylmannosamine kinase (NanK), and N-acetylglucosamine kinase (NagK). The catalytic efficiencies of these enzymes are 103- to 105-fold lower than native glucokinases, ranging from 2,400 M-1 s-1 for the most active candidate, NagK, to 15 M-1 s-1 for the least active candidate, AlsK. To investigate the relationship between catalytic activities of promiscuous enzymes and their recruitment, we performed adaptive evolution of a glucokinase-deficient E. coli strain to restore glycolytic metabolism. We observed preferential recruitment of NanK via a trajectory involving early mutations that facilitate glucose uptake and amplify nanK transcription, followed by nonsynonymous substitutions in NanK that enhance the enzyme's promiscuous glucokinase activity. These substitutions reduced the native activity of NanK and reduced organismal fitness during growth on an N-acetylated carbon source, indicating that enzyme recruitment comes at a cost for growth on other substrates. Notably, the two most active candidates, NagK and Mak, were not recruited, suggesting that catalytic activity alone does not dictate evolutionary outcomes. The results highlight our lack of knowledge regarding biological drivers of enzyme recruitment and emphasize the need for a systems-wide approach to identify factors facilitating or constraining this important adaptive process.

Molecular Typing of fumC, icd, and mdh Genes in Serratia Marcescens

Aim

Serratia marcescens genes fumC, icd, and mdh were molecularly typed in various groups of 200 clinical samples.

Results

According to the findings, 38 (19%) of the isolates are Serratia marcescens. All these bacterial isolates had their DNA extracted. Then, using particular primers, the genes fumC, icd, and mdh are detected and amplified. These genes were sequenced, and the results were aligned with NCBI sequences. Using the Geneious version 9 software, gene sequences were analyzed. Sequencing of these genes revealed variant regions when compared to global isolates in NCBI. Energy levels in bacterial cells may be impacted by TCA cycle enzyme variant sequence genes.

Conclusion

The bacterial sequences from Iraq that were listed in NCBI with an accession number were LC735551 is a gene bank. (Gene Bank: LC735550, 1 Iraqi 40 fumC gene. One Iraq 41 icd gene; gene accession number: LC735549.42 mdh gene in Iraq.).

Supramolecular enzyme-mimicking catalysts self-assembled from peptides

Natural enzymes catalyze biochemical transformations in superior catalytic efficiency and remarkable substrate specificity. The excellent catalytic repertoire of enzymes is attributed to the sophisticated chemical structures of their active sites, as a result of billions-of-years natural evolution. However, large-scale practical applications of natural enzymes are restricted due to their poor stability, difficulty in modification, and high costs of production. One viable solution is to fabricate supramolecular catalysts with enzyme-mimetic active sites. In this review, we introduce the principles and strategies of designing peptide-based artificial enzymes which display catalytic activities similar to those of natural enzymes, such as aldolases, laccases, peroxidases, and hydrolases (mainly the esterases and phosphatases). We also discuss some multifunctional enzyme-mimicking systems which are capable of catalyzing orthogonal or cascade reactions. We highlight the relationship between structures of enzyme-like active sites and the catalytic properties, as well as the significance of these studies from an evolutionary point of view.

Comparative genomics of Leuconostoc lactis strains isolated from human gastrointestinal system and fermented foods microbiomes

Background

Leuconostoc lactis forms a crucial member of the genus Leuconostoc and has been widely used in the fermentation industry to convert raw material into acidified and flavored products in dairy and plant-based food systems. Since the ecological niches that strains of Ln. lactis being isolated from were truly diverse such as the human gut, dairy, and plant environments, comparative genome analysis studies are needed to better understand the strain differences from a metabolic adaptation point of view across diverse sources of origin. We compared eight Ln. lactis strains of 1.2.28, aa_0143, BIOML-A1, CBA3625, LN19, LN24, WIKIM21, and WiKim40 using bioinformatics to elucidate genomic level characteristics of each strain for better utilization of this species in a broad range of applications in food industry.

Results

Phylogenomic analysis of twenty-nine Ln. lactis strains resulted in nine clades. Whole-genome sequence analysis was performed on eight Ln. lactis strains representing human gastrointestinal tract and fermented foods microbiomes. The findings of the present study are based on comparative genome analysis against the reference Ln. lactis CBA3625 genome. Overall, a ~ 41% of all CDS were conserved between all strains. When the coding sequences were assigned to a function, mobile genetic elements, mainly insertion sequences were carried by all eight strains. All strains except LN24 and WiKim40 harbor at least one intact putative prophage region, and two of the strains contained CRISPR-Cas system. All strains encoded Lactococcin 972 bacteriocin biosynthesis gene clusters except for CBA3625.

Conclusions

The findings in the present study put forth new perspectives on genomics of Ln. lactis via complete genome sequence based comparative analysis and further determination of genomic characteristics. The outcomes of this work could potentially pave the way for developing elements for future strain engineering applications.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01074-6.

Human cytosolic transaminases: side activities and patterns of discrimination towards physiologically available alternative substrates

Transaminases play key roles in central metabolism, transferring the amino group from a donor substrate to an acceptor. These enzymes can often act, with low efficiency, on compounds different from the preferred substrates. To understand what might have shaped the substrate specificity of this class of enzymes, we examined the reactivity of six human cytosolic transaminases towards amino acids whose main degradative pathways do not include any transamination. We also tested whether sugars and sugar phosphates could serve as alternative amino group acceptors for these cytosolic enzymes. Each of the six aminotransferases reacted appreciably with at least three of the alternative amino acid substrates in vitro, albeit at usually feeble rates. Reactions with L-Thr, L-Arg, L-Lys and L-Asn were consistently very slow-a bias explained in part by the structural differences between these amino acids and the preferred substrates of the transaminases. On the other hand, L-His and L-Trp reacted more efficiently, particularly with GTK (glutamine transaminase K; also known as KYAT1). This points towards a role of GTK in the salvage of L-Trp (in cooperation with ω-amidase and possibly with the cytosolic malate dehydrogenase, MDH1, which efficiently reduced the product of L-Trp transamination). Finally, the transaminases were extremely ineffective at utilizing sugars and sugar derivatives, with the exception of the glycolytic intermediate dihydroxyacetone phosphate, which was slowly but appreciably transaminated by some of the enzymes to yield serinol phosphate. Evidence for the formation of this compound in a human cell line was also obtained. We discuss the biological and evolutionary implications of our results.

Ancestral sequence reconstruction of a cytochrome P450 family involved in chemical defence reveals the functional evolution of a promiscuous, xenobiotic-metabolizing enzyme in vertebrates

The cytochrome P450 family 1 enzymes (CYP1s) are a diverse family of hemoprotein monooxygenases, which metabolize many xenobiotics including numerous environmental carcinogens. However, their historical function and evolution remain largely unstudied. Here we investigate CYP1 evolution via the reconstruction and characterization of the vertebrate CYP1 ancestors. Younger ancestors and extant forms generally demonstrated higher activity toward typical CYP1 xenobiotic and steroid substrates than older ancestors, suggesting significant diversification away from the original CYP1 function. Caffeine metabolism appears to be a recently evolved trait of the CYP1A subfamily, observed in the mammalian CYP1A lineage, and may parallel the recent evolution of caffeine synthesis in multiple separate plant species. Likewise, the aryl hydrocarbon receptor agonist, 6-formylindolo[3,2-b]carbazole (FICZ) was metabolized to a greater extent by certain younger ancestors and extant forms, suggesting that activity toward FICZ increased in specific CYP1 evolutionary branches, a process that may have occurred in parallel to the exploitation of land where UV-exposure was higher than in aquatic environments. As observed with previous reconstructions of P450 enzymes, thermostability correlated with evolutionary age; the oldest ancestor was up to 35 °C more thermostable than the extant forms, with a ¹⁰T₅₀ (temperature at which 50% of the hemoprotein remains intact after 10 min) of 71 °C. This robustness may have facilitated evolutionary diversification of the CYP1s by buffering the destabilizing effects of mutations that conferred novel functions, a phenomenon which may also be useful in exploiting the catalytic versatility of these ancestral enzymes for commercial application as biocatalysts.

Multiple genetic paths including massive gene amplification allow Mycobacterium tuberculosis to overcome loss of ESX-3 secretion system substrates

The Mycobacterium tuberculosis (Mtb) ESX-3 type VII secretion system plays a critical role in iron acquisition. Infection of mice with highly attenuated Mtb deletion mutants lacking esxG or esxH, genes encoding key ESX-3 substrates, unexpectedly yielded suppressor mutants with restored capacity to grow in vivo and in vitro in the absence of iron supplementation. Whole-genome sequencing identified two mechanisms of suppression, the disruption of a transcriptional repressor that regulates expression of an ESX-3 paralogous region encoding EsxR and EsxS, and a massive 38- to 60-fold gene amplification of this same region. These data are significant because they reveal a previously unrecognized iron acquisition regulon and inform mechanisms of Mtb chromosome evolution.

Mycobacterium tuberculosis (Mtb) possesses five type VII secretion systems (T7SS), virulence determinants that include the secretion apparatus and associated secretion substrates. Mtb strains deleted for the genes encoding substrates of the ESX-3 T7SS, esxG or esxH, require iron supplementation for in vitro growth and are highly attenuated in vivo. In a subset of infected mice, suppressor mutants of esxG or esxH deletions were isolated, which enabled growth to high titers or restored virulence. Suppression was conferred by mechanisms that cause overexpression of an ESX-3 paralogous region that lacks genes for the secretion apparatus but encodes EsxR and EsxS, apparent ESX-3 orphan substrates that functionally compensate for the lack of EsxG or EsxH. The mechanisms include the disruption of a transcriptional repressor and a massive 38- to 60-fold gene amplification. These data identify an iron acquisition regulon, provide insight into T7SS, and reveal a mechanism of Mtb chromosome evolution involving “accordion-type” amplification.

Using Steady-State Kinetics to Quantitate Substrate Selectivity and Specificity: A Case Study with Two Human Transaminases

We examined the ability of two human cytosolic transaminases, aspartate aminotransferase (GOT1) and alanine aminotransferase (GPT), to transform their preferred substrates whilst discriminating against similar metabolites. This offers an opportunity to survey our current understanding of enzyme selectivity and specificity in a biological context. Substrate selectivity can be quantitated based on the ratio of the k_cat/K_M values for two alternative substrates (the 'discrimination index'). After assessing the advantages, implications and limits of this index, we analyzed the reactions of GOT1 and GPT with alternative substrates that are metabolically available and show limited structural differences with respect to the preferred substrates. The transaminases' observed selectivities were remarkably high. In particular, GOT1 reacted ~10⁶-fold less efficiently when the side-chain carboxylate of the 'physiological' substrates (aspartate and glutamate) was replaced by an amido group (asparagine and glutamine). This represents a current empirical limit of discrimination associated with this chemical difference. The structural basis of GOT1 selectivity was addressed through substrate docking simulations, which highlighted the importance of electrostatic interactions and proper substrate positioning in the active site. We briefly discuss the biological implications of these results and the possibility of using k_cat/K_M values to derive a global measure of enzyme specificity.

Highly thermostable and promiscuous β-1,3-xylanasen designed by optimized ancestral sequence reconstruction

The ancestor of β-1,3-xylanases (AncXyl09) were reconstructed by the optimized ancestral sequences reconstruction strategy to solve the poor catalytic performances of existing β-1,3-xylanases. The results showed that the half-life at 50 °C was 65.08 h, indicating good thermostability. The large number of hydrogen bonds and the disulfide bonds were the major attributes related with the thermal stability of Anxyl09. Interestingly, AncXyl09 could hydrolyze lichen besides the original substrate of β-1, 3-xylan, which is the first reported β-1,3-xylanase with substrate promiscuity. Moreover, the hydrolytic products are mainly disaccharides, the content of β-1,3-xylobiose and lichoridiose more than 70% as determined by high performance liquid chromatography (HPLC), which could significantly facilitate the separation and purification of oligosaccharides. The successful design of AncXyl09 was the representative of the semi-rationally engineered β-1, 3-xylanase, which will shield a new light on the β-1,3-xylanase engineering, active oligosaccharide preparation and marine algae resource utilization.

Were Ancestral Proteins Less Specific?

Some have hypothesized that ancestral proteins were, on average, less specific than their descendants. If true, this would provide a universal axis along which to organize protein evolution and suggests that reconstructed ancestral proteins may be uniquely powerful tools for protein engineering. Ancestral sequence reconstruction studies are one line of evidence used to support this hypothesis. Previously, we performed such a study, investigating the evolution of peptide-binding specificity for the paralogs S100A5 and S100A6. The modern proteins appeared more specific than their last common ancestor (ancA5/A6), as each paralog bound a subset of the peptides bound by ancA5/A6. In this study, we revisit this transition, using quantitative phage display to measure the interactions of 30,533 random peptides with human S100A5, S100A6, and ancA5/A6. This unbiased screen reveals a different picture. While S100A5 and S100A6 do indeed bind to a subset of the peptides recognized by ancA5/A6, they also acquired new peptide partners outside of the set recognized by ancA5/A6. Our previous work showed that ancA5/A6 had lower specificity than its descendants when measured against biological targets; our new work shows that ancA5/A6 has similar specificity to the modern proteins when measured against a random set of peptide targets. This demonstrates that altered biological specificity does not necessarily indicate altered intrinsic specificity, and sounds a cautionary note for using ancestral reconstruction studies with biological targets as a means to infer global evolutionary trends in specificity.

Comparative genomics of eight Lactobacillus buchneri strains isolated from food spoilage

Abstract

Background

Lactobacillus buchneri is a lactic acid bacterium frequently associated with food bioprocessing and fermentation and has been found to be either beneficial or detrimental to industrial food processes depending on the application. The ability to metabolize lactic acid into acetic acid and 1,2-propandiol makes L. buchneri invaluable to the ensiling process, however, this metabolic activity leads to spoilage in other applications, and is especially damaging to the cucumber fermentation industry. This study aims to augment our genomic understanding of L. buchneri in order to make better use of the species in a wide range of applicable industrial settings.

Results

Whole-genome sequencing (WGS) was performed on seven phenotypically diverse strains isolated from spoiled, fermented cucumber and the ATCC type strain for L. buchneri, ATCC 4005. Here, we present our findings from the comparison of eight newly-sequenced and assembled genomes against two publicly available closed reference genomes, L. buchneri CD034 and NRRL B-30929. Overall, we see ~ 50% of all coding sequences are conserved across these ten strains. When these coding sequences are clustered by functional description, the strains appear to be enriched in mobile genetic elements, namely transposons. All isolates harbor at least one CRISPR-Cas system, and many contain putative prophage regions, some of which are targeted by the host’s own DNA-encoded spacer sequences.

Conclusions

Our findings provide new insights into the genomics of L. buchneri through whole genome sequencing and subsequent characterization of genomic features, building a platform for future studies and identifying elements for potential strain manipulation or engineering.

An efficient expression tag library based on self-assembling amphipathic peptides

Background

Self-assembling amphipathic peptides (SAPs) may improve protein production or induce the formation of inclusion bodies by fusing them to the N-terminus of proteins. However, they do not function uniformly well with all target enzymes and systematic research on how the composition of SAPs influence the production of fusion protein is still limited.

Results

To improve the efficiency of SAPs, we studied factors that might be involved in SAP-mediated protein production using S1 (AEAEAKAK)₂ as the original SAP and green fluorescent protein (GFP) as the reporter. The results indicate that hydrophobicity and net charges of SAPs play a key role in protein expression. As hydrophobicity regulation tend to cause the formation of insoluble inclusion bodies of protein, an expression tag library composed of SAPs, which varied in net charge (from + 1 to + 20), was constructed based on the random amplification of S1nv1 (ANANARAR)₁₀. The efficiency of the library was validated by polygalacturonate lyase (PGL), lipoxygenase (LOX), l-asparaginase (ASN) and transglutaminase (MTG). To accelerate preliminary screening, each enzyme was fused at the C-terminus with GFP. Among the four enzyme fusions, the SAPs with + 2 – + 6 net charges were optimal for protein expression. Finally, application of the library improved the expression of PGL, LOX, ASN, and MTG by 8.3, 3.5, 2.64, and 3.68-fold relative to that of the corresponding wild-type enzyme, respectively.

Conclusions

This is the first report to study key factors of SAPs as an expression tag to enhance recombinant enzyme production. The SAP library could be used as a novel plug-and-play protein-engineering method to screen for enzymes or proteins with enhanced production.

Electronic supplementary material

The online version of this article (10.1186/s12934-019-1142-9) contains supplementary material, which is available to authorized users.

Selection for novel metabolic capabilities in Salmonella enterica

Bacteria are known to display extensive metabolic diversity and many studies have shown that they can use an extensive repertoire of small molecules as carbon‐ and energy sources. However, it is less clear to what extent a bacterium can expand its existing metabolic capabilities by acquiring mutations that, for example, rewire its metabolic pathways. To investigate this capability and potential for evolution of novel phenotypes, we sampled large populations of mutagenized Salmonella enterica to select very rare mutants that can grow on minimal media containing 124 low molecular weight compounds as sole carbon sources. We found mutants growing on 18 of these novel carbon sources, and identified the causal mutations that allowed growth for four of them. Mutations that relieve physiological constraints or increase expression of existing pathways were found to be important contributors to the novel phenotypes. For the remaining 14 novel phenotypes, whole genome sequencing of independent mutants and genetic analysis suggested that these novel metabolic phenotypes result from a combination of multiple mutations. This work, by virtue of identifying the genetic and mechanistic basis for new metabolic capabilities, sheds light on the properties of adaptive landscapes underlying the evolution of novel phenotypes.

Local applications but global implications: Can pesticides drive microorganisms to develop antimicrobial resistance?

Pesticides are an important agricultural input, and the introduction of new active ingredients with increased efficiencies drives their higher production and consumption worldwide. Inappropriate application and storage of these chemicals often contaminate plant tissues, air, water, or soil environments. The presence of pesticides can lead to developing tolerance, resistance or persistence and even the capabilities to degrade them by the microbiomes of theses environments. The pesticide-degrading microorganisms gain and employ several mechanisms for attraction (chemotaxis), membrane transport systems, efflux pumps, enzymes and genetical make-up with plasmid and chromosome encoded catabolic genes for degradation. Even the evolution and the mechanisms of inheritance for pesticide-degradation as a functional trait in several microorganisms are beginning to be understood. Because of the commonalities in the microbial responses of sensing and uptake, and adaptation due to the selection pressures of pesticides and antimicrobial substances including antibiotics, the pesticide-degraders have higher chances of possessing antimicrobial resistance as a surplus functional trait. This review critically examines the probabilities of pesticide contamination of soil and foliage, the knowledge gaps in the regulation and storage of pesticide chemicals, and the human implications of pesticide-degrading microorganisms with antimicrobial resistance in the global strategy of 'One Health'.

CoMetGeNe: mining conserved neighborhood patterns in metabolic and genomic contexts

Background

In systems biology, there is an acute need for integrative approaches in heterogeneous network mining in order to exploit the continuous flux of genomic data. Simultaneous analysis of the metabolic pathways and genomic context of a given species leads to the identification of patterns consisting in reaction chains catalyzed by products of neighboring genes. Similar such patterns across several species can reveal their mode of conservation throughout the tree of life.

Results

We present CoMetGeNe (COnserved METabolic and GEnomic NEighborhoods), a novel method that identifies metabolic and genomic patterns consisting in maximal trails of reactions being catalyzed by products of neighboring genes. Patterns determined by CoMetGeNe in one species are subsequently employed in order to reflect their degree of conservation across multiple prokaryotic species. These interspecies comparisons help to improve genome annotation and can reveal putative alternative metabolic routes as well as unexpected gene ordering occurrences.

Conclusions

CoMetGeNe is an exploratory tool at both the genomic and the metabolic levels, leading to insights into the conservation of functionally related clusters of neighboring enzyme-coding genes. The open-source CoMetGeNe pipeline is freely available at https://cometgene.lri.fr.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2542-2) contains supplementary material, which is available to authorized users.

Balancing Specificity and Promiscuity in Enzyme Evolution: Multidimensional Activity Transitions in the Alkaline Phosphatase Superfamily

Highly proficient, promiscuous enzymes can be springboards for functional evolution, able to avoid loss of function during adaptation by their capacity to promote multiple reactions. We employ a systematic comparative study of structure, sequence, and substrate specificity to track the evolution of specificity and reactivity between promiscuous members of clades of the alkaline phosphatase (AP) superfamily. Construction of a phylogenetic tree of protein sequences maps out the likely transition zone between arylsulfatases (ASs) and phosphonate monoester hydrolases (PMHs). Kinetic analysis shows that all enzymes characterized have four chemically distinct phospho- and sulfoesterase activities, with rate accelerations ranging from 10¹¹- to 10¹⁷-fold for their primary and 10⁹- to 10¹²-fold for their promiscuous reactions, suggesting that catalytic promiscuity is widespread in the AP-superfamily. This functional characterization and crystallography reveal a novel class of ASs that is so similar in sequence to known PMHs that it had not been recognized as having diverged in function. Based on analysis of snapshots of catalytic promiscuity "in transition", we develop possible models that would allow functional evolution and determine scenarios for trade-off between multiple activities. For the new ASs, we observe largely invariant substrate specificity that would facilitate the transition from ASs to PMHs via trade-off-free molecular exaptation, that is, evolution without initial loss of primary activity and specificity toward the original substrate. This ability to bypass low activity generalists provides a molecular solution to avoid adaptive conflict.

Protein moonlighting elucidates the essential human pathway catalyzing lipoic acid assembly on its cognate enzymes

The lack of attachment of lipoic acid to its cognate enzyme proteins results in devastating human metabolic disorders. These mitochondrial disorders are evident soon after birth and generally result in early death. The mutations causing specific defects in lipoyl assembly map in three genes, LIAS, LIPT1, and LIPT2 Although physiological roles have been proposed for the encoded proteins, only the LIPT1 protein had been studied at the enzyme level. LIPT1 was reported to catalyze only the second partial reaction of the classical lipoate ligase mechanism. We report that the physiologically relevant LIPT1 enzyme activity is transfer of lipoyl moieties from the H protein of the glycine cleavage system to the E2 subunits of the 2-oxoacid dehydrogenases required for respiration (e.g., pyruvate dehydrogenase) and amino acid degradation. We also report that LIPT2 encodes an octanoyl transferase that initiates lipoyl group assembly. The human pathway is now biochemically defined.

Systematic analysis reveals long noncoding RNAs regulating neighboring transcription factors in human cancers

Long noncoding RNAs (lncRNAs) are proposed to play essential roles in regulating gene transcription. Moreover, a subset has been implicated in modulating the expression of the nearby loci. Here we systematically evaluated the relationship between lncRNAs and their neighboring genes based on transcriptome expression profiles from 4900 samples across 12 cancer types. Our findings reveal that lncRNAs, especially those of high syntenic conservation across species, are spatially correlated with transcription factors across the genome. Combining the methods of conservation, co-expression, and causal inference test, we identified a list of 28 lncRNA/TF regulatory pairs across 12 TCGA cancer types, and 19 of which were further confirmed in additional cancer cell lines. Several of these pairs, including PTV1/MYC and GATA6-AS1/GATA6, show prior evidence of regulatory relationships. Other candidates such as LINC00261/FOXA2 and PITRM1-AS1/KLF6 were novel. Our study highlights the significant roles of lncRNAs in tumorigenesis and provides a comprehensive overview of lncRNA regulation on its neighboring TF genes in human cancers.

Evolution of enzymes with new specificity by high-throughput screening using DmpR-based genetic circuits and multiple flow cytometry rounds

Genetic circuit-based biosensors are useful in detecting target metabolites or in vivo enzymes using transcription factors (Tx) as a molecular switch to express reporter signals, such as cellular fluorescence and antibiotic resistance. Herein, a phenol-detecting Tx (DmpR) was employed as a critical tool for enzyme engineering, specifically for the rapid analysis of numerous mutants with multiple mutations at the active site of tryptophan-indole lyase (TIL, EC 4.1.99.1). Cellular fluorescence was monitored cell-by-cell using flow cytometry to detect the creation of phenolic compounds by a new tyrosine-phenol-lyase (TPL, EC 4.1.99.2). In the TIL scaffold, target amino acids near the indole ring (Asp¹³⁷, Phe³⁰⁴, Val³⁹⁴, Ile³⁹⁶ and His⁴⁶³) were mutated randomly to construct a large diversity of specificity variations. Collection of candidate positives by cell sorting using flow cytometry and subsequent shuffling of beneficial mutations identified a critical hit with four mutations (D137P, F304D, V394L, and I396R) in the TIL sequence. The variant displayed one-thirteenth the level of TPL activity, compared with native TPLs, and completely lost the original TIL activity. The findings demonstrate that hypersensitive, Tx-based biosensors could be useful critically to generate new activity from a related template, which would alleviate the current burden to high-throughput screening.

Conservation of Specificity in Two Low-Specificity Proteins

Many regulatory proteins bind peptide regions of target proteins and modulate their activity. Such regulatory proteins can often interact with highly diverse target peptides. In many instances, it is not known if the peptide-binding interface discriminates targets in a biological context, or whether biological specificity is achieved exclusively through external factors such as subcellular localization. We used an evolutionary biochemical approach to distinguish these possibilities for two such low-specificity proteins: S100A5 and S100A6. We used isothermal titration calorimetry to study the binding of peptides with diverse sequence and biochemistry to human S100A5 and S100A6. These proteins bound distinct, but overlapping, sets of peptide targets. We then studied the peptide binding properties of orthologs sampled from across five amniote species. Binding specificity was conserved along all lineages, for the last 320 million years, despite the low specificity of each protein. We used ancestral sequence reconstruction to determine the binding specificity of the last common ancestor of the paralogs. The ancestor bound the entire set of peptides bound by modern S100A5 and S100A6 proteins, suggesting that paralog specificity evolved via subfunctionalization. To rule out the possibility that specificity is conserved because it is difficult to modify, we identified a single historical mutation that, when reverted in human S100A5, gave it the ability to bind an S100A6-specific peptide. These results reveal strong evolutionary constraints on peptide binding specificity. Despite being able to bind a large number of targets, the specificity of S100 peptide interfaces is likely important for the biology of these proteins.

How the chemical features of molecules may have addressed the settlement of metabolic steps

INTRODUCTION

While the evolutionary adaptation of enzymes to their own substrates is a well assessed and rationalized field, how molecules have been originally selected in order to initiate and assemble convenient metabolic pathways is a fascinating, but still debated argument.

OBJECTIVES

Aim of the present study is to give a rationale for the preferential selection of specific molecules to generate metabolic pathways.

METHODS

The comparison of structural features of molecules, through an inductive methodological approach, offer a reading key to cautiously propose a determining factor for their metabolic recruitment.

RESULTS

Starting with some commonplaces occurring in the structural representation of relevant carbohydrates, such as glucose, fructose and ribose, arguments are presented in associating stable structural determinants of these molecules and their peculiar occurrence in metabolic pathways.

CONCLUSIONS

Among other possible factors, the reliability of the structural asset of a molecule may be relevant or its selection among structurally and, a priori, functionally similar molecules.

A two-step evolutionary process establishes a non-native vitamin B6 pathway in Bacillus subtilis

Pyridoxal 5'-phosphate (PLP), the most important form of vitamin B6 serves as a cofactor for many proteins. Two alternative pathways for de novo PLP biosynthesis are known: the short deoxy-xylulose-5-phosphate (DXP)-independent pathway, which is present in the Gram-positive model bacterium Bacillus subtilis and the longer DXP-dependent pathway, which has been intensively studied in the Gram-negative model bacterium Escherichia coli. Previous studies revealed that bacteria contain many promiscuous enzymes causing a so-called 'underground metabolism', which can be important for the evolution of novel pathways. Here, we evaluated the potential of B. subtilis to use a truncated non-native DXP-dependent PLP pathway from E. coli for PLP synthesis. Adaptive laboratory evolution experiments revealed that two non-native enzymes catalysing the last steps of the DXP-dependent PLP pathway and two genomic alterations are sufficient to allow growth of vitamin B6 auxotrophic bacteria as rapid as the wild type. Thus, the existence of an underground metabolism in B. subtilis facilitates the generation of a pathway for synthesis of PLP using parts of a non-native vitamin B6 pathway. The introduction of non-native enzymes into a metabolic network and rewiring of native metabolism could be helpful to generate pathways that might be optimized for producing valuable substances.

Evolution of biosynthetic diversity

Since the emergence of the last common ancestor from which all extant life evolved, the metabolite repertoire of cells has increased and diversified. Not only has the metabolite cosmos expanded, but the ways in which the same metabolites are made have diversified. Enzymes catalyzing the same reaction have evolved independently from different protein folds; the same protein fold can produce enzymes recognizing different substrates, and enzymes performing different chemistries. Genes encoding useful enzymes can be transferred between organisms and even between the major domains of life. Organisms that live in metabolite-rich environments sometimes lose the pathways that produce those same metabolites. Fusion of different protein domains results in enzymes with novel properties. This review will consider the major evolutionary mechanisms that generate biosynthetic diversity: gene duplication (and gene loss), horizontal and endosymbiotic gene transfer, and gene fusion. It will also discuss mechanisms that lead to convergence as well as divergence. To illustrate these mechanisms, one of the original metabolisms present in the last universal common ancestor will be employed: polyamine metabolism, which is essential for the growth and cell proliferation of archaea and eukaryotes, and many bacteria.

One Enzyme, Three Metabolites: Shewanella algae Controls Siderophore Production via the Cellular Substrate Pool

Shewanella algae B516 produces avaroferrin, an asymmetric hydroxamate siderophore, which has been shown to inhibit swarming motility of Vibrio alginolyticus. We aimed to elucidate the biosynthesis of this siderophore and to investigate how S. algae coordinates the production of avaroferrin and its two symmetric counterparts. We reconstituted the reaction in vitro with the main enzyme AvbD and the putative biosynthetic precursors, and demonstrate that multispecificity of this enzyme results in the production of all three cyclic hydroxamate siderophores that were previously isolated as natural products from S. algae. Surprisingly, purified AvbD exhibited a clear preference for the larger cadaverine-derived substrate. In live cells, however, siderophore ratios are maximized toward avaroferrin production, and we demonstrate that these siderophore ratios are the result of a regulation on substrate pool level, which may allow rapid evolutionary adaptation to environmental changes. Our results thereby give insights into a unique evolutionary strategy toward metabolite diversity.

Experimental Evolution of Escherichia coli Harboring an Ancient Translation Protein

The ability to design synthetic genes and engineer biological systems at the genome scale opens new means by which to characterize phenotypic states and the responses of biological systems to perturbations. One emerging method involves inserting artificial genes into bacterial genomes and examining how the genome and its new genes adapt to each other. Here we report the development and implementation of a modified approach to this method, in which phylogenetically inferred genes are inserted into a microbial genome, and laboratory evolution is then used to examine the adaptive potential of the resulting hybrid genome. Specifically, we engineered an approximately 700-million-year-old inferred ancestral variant of tufB, an essential gene encoding elongation factor Tu, and inserted it in a modern Escherichia coli genome in place of the native tufB gene. While the ancient homolog was not lethal to the cell, it did cause a twofold decrease in organismal fitness, mainly due to reduced protein dosage. We subsequently evolved replicate hybrid bacterial populations for 2000 generations in the laboratory and examined the adaptive response via fitness assays, whole genome sequencing, proteomics, and biochemical assays. Hybrid lineages exhibit a general adaptive strategy in which the fitness cost of the ancient gene was ameliorated in part by upregulation of protein production. Our results suggest that an ancient-modern recombinant method may pave the way for the synthesis of organisms that exhibit ancient phenotypes, and that laboratory evolution of these organisms may prove useful in elucidating insights into historical adaptive processes.

Vitamin B6 metabolism in microbes and approaches for fermentative production

Vitamin B6 is a designation for the six vitamers pyridoxal, pyridoxine, pyridoxamine, pyridoxal 5'-phosphate (PLP), pyridoxine 5'-phosphate, and pyridoxamine. PLP, being the most important B6 vitamer, serves as a cofactor for many proteins and enzymes. In contrast to other organisms, animals and humans have to ingest vitamin B6 with their food. Several disorders are associated with vitamin B6 deficiency. Moreover, pharmaceuticals interfere with metabolism of the cofactor, which also results in vitamin B6 deficiency. Therefore, vitamin B6 is a valuable compound for the pharmaceutical and the food industry. Although vitamin B6 is currently chemically synthesized, there is considerable interest on the industrial side to shift from chemical processes to sustainable fermentation technologies. Here, we review recent findings regarding biosynthesis and homeostasis of vitamin B6 and describe the approaches that have been made in the past to develop microbial production processes. Moreover, we will describe novel routes for vitamin B6 biosynthesis and discuss their potential for engineering bacteria that overproduce the commercially valuable substance. We also highlight bottlenecks of the vitamin B6 biosynthetic pathways and propose strategies to circumvent these limitations.

Computational Identification of the Paralogs and Orthologs of Human Cytochrome P450 Superfamily and the Implication in Drug Discovery

The human cytochrome P450 (CYP) superfamily consisting of 57 functional genes is the most important group of Phase I drug metabolizing enzymes that oxidize a large number of xenobiotics and endogenous compounds, including therapeutic drugs and environmental toxicants. The CYP superfamily has been shown to expand itself through gene duplication, and some of them become pseudogenes due to gene mutations. Orthologs and paralogs are homologous genes resulting from speciation or duplication, respectively. To explore the evolutionary and functional relationships of human CYPs, we conducted this bioinformatic study to identify their corresponding paralogs, homologs, and orthologs. The functional implications and implications in drug discovery and evolutionary biology were then discussed. GeneCards and Ensembl were used to identify the paralogs of human CYPs. We have used a panel of online databases to identify the orthologs of human CYP genes: NCBI, Ensembl Compara, GeneCards, OMA ("Orthologous MAtrix") Browser, PATHER, TreeFam, EggNOG, and Roundup. The results show that each human CYP has various numbers of paralogs and orthologs using GeneCards and Ensembl. For example, the paralogs of CYP2A6 include CYP2A7, 2A13, 2B6, 2C8, 2C9, 2C18, 2C19, 2D6, 2E1, 2F1, 2J2, 2R1, 2S1, 2U1, and 2W1; CYP11A1 has 6 paralogs including CYP11B1, 11B2, 24A1, 27A1, 27B1, and 27C1; CYP51A1 has only three paralogs: CYP26A1, 26B1, and 26C1; while CYP20A1 has no paralog. The majority of human CYPs are well conserved from plants, amphibians, fishes, or mammals to humans due to their important functions in physiology and xenobiotic disposition. The data from different approaches are also cross-validated and validated when experimental data are available. These findings facilitate our understanding of the evolutionary relationships and functional implications of the human CYP superfamily in drug discovery.

Cupincin: A Unique Protease Purified from Rice (Oryza sativa L.) Bran Is a New Member of the Cupin Superfamily

Cupin superfamily is one of the most diverse super families. This study reports the purification and characterization of a novel cupin domain containing protease from rice bran for the first time. Hypothetical protein OsI_13867 was identified and named as cupincin. Cupincin was purified to 4.4 folds with a recovery of 4.9%. Cupincin had an optimum pH and temperature of pH 4.0 and 60°C respectively. Cupincin was found to be a homotrimer, consisting of three distinct subunits with apparent molecular masses of 33.45 kDa, 22.35 kDa and 16.67 kDa as determined by MALDI-TOF, whereas it eluted as a single unit with an apparent molecular mass of 135.33 ± 3.52 kDa in analytical gel filtration and migrated as a single band in native page, suggesting its homogeneity. Sequence identity of cupincin was deduced by determining the amino-terminal sequence of the polypeptide chains and by and de novo sequencing. For understanding the hydrolysing mechanism of cupincin, its three-dimensional model was developed. Structural analysis indicated that cupincin contains His313, His326 and Glu318 with zinc ion as the putative active site residues, inhibition of enzyme activity by 1,10-phenanthroline and atomic absorption spectroscopy confirmed the presence of zinc ion. The cleavage specificity of cupincin towards oxidized B-chain of insulin was highly specific; cleaving at the Leu₁₅-Tyr₁₆ position, the specificity was also determined using neurotensin as a substrate, where it cleaved only at the Glu₁-Tyr₂ position. Limited proteolysis of the protease suggests a specific function for cupincin. These results demonstrated cupincin as a completely new protease.

Ancestral Protein Reconstruction Yields Insights into Adaptive Evolution of Binding Specificity in Solute-Binding Proteins

The promiscuous functions of proteins are an important reservoir of functional novelty in protein evolution, but the molecular basis for binding promiscuity remains elusive. We used ancestral protein reconstruction to experimentally characterize evolutionary intermediates in the functional expansion of the polar amino acid-binding protein family, which has evolved to bind a variety of amino acids with high affinity and specificity. High-resolution crystal structures of an ancestral arginine-binding protein in complex with l-arginine and l-glutamine show that the promiscuous binding of l-glutamine is enabled by multi-scale conformational plasticity, water-mediated interactions, and selection of an alternative conformational substate productive for l-glutamine binding. Evolution of specialized glutamine-binding proteins from this ancestral protein was achieved by displacement of water molecules from the protein-ligand interface, reducing the entropic penalty associated with the promiscuous interaction. These results provide a structural and thermodynamic basis for the co-option of a promiscuous interaction in the evolution of binding specificity.

Targeted metagenomics unveils the molecular basis for adaptive evolution of enzymes to their environment

Microorganisms have a wonderful ability to adapt rapidly to new or altered environmental conditions. Enzymes are the basis of metabolism in all living organisms and, therefore, enzyme adaptation plays a crucial role in the adaptation of microorganisms. Comparisons of homology and parallel beneficial mutations in an enzyme family provide valuable hints of how an enzyme adapted to an ecological system; consequently, a series of enzyme collections is required to investigate enzyme evolution. Targeted metagenomics is a promising tool for the construction of enzyme pools and for studying the adaptive evolution of enzymes. This perspective article presents a summary of targeted metagenomic approaches useful for this purpose.

The Impact of Non-Enzymatic Reactions and Enzyme Promiscuity on Cellular Metabolism during (Oxidative) Stress Conditions

Cellular metabolism assembles in a structurally highly conserved, but functionally dynamic system, known as the metabolic network. This network involves highly active, enzyme-catalyzed metabolic pathways that provide the building blocks for cell growth. In parallel, however, chemical reactivity of metabolites and unspecific enzyme function give rise to a number of side products that are not part of canonical metabolic pathways. It is increasingly acknowledged that these molecules are important for the evolution of metabolism, affect metabolic efficiency, and that they play a potential role in human disease—age-related disorders and cancer in particular. In this review we discuss the impact of oxidative and other cellular stressors on the formation of metabolic side products, which originate as a consequence of: (i) chemical reactivity or modification of regular metabolites; (ii) through modifications in substrate specificity of damaged enzymes; and (iii) through altered metabolic flux that protects cells in stress conditions. In particular, oxidative and heat stress conditions are causative of metabolite and enzymatic damage and thus promote the non-canonical metabolic activity of the cells through an increased repertoire of side products. On the basis of selected examples, we discuss the consequences of non-canonical metabolic reactivity on evolution, function and repair of the metabolic network.

The Classification and Evolution of Enzyme Function

Enzymes are the proteins responsible for the catalysis of life. Enzymes sharing a common ancestor as defined by sequence and structure similarity are grouped into families and superfamilies. The molecular function of enzymes is defined as their ability to catalyze biochemical reactions; it is manually classified by the Enzyme Commission and robust approaches to quantitatively compare catalytic reactions are just beginning to appear. Here, we present an overview of studies at the interface of the evolution and function of enzymes.

From microbiology to cancer biology: the Rid protein family prevents cellular damage caused by endogenously generated reactive nitrogen species

The Rid family of proteins is highly conserved and broadly distributed throughout the domains of life. Genetic and biochemical studies, primarily in Salmonella enterica, have defined a role for RidA in responding to endogenously generated reactive metabolites. The data show that 2-aminoacrylate (2AA), a reactive enamine intermediate generated by some pyridoxal 5'-phosphate-dependent enzymes, accumulates in the absence of RidA. The accumulation of 2AA leads to covalent modification and inactivation of several enzymes involved in essential metabolic processes. This review describes the 2AA hydrolyzing activity of RidA and the effect of this biochemical activity on the metabolic network, which impacts organism fitness. The reported activity of RidA and the consequences encountered in vivo when RidA is absent have challenged fundamental assumptions in enzymology, biochemistry and cell metabolism regarding the fate of transiently generated reactive enamine intermediates. The current understanding of RidA in Salmonella and the broad distribution of Rid family proteins provide exciting opportunities for future studies to define metabolic roles of Rid family members from microbes to man.

Repurposing TRASH: emergence of the enzyme organomercurial lyase from a non-catalytic zinc finger scaffold

The mercury resistance pathway enzyme organomercurial lyase (MerB) catalyzes the conversion of organomercurials to ionic mercury (Hg(2+)). Here, we provide evidence for the emergence of this enzyme from a TRASH-like, non-enzymatic, treble-clef zinc finger ancestor by domain duplication and fusion. Surprisingly, the structure-stabilizing metal-binding core of the treble-clef appears to have been repurposed in evolution to serve a catalytic role. Novel enzymatic functions are believed to have evolved from ancestral generalist catalytic scaffolds or from already specialized enzymes with catalytic promiscuity. The emergence of MerB from a zinc finger ancestor serves as a rare example of how a novel enzyme may emerge from a non-catalytic scaffold with a related binding function.

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.

The evolution of enzyme function in the isomerases

The advent of computational approaches to measure functional similarity between enzymes adds a new dimension to existing evolutionary studies based on sequence and structure. This paper reviews research efforts aiming to understand the evolution of enzyme function in superfamilies, presenting a novel strategy to provide an overview of the evolution of enzymes belonging to an individual EC class, using the isomerases as an exemplar.

The future of drug discovery: enabling technologies for enhancing lead characterization and profiling therapeutic potential

Technology often serves as a handmaiden and catalyst of invention. The discovery of safe, effective medications depends critically upon experimental approaches capable of providing high-impact information on the biological effects of drug candidates early in the discovery pipeline. This information can enable reliable lead identification, pharmacological compound differentiation and successful translation of research output into clinically useful therapeutics. The shallow preclinical profiling of candidate compounds promulgates a minimalistic understanding of their biological effects and undermines the level of value creation necessary for finding quality leads worth moving forward within the development pipeline with efficiency and prognostic reliability sufficient to help remediate the current pharma-industry productivity drought. Three specific technologies discussed herein, in addition to experimental areas intimately associated with contemporary drug discovery, appear to hold particular promise for strengthening the preclinical valuation of drug candidates by deepening lead characterization. These are: i) hydrogen-deuterium exchange mass spectrometry for characterizing structural and ligand-interaction dynamics of disease-relevant proteins; ii) activity-based chemoproteomics for profiling the functional diversity of mammalian proteomes; and iii) nuclease-mediated precision gene editing for developing more translatable cellular and in vivo models of human diseases. When applied in an informed manner congruent with the clinical understanding of disease processes, technologies such as these that span levels of biological organization can serve as valuable enablers of drug discovery and potentially contribute to reducing the current, unacceptably high rates of compound clinical failure.

Enzyme informatics

Over the last 50 years, sequencing, structural biology and bioinformatics have completely revolutionised biomolecular science, with millions of sequences and tens of thousands of three dimensional structures becoming available. The bioinformatics of enzymes is well served by, mostly free, online databases. BRENDA describes the chemistry, substrate specificity, kinetics, preparation and biological sources of enzymes, while KEGG is valuable for understanding enzymes and metabolic pathways. EzCatDB, SFLD and MACiE are key repositories for data on the chemical mechanisms by which enzymes operate. At the current rate of genome sequencing and manual annotation, human curation will never finish the functional annotation of the ever-expanding list of known enzymes. Hence there is an increasing need for automated annotation, though it is not yet widespread for enzyme data. In contrast, functional ontologies such as the Gene Ontology already profit from automation. Despite our growing understanding of enzyme structure and dynamics, we are only beginning to be able to design novel enzymes. One can now begin to trace the functional evolution of enzymes using phylogenetics. The ability of enzymes to perform secondary functions, albeit relatively inefficiently, gives clues as to how enzyme function evolves. Substrate promiscuity in enzymes is one example of imperfect specificity in protein-ligand interactions. Similarly, most drugs bind to more than one protein target. This may sometimes result in helpful polypharmacology as a drug modulates plural targets, but also often leads to adverse side-effects. Many chemoinformatics approaches can be used to model the interactions between druglike molecules and proteins in silico. We can even use quantum chemical techniques like DFT and QM/MM to compute the structural and energetic course of enzyme catalysed chemical reaction mechanisms, including a full description of bond making and breaking.

Evolution of a new chlorophyll metabolic pathway driven by the dynamic changes in enzyme promiscuous activity

Organisms generate an enormous number of metabolites; however, the mechanisms by which a new metabolic pathway is acquired are unknown. To elucidate the importance of promiscuous enzyme activity for pathway evolution, the catalytic and substrate specificities of Chl biosynthetic enzymes were examined. In green plants, Chl a and Chl b are interconverted by the Chl cycle: Chl a is hydroxylated to 7-hydroxymethyl chlorophyll a followed by the conversion to Chl b, and both reactions are catalyzed by chlorophyllide a oxygenase. Chl b is reduced to 7-hydroxymethyl chlorophyll a by Chl b reductase and then converted to Chl a by 7-hydroxymethyl chlorophyll a reductase (HCAR). A phylogenetic analysis indicated that HCAR evolved from cyanobacterial 3,8-divinyl chlorophyllide reductase (DVR), which is responsible for the reduction of an 8-vinyl group in the Chl biosynthetic pathway. In addition to vinyl reductase activity, cyanobacterial DVR also has Chl b reductase and HCAR activities; consequently, three of the four reactions of the Chl cycle already existed in cyanobacteria, the progenitor of the chloroplast. During the evolution of cyanobacterial DVR to HCAR, the HCAR activity, a promiscuous reaction of cyanobacterial DVR, became the primary reaction. Moreover, the primary reaction (vinyl reductase activity) and some disadvantageous reactions were lost, but the neutral promiscuous reaction (NADH dehydrogenase) was retained in both DVR and HCAR. We also show that a portion of the Chl c biosynthetic pathway already existed in cyanobacteria. We discuss the importance of dynamic changes in promiscuous activity and of the latent pathways for metabolic evolution.

Highly diverse protein library based on the ubiquitous (β/α)₈ enzyme fold yields well-structured proteins through in vitro folding selection

Proper protein folding is a prerequisite for protein stability and enzymatic activity. Although directed evolution can be a powerful tool to investigate enzymatic function and to isolate novel activities, well-designed libraries of folded proteins are essential. In vitro selection methods are particularly capable of searching for enzymatic activities in libraries of trillions of protein variants, yet high-quality libraries of well-folded enzymes with such high diversity are lacking. We describe the construction and detailed characterization of a folding-enriched protein library based on the ubiquitous (β/α)₈ barrel fold, which is found in five of the six enzyme classes. We introduced seven randomized loops on the catalytic face of the monomeric, thermostable (β/α)₈ barrel of glycerophosphodiester phosphodiesterase (GDPD) from Thermotoga maritima. We employed in vitro folding selection based on protease digestion to enrich intermediate libraries containing three to four randomized loops for folded variants, and then combined them to assemble the final library (10¹⁴ DNA sequences). The resulting library was analyzed by using the in vitro protease assay and an in vivo GFP-folding assay; it contains ∼10¹² soluble monomeric protein variants. We isolated six library members and demonstrated that these proteins are soluble, monomeric and show (β/α)₈-barrel fold-like secondary and tertiary structure. The quality of the folding-enriched library improved up to 50-fold compared to a control library that was assembled without the folding selection. To the best of our knowledge, this work is the first example of combining the ultra-high throughput mRNA display method with selection for folding. The resulting (β/α)₈ barrel libraries provide a valuable starting point to study the unique catalytic capabilities of the (β/α)₈ fold, and to isolate novel enzymes.