Positional Information Generated by Spatially Distributed Signaling Cascades

The temporal and stationary behavior of protein modification cascades has been extensively studied, yet little is known about the spatial aspects of signal propagation. We have previously shown that the spatial separation of opposing enzymes, such as a kinase and a phosphatase, creates signaling activity gradients. Here we show under what conditions signals stall in the space or robustly propagate through spatially distributed signaling cascades. Robust signal propagation results in activity gradients with long plateaus, which abruptly decay at successive spatial locations. We derive an approximate analytical solution that relates the maximal amplitude and propagation length of each activation profile with the cascade level, protein diffusivity, and the ratio of the opposing enzyme activities. The control of the spatial signal propagation appears to be very different from the control of transient temporal responses for spatially homogenous cascades. For spatially distributed cascades where activating and deactivating enzymes operate far from saturation, the ratio of the opposing enzyme activities is shown to be a key parameter controlling signal propagation. The signaling gradients characteristic for robust signal propagation exemplify a pattern formation mechanism that generates precise spatial guidance for multiple cellular processes and conveys information about the cell size to the nucleus

Signaling Cascades Modulate the Speed of Signal Propagation through Space

Cells are not mixed bags of signaling molecules. As a consequence, signals must travel from their origin to distal locations. Much is understood about the purely diffusive propagation of signals through space. Many signals, however, propagate via signaling cascades. Here, we show that, depending on their kinetics, cascades speed up or slow down the propagation of signals through space, relative to pure diffusion.We modeled simple cascades operating under different limits of Michaelis-Menten kinetics using deterministic reaction-diffusion equations. Cascades operating far from enzyme saturation speed up signal propagation; the second mobile species moves more quickly than the first through space, on average. The enhanced speed is due to more efficient serial activation of a downstream signaling module (by the signaling molecule immediately upstream in the cascade) at points distal from the signaling origin, compared to locations closer to the source. Conversely, cascades operating under saturated kinetics, which exhibit zero-order ultrasensitivity, can slow down signals, ultimately localizing them to regions around the origin.Signal speed modulation may be a fundamental function of cascades, affecting the ability of signals to penetrate within a cell, to cross-react with other signals, and to activate distant targets. In particular, enhanced speeds provide a way to increase signal penetration into a cell without needing to flood the cell with large numbers of active signaling molecules; conversely, diminished speeds in zero-order ultrasensitive cascades facilitate strong, but localized, signaling

A Hidden Feedback in Signaling Cascades Is Revealed

Cycles involving covalent modification of proteins are key components of the intracellular signaling machinery. Each cycle is comprised of two interconvertable forms of a particular protein. A classic signaling pathway is structured by a chain or cascade of basic cycle units in such a way that the activated protein in one cycle promotes the activation of the next protein in the chain, and so on. Starting from a mechanistic kinetic description and using a careful perturbation analysis, we have derived, to our knowledge for the first time, a consistent approximation of the chain with one variable per cycle. The model we derive is distinct from the one that has been in use in the literature for several years, which is a phenomenological extension of the Goldbeter-Koshland biochemical switch. Even though much has been done regarding the mathematical modeling of these systems, our contribution fills a gap between existing models and, in doing so, we have unveiled critical new properties of this type of signaling cascades. A key feature of our new model is that a negative feedback emerges naturally, exerted between each cycle and its predecessor. Due to this negative feedback, the system displays damped temporal oscillations under constant stimulation and, most important, propagates perturbations both forwards and backwards. This last attribute challenges the widespread notion of unidirectionality in signaling cascades. Concrete examples of applications to MAPK cascades are discussed. All these properties are shared by the complete mechanistic description and our simplified model, but not by previously derived phenomenological models of signaling cascades

Pathwise Sensitivity Analysis in Transient Regimes

The instantaneous relative entropy (IRE) and the corresponding instanta- neous Fisher information matrix (IFIM) for transient stochastic processes are pre- sented in this paper. These novel tools for sensitivity analysis of stochastic models serve as an extension of the well known relative entropy rate (RER) and the corre- sponding Fisher information matrix (FIM) that apply to stationary processes. Three cases are studied here, discrete-time Markov chains, continuous-time Markov chains and stochastic differential equations. A biological reaction network is presented as a demonstration numerical example

Coulomb Interactions between Cytoplasmic Electric Fields and Phosphorylated Messenger Proteins Optimize Information Flow in Cells

