Homology among (betaalpha)(8) barrels: implications for the evolution of metabolic pathways

Lysophosphatidylglycerol (LPG) phospholipase D maintains membrane homeostasis in Staphylococcus aureus by converting LPG to lysophosphatidic acid

Lysophospholipids are deacylated derivatives of their bilayer forming phospholipid counterparts that are present at low concentrations in cells. Phosphatidylglycerol (PG) is the principal membrane phospholipid in Staphylococcus aureus and lysophosphatidylglycerol (LPG) is detected in low abundance. Here, we used a mass spectrometry screen to identify locus SAUSA300_1020 as the gene responsible for maintaining low concentrations of 1-acyl-LPG in S. aureus. The SAUSA300_1020 gene encodes a protein with a predicted amino terminal transmembrane α-helix attached to a globular glycerophosphodiester phosphodiesterase (GDPD) domain. We determined that the purified protein lacking the hydrophobic helix (LpgDΔN) possesses cation-dependent lysophosphatidylglycerol phospholipase D activity that generates both lysophosphatidic acid (LPA) and cyclic-LPA products and hydrolyzes cyclic-LPA to LPA. Mn2+ was the highest affinity cation and stabilized LpgDΔN to thermal denaturation. LpgDΔN was not specific for the phospholipid headgroup and degraded 1-acyl-LPG, but not 2-acyl-LPG. Furthermore, a 2.1 Å crystal structure shows that LpgDΔN adopts the GDPD variation of the TIM barrel architecture except for the length and positioning of helix α6 and sheet β7. These alterations create a hydrophobic diffusion path for LPG to access the active site. The LpgD active site has the canonical GDPD metal binding and catalytic residues, and our biochemical characterization of site-directed mutants support a two-step mechanism involving a cyclic-LPA intermediate. Thus, the physiological function of LpgD in S. aureus is to convert LPG to LPA, which is re-cycled into the PG biosynthetic pathway at the LPA acyltransferase step to maintain membrane PG molecular species homeostasis.

Uniform binding and negative catalysis at the origin of enzymes

Enzymes are well known for their catalytic abilities, some even reaching "catalytic perfection" in the sense that the reaction they catalyze has reached the physical bound of the diffusion rate. However, our growing understanding of enzyme superfamilies has revealed that only some share a catalytic chemistry while others share a substrate-handle binding motif, for example, for a particular phosphate group. This suggests that some families emerged through a "substrate-handle-binding-first" mechanism ("binding-first" for brevity) instead of "chemistry-first" and we are, therefore, left to wonder what the role of non-catalytic binders might have been during enzyme evolution. In the last of their eight seminal, back-to-back articles from 1976, John Albery and Jeremy Knowles addressed the question of enzyme evolution by arguing that the simplest mode of enzyme evolution is what they defined as "uniform binding" (parallel stabilization of all enzyme-bound states to the same degree). Indeed, we show that a uniform-binding proto-catalyst can accelerate a reaction, but only when catalysis is already present, that is, when the transition state is already stabilized to some degree. Thus, we sought an alternative explanation for the cases where substrate-handle-binding preceded any involvement of a catalyst. We find that evolutionary starting points that exhibit negative catalysis can redirect the reaction's course to a preferred product without need for rate acceleration or product release; that is, if they do not stabilize, or even destabilize, the transition state corresponding to an undesired product. Such a mechanism might explain the emergence of "binding-first" enzyme families like the aldolase superfamily.

Enabling Role of Ligand-Driven Conformational Changes in Enzyme Evolution

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

Effect of Non-Lethal Selection on Spontaneous Revertants of Frameshift Mutations: The Escherichia coli hisF Case

Microorganisms possess the potential to adapt to fluctuations in environmental parameters, and their evolution is driven by the continuous generation of mutations. The reversion of auxotrophic mutations has been widely studied; however, little is known about the reversion of frameshift mutations resulting in amino acid auxotrophy and on the structure and functioning of the protein encoded by the revertant mutated gene. The aims of this work were to analyze the appearance of reverse mutations over time and under different selective pressures and to investigate revertant enzymes' three-dimensional structures and their correlation with a different growth ability. Escherichia coli FB182 strain, carrying the hisF892 single nucleotide deletion resulting in histidine auxotrophy, was subjected to different selective pressures, and revertant mutants were isolated and characterized. The obtained results allowed us to identify different indels of different lengths located in different positions in the hisF gene, and relations with the incubation time and the selective pressure applied were observed. Moreover, the structure of the different mutant proteins was consistent with the respective revertant ability to grow in absence of histidine, highlighting a correlation between the mutations and the catalytic activity of the mutated HisF enzyme.

Protein sequence design with a learned potential

The task of protein sequence design is central to nearly all rational protein engineering problems, and enormous effort has gone into the development of energy functions to guide design. Here, we investigate the capability of a deep neural network model to automate design of sequences onto protein backbones, having learned directly from crystal structure data and without any human-specified priors. The model generalizes to native topologies not seen during training, producing experimentally stable designs. We evaluate the generalizability of our method to a de novo TIM-barrel scaffold. The model produces novel sequences, and high-resolution crystal structures of two designs show excellent agreement with in silico models. Our findings demonstrate the tractability of an entirely learned method for protein sequence design.

Rational protein design to achieve a given protein backbone conformation is needed to engineer specific functions. Here Anand et al. describe a machine learning method using a learned neural network potential for fixed-backbone protein design.

Novel protein from larval sponge cells, ilborin, is related to energy turnover and calcium binding and is conserved among marine invertebrates

