The compositional and evolutionary logic of metabolism

Metabolism displays striking and robust regularities in the forms of modularity and hierarchy, whose composition may be compactly described. This renders metabolic architecture comprehensible as a system, and suggests the order in which layers of that system emerged. Metabolism also serves as the foundation in other hierarchies, at least up to cellular integration including bioenergetics and molecular replication, and trophic ecology. The recapitulation of patterns first seen in metabolism, in these higher levels, suggests metabolism as a source of causation or constraint on many forms of organization in the biosphere. We identify as modules widely reused subsets of chemicals, reactions, or functions, each with a conserved internal structure. At the small molecule substrate level, module boundaries are generally associated with the most complex reaction mechanisms and the most conserved enzymes. Cofactors form a structurally and functionally distinctive control layer over the small-molecule substrate. Complex cofactors are often used at module boundaries of the substrate level, while simpler ones participate in widely used reactions. Cofactor functions thus act as "keys" that incorporate classes of organic reactions within biochemistry. The same modules that organize the compositional diversity of metabolism are argued to have governed long-term evolution. Early evolution of core metabolism, especially carbon-fixation, appears to have required few innovations among a small number of conserved modules, to produce adaptations to simple biogeochemical changes of environment. We demonstrate these features of metabolism at several levels of hierarchy, beginning with the small-molecule substrate and network architecture, continuing with cofactors and key conserved reactions, and culminating in the aggregation of multiple diverse physical and biochemical processes in cells.Comment: 56 pages, 28 figure

Analytical study of non Gaussian fluctuations in a stochastic scheme of autocatalytic reactions

A stochastic model of autocatalytic chemical reactions is studied both numerically and analytically. The van Kampen perturbative scheme is implemented, beyond the second order approximation, so to capture the non Gaussianity traits as displayed by the simulations. The method is targeted to the characterization of the third moments of the distribution of fluctuations, originating from a system of four populations in mutual interaction. The theory predictions agree well with the simulations, pointing to the validity of the van Kampen expansion beyond the conventional Gaussian solution.Comment: 15 pages, 8 figures, submitted to Phys. Rev.

Nonintegrability, Chaos, and Complexity

Two-dimensional driven dissipative flows are generally integrable via a conservation law that is singular at equilibria. Nonintegrable dynamical systems are confined to n*3 dimensions. Even driven-dissipative deterministic dynamical systems that are critical, chaotic or complex have n-1 local time-independent conservation laws that can be used to simplify the geometric picture of the flow over as many consecutive time intervals as one likes. Those conserevation laws generally have either branch cuts, phase singularities, or both. The consequence of the existence of singular conservation laws for experimental data analysis, and also for the search for scale-invariant critical states via uncontrolled approximations in deterministic dynamical systems, is discussed. Finally, the expectation of ubiquity of scaling laws and universality classes in dynamics is contrasted with the possibility that the most interesting dynamics in nature may be nonscaling, nonuniversal, and to some degree computationally complex

Unified analysis of terminal-time control in classical and quantum systems

Many phenomena in physics, chemistry, and biology involve seeking an optimal control to maximize an objective for a classical or quantum system which is open and interacting with its environment. The complexity of finding an optimal control for maximizing an objective is strongly affected by the possible existence of sub-optimal maxima. Within a unified framework under specified conditions, control objectives for maximizing at a terminal time physical observables of open classical and quantum systems are shown to be inherently free of sub-optimal maxima. This attractive feature is of central importance for enabling the discovery of controls in a seamless fashion in a wide range of phenomena transcending the quantum and classical regimes.Comment: 10 page

Signatures of arithmetic simplicity in metabolic network architecture

Metabolic networks perform some of the most fundamental functions in living cells, including energy transduction and building block biosynthesis. While these are the best characterized networks in living systems, understanding their evolutionary history and complex wiring constitutes one of the most fascinating open questions in biology, intimately related to the enigma of life's origin itself. Is the evolution of metabolism subject to general principles, beyond the unpredictable accumulation of multiple historical accidents? Here we search for such principles by applying to an artificial chemical universe some of the methodologies developed for the study of genome scale models of cellular metabolism. In particular, we use metabolic flux constraint-based models to exhaustively search for artificial chemistry pathways that can optimally perform an array of elementary metabolic functions. Despite the simplicity of the model employed, we find that the ensuing pathways display a surprisingly rich set of properties, including the existence of autocatalytic cycles and hierarchical modules, the appearance of universally preferable metabolites and reactions, and a logarithmic trend of pathway length as a function of input/output molecule size. Some of these properties can be derived analytically, borrowing methods previously used in cryptography. In addition, by mapping biochemical networks onto a simplified carbon atom reaction backbone, we find that several of the properties predicted by the artificial chemistry model hold for real metabolic networks. These findings suggest that optimality principles and arithmetic simplicity might lie beneath some aspects of biochemical complexity

Microbial catabolic activities are naturally selected by metabolic energy harvest rate

