Genotype networks in metabolic reaction spaces

Background: A metabolic genotype comprises all chemical reactions an organism can catalyze via enzymes encoded in its genome. A genotype is viable in a given environment if it is capable of producing all biomass components the organism needs to survive and reproduce. Previous work has focused on the properties of individual genotypes while little is known about how genome-scale metabolic networks with a given function can vary in their reaction content. Results: We here characterize spaces of such genotypes. Specifically, we study metabolic genotypes whose phenotype is viability in minimal chemical environments that differ in their sole carbon sources. We show that regardless of the number of reactions in a metabolic genotype, the genotypes of a given phenotype typically form vast, connected, and unstructured sets -- genotype networks -- that nearly span the whole of genotype space. The robustness of metabolic phenotypes to random reaction removal in such spaces has a narrow distribution with a high mean. Different carbon sources differ in the number of metabolic genotypes in their genotype network; this number decreases as a genotype is required to be viable on increasing numbers of carbon sources, but much less than if metabolic reactions were used independently across different chemical environments. Conclusions: Our work shows that phenotype-preserving genotype networks have generic organizational properties and that these properties are insensitive to the number of reactions in metabolic genotypes.Comment: 48 pages, 10 main figures, 14 supplementary figure

The organization of metabolic genotype space facilitates adaptive evolution in nitrogen metabolism

A metabolism is a complex chemical reaction system, whose metabolic genotype – the DNA encoding the enzymes catalyzing these reactions – can be compactly represented by its complement of metabolic reactions. Here, we analyze a space of such metabolic genotypes. Specifically, we study nitrogen metabolism and focus on nitrogen utilization phenotypes that are defined through the viability of a metabolism – its ability to synthesize all essential small biomass precursors – on a given combination of sole nitrogen sources. We randomly sample metabolisms with known phenotypes from metabolic genotype space with the aid of a method based on Markov Chain Monte Carlo sampling. We find that metabolisms viable on a given nitrogen source or a combination of nitrogen sources can differ in as much as 80 percent of their reactions, but can form networks of genotypes that are connected to one another through sequences of single reaction changes. The reactions that cannot vary in any one metabolism differ among metabolisms, and include a small core of “absolutely superessential” reactions that are required in all metabolisms we study. Only a small number of reaction changes are needed to reach the genotype network of one metabolic phenotype from the genotype network of another metabolic phenotype. Our observations indicate deep similarities between the genotype spaces of macromolecules, regulatory circuits, and metabolism that can facilitate the origin of novel phenotypes in evolution.  

Challenges in experimental data integration within genome-scale metabolic models

A report of the meeting "Challenges in experimental data integration within genome-scale metabolic models", Institut Henri Poincar\'e, Paris, October 10-11 2009, organized by the CNRS-MPG joint program in Systems Biology.Comment: 5 page

Genotype networks, innovation, and robustness in sulfur metabolism

Background: A metabolism is a complex network of chemical reactions. This network synthesizes multiple small precursor molecules of biomass from chemicals that occur in the environment. The metabolic network of any one organism is encoded by a metabolic genotype, defined as the set of enzyme-coding genes whose products catalyze the network's reactions. Each metabolic genotype has a metabolic phenotype. We define this metabolic phenotype as the spectrum of different sources of a chemical element that a metabolism can use to synthesize biomass. We here focus on the element sulfur. We study properties of the space of all possible metabolic genotypes in sulfur metabolism by analyzing random metabolic genotypes that are viable on different numbers of sulfur sources. Results: We show that metabolic genotypes with the same phenotype form large connected genotype networks - networks of metabolic networks - that extend far through metabolic genotype space. How far they reach through this space depends linearly on the number of super-essential reactions. A super-essential reaction is an essential reaction that occurs in all networks viable in a given environment. Metabolic networks can differ in how robust their phenotype is to the removal of individual reactions. We find that this robustness depends on metabolic network size, and on other variables, such as the size of minimal metabolic networks whose reactions are all essential in a specific environment. We show that different neighborhoods of any genotype network harbor very different novel phenotypes, metabolic innovations that can sustain life on novel sulfur sources. We also analyze the ability of evolving populations of metabolic networks to explore novel metabolic phenotypes. This ability is facilitated by the existence of genotype networks, because different neighborhoods of these networks contain very different novel phenotypes. Conclusions: We show that the space of metabolic genotypes involved in sulfur metabolism is organized similarly to that of carbon metabolism. We demonstrate that the maximum genotype distance and robustness of metabolic networks can be explained by the number of superessential reactions and by the sizes of minimal metabolic networks viable in an environment. In contrast to the genotype space of macromolecules, where phenotypic robustness may facilitate phenotypic innovation, we show that here the ability to access novel phenotypes does not monotonically increase with robustness