Sponges (phylum Porifera) are early-branching animals, whose outwardly simple body plan is underlain by a complex genetic repertoire. The transition from a mobile larva to an attached filter-feeding organism occurs by metamorphosis, a process accompanied by a radical change of the body plan and cell transdifferentiation. The continuity between larval cells and adult tissues is still obscure. In a previous study, we have produced polyclonal antibodies against the major protein of the flagellated cells covering the larva of the sponge Halisarca dujardini, used them to trace the fate of these cells and shown that the larval flagellated cells transdifferentiate into the choanocytes. In the present work, we identified the sequence of this novel protein, which we named ilborin. A search in the open databases showed that multiple orthologues of the newly identified protein are present in sponges, cnidarians, flatworms, ctenophores and echinoderms, but none of them has been described yet. Ilborin has two conserved domains: triosephosphate isomerase-barrel, which has enzymatic activity against macroergic compounds, and canonical EF-hand, which binds calcium. mRNA of ilborin is expressed in the larval flagellated cells. We suggest that the new protein is involved in the calcium-mediated regulation of energy metabolism, whose activation precedes metamorphosis.

Glycosyl hydrolases family 5, subfamily 5: Relevance and structural insights for designing improved biomass degrading cocktails

Endoglucanases are carbohydrate-degrading enzymes widely used for bioethanol production as part of the enzymatic cocktail. However, family 5 subfamily 5 (GH5_5) endoglucanases are still poorly explored in depth. The Trichoderma reesei representative is the most studied enzyme, presenting catalytic activity in acidic media and mild temperature conditions. Though biochemically similar, its modular structure and synergy with other components vary greatly compared to other GH5_5 members and there is still a lack of specific studies regarding their interaction with other cellulases and application on novel and better mixtures. In this regard, the threedimensional structure elucidation is a highly valuable tool to both uncover basic catalytic mechanisms and implement engineering techniques, proved by the high success rate GH5_5 endoglucanases show. GH5_5 enzymes must be carefully evaluated to fully uncover their potential in biomass-degrading cocktails: the optimal industrial conditions, synergy with other cellulases, structural studies, and enzyme engineering approaches. We aimed to provide the current understanding of these main topics, collecting all available information about characterized GH5_5 endoglucanases function, structure, and bench experiments, in order to suggest future directions to a better application of these enzymes in the industry.

Expanded analyses of the functional correlations within structural classifications of glycoside hydrolases

Glycoside hydrolases (GHs) are greatly diverse in sequences and functions, but systematic studies of GH relationships based on structural information are lacking. Here, we report that GHs have multiple evolutionary origins and are structurally derived from 27 homologous superfamilies and 16 folds, but GHs are highly biased to distribute in a few superfamilies and folds. Six of these superfamilies are widely encoded by archaea, bacteria, and eukaryotes, indicating that they may be the most ancient in origin. Most superfamilies vary in enzyme function, and some, such as the superfamilies of (β/α)₈-barrel and (α/α)₆-barrel structures, exhibit extreme functional diversity; this is highly positively correlated with sequence diversity. More than one-third of glycosidase activities show a phenomenon of convergent evolution, especially the degradation functions of GHs on polysaccharides. The GHs of most superfamilies have relatively narrow environmental distributions, normally with the highest abundance in host-associated environments and a distribution preference for moderate low-temperature and acidic environments. Overall, our expanded analysis facilitates an understanding of complex GH sequence-structure-function relationships and may guide our screening and engineering of GHs.

Short non-coded peptides interacting with cofactors facilitated the integration of early chemical networks

Independently developed iron-sulphur/thioester- and phosphate-driven chemical reactions would have set up two distinct reaction networks prior to coupling in a proto-metabolic system supporting a minimal organisation closure. Each chemical system assisted initially by simple catalysts and then by more complex cofactors would have provided the precursors of the small metabolites and monomer units along with their respective polymers through dehydrating template-independent assemblies. For example, acylation reactions mediated by activated thioester groups produced peptides, fatty acids and polyhydroxyalkanoates, while phosphorylation reactions by phosphorylating agents allowed the synthesis of polysaccharides, polyribonucleotides and polyphosphates. Here, we address how these independent chemical systems might fit together and shaped a proto-metabolic system, focusing specifically on cofactors as molecular fossils of metabolism. As a result, the proposed overview suggests that non-coded peptides capable of binding a variety of ligands, but in particular with a redox active versatility and/or group transfer potential could have facilitated the chemical connections that led to a minimal closure with a proto-metabolism. Later developments would have made it possible to establish a cellular organisation with more complex and interdependent metabolic pathways.

Evolutionary Origins of DNA Repair Pathways: Role of Oxygen Catastrophe in the Emergence of DNA Glycosylases

It was proposed that the last universal common ancestor (LUCA) evolved under high temperatures in an oxygen-free environment, similar to those found in deep-sea vents and on volcanic slopes. Therefore, spontaneous DNA decay, such as base loss and cytosine deamination, was the major factor affecting LUCA’s genome integrity. Cosmic radiation due to Earth’s weak magnetic field and alkylating metabolic radicals added to these threats. Here, we propose that ancient forms of life had only two distinct repair mechanisms: versatile apurinic/apyrimidinic (AP) endonucleases to cope with both AP sites and deaminated residues, and enzymes catalyzing the direct reversal of UV and alkylation damage. The absence of uracil–DNA N-glycosylases in some Archaea, together with the presence of an AP endonuclease, which can cleave uracil-containing DNA, suggests that the AP endonuclease-initiated nucleotide incision repair (NIR) pathway evolved independently from DNA glycosylase-mediated base excision repair. NIR may be a relic that appeared in an early thermophilic ancestor to counteract spontaneous DNA damage. We hypothesize that a rise in the oxygen level in the Earth’s atmosphere ~2 Ga triggered the narrow specialization of AP endonucleases and DNA glycosylases to cope efficiently with a widened array of oxidative base damage and complex DNA lesions.

The enzymology of oxazolone and thioamide synthesis in methanobactin

Methanobactins are ribosomally synthesized and post-translationally modified peptidic (RiPP) natural products that are known for their ability to chelate copper ions. Crucial for their high copper affinity is a pair of bidentate ligands comprising a nitrogen-containing heterocycle and an adjacent thioamide or enethiol group. The previously uncharacterized proteins MbnB and MbnC were recently shown to synthesize these groups. In this chapter, we describe the methods that were used to determine that MbnB and MbnC are the core biosynthetic enzymes in methanobactin biosynthesis. The two proteins form a heterodimeric complex (MbnBC) which, through a dioxygen-dependent four-electron oxidation of the precursor peptide (MbnA), modifies a cysteine residue in order to install the oxazolone and thioamide moieties. This overview covers the heterologous expression and purification of MbnBC, characterization of the iron cluster found in MbnB, and characterization of the modification installed on MbnA. While this chapter is specific to MbnBC, the methods outlined here can be broadly applied to the enzymology of other proteins that install similar groups as well as enzyme pairs related to MbnB and MbnC.

