133,307 research outputs found
Environmental boundary conditions for the origin of life converge to an organo-sulfur metabolism
Published in final edited form as:
Nat Ecol Evol. 2019 December ; 3(12): 1715–1724. doi:10.1038/s41559-019-1018-8.It has been suggested that a deep memory of early life is hidden in the architecture of metabolic networks, whose reactions could have been catalyzed by small molecules or minerals before genetically encoded enzymes. A major challenge in unravelling these early steps is assessing the plausibility of a connected, thermodynamically consistent proto-metabolism under different geochemical conditions, which are still surrounded by high uncertainty. Here we combine network-based algorithms with physico-chemical constraints on chemical reaction networks to systematically show how different combinations of parameters (temperature, pH, redox potential and availability of molecular precursors) could have affected the evolution of a proto-metabolism. Our analysis of possible trajectories indicates that a subset of boundary conditions converges to an organo-sulfur-based proto-metabolic network fuelled by a thioester- and redox-driven variant of the reductive tricarboxylic acid cycle that is capable of producing lipids and keto acids. Surprisingly, environmental sources of fixed nitrogen and low-potential electron donors are not necessary for the earliest phases of biochemical evolution. We use one of these networks to build a steady-state dynamical metabolic model of a protocell, and find that different combinations of carbon sources and electron donors can support the continuous production of a minimal ancient 'biomass' composed of putative early biopolymers and fatty acids.80NSSC17K0295 - Intramural NASA; 80NSSC17K0296 - Intramural NASA; T32 GM100842 - NIGMS NIH HHSAccepted manuscrip
The Nondeterministic Waiting Time Algorithm: A Review
We present briefly the Nondeterministic Waiting Time algorithm. Our technique
for the simulation of biochemical reaction networks has the ability to mimic
the Gillespie Algorithm for some networks and solutions to ordinary
differential equations for other networks, depending on the rules of the
system, the kinetic rates and numbers of molecules. We provide a full
description of the algorithm as well as specifics on its implementation. Some
results for two well-known models are reported. We have used the algorithm to
explore Fas-mediated apoptosis models in cancerous and HIV-1 infected T cells
Modeling the complex dynamics of enzyme-pathway coevolution.
Peer reviewedPostprin
Modelling Cell Cycle using Different Levels of Representation
Understanding the behaviour of biological systems requires a complex setting
of in vitro and in vivo experiments, which attracts high costs in terms of time
and resources. The use of mathematical models allows researchers to perform
computerised simulations of biological systems, which are called in silico
experiments, to attain important insights and predictions about the system
behaviour with a considerably lower cost. Computer visualisation is an
important part of this approach, since it provides a realistic representation
of the system behaviour. We define a formal methodology to model biological
systems using different levels of representation: a purely formal
representation, which we call molecular level, models the biochemical dynamics
of the system; visualisation-oriented representations, which we call visual
levels, provide views of the biological system at a higher level of
organisation and are equipped with the necessary spatial information to
generate the appropriate visualisation. We choose Spatial CLS, a formal
language belonging to the class of Calculi of Looping Sequences, as the
formalism for modelling all representation levels. We illustrate our approach
using the budding yeast cell cycle as a case study
Autocatalytic sets in a partitioned biochemical network
In previous work, RAF theory has been developed as a tool for making
theoretical progress on the origin of life question, providing insight into the
structure and occurrence of self-sustaining and collectively autocatalytic sets
within catalytic polymer networks. We present here an extension in which there
are two "independent" polymer sets, where catalysis occurs within and between
the sets, but there are no reactions combining polymers from both sets. Such an
extension reflects the interaction between nucleic acids and peptides observed
in modern cells and proposed forms of early life.Comment: 28 pages, 8 figure
Split histidine kinases enable ultrasensitivity and bistability in two-component signaling networks
Bacteria sense and respond to their environment through signaling cascades generally referred to as two-component signaling networks. These networks comprise histidine kinases and their cognate response regulators. Histidine kinases have a number of biochemical activities: ATP binding, autophosphorylation, the ability to act as a phosphodonor for their response regulators, and in many cases the ability to catalyze the hydrolytic dephosphorylation of their response regulator. Here, we explore the functional role of “split kinases” where the ATP binding and phosphotransfer activities of a conventional histidine kinase are split onto two distinct proteins that form a complex. We find that this unusual configuration can enable ultrasensitivity and bistability in the signal-response relationship of the resulting system. These dynamics are displayed under a wide parameter range but only when specific biochemical requirements are met. We experimentally show that one of these requirements, namely segregation of the phosphatase activity predominantly onto the free form of one of the proteins making up the split kinase, is met in Rhodobacter sphaeroides. These findings indicate split kinases as a bacterial alternative for enabling ultrasensitivity and bistability in signaling networks. Genomic analyses reveal that up 1.7% of all identified histidine kinases have the potential to be split and bifunctional
Emergence of switch-like behavior in a large family of simple biochemical networks
Bistability plays a central role in the gene regulatory networks (GRNs)
controlling many essential biological functions, including cellular
differentiation and cell cycle control. However, establishing the network
topologies that can exhibit bistability remains a challenge, in part due to the
exceedingly large variety of GRNs that exist for even a small number of
components. We begin to address this problem by employing chemical reaction
network theory in a comprehensive in silico survey to determine the capacity
for bistability of more than 40,000 simple networks that can be formed by two
transcription factor-coding genes and their associated proteins (assuming only
the most elementary biochemical processes). We find that there exist reaction
rate constants leading to bistability in ~90% of these GRN models, including
several circuits that do not contain any of the TF cooperativity commonly
associated with bistable systems, and the majority of which could only be
identified as bistable through an original subnetwork-based analysis. A
topological sorting of the two-gene family of networks based on the presence or
absence of biochemical reactions reveals eleven minimal bistable networks
(i.e., bistable networks that do not contain within them a smaller bistable
subnetwork). The large number of previously unknown bistable network topologies
suggests that the capacity for switch-like behavior in GRNs arises with
relative ease and is not easily lost through network evolution. To highlight
the relevance of the systematic application of CRNT to bistable network
identification in real biological systems, we integrated publicly available
protein-protein interaction, protein-DNA interaction, and gene expression data
from Saccharomyces cerevisiae, and identified several GRNs predicted to behave
in a bistable fashion.Comment: accepted to PLoS Computational Biolog
- …