Networks in molecular evolution

Networks are a common theme at all levels of molecular evolution: Networks of metastable states and their connecting saddle points determine structure and folding kinetics of biopolymers. Neutral networks in sequence space explain the evolvability of both nucleic acids and polypeptides by linking Darwinian selection with neutral drift. Interacting replicators, be they simple molecules or highly complex mammals, form intricate ecological networks that are crucial for their survival. Chemical reactions are collected in extensive metabolic networks by means of specific enzymes; both the enzymes and the chemical reaction network that they govern undergoes evolutionary changes. Networks of gene regulation, protein-protein interaction, and cell signaling form the physico-chemical basis for morphogenesis and development. The nervous systems of higher animals form another distinct level of network architecture. We are beginning to understand the structure and function of each of the individual levels in some detail. Yet, their interplay largely remains still in the dark

A Sequence-to-Function Map for Ribozyme-catalyzed Metabolisms

We introduce a novel genotype-phenotype mapping based on the relation between RNA sequence and its secondary structure for the use in evolutionary studies. Various extensive studies concerning RNA folding in the context of neutral theory yielded insights about properties of the structure space and the mapping itself. We intend to get a better understanding of some of these properties and especially of the evolution of RNA-molecules as well as their effect on the evolution of the entire molecular system. We investigate the constitution of the neutral network and compare our mapping with other artificial approaches using cellular automatons, random boolean networks and others also based on RNA folding. We yield the highest extent, connectivity and evolvability of the underlying neutral network. Further, we successfully apply the mapping in an existing model for the evolution of a ribozyme-catalyzed metabolism

The cultural epigenetics of psychopathology: The missing heritability of complex diseases found?

We extend a cognitive paradigm for gene expression based on the asymptotic limit theorems of information theory to the epigenetic epidemiology of mental disorders. In particular, we recognize the fundamental role culture plays in human biology, another heritage mechanism parallel to, and interacting with, the more familiar genetic and epigenetic systems. We do this via a model through which culture acts as another tunable epigenetic catalyst that both directs developmental trajectories, and becomes convoluted with individual ontology, via a mutually-interacting crosstalk mediated by a social interaction that is itself culturally driven. We call for the incorporation of embedding culture as an essential component of the epigenetic regulation of human mental development and its dysfunctions, bringing what is perhaps the central reality of human biology into the center of biological psychiatry. Current US work on gene-environment interactions in psychiatry must be extended to a model of gene-environment-culture interaction to avoid becoming victim of an extreme American individualism that threatens to create paradigms particular to that culture and that are, indeed, peculiar in the context of the world's cultures. The cultural and epigenetic systems of heritage may well provide the 'missing' heritability of complex diseases now under so much intense discussion

Irreversible thermodynamics of open chemical networks I: Emergent cycles and broken conservation laws

In this and a companion paper we outline a general framework for the thermodynamic description of open chemical reaction networks, with special regard to metabolic networks regulating cellular physiology and biochemical functions. We first introduce closed networks "in a box", whose thermodynamics is subjected to strict physical constraints: the mass-action law, elementarity of processes, and detailed balance. We further digress on the role of solvents and on the seemingly unacknowledged property of network independence of free energy landscapes. We then open the system by assuming that the concentrations of certain substrate species (the chemostats) are fixed, whether because promptly regulated by the environment via contact with reservoirs, or because nearly constant in a time window. As a result, the system is driven out of equilibrium. A rich algebraic and topological structure ensues in the network of internal species: Emergent irreversible cycles are associated to nonvanishing affinities, whose symmetries are dictated by the breakage of conservation laws. These central results are resumed in the relation $a + b = s^Y$ between the number of fundamental affinities $a$, that of broken conservation laws $b$ and the number of chemostats $s^Y$. We decompose the steady state entropy production rate in terms of fundamental fluxes and affinities in the spirit of Schnakenberg's theory of network thermodynamics, paving the way for the forthcoming treatment of the linear regime, of efficiency and tight coupling, of free energy transduction and of thermodynamic constraints for network reconstruction.Comment: 18 page