The evolution of metabolism: How to test evolutionary hypotheses at the genomic level

The origin of primordial metabolism and its expansion to form the metabolic networks extant today represent excellent systems to study the impact of natural selection and the potential adaptive role of novel compounds. Here we present the current hypotheses made on the origin of life and ancestral metabolism and present the theories and mechanisms by which the large chemical diversity of plants might have emerged along evolution. In particular, we provide a survey of statistical methods that can be used to detect signatures of selection at the gene and population level, and discuss potential and limits of these methods for investigating patterns of molecular adaptation in plant metabolism.

The ammonia-lyases: enzymes that use a wide range of approaches to catalyze the same type of reaction

The paradigm that protein structure determines protein function has been clearly established. What is less clear is whether a specific protein structure is always required to carry out a specific function. Numerous cases are now known where there is no apparent connection between the biological function of a protein and the other members of its structural class, and where functionally related proteins can have quite diverse structures. A set of enzymes with these diverse properties, the ammonia-lyases, will be examined in this review. These are a class of enzymes that catalyze a relatively straightforward deamination reaction. However, the individual enzymes of this class possess a wide variety of different structures, utilize a diverse set of cofactors, and appear to catalyze this related reaction through a range of different mechanisms. This review aims to address a basic question: if there is not a specific protein structure and active site architecture that is both required and sufficient to define a catalyst for a given chemical reaction, then what factor(s) determine the structure and the mechanism that is selected to catalyze a particular reaction?

Insights Into the Origin of Life: Did It Begin from HCN and H2O?

The seminal Urey-Miller experiments showed that molecules crucial to life such as HCN could have formed in the reducing atmosphere of the Hadean Earth and then dissolved in the oceans. Subsequent proponents of the "RNA World" hypothesis have shown aqueous HCN to be the starting point for the formation of the precursors of RNA and proteins. However, the conditions of early Earth suggest that aqueous HCN would have had to react under a significant number of constraints. Therefore, given the limiting conditions, could RNA and protein precursors still have formed from aqueous HCN? If so, what mechanistic routes would have been followed? The current computational study, with the aid of the ab initio nanoreactor (AINR), a powerful new tool in computational chemistry, addresses these crucial questions. Gratifyingly, not only do the results from the AINR approach show that aqueous HCN could indeed have been the source of RNA and protein precursors, but they also indicate that just the interaction of HCN with water would have sufficed to begin a series of reactions leading to the precursors. The current work therefore provides important missing links in the story of prebiotic chemistry and charts the road from aqueous HCN to the precursors of RNA and proteins.

Graph Splitting: A Graph-Based Approach for Superfamily-Scale Phylogenetic Tree Reconstruction

A protein superfamily contains distantly related proteins that have acquired diverse biological functions through a long evolutionary history. Phylogenetic analysis of the early evolution of protein superfamilies is a key challenge because existing phylogenetic methods show poor performance when protein sequences are too diverged to construct an informative multiple sequence alignment (MSA). Here, we propose the Graph Splitting (GS) method, which rapidly reconstructs a protein superfamily-scale phylogenetic tree using a graph-based approach. Evolutionary simulation showed that the GS method can accurately reconstruct phylogenetic trees and be robust to major problems in phylogenetic estimation, such as biased taxon sampling, heterogeneous evolutionary rates, and long-branch attraction when sequences are substantially diverge. Its application to an empirical data set of the triosephosphate isomerase (TIM)-barrel superfamily suggests rapid evolution of protein-mediated pyrimidine biosynthesis, likely taking place after the RNA world. Furthermore, the GS method can also substantially improve performance of widely used MSA methods by providing accurate guide trees.

Protein Flexibility and Stiffness Enable Efficient Enzymatic Catalysis

The enormous rate accelerations observed for many enzyme catalysts are due to strong stabilizing interactions between the protein and reaction transition state. The defining property of these catalysts is their specificity for binding the transition state with a much higher affinity than substrate. Experimental results are presented which show that the phosphodianion-binding energy of phosphate monoester substrates is used to drive conversion of their protein catalysts from flexible and entropically rich ground states to stiff and catalytically active Michaelis complexes. These results are generalized to other enzyme-catalyzed reactions. The existence of many enzymes in flexible, entropically rich, and inactive ground states provides a mechanism for utilization of ligand-binding energy to mold these catalysts into stiff and active forms. This reduces the substrate-binding energy expressed at the Michaelis complex, while enabling the full and specific expression of large transition-state binding energies. Evidence is presented that the complexity of enzyme conformational changes increases with increases in the enzymatic rate acceleration. The requirement that a large fraction of the total substrate-binding energy be utilized to drive conformational changes of floppy enzymes is proposed to favor the selection and evolution of protein folds with multiple flexible unstructured loops, such as the TIM-barrel fold. The effect of protein motions on the kinetic parameters for enzymes that undergo ligand-driven conformational changes is considered. The results of computational studies to model the complex ligand-driven conformational change in catalysis by triosephosphate isomerase are presented.

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

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

Methanobactins: from genome to function

Methanobactins (Mbns) are ribosomally produced, post-translationally modified peptide (RiPP) natural products that bind copper with high affinity using nitrogen-containing heterocycles and thioamide groups. In some methanotrophic bacteria, Mbns are secreted under conditions of copper starvation and then re-internalized as a copper source for the enzyme particulate methane monooxygenase (pMMO). Genome mining studies have led to the identification and classification of operons encoding the Mbn precursor peptide (MbnA) as well as a number of putative transport, regulatory, and biosynthetic proteins. These Mbn operons are present in non-methanotrophic bacteria as well, suggesting a broader role in and perhaps beyond copper acquisition. Genetic and biochemical studies indicate that specific operon-encoded proteins are involved in Mbn transport and provide insight into copper-responsive gene regulation in methanotrophs. Mbn biosynthesis is not yet understood, but combined analysis of Mbn structures, MbnA sequences, and operon content represents a powerful approach to elucidating the roles of specific biosynthetic enzymes. Future work will likely lead to the discovery of unique pathways for natural product biosynthesis and new mechanisms of microbial metal homeostasis.