Normal cell function requires timely and accurate transmission of information from receptors on the cell membrane (CM) to the nucleus. Movement of messenger proteins in the cytoplasm is thought to be dependent on random walk. However, Brownian motion will disperse messenger proteins throughout the cytosol resulting in slow and highly variable transit times. We propose that a critical component of information transfer is an intracellular electric field generated by distribution of charge on the nuclear membrane (NM). While the latter has been demonstrated experimentally for decades, the role of the consequent electric field has been assumed to be minimal due to a Debye length of about 1 nanometer that results from screening by intracellular Cl- and K+. We propose inclusion of these inorganic ions in the Debye-Huckel equation is incorrect because nuclear pores allow transit through the membrane at a rate far faster than the time to thermodynamic equilibrium. In our model, only the charged, mobile messenger proteins contribute to the Debye length.Using this revised model and published data, we estimate the NM possesses a Debye-Huckel length of a few microns and find this is consistent with recent measurement using intracellular nano-voltmeters. We demonstrate the field will accelerate isolated messenger proteins toward the nucleus through Coulomb interactions with negative charges added by phosphorylation. We calculate transit times as short as 0.01 sec. When large numbers of phosphorylated messenger proteins are generated by increasing concentrations of extracellular ligands, we demonstrate they generate a self-screening environment that regionally attenuates the cytoplasmic field, slowing movement but permitting greater cross talk among pathways. Preliminary experimental results with phosphorylated RAF are consistent with model predictions.This work demonstrates that previously unrecognized Coulomb interactions between phosphorylated messenger proteins and intracellular electric fields will optimize information transfer from the CM to the NM in cells

Modular cell biology: retroactivity and insulation

Modularity plays a fundamental role in the prediction of the behavior of a system from the behavior of its components, guaranteeing that the properties of individual components do not change upon interconnection. Just as electrical, hydraulic, and other physical systems often do not display modularity, nor do many biochemical systems, and specifically, genetic networks. Here, we study the effect of interconnections on the input–output dynamic characteristics of transcriptional components, focusing on a property, which we call ‘retroactivity', that plays a role analogous to non-zero output impedance in electrical systems. In transcriptional networks, retroactivity is large when the amount of transcription factor is comparable to, or smaller than, the amount of promoter-binding sites, or when the affinity of such binding sites is high. To attenuate the effect of retroactivity, we propose a feedback mechanism inspired by the design of amplifiers in electronics. We introduce, in particular, a mechanism based on a phosphorylation–dephosphorylation cycle. This mechanism enjoys a remarkable insulation property, due to the fast timescales of the phosphorylation and dephosphorylation reactions

Computational and Mathematical Modelling of the EGF Receptor System

This chapter gives an overview of computational and mathematical modelling of the EGF receptor system. It begins with a survey of motivations for producing such models, then describes the main approaches that are taken to carrying out such modelling, viz. differential equations and individual-based modelling. Finally, a number of projects that applying modelling and simulation techniques to various aspects of the EGF receptor system are described

A Pycellerator Tutorial.

We present a tutorial on using Pycellerator for biomolecular simulations. Models are described in human readable (and editable) text files (UTF8 or ASCII) containing collections of reactions, assignments, initial conditions, function definitions, and rate constants. These models are then converted into a Python program that can optionally solve the system, e.g., as a system of differential equations using ODEINT, or be run by another program. The input language implements an extended version of the Cellerator arrow notation, including mass action, Hill functions, S-Systems, MWC, and reactions with user-defined kinetic laws. Simple flux balance analysis is also implemented. We will demonstrate the implementation and analysis of progressively more complex models, starting from simple mass action through indexed cascades. Pycellerator can be used as a library that is integrated into other programs, run as a command line program, or in iPython notebooks. It is implemented in Python 2.7 and available under an open source GPL license

Unraveling the Design Principle for Motif Organization in Signaling Networks

Cellular signaling networks display complex architecture. Defining the design principle of this architecture is crucial for our understanding of various biological processes. Using a mathematical model for three-node feed-forward loops, we identify that the organization of motifs in specific manner within the network serves as an important regulator of signal processing. Further, incorporating a systemic stochastic perturbation to the model we could propose a possible design principle, for higher-order organization of motifs into larger networks in order to achieve specific biological output. The design principle was then verified in a large, complex human cancer signaling network. Further analysis permitted us to classify signaling nodes of the network into robust and vulnerable nodes as a result of higher order motif organization. We show that distribution of these nodes within the network at strategic locations then provides for the range of features displayed by the signaling network