Emergence and maintenance of functional modules in signaling pathways

© 2007 Soyer; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.Background: While detection and analysis of functional modules in biological systems have received great attention in recent years, we still lack a complete understanding of how such modules emerge. One theory is that systems must encounter a varying selection (i.e. environment) in order for modularity to emerge. Here, we provide an alternative and simpler explanation using a realistic model of biological signaling pathways and simulating their evolution.Results: These evolutionary simulations start with a homogenous population of a minimal pathway containing two effectors coupled to two signals via a single receptor. This population is allowed to evolve under a constant selection pressure for mediating two separate responses. Results of these evolutionary simulations show that under such a selective pressure, mutational processes easily lead to the emergence of pathways with two separate sub-pathways (i.e. modules) each mediating a distinct response only to one of the signals. Such functional modules are maintained as long as mutations leading to new interactions among existing proteins in the pathway are rare.Conclusion: While supporting a neutralistic view for the emergence of modularity in biological systems, these findings highlight the relevant rate of different mutational processes and the distribution of functional pathways in the topology space as key factors for its maintenance

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

Thermodynamic modelling of synthetic communities predicts minimum free energy requirements for sulfate reduction and methanogenesis

Microbial communities are complex dynamical systems harbouring many species interacting together to implement higher-level functions. Among these higher-level functions, conversion of organic matter into simpler building blocks by microbial communities underpins biogeochemical cycles and animal and plant nutrition, and is exploited in biotechnology. A prerequisite to predicting the dynamics and stability of community-mediated metabolic conversions is the development and calibration of appropriate mathematical models. Here, we present a generic, extendable thermodynamic model for community dynamics and calibrate a key parameter of this thermodynamic model, the minimum energy requirement associated with growth-supporting metabolic pathways, using experimental population dynamics data from synthetic communities composed of a sulfate reducer and two methanogens. Our findings show that accounting for thermodynamics is necessary in capturing the experimental population dynamics of these synthetic communities that feature relevant species using low energy growth pathways. Furthermore, they provide the first estimates for minimum energy requirements of methanogenesis (in the range of −30 kJ mol−1) and elaborate on previous estimates of lactate fermentation by sulfate reducers (in the range of −30 to −17 kJ mol−1 depending on the culture conditions). The open-source nature of the developed model and demonstration of its use for estimating a key thermodynamic parameter should facilitate further thermodynamic modelling of microbial communities

Interrogating metabolism as an electron flow system

Metabolism is generally considered as a neatly organised system of modular pathways, shaped by evolution under selection for optimal cellular growth. This view falls short of explaining and predicting a number of key observations about the structure and dynamics of metabolism. We highlight these limitations of a pathway-centric view on metabolism and summarise studies suggesting how these could be overcome by viewing metabolism as a thermodynamically and kinetically constrained, dynamical flow system. Such a systems-level, first-principles based view of metabolism can open up new avenues of metabolic engineering and cures for metabolic diseases and allow better insights to a myriad of physiological processes that are ultimately linked to metabolism. Towards further developing this view, we call for a closer interaction among physical and biological disciplines and an increased use of electrochemical and biophysical approaches to interrogate cellular metabolism together with the microenvironment in which it exists

Phosphorelays provide tunable signal processing capabilities for the cell

Achieving a complete understanding of cellular signal transduction requires deciphering the relation between structural and biochemical features of a signaling system and the shape of the signal-response relationship it embeds. Using explicit analytical expressions and numerical simulations, we present here this relation for four-layered phosphorelays, which are signaling systems that are ubiquitous in prokaryotes and also found in lower eukaryotes and plants. We derive an analytical expression that relates the shape of the signal-response relationship in a relay to the kinetic rates of forward, reverse phosphorylation and hydrolysis reactions. This reveals a set of mathematical conditions which, when satisfied, dictate the shape of the signal-response relationship. We find that a specific topology also observed in nature can satisfy these conditions in such a way to allow plasticity among hyperbolic and sigmoidal signal-response relationships. Particularly, the shape of the signal-response relationship of this relay topology can be tuned by altering kinetic rates and total protein levels at different parts of the relay. These findings provide an important step towards predicting response dynamics of phosphorelays, and the nature of subsequent physiological responses that they mediate, solely from topological features and few composite measurements; measuring the ratio of reverse and forward phosphorylation rate constants could be sufficient to determine the shape of the signal-response relationship the relay exhibits. Furthermore, they highlight the potential ways in which selective pressures on signal processing could have played a role in the evolution of the observed structural and biochemical characteristic in phosphorelays

Response dynamics of phosphorelays suggest their potential utility in cell signalling

Phosphorelays are extended two-component signalling systems found in diverse bacteria, lower eukaryotes and plants. Only few of these systems are characterized, and we still lack a full understanding of their signalling abilities. Here, we aim to achieve a global understanding of phosphorelay signalling and its dynamical properties. We develop a generic model, allowing us to systematically analyse response dynamics under different assumptions. Using this model, we find that the steady-state concentration of phosphorylated protein at the final layer of a phosphorelay is a linearly increasing, but eventually saturating function of the input. In contrast, the intermediate layers can display ultrasensitivity. We find that such ultrasensitivity is a direct result of the phosphorelay biochemistry; shuttling of a single phosphate group from the first to the last layer. The response dynamics of the phosphorelay results in tolerance of cross-talk, especially when it occurs as cross-deactivation. Further, it leads to a high signal-to-noise ratio for the final layer. We find that a relay length of four, which is most commonly observed, acts as a saturating point for these dynamic properties. These findings suggest that phosphorelays could act as a mechanism to reduce noise and effects of cross-talk on the final layer of the relay and enforce its input–response relation to be linear. In addition, our analysis suggests that middle layers of phosphorelays could embed thresholds. We discuss the consequence of these findings in relation to why cells might use phosphorelays along with enzymatic kinase cascades