A Paradigm for CH Bond Cleavage: Structural and Functional Aspects of Transition State Stabilization by Mandelate Racemase

Mandelate racemase (MR) from Pseudomonas putida catalyzes the Mg²⁺-dependent, 1,1-proton transfer reaction that racemizes (R)- and (S)-mandelate. MR shares a partial reaction (i.e., the metal ion-assisted, Brønsted base-catalyzed proton abstraction of the α-proton of carboxylic acid substrates) and structural features ((β/α)₇β-barrel and N-terminal α + β capping domains) with a vast group of homologous, yet functionally diverse, enzymes in the enolase superfamily. Mechanistic and structural studies have developed this enzyme into a paradigm for understanding how enzymes such as those of the enolase superfamily overcome kinetic and thermodynamic barriers to catalyze the abstraction of an α-proton from a carbon acid substrate with a relatively high pK_a value. Structural studies on MR bound to intermediate/transition state analogues have delineated those structural features that MR uses to stabilize transition states and enhance reaction rates of proton abstraction. Kinetic, site-directed mutagenesis, and structural studies have also revealed that the phenyl ring of the substrate migrates through the hydrophobic cavity within the active site during catalysis and that the Brønsted acid-base catalysts (Lys 166 and His 297) may be utilized as binding determinants for inhibitor recognition. In addition, structural studies on the adduct formed from the irreversible inhibition of MR by 3-hydroxypyruvate revealed that MR can form and deprotonate a Schiff-base with 3-hydroxypyruvate to yield an enol(ate)-aldehyde adduct, suggesting a possible evolutionary link between MR and the Schiff-base forming aldolases. As the archetype of the enolase superfamily, mechanistic and structural studies on MR will continue to enhance our understanding of enzyme catalysis and furnish insights into the evolution of enzyme function.

Specificity Effects of Amino Acid Substitutions in Promiscuous Hydrolases: Context-Dependence of Catalytic Residue Contributions to Local Fitness Landscapes in Nearby Sequence Space

Catalytic promiscuity can facilitate evolution of enzyme functions-a multifunctional catalyst may act as a springboard for efficient functional adaptation. We test the effect of single mutations on multiple activities in two groups of promiscuous AP superfamily members to probe this hypothesis. We quantify the effect of site-saturating mutagenesis of an analogous, nucleophile-flanking residue in two superfamily members: an arylsulfatase (AS) and a phosphonate monoester hydrolase (PMH). Statistical analysis suggests that no one physicochemical characteristic alone explains the mutational effects. Instead, these effects appear to be dominated by their structural context. Likewise, the effect of changing the catalytic nucleophile itself is not reaction-type-specific. Mapping of "fitness landscapes" of four activities onto the possible variation of a chosen sequence position revealed tremendous potential for respecialization of AP superfamily members through single-point mutations, highlighting catalytic promiscuity as a powerful predictor of adaptive potential.

Lysine Biosynthesis of Thermococcus kodakarensis with the Capacity to Function as an Ornithine Biosynthetic System

We recently discovered a biosynthetic system using a novel amino group carrier protein called LysW for lysine biosynthesis via α-aminoadipate (AAA), and revealed that this system is also utilized in the biosynthesis of arginine by Sulfolobus In the present study, we focused on the biosynthesis of lysine and ornithine in the hyperthermophilic archaeon Thermococcus kodakarensis, and showed that their biosynthesis is accomplished by a single set of metabolic enzymes. We also determined the crystal structure of the LysX family protein from T. kodakarensis, which catalyzes the conjugation of LysW with either AAA or glutamate, in a complex with LysW-γ-AAA. This crystal structure is the first example to show how LysX recognizes AAA as a substrate and provides a structural basis for the bifunctionality of the LysX family protein from T. kodakarensis Based on comparisons with other LysX family proteins, we propose a mechanism for substrate recognition and its relationship with molecular evolution among LysX family proteins, which have different substrate specificities.

Structural and mechanistic insights into a Bacteroides vulgatus retaining N-acetyl-β-galactosaminidase that uses neighbouring group participation

Bacteroides vulgatus is a member of the human microbiota whose abundance is increased in patients with Crohn's disease. We show that a B. vulgatus glycoside hydrolase from the carbohydrate active enzyme family GH123, BvGH123, is an N-acetyl-β-galactosaminidase that acts with retention of stereochemistry, and, through a 3-D structure in complex with Gal-thiazoline, provide evidence in support of a neighbouring group participation mechanism.

The TIM Barrel Architecture Facilitated the Early Evolution of Protein-Mediated Metabolism

The triosephosphate isomerase (TIM) barrel protein fold is a structurally repetitive architecture that is present in approximately 10 % of all enzymes. It is generally assumed that this ubiquity in modern proteomes reflects an essential historical role in early protein-mediated metabolism. Here, we provide quantitative and comparative analyses to support several hypotheses about the early importance of the TIM barrel architecture. An information theoretical analysis of protein structures supports the hypothesis that the TIM barrel architecture could arise more easily by duplication and recombination compared to other mixed α/β structures. We show that TIM barrel enzymes corresponding to the most taxonomically broad superfamilies also have the broadest range of functions, often aided by metal and nucleotide-derived cofactors that are thought to reflect an earlier stage of metabolic evolution. By comparison to other putatively ancient protein architectures, we find that the functional diversity of TIM barrel proteins cannot be explained simply by their antiquity. Instead, the breadth of TIM barrel functions can be explained, in part, by the incorporation of a broad range of cofactors, a trend that does not appear to be shared by proteins in general. These results support the hypothesis that the simple and functionally general TIM barrel architecture may have arisen early in the evolution of protein biosynthesis and provided an ideal scaffold to facilitate the metabolic transition from ribozymes, peptides, and geochemical catalysts to modern protein enzymes.