The fundamental trade-off between yield and rate of energy harvest per unit of substrate has been largely discussed as a main characteristic for microbial established cooperation or competition. In this study, this point is addressed by developing a generalized model that simulates competition between existing and not experimentally reported microbial catabolic activities defined only based on well-known biochemical pathways. No specific microbial physiological adaptations are considered, growth yield is calculated coupled to catabolism energetics and a common maximum biomass-specific catabolism rate (expressed as electron transfer rate) is assumed for all microbial groups. Under this approach, successful microbial metabolisms are predicted in line with experimental observations under the hypothesis of maximum energy harvest rate. Two microbial ecosystems, typically found in wastewater treatment plants, are simulated, namely: (i) the anaerobic fermentation of glucose and (ii) the oxidation and reduction of nitrogen under aerobic autotrophic (nitrification) and anoxic heterotrophic and autotrophic (denitrification) conditions. The experimentally observed cross feeding in glucose fermentation, through multiple intermediate fermentation pathways, towards ultimately methane and carbon dioxide is predicted. Analogously, two-stage nitrification (by ammonium and nitrite oxidizers) is predicted as prevailing over nitrification in one stage. Conversely, denitrification is predicted in one stage (by denitrifiers) as well as anammox (anaerobic ammonium oxidation). The model results suggest that these observations are a direct consequence of the different energy yields per electron transferred at the different steps of the pathways. Overall, our results theoretically support the hypothesis that successful microbial catabolic activities are selected by an overall maximum energy harvest rate

Toward homochiral protocells in noncatalytic peptide systems

The activation-polymerization-epimerization-depolymerization (APED) model of Plasson et al. has recently been proposed as a mechanism for the evolution of homochirality on prebiotic Earth. The dynamics of the APED model in two-dimensional spatially-extended systems is investigated for various realistic reaction parameters. It is found that the APED system allows for the formation of isolated homochiral proto-domains surrounded by a racemate. A diffusive slowdown of the APED network such as induced through tidal motion or evaporating pools and lagoons leads to the stabilization of homochiral bounded structures as expected in the first self-assembled protocells.Comment: 10 pages, 5 figure

Clades of huge phages from across Earth's ecosystems.

Bacteriophages typically have small genomes1 and depend on their bacterial hosts for replication2. Here we sequenced DNA from diverse ecosystems and found hundreds of phage genomes with lengths of more than 200 kilobases (kb), including a genome of 735 kb, which is-to our knowledge-the largest phage genome to be described to date. Thirty-five genomes were manually curated to completion (circular and no gaps). Expanded genetic repertoires include diverse and previously undescribed CRISPR-Cas systems, transfer RNAs (tRNAs), tRNA synthetases, tRNA-modification enzymes, translation-initiation and elongation factors, and ribosomal proteins. The CRISPR-Cas systems of phages have the capacity to silence host transcription factors and translational genes, potentially as part of a larger interaction network that intercepts translation to redirect biosynthesis to phage-encoded functions. In addition, some phages may repurpose bacterial CRISPR-Cas systems to eliminate competing phages. We phylogenetically define the major clades of huge phages from human and other animal microbiomes, as well as from oceans, lakes, sediments, soils and the built environment. We conclude that the large gene inventories of huge phages reflect a conserved biological strategy, and that the phages are distributed across a broad bacterial host range and across Earth's ecosystems

Bootstrapping the energy flow in the beginning of life.

This paper suggests that the energy flow on which all living structures depend only started up slowly, the low-energy, initial phase starting up a second, slightly more energetic phase, and so on. In this way, the build up of the energy flow follows a bootstrapping process similar to that found in the development of computers, the first generation making possible the calculations necessary for constructing the second one, etc. In the biogenetic upstart of an energy flow, non-metals in the lower periods of the Periodic Table of Elements would have constituted the most primitive systems, their operation being enhanced and later supplanted by elements in the higher periods that demand more energy. This bootstrapping process would put the development of the metabolisms based on the second period elements carbon, nitrogen and oxygen at the end of the evolutionary process rather than at, or even before, the biogenetic even

Genomic variation landscape of the human gut microbiome

While large-scale efforts have rapidly advanced the understanding and practical impact of human genomic variation, the latter is largely unexplored in the human microbiome. We therefore developed a framework for metagenomic variation analysis and applied it to 252 fecal metagenomes of 207 individuals from Europe and North America. Using 7.4 billion reads aligned to 101 reference species, we detected 10.3 million single nucleotide polymorphisms (SNPs), 107,991 short indels, and 1,051 structural variants. The average ratio of non-synonymous to synonymous polymorphism rates of 0.11 was more variable between gut microbial species than across human hosts. Subjects sampled at varying time intervals exhibited individuality and temporal stability of SNP variation patterns, despite considerable composition changes of their gut microbiota. This implies that individual-specific strains are not easily replaced and that an individual might have a unique metagenomic genotype, which may be exploitable for personalized diet or drug intake