Unraveling Adaptation in Eukaryotic Pathways: Lessons from Protocells

Eukaryotic adaptation pathways operate within wide-ranging environmental conditions without stimulus saturation. Despite numerous differences in the adaptation mechanisms employed by bacteria and eukaryotes, all require energy consumption. Here, we present two minimal models showing that expenditure of energy by the cell is not essential for adaptation. Both models share important features with large eukaryotic cells: they employ small diffusible molecules and involve receptor subunits resembling highly conserved G-protein cascades. Analyzing the drawbacks of these models helps us understand the benefits of energy consumption, in terms of adjustability of response and adaptation times as well as separation of cell-external sensing and cell-internal signaling. Our work thus sheds new light on the evolution of adaptation mechanisms in complex systems.Comment: accepted for publication in PLoS Computational Biology; 19 pages, 8 figure

Narrative-based computational modelling of the Gp130/JAK/STAT signalling pathway.

BACKGROUND: Appropriately formulated quantitative computational models can support researchers in understanding the dynamic behaviour of biological pathways and support hypothesis formulation and selection by "in silico" experimentation. An obstacle to widespread adoption of this approach is the requirement to formulate a biological pathway as machine executable computer code. We have recently proposed a novel, biologically intuitive, narrative-style modelling language for biologists to formulate the pathway which is then automatically translated into an executable format and is, thus, usable for analysis via existing simulation techniques. RESULTS: Here we use a high-level narrative language in designing a computational model of the gp130/JAK/STAT signalling pathway and show that the model reproduces the dynamic behaviour of the pathway derived by biological observation. We then "experiment" on the model by simulation and sensitivity analysis to define those parameters which dominate the dynamic behaviour of the pathway. The model predicts that nuclear compartmentalisation and phosphorylation status of STAT are key determinants of the pathway and that alternative mechanisms of signal attenuation exert their influence on different timescales. CONCLUSION: The described narrative model of the gp130/JAK/STAT pathway represents an interesting case study showing how, by using this approach, researchers can model biological systems without explicitly dealing with formal notations and mathematical expressions (typically used for biochemical modelling), nevertheless being able to obtain simulation and analysis results. We present the model and the sensitivity analysis results we have obtained, that allow us to identify the parameters which are most sensitive to perturbations. The results, which are shown to be in agreement with existing mathematical models of the gp130/JAK/STAT pathway, serve us as a form of validation of the model and of the approach itself

Bioengineering models of cell signaling

Strategies for rationally manipulating cell behavior in cell-based technologies and molecular therapeutics and understanding effects of environmental agents on physiological systems may be derived from a mechanistic understanding of underlying signaling mechanisms that regulate cell functions. Three crucial attributes of signal transduction necessitate modeling approaches for analyzing these systems: an ever-expanding plethora of signaling molecules and interactions, a highly interconnected biochemical scheme, and concurrent biophysical regulation. Because signal flow is tightly regulated with positive and negative feedbacks and is bidirectional with commands traveling both from outside-in and inside-out, dynamic models that couple biophysical and biochemical elements are required to consider information processing both during transient and steady-state conditions. Unique mathematical frameworks will be needed to obtain an integrated perspective on these complex systems, which operate over wide length and time scales. These may involve a two-level hierarchical approach wherein the overall signaling network is modeled in terms of effective "circuit" or "algorithm" modules, and then each module is correspondingly modeled with more detailed incorporation of its actual underlying biochemical/biophysical molecular interactions

Adaptive response and enlargement of dynamic range

Many membrane channels and receptors exhibit adaptive, or desensitized, response to a strong sustained input stimulus, often supported by protein activity-dependent inactivation. Adaptive response is thought to be related to various cellular functions such as homeostasis and enlargement of dynamic range by background compensation. Here we study the quantitative relation between adaptive response and background compensation within a modeling framework. We show that any particular type of adaptive response is neither sufficient nor necessary for adaptive enlargement of dynamic range. In particular a precise adaptive response, where system activity is maintained at a constant level at steady state, does not ensure a large dynamic range neither in input signal nor in system output. A general mechanism for input dynamic range enlargement can come about from the activity-dependent modulation of protein responsiveness by multiple biochemical modification, regardless of the type of adaptive response it induces. Therefore hierarchical biochemical processes such as methylation and phosphorylation are natural candidates to induce this property in signaling systems.Comment: Corrected typos, minor text revision

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

A Characterization of Scale Invariant Responses in Enzymatic Networks

An ubiquitous property of biological sensory systems is adaptation: a step increase in stimulus triggers an initial change in a biochemical or physiological response, followed by a more gradual relaxation toward a basal, pre-stimulus level. Adaptation helps maintain essential variables within acceptable bounds and allows organisms to readjust themselves to an optimum and non-saturating sensitivity range when faced with a prolonged change in their environment. Recently, it was shown theoretically and experimentally that many adapting systems, both at the organism and single-cell level, enjoy a remarkable additional feature: scale invariance, meaning that the initial, transient behavior remains (approximately) the same even when the background signal level is scaled. In this work, we set out to investigate under what conditions a broadly used model of biochemical enzymatic networks will exhibit scale-invariant behavior. An exhaustive computational study led us to discover a new property of surprising simplicity and generality, uniform linearizations with fast output (ULFO), whose validity we show is both necessary and sufficient for scale invariance of enzymatic networks. Based on this study, we go on to develop a mathematical explanation of how ULFO results in scale invariance. Our work provides a surprisingly consistent, simple, and general framework for understanding this phenomenon, and results in concrete experimental predictions