Inactivation of Mandelate Racemase by 3-Hydroxypyruvate Reveals a Potential Mechanistic Link between Enzyme Superfamilies

Mandelate racemase (MR), a member of the enolase superfamily, catalyzes the Mg(2+)-dependent interconversion of the enantiomers of mandelate. Several α-keto acids are modest competitive inhibitors of MR [e.g., mesoxalate (Ki = 1.8 ± 0.3 mM) and 3-fluoropyruvate (Ki = 1.3 ± 0.1 mM)], but, surprisingly, 3-hydroxypyruvate (3-HP) is an irreversible, time-dependent inhibitor (kinact/KI = 83 ± 8 M(-1) s(-1)). Protection from inactivation by the competitive inhibitor benzohydroxamate, trypsinolysis and electrospray ionization tandem mass spectrometry analyses, and X-ray crystallographic studies reveal that 3-HP undergoes Schiff-base formation with Lys 166 at the active site, followed by formation of an aldehyde/enol(ate) adduct. Such a reaction is unprecedented in the enolase superfamily and may be a relic of an activity possessed by a promiscuous progenitor enzyme. The ability of MR to form and deprotonate a Schiff-base intermediate furnishes a previously unrecognized mechanistic link to other α/β-barrel enzymes utilizing Schiff-base chemistry and is in accord with the sequence- and structure-based hypothesis that members of the metal-dependent enolase superfamily and the Schiff-base-forming N-acetylneuraminate lyase superfamily and aldolases share a common ancestor.

ECOD: an evolutionary classification of protein domains

Understanding the evolution of a protein, including both close and distant relationships, often reveals insight into its structure and function. Fast and easy access to such up-to-date information facilitates research. We have developed a hierarchical evolutionary classification of all proteins with experimentally determined spatial structures, and presented it as an interactive and updatable online database. ECOD (Evolutionary Classification of protein Domains) is distinct from other structural classifications in that it groups domains primarily by evolutionary relationships (homology), rather than topology (or "fold"). This distinction highlights cases of homology between domains of differing topology to aid in understanding of protein structure evolution. ECOD uniquely emphasizes distantly related homologs that are difficult to detect, and thus catalogs the largest number of evolutionary links among structural domain classifications. Placing distant homologs together underscores the ancestral similarities of these proteins and draws attention to the most important regions of sequence and structure, as well as conserved functional sites. ECOD also recognizes closer sequence-based relationships between protein domains. Currently, approximately 100,000 protein structures are classified in ECOD into 9,000 sequence families clustered into close to 2,000 evolutionary groups. The classification is assisted by an automated pipeline that quickly and consistently classifies weekly releases of PDB structures and allows for continual updates. This synchronization with PDB uniquely distinguishes ECOD among all protein classifications. Finally, we present several case studies of homologous proteins not recorded in other classifications, illustrating the potential of how ECOD can be used to further biological and evolutionary studies.

Evolutionary relationship of two ancient protein superfolds

Proteins are the molecular machines of the cell that fold into specific three-dimensional structures to fulfill their functions. To improve our understanding of how the structure and function of proteins arises, it is crucial to understand how evolution has generated the structural diversity we observe today. Classically, proteins that adopt different folds are considered to be nonhomologous. However, using state-of-the-art tools for homology detection, we found evidence of homology between proteins of two ancient and highly populated protein folds, the (βα)8-barrel and the flavodoxin-like fold. We detected a family of sequences that show intermediate features between both folds and determined what is to our knowledge the first representative crystal structure of one of its members, giving new insights into the evolutionary link of two of the earliest folds. Our findings contribute to an emergent vision where protein superfolds share common ancestry and encourage further approaches to complete the mapping of structure space onto sequence space.

Biochemical fossils of the ancient transition from geoenergetics to bioenergetics in prokaryotic one carbon compound metabolism

The deep dichotomy of archaea and bacteria is evident in many basic traits including ribosomal protein composition, membrane lipid synthesis, cell wall constituents, and flagellar composition. Here we explore that deep dichotomy further by examining the distribution of genes for the synthesis of the central carriers of one carbon units, tetrahydrofolate (H4F) and tetrahydromethanopterin (H4MPT), in bacteria and archaea. The enzymes underlying those distinct biosynthetic routes are broadly unrelated across the bacterial-archaeal divide, indicating that the corresponding pathways arose independently. That deep divergence in one carbon metabolism is mirrored in the structurally unrelated enzymes and different organic cofactors that methanogens (archaea) and acetogens (bacteria) use to perform methyl synthesis in their H4F- and H4MPT-dependent versions, respectively, of the acetyl-CoA pathway. By contrast, acetyl synthesis in the acetyl-CoA pathway - from a methyl group, CO2 and reduced ferredoxin - is simpler, uniform and conserved across acetogens and methanogens, and involves only transition metals as catalysts. The data suggest that the acetyl-CoA pathway, while being the most ancient of known CO2 assimilation pathways, reflects two phases in early evolution: an ancient phase in a geochemically confined and non-free-living universal common ancestor, in which acetyl thioester synthesis proceeded spontaneously with the help of geochemically supplied methyl groups, and a later phase that reflects the primordial divergence of the bacterial and archaeal stem groups, which independently invented genetically-encoded means to synthesize methyl groups via enzymatic reactions. This article is part of a Special Issue entitled: 18th European Bioenergetic Conference.

Bioorganometallic chemistry with IspG and IspH: structure, function, and inhibition of the [Fe(4)S(4)] proteins involved in isoprenoid biosynthesis