MetQy – an R package to query metabolic functions of genes and genomes

Summary: With the rapid accumulation of sequencing data from genomic and metagenomic studies, there is an acute need for better tools that facilitate their analyses against biological functions. To this end, we developed MetQy, an open–source R package designed for query–based analysis of functional units in [meta]genomes and/or sets of genes using the The Kyoto Encyclopedia of Genes and Genomes (KEGG). Furthermore, MetQy contains visualisation and analysis tools and facilitates KEGG’s flat file manipulation. Thus, MetQy enables better understanding of metabolic capabilities of known genomes or user–specified [meta]genomes by using the available information and can help guide studies in microbial ecology, metabolic engineering and synthetic biology. Availability and Implementation: The MetQy R package is freely available and can be downloaded from our group’s website (http://osslab.lifesci.warwick.ac.uk) or GitHub (https://github.com/OSS-Lab/MetQy)

Manganese oxide biomineralization is a social trait protecting

Manganese bio-mineralization is a widespread process among bacteria and fungi. To date there is no conclusive experimental evidence for, how and if this process impacts microbial fitness in the environment. Here we show how a model organism for manganese oxidation is growth-inhibited by nitrite, and that this inhibition is mitigated in presence of manganese. We show that such manganese-mediated mitigation of nitrite-inhibition is dependent on the culture inoculum size and that manganese oxide (MnOX) forms granular precipitates in the culture, rather than sheaths around individual cells. We provide evidence that MnOX protection involves both its ability to catalyze nitrite oxidation into (non-toxic) nitrate under physiological conditions, and its potential role in influencing processes involving reactive oxygen species (ROS). Taken together, these results demonstrate improved microbial fitness through MnOX deposition in an ecological setting, i.e. mitigation of nitrite toxicity, and point to a key role of MnOX in handling stresses arising from ROS

BioJazz : In silico evolution of cellular networks with unbounded complexity using rule-based modeling

Systems biologists aim to decipher the structure and dynamics of signaling and regulatory networks underpinning cellular responses; synthetic biologists can use this insight to alter existing networks or engineer de novo ones. Both tasks will benefit from an understanding of which structural and dynamic features of networks can emerge from evolutionary processes, through which intermediary steps these arise, and whether they embody general design principles. As natural evolution at the level of network dynamics is difficult to study, in silico evolution of network models can provide important insights. However, current tools used for in silico evolution of network dynamics are limited to ad hoc computer simulations and models. Here we introduce BioJazz, an extendable, user-friendly tool for simulating the evolution of dynamic biochemical networks. Unlike previous tools for in silico evolution, BioJazz allows for the evolution of cellular networks with unbounded complexity by combining rule-based modeling with an encoding of networks that is akin to a genome. We show that BioJazz can be used to implement biologically realistic selective pressures and allows exploration of the space of network architectures and dynamics that implement prescribed physiological functions. BioJazz is provided as an open-source tool to facilitate its further development and use. Source code and user manuals are available at: http://oss-lab.github.io/biojazz and http://osslab.lifesci.warwick.ac.uk/BioJazz.aspx

Nonlinear dynamics in gene regulation promote robustness and evolvability of gene expression levels

Cellular phenotypes underpinned by regulatory networks need to respond to evolutionary pressures to allow adaptation, but at the same time be robust to perturbations. This creates a conflict in which mutations affecting regulatory networks must both generate variance but also be tolerated at the phenotype level. Here, we perform mathematical analyses and simulations of regulatory networks to better understand the potential trade-off between robustness and evolvability. Examining the phenotypic effects of mutations, we find an inverse correlation between robustness and evolvability that breaks only with nonlinearity in the network dynamics, through the creation of regions presenting sudden changes in phenotype with small changes in genotype. For genotypes embedding low levels of nonlinearity, robustness and evolvability correlate negatively and almost perfectly. By contrast, genotypes embedding nonlinear dynamics allow expression levels to be robust to small perturbations, while generating high diversity (evolvability) under larger perturbations. Thus, nonlinearity breaks the robustness-evolvability trade-off in gene expression levels by allowing disparate responses to different mutations. Using analytical derivations of robustness and system sensitivity, we show that these findings extend to a large class of gene regulatory network architectures and also hold for experimentally observed parameter regimes. Further, the effect of nonlinearity on the robustness-evolvability trade-off is ensured as long as key parameters of the system display specific relations irrespective of their absolute values. We find that within this parameter regime genotypes display low and noisy expression levels. Examining the phenotypic effects of mutations, we find an inverse correlation between robustness and evolvability that breaks only with nonlinearity in the network dynamics. Our results provide a possible solution to the robustness-evolvability trade-off, suggest an explanation for the ubiquity of nonlinear dynamics in gene expression networks, and generate useful guidelines for the design of synthetic gene circuits