Enzymes of the methylerythritol phosphate pathway of isoprenoid biosynthesis are attractive anti-infective drug targets. The last two enzymes of this pathway, IspG and IspH, are [Fe4 S4 ] proteins that are not produced by humans and catalyze 2 H(+) / 2 e(-) reductions with novel mechanisms. In this Review, we summarize recent advances in structural, mechanistic, and inhibitory studies of these two enzymes. In particular, mechanistic proposals involving bioorganometallic intermediates are presented, and compared with other mechanistic possibilities. In addition, inhibitors based on substrate analogues as well as developed by rational design and compound-library screening, are discussed. The results presented support bioorganometallic catalytic mechanisms for IspG and IspH, and open up new routes to anti-infective drug design targeting [Fe4 S4 ] clusters in proteins.

A common structural scaffold in CTD phosphatases that supports distinct catalytic mechanisms

The phosphorylation and dephosphorylation of the carboxyl-terminal domain (CTD) of the largest RNA polymerase II (RNAPII) subunit is a critical regulatory checkpoint for transcription and mRNA processing. This CTD is unique to eukaryotic organisms and it contains multiple tandem-repeats with the consensus sequence Tyr(1) -Ser(2) -Pro(3) -Thr(4) -Ser(5) -Pro(6) -Ser(7) . Traditionally, CTD phosphatases that use metal-ion-independent (cysteine-based) and metal-ion-assisted (aspartate-based) catalytic mechanisms have been considered to belong to two independent groups. However, using structural comparisons we have identified a common structural scaffold in these two groups of CTD phosphatases. This common scaffold accommodates different catalytic processes with the same substrate specificity, in this case phospho-serine/threonine residues flanked by prolines. Furthermore, this scaffold provides a structural connection between two groups of protein tyrosine phosphatases (PTPs): Cys-based (classes I, II, and III) and Asp-based (class IV) PTPs. Redundancy in catalytic mechanisms is not infrequent and may arise in specific biological settings. To better understand the activity of the CTD phosphatases, we combined our structural analyses with data on CTD phosphatase expression in different human and mouse tissues. The results suggest that aspartate- and cysteine-based CTD-dephosphorylation acts in concert during cellular stress, when high levels of reactive oxygen species can inhibit the nucleophilic function of the catalytic cysteine, as occurs in mental and neurodegenerative disorders like schizophrenia, Alzheimer's and Parkinson's diseases. Moreover, these findings have significant implications for the study of the RNAPII-CTD dephosphorylation in eukaryotes.

Genome mining for methanobactins

Background

Methanobactins (Mbns) are a family of copper-binding natural products involved in copper uptake by methanotrophic bacteria. The few Mbns that have been structurally characterized feature copper coordination by two nitrogen-containing heterocycles next to thioamide groups embedded in a peptidic backbone of varying composition. Mbns are proposed to derive from post-translational modification of ribosomally synthesized peptides, but only a few genes encoding potential precursor peptides have been identified. Moreover, the relevance of neighboring genes in these genomes has been unclear.

Results

The potential for Mbn production in a wider range of bacterial species was assessed by mining microbial genomes. Operons encoding Mbn-like precursor peptides, MbnAs, were identified in 16 new species, including both methanotrophs and, surprisingly, non-methanotrophs. Along with MbnA, the core of the operon is formed by two putative biosynthetic genes denoted MbnB and MbnC. The species can be divided into five groups on the basis of their MbnA and MbnB sequences and their operon compositions. Additional biosynthetic proteins, including aminotransferases, sulfotransferases and flavin adenine dinucleotide (FAD)-dependent oxidoreductases were also identified in some families. Beyond biosynthetic machinery, a conserved set of transporters was identified, including MATE multidrug exporters and TonB-dependent transporters. Additional proteins of interest include a di-heme cytochrome c peroxidase and a partner protein, the roles of which remain a mystery.

Conclusions

This study indicates that Mbn-like compounds may be more widespread than previously thought, but are not present in all methanotrophs. This distribution of species suggests a broader role in metal homeostasis. These data provide a link between precursor peptide sequence and Mbn structure, facilitating predictions of new Mbn structures and supporting a post-translational modification biosynthetic pathway. In addition, testable models for Mbn transport and for methanotrophic copper regulation have emerged. Given the unusual modifications observed in Mbns characterized thus far, understanding the roles of the putative biosynthetic proteins is likely to reveal novel pathways and chemistry.

Accurate prediction of protein enzymatic class by N-to-1 Neural Networks

We present a novel ab initio predictor of protein enzymatic class. The predictor can classify proteins, solely based on their sequences, into one of six classes extracted from the enzyme commission (EC) classification scheme and is trained on a large, curated database of over 6,000 non-redundant proteins which we have assembled in this work. The predictor is powered by an ensemble of N-to-1 Neural Network, a novel architecture which we have recently developed. N-to-1 Neural Networks operate on the full sequence and not on predefined features. All motifs of a predefined length (31 residues in this work) are considered and are compressed by an N-to-1 Neural Network into a feature vector which is automatically determined during training. We test our predictor in 10-fold cross-validation and obtain state of the art results, with a 96% correct classification and 86% generalized correlation. All six classes are predicted with a specificity of at least 80% and false positive rates never exceeding 7%. We are currently investigating enhanced input encoding schemes which include structural information, and are analyzing trained networks to mine motifs that are most informative for the prediction, hence, likely, functionally relevant.

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

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

Polymorphic toxin systems: Comprehensive characterization of trafficking modes, processing, mechanisms of action, immunity and ecology using comparative genomics

BACKGROUND

Proteinaceous toxins are observed across all levels of inter-organismal and intra-genomic conflicts. These include recently discovered prokaryotic polymorphic toxin systems implicated in intra-specific conflicts. They are characterized by a remarkable diversity of C-terminal toxin domains generated by recombination with standalone toxin-coding cassettes. Prior analysis revealed a striking diversity of nuclease and deaminase domains among the toxin modules. We systematically investigated polymorphic toxin systems using comparative genomics, sequence and structure analysis.

RESULTS

Polymorphic toxin systems are distributed across all major bacterial lineages and are delivered by at least eight distinct secretory systems. In addition to type-II, these include type-V, VI, VII (ESX), and the poorly characterized "Photorhabdus virulence cassettes (PVC)", PrsW-dependent and MuF phage-capsid-like systems. We present evidence that trafficking of these toxins is often accompanied by autoproteolytic processing catalyzed by HINT, ZU5, PrsW, caspase-like, papain-like, and a novel metallopeptidase associated with the PVC system. We identified over 150 distinct toxin domains in these systems. These span an extraordinary catalytic spectrum to include 23 distinct clades of peptidases, numerous previously unrecognized versions of nucleases and deaminases, ADP-ribosyltransferases, ADP ribosyl cyclases, RelA/SpoT-like nucleotidyltransferases, glycosyltranferases and other enzymes predicted to modify lipids and carbohydrates, and a pore-forming toxin domain. Several of these toxin domains are shared with host-directed effectors of pathogenic bacteria. Over 90 families of immunity proteins might neutralize anywhere between a single to at least 27 distinct types of toxin domains. In some organisms multiple tandem immunity genes or immunity protein domains are organized into polyimmunity loci or polyimmunity proteins. Gene-neighborhood-analysis of polymorphic toxin systems predicts the presence of novel trafficking-related components, and also the organizational logic that allows toxin diversification through recombination. Domain architecture and protein-length analysis revealed that these toxins might be deployed as secreted factors, through directed injection, or via inter-cellular contact facilitated by filamentous structures formed by RHS/YD, filamentous hemagglutinin and other repeats. Phyletic pattern and life-style analysis indicate that polymorphic toxins and polyimmunity loci participate in cooperative behavior and facultative 'cheating' in several ecosystems such as the human oral cavity and soil. Multiple domains from these systems have also been repeatedly transferred to eukaryotes and their viruses, such as the nucleo-cytoplasmic large DNA viruses.

CONCLUSIONS

Along with a comprehensive inventory of toxins and immunity proteins, we present several testable predictions regarding active sites and catalytic mechanisms of toxins, their processing and trafficking and their role in intra-specific and inter-specific interactions between bacteria. These systems provide insights regarding the emergence of key systems at different points in eukaryotic evolution, such as ADP ribosylation, interaction of myosin VI with cargo proteins, mediation of apoptosis, hyphal heteroincompatibility, hedgehog signaling, arthropod toxins, cell-cell interaction molecules like teneurins and different signaling messengers.

Out of fuzzy chemistry: from prebiotic chemistry to metabolic networks

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

Folding mechanism of an extremely thermostable (βα)(8)-barrel enzyme: a high kinetic barrier protects the protein from denaturation

HisF, the cyclase subunit of imidazole glycerol phosphate synthase (ImGPS) from Thermotoga maritima, is an extremely thermostable (βα)(8)-barrel protein. We elucidated the unfolding and refolding mechanism of HisF. Its unfolding transition is reversible and adequately described by the two-state model, but 6 weeks is necessary to reach equilibrium (at 25 °C). During refolding, initially a burst-phase off-pathway intermediate is formed. The subsequent productive folding occurs in two kinetic phases with time constants of ~3 and ~20 s. They reflect a sequential process via an on-pathway intermediate, as revealed by stopped-flow double-mixing experiments. The final step leads to native HisF, which associates with the glutaminase subunit HisH to form the functional ImGPS complex. The conversion of the on-pathway intermediate to the native protein results in a 10(6)-fold increase of the time constant for unfolding from 89 ms to 35 h (at 4.0 M GdmCl) and thus establishes a high energy barrier to denaturation. We conclude that the extra stability of HisF is used for kinetic protection against unfolding. In its refolding mechanism, HisF resembles other (βα)(8)-barrel proteins.

Evolutionary relationships of microbial aromatic prenyltransferases

The linkage of isoprenoid and aromatic moieties, catalyzed by aromatic prenyltransferases (PTases), leads to an impressive diversity of primary and secondary metabolites, including important pharmaceuticals and toxins. A few years ago, a hydroxynaphthalene PTase, NphB, featuring a novel ten-stranded β-barrel fold was identified in Streptomyces sp. strain CL190. This fold, termed the PT-barrel, is formed of five tandem ααββ structural repeats and remained exclusive to the NphB family until its recent discovery in the DMATS family of indole PTases. Members of these two families exist only in fungi and bacteria, and all of them appear to catalyze the prenylation of aromatic substrates involved in secondary metabolism. Sequence comparisons using PSI-BLAST do not yield matches between these two families, suggesting that they may have converged upon the same fold independently. However, we now provide evidence for a common ancestry for the NphB and DMATS families of PTases. We also identify sequence repeats that coincide with the structural repeats in proteins belonging to these two families. Therefore we propose that the PT-barrel arose by amplification of an ancestral ααββ module. In view of their homology and their similarities in structure and function, we propose to group the NphB and DMATS families together into a single superfamily, the PT-barrel superfamily.

Horizontal and vertical growth of S. cerevisiae metabolic network

BACKGROUND

The growth and development of a biological organism is reflected by its metabolic network, the evolution of which relies on the essential gene duplication mechanism. There are two current views about the evolution of metabolic networks. The retrograde model hypothesizes that a pathway evolves by recruiting novel enzymes in a direction opposite to the metabolic flow. The patchwork model is instead based on the assumption that the evolution is based on the exploitation of broad-specificity enzymes capable of catalysing a variety of metabolic reactions.

RESULTS

We analysed a well-studied unicellular eukaryotic organism, S. cerevisiae, and studied the effect of the removal of paralogous gene products on its metabolic network. Our results, obtained using different paralog and network definitions, show that, after an initial period when gene duplication was indeed instrumental in expanding the metabolic space, the latter reached an equilibrium and subsequent gene duplications were used as a source of more specialized enzymes rather than as a source of novel reactions. We also show that the switch between the two evolutionary strategies in S. cerevisiae can be dated to about 350 million years ago.

CONCLUSIONS

Our data, obtained through a novel analysis methodology, strongly supports the hypothesis that the patchwork model better explains the more recent evolution of the S. cerevisiae metabolic network. Interestingly, the effects of a patchwork strategy acting before the Euascomycete-Hemiascomycete divergence are still detectable today.

The molecular origins of evolutionary innovations

The history of life is a history of evolutionary innovations, qualitatively new phenotypic traits that endow their bearers with new, often game-changing abilities. We know many individual examples of innovations and their natural history, but we know little about the fundamental principles of phenotypic variability that permit new phenotypes to arise. Most phenotypic innovations result from changes in three classes of systems: metabolic networks, regulatory circuits, and macromolecules. I here highlight two important features that these classes of systems share. The first is the ubiquity of vast genotype networks - connected sets of genotypes with the same phenotype. The second is the great phenotypic diversity of small neighborhoods around different genotypes in genotype space. I here explain that both features are essential for the phenotypic variability that can bring forth qualitatively new phenotypes. Both features emerge from a common cause, the robustness of phenotypes to perturbations, whose origins are linked to life in changing environments.

Phylogenetics in the bioinformatics culture of understanding

Bioinformatics, as a relatively young discipline, has grown up in a world of high-throughput large volume data that requires automatic analysis to enable us to stay on top of it all. As a response, the bioinformatics discipline has developed strategies to find patterns in a ‘low signal : noise ratio’ environment. While the need to process large amounts of information and extract hypotheses is both laudable and inescapable, the pressures that such requirements have introduced can lead to short cuts and misapprehensions. This is particularly the case with reference to assumptions about the underlying evolutionary theories that are implicitly invoked by the algorithms utilised in the analysis pipelines. The classic example is the misuse of the term ‘homologous’ to mean ‘similar’ or even ‘functionally similar’, rather than the correct definition of ‘having the same evolutionary origin’, which may or may not imply similarity of function. In this review, we outline some of the common phylogenetic questions from a bioinformatics perspective that can be better addressed with a deeper understanding of evolutionary principles and show, with examples from the amidohydrolase and Toll families, that quite different conclusions can be drawn if such approaches are taken. This review focuses on the importance of the underlying evolutionary biology, rather than assessing the merits of different phylogenetic techniques. The relative merits of a priori and a posteriori inclusion of biological information are discussed.

The low cost of recombination in creating novel phenotypes: Recombination can create new phenotypes while disrupting well-adapted phenotypes much less than mutation

Recombination is often considered a disruptive force for well-adapted phenotypes, but recent evidence suggests that this cost of recombination can be small. A key benefit of recombination is that it can help create proteins and regulatory circuits with novel and useful phenotypes more efficiently than point mutation. Its effectiveness stems from the large-scale reorganization of genotypes that it causes, which can help explore far-flung regions in genotype space. Recent work on complex phenotypes in model gene regulatory circuits and proteins shows that the disruptive effects of recombination can be very mild compared to the effects of mutation. Recombination thus can have great benefits at a modest cost, but we do not understand the reasons well. A better understanding might shed light on the evolution of recombination and help improve evolutionary strategies in biochemical engineering.

Structural and functional energetic linkages in allosteric regulation of muscle pyruvate kinase

The understanding of the molecular mechanisms of allostery in rabbit muscle pyruvate kinase (RMPK) is still in its infancy. Although, there is a paucity of knowledge on the ground rules on how its functions are regulated, RMPK is an ideal system to address basic questions regarding the fundamental chemical principles governing the regulatory mechanisms about this enzyme which has a TIM (α/β)(8) barrel structural motif [Copley, R. R., and Bork, P. (2000). Homology among (βα)8 barrels: Implications for the evolution of metabolic pathways. J. Mol. Biol.303, 627-640; Farber, G. K., and Petsko, G. A. (1990). The evolution of α/ß barrel enzymes. Trends Biochem.15, 228-234; Gerlt, J. A., and Babbitt, P. C. (2001). Divergent evolution of enzymatic function: Mechanistically diverse superfamilies and functionally distinct superfamilies. Annu. Rev. Biochem.70, 209-246; Heggi, H., and Gerstein, M. (1999). The relationship between protein structure and function: A comprehensive survey with application to the yeast genome. J. Mol. Biol.288, 147-164; Wierenga, R. K. (2001). The TIM-barrel fold: A versatile framework for efficient enzymes. FEB Lett.492, 193-198]. RMPK is a homotetramer. Each subunit consists of 530 amino acids and multiple domains. The active site resides between the A and B domains. Besides the basic TIM-barrel motif, RMPK also exhibits looped-out regions in the α/β barrel of each monomer forming the B- and C-domains. The two isozymes of PK, namely, the kidney and muscle isozymes, exhibit very different allosteric behaviors under the same experimental condition. The only amino acid sequence differences between the mammalian kidney and muscle PK isozymes are located in the C-domain and are involved in intersubunit interactions. Thus, embedded in these two isozymes of PK are the rules involved in engineering the popular TIM (α/β)(8) motif to modulate its allosteric properties. The PK system exhibits a lot of the properties that will allow mining of the ground rules governing the correlative linkages between sequence-fold-function. In this chapter, we review the approaches to acquire the fundamental functional and structural energetics that establish the linkages among this intricate network of linked multiequilibria. Results from these diverse approaches are integrated to establish a working model to represent the complex network of multiple linked reactions which ultimately leads to the observation of allosteric regulation of PK.

In silico evolution of early metabolism

We developed a simulation tool for investigating the evolution of early metabolism, allowing us to speculate on the formation of metabolic pathways from catalyzed chemical reactions and on the development of their characteristic properties. Our model consists of a protocellular entity with a simple RNA-based genetic system and an evolving metabolism of catalytically active ribozymes that manipulate a rich underlying chemistry. Ensuring an almost open-ended and fairly realistic simulation is crucial for understanding the first steps in metabolic evolution. We show here how our simulation tool can be helpful in arguing for or against hypotheses on the evolution of metabolic pathways. We demonstrate that seemingly mutually exclusive hypotheses may well be compatible when we take into account that different processes dominate different phases in the evolution of a metabolic system. Our results suggest that forward evolution shapes metabolic network in the very early steps of evolution. In later and more complex stages, enzyme recruitment supersedes forward evolution, keeping a core set of pathways from the early phase.