Robustness of MEK-ERK Dynamics and Origins of Cell-to-Cell Variability in MAPK Signaling

Liebig's law of the minimum in the TGF-β/SMAD pathway

Cells use signaling pathways to sense and respond to their environments. The transforming growth factor-β (TGF-β) pathway produces context-specific responses. Here, we combined modeling and experimental analysis to study the dependence of the output of the TGF-β pathway on the abundance of signaling molecules in the pathway. We showed that the TGF-β pathway processes the variation of TGF-β receptor abundance using Liebig's law of the minimum, meaning that the output-modifying factor is the signaling protein that is most limited, to determine signaling responses across cell types and in single cells. We found that the abundance of either the type I (TGFBR1) or type II (TGFBR2) TGF-β receptor determined the responses of cancer cell lines, such that the receptor with relatively low abundance dictates the response. Furthermore, nuclear SMAD2 signaling correlated with the abundance of TGF-β receptor in single cells depending on the relative expression levels of TGFBR1 and TGFBR2. A similar control principle could govern the heterogeneity of signaling responses in other signaling pathways.

A guide to ERK dynamics, part 1: mechanisms and models

Extracellular signal-regulated kinase (ERK) has long been studied as a key driver of both essential cellular processes and disease. A persistent question has been how this single pathway is able to direct multiple cell behaviors, including growth, proliferation, and death. Modern biosensor studies have revealed that the temporal pattern of ERK activity is highly variable and heterogeneous, and critically, that these dynamic differences modulate cell fate. This two-part review discusses the current understanding of dynamic activity in the ERK pathway, how it regulates cellular decisions, and how these cell fates lead to tissue regulation and pathology. In part 1, we cover the optogenetic and live-cell imaging technologies that first revealed the dynamic nature of ERK, as well as current challenges in biosensor data analysis. We also discuss advances in mathematical models for the mechanisms of ERK dynamics, including receptor-level regulation, negative feedback, cooperativity, and paracrine signaling. While hurdles still remain, it is clear that higher temporal and spatial resolution provide mechanistic insights into pathway circuitry. Exciting new algorithms and advanced computational tools enable quantitative measurements of single-cell ERK activation, which in turn inform better models of pathway behavior. However, the fact that current models still cannot fully recapitulate the diversity of ERK responses calls for a deeper understanding of network structure and signal transduction in general.

Julia for biologists

Major computational challenges exist in relation to the collection, curation, processing and analysis of large genomic and imaging datasets, as well as the simulation of larger and more realistic models in systems biology. Here we discuss how a relative newcomer among programming languages-Julia-is poised to meet the current and emerging demands in the computational biosciences and beyond. Speed, flexibility, a thriving package ecosystem and readability are major factors that make high-performance computing and data analysis available to an unprecedented degree. We highlight how Julia's design is already enabling new ways of analyzing biological data and systems, and we provide a list of resources that can facilitate the transition into Julian computing.

The Mechanisms of BDNF Promoting the Proliferation of Porcine Follicular Granulosa Cells: Role of miR-127 and Involvement of the MAPK-ERK1/2 Pathway

As a member of the neurotrophic family, brain-derived neurotrophic factor (BDNF) provides a key link in the physiological process of mammalian ovarian follicle development, in addition to its functions in the nervous system. The emphasis of this study lay in the impact of BDNF on the proliferation of porcine follicular granulosa cells (GCs) in vitro. BDNF and tyrosine kinase B (TrkB, receptor of BDNF) were detected in porcine follicular GCs. Additionally, cell viability significantly increased during the culture of porcine GCs with BDNF (100 ng/mL) in vitro. However, BDNF knockdown in GCs decreased cell viability and S-phase cells proportion-and BDNF simultaneously regulated the expression of genes linked with cell proliferation (CCND1, p21 and Bcl2) and apoptosis (Bax). Then, the results of the receptor blocking experiment showed that BDNF promoted GC proliferation via TrkB. The high-throughput sequencing showed that BDNF also regulated the expression profiles of miRNAs in GCs. The differential expression profiles were obtained by miRNA sequencing after BDNF (100 ng/mL) treatment with GCs. The sequencing results showed that, after BDNF treatment, 72 significant differentially-expressed miRNAs were detected-5 of which were related to cell process and proliferation signaling pathways confirmed by RT-PCR. Furthermore, studies showed that BDNF promoted GCs' proliferation by increasing the expression of CCND1, downregulating miR-127 and activating the ERK1/2 signal pathway. Moreover, BDNF indirectly activated the ERK1/2 signal pathway by downregulating miR-127. In conclusion, BDNF promoted porcine GC proliferation by increasing CCND1 expression, downregulating miR-127 and stimulating the MAPK-ERK1/2 signaling cascade.

The role of mitosis in generating fitness heterogeneity

Cancer cells have heterogeneous fitness, and this heterogeneity stems from genetic and epigenetic sources. Here, we sought to assess the contribution of asymmetric mitosis (AM) and time on the variability of fitness in sister cells. Around one quarter of sisters had differences in fitness, assessed as the intermitotic time (IMT), from 330 to 510 min. Phenotypes related to fitness, such as ERK activity (herein referring to ERK1 and ERK2, also known as MAPK3 and MAPK1, respectively), DNA damage and nuclear morphological phenotypes were also asymmetric at mitosis or turned asymmetric over the course of the cell cycle. The ERK activity of mother cell was found to influence the ERK activity and the IMT of the daughter cells, and cells with ERK asymmetry at mitosis produced more offspring with AMs, suggesting heritability of the AM phenotype for ERK activity. Our findings demonstrate how variabilities in sister cells can be generated, contributing to the phenotype heterogeneities in tumor cells.

Modeling Cellular Signaling Variability Based on Single-Cell Data: The TGFβ-SMAD Signaling Pathway

Nongenetic heterogeneity is key to cellular decisions, as even genetically identical cells respond in very different ways to the same external stimulus, e.g., during cell differentiation or therapeutic treatment of disease. Strong heterogeneity is typically already observed at the level of signaling pathways that are the first sensors of external inputs and transmit information to the nucleus where decisions are made. Since heterogeneity arises from random fluctuations of cellular components, mathematical models are required to fully describe the phenomenon and to understand the dynamics of heterogeneous cell populations. Here, we review the experimental and theoretical literature on cellular signaling heterogeneity, with special focus on the TGFβ/SMAD signaling pathway.

Open problems in mathematical biology

Biology is data-rich, and it is equally rich in concepts and hypotheses. Part of trying to understand biological processes and systems is therefore to confront our ideas and hypotheses with data using statistical methods to determine the extent to which our hypotheses agree with reality. But doing so in a systematic way is becoming increasingly challenging as our hypotheses become more detailed, and our data becomes more complex. Mathematical methods are therefore gaining in importance across the life- and biomedical sciences. Mathematical models allow us to test our understanding, make testable predictions about future behaviour, and gain insights into how we can control the behaviour of biological systems. It has been argued that mathematical methods can be of great benefit to biologists to make sense of data. But mathematics and mathematicians are set to benefit equally from considering the often bewildering complexity inherent to living systems. Here we present a small selection of open problems and challenges in mathematical biology. We have chosen these open problems because they are of both biological and mathematical interest.

Towards 'end-to-end' analysis and understanding of biological timecourse data

Petabytes of increasingly complex and multidimensional live cell and tissue imaging data are generated every year. These videos hold large promise for understanding biology at a deep and fundamental level, as they capture single-cell and multicellular events occurring over time and space. However, the current modalities for analysis and mining of these data are scattered and user-specific, preventing more unified analyses from being performed over different datasets and obscuring possible scientific insights. Here, we propose a unified pipeline for storage, segmentation, analysis, and statistical parametrization of live cell imaging datasets.

HyperGraphs.jl - representing high-order relationships in Julia

Summary

HyperGraphs.jl is a Julia package that implements hypergraphs. These are a generalization of graphs that allow us to represent n-ary relationships and not just binary, pairwise relationships. High-order interactions are commonplace in biological systems and are of critical importance to their dynamics; hypergraphs thus offer a natural way to accurately describe and model these systems.

Availability and implementation

HyperGraphs.jl is freely available under the MIT license. Source code and documentation can be found at https://github.com/lpmdiaz/HyperGraphs.jl.

Supplementary information

Supplementary data are available at Bioinformatics online.

microRNA-Mediated Encoding and Decoding of Time-Dependent Signals in Tumorigenesis

microRNAs, pivotal post-transcriptional regulators of gene expression, in the past decades have caught the attention of researchers for their involvement in different biological processes, ranging from cell development to cancer. Although lots of effort has been devoted to elucidate the topological features and the equilibrium properties of microRNA-mediated motifs, little is known about how the information encoded in frequency, amplitude, duration, and other features of their regulatory signals can affect the resulting gene expression patterns. Here, we review the current knowledge about microRNA-mediated gene regulatory networks characterized by time-dependent input signals, such as pulses, transient inputs, and oscillations. First, we identify the general characteristic of the main motifs underlying temporal patterns. Then, we analyze their impact on two commonly studied oncogenic networks, showing how their dysfunction can lead to tumorigenesis.

A single-cell resolved cell-cell communication model explains lineage commitment in hematopoiesis

Cells do not make fate decisions independently. Arguably, every cell-fate decision occurs in response to environmental signals. In many cases, cell-cell communication alters the dynamics of the internal gene regulatory network of a cell to initiate cell-fate transitions, yet models rarely take this into account. Here, we have developed a multiscale perspective to study the granulocyte-monocyte versus megakaryocyte-erythrocyte fate decisions. This transition is dictated by the GATA1-PU.1 network: a classical example of a bistable cell-fate system. We show that, for a wide range of cell communication topologies, even subtle changes in signaling can have pronounced effects on cell-fate decisions. We go on to show how cell-cell coupling through signaling can spontaneously break the symmetry of a homogenous cell population. Noise, both intrinsic and extrinsic, shapes the decision landscape profoundly, and affects the transcriptional dynamics underlying this important hematopoietic cell-fate decision-making system.

This article has an associated ‘The people behind the papers’ interview.

Summary: Through theory and computational modeling, cell-cell communication is revealed to be a crucial and under-appreciated determinant of cell-fate decision-making during hematopoiesis.

Pathway dynamics can delineate the sources of transcriptional noise in gene expression

Single-cell expression profiling opens up new vistas on cellular processes. Extensive cell-to-cell variability at the transcriptomic and proteomic level has been one of the stand-out observations. Because most experimental analyses are destructive we only have access to snapshot data of cellular states. This loss of temporal information presents significant challenges for inferring dynamics, as well as causes of cell-to-cell variability. In particular, we typically cannot separate dynamic variability from within cells ('intrinsic noise') from variability across the population ('extrinsic noise'). Here, we make this non-identifiability mathematically precise, allowing us to identify new experimental set-ups that can assist in resolving this non-identifiability. We show that multiple generic reporters from the same biochemical pathways (e.g. mRNA and protein) can infer magnitudes of intrinsic and extrinsic transcriptional noise, identifying sources of heterogeneity. Stochastic simulations support our theory, and demonstrate that 'pathway-reporters' compare favourably to the well-known, but often difficult to implement, dual-reporter method.

Somatic CG6015 mediates cyst stem cell maintenance and germline stem cell differentiation via EGFR signaling in Drosophila testes

Stem cell niche is regulated by intrinsic and extrinsic factors. In the Drosophila testis, cyst stem cells (CySCs) support the differentiation of germline stem cells (GSCs). However, the underlying mechanisms remain unclear. In this study, we found that somatic CG6015 is required for CySC maintenance and GSC differentiation in a Drosophila model. Knockdown of CG6015 in CySCs caused aberrant activation of dpERK in undifferentiated germ cells in the Drosophila testis, and disruption of key downstream targets of EGFR signaling (Dsor1 and rl) in CySCs results in a phenotype resembling that of CG6015 knockdown. CG6015, Dsor1, and rl are essential for the survival of Drosophila cell line Schneider 2 (S2) cells. Our data showed that somatic CG6015 regulates CySC maintenance and GSC differentiation via EGFR signaling, and inhibits aberrant activation of germline dpERK signals. These findings indicate regulatory mechanisms of stem cell niche homeostasis in the Drosophila testis.

Phenotypic Heterogeneity in Tumor Progression, and Its Possible Role in the Onset of Cancer

Heterogeneity among isogenic cells/individuals has been known for at least 150 years. Even Mendel, working on pea plants, realized that not all tall plants were identical. However, Mendel was more interested in the discontinuous variation between genetically distinct individuals. The concept of environment dictating distinct phenotypes among isogenic individuals has since been shown to impact the evolution of populations in numerous examples at different scales of life. In this review, we discuss how phenotypic heterogeneity and its evolutionary implications exist at all levels of life, from viruses to mammals. In particular, we discuss how a particular disease condition (cancer) is impacted by heterogeneity among isogenic cells, and propose a potential role that phenotypic heterogeneity might play toward the onset of the disease.

Cell-to-cell diversification in ERBB-RAS-MAPK signal transduction that produces cell-type specific growth factor responses

Growth factors regulate cell fates, including their proliferation, differentiation, survival, and death, according to the cell type. Even when the response to a specific growth factor is deterministic for collective cell behavior, significant levels of fluctuation are often observed between single cells. Statistical analyses of single-cell responses provide insights into the mechanism of cell fate decisions but very little is known about the distributions of the internal states of cells responding to growth factors. Using multi-color immunofluorescent staining, we have here detected the phosphorylation of seven elements in the early response of the ERBB-RAS-MAPK system to two growth factors. Among these seven elements, five were analyzed simultaneously in distinct combinations in the same single cells. Although principle component analysis suggested cell-type and input specific phosphorylation patterns, cell-to-cell fluctuation was large. Mutual information analysis suggested that each cell type uses multitrack (bush-like) signal transduction pathways under conditions in which clear fate changes have been reported. The clustering of single-cell response patterns indicated that the fate change in a cell population correlates with the large entropy of the response, suggesting a bet-hedging strategy is used in decision making. A comparison of true and randomized datasets further indicated that this large variation is not produced by simple reaction noise, but is defined by the properties of the signal-processing network.

Inference of Multisite Phosphorylation Rate Constants and Their Modulation by Pathogenic Mutations

Multisite protein phosphorylation plays a critical role in cell regulation [1-3]. It is widely appreciated that the functional capabilities of multisite phosphorylation depend on the order and kinetics of phosphorylation steps, but kinetic aspects of multisite phosphorylation remain poorly understood [4-6]. Here, we focus on what appears to be the simplest scenario, when a protein is phosphorylated on only two sites in a strict, well-defined order. This scenario describes the activation of ERK, a highly conserved cell-signaling enzyme. We use Bayesian parameter inference in a structurally identifiable kinetic model to dissect dual phosphorylation of ERK by MEK, a kinase that is mutated in a large number of human diseases [7-12]. Our results reveal how enzyme processivity and efficiencies of individual phosphorylation steps are altered by pathogenic mutations. The presented approach, which connects specific mutations to kinetic parameters of multisite phosphorylation mechanisms, provides a systematic framework for closing the gap between studies with purified enzymes and their effects in the living organism.

Information-theoretic analysis of multivariate single-cell signaling responses

Mathematical methods of information theory appear to provide a useful language to describe how stimuli are encoded in activities of signaling effectors. Exploring the information-theoretic perspective, however, remains conceptually, experimentally and computationally challenging. Specifically, existing computational tools enable efficient analysis of relatively simple systems, usually with one input and output only. Moreover, their robust and readily applicable implementations are missing. Here, we propose a novel algorithm, SLEMI—statistical learning based estimation of mutual information, to analyze signaling systems with high-dimensional outputs and a large number of input values. Our approach is efficient in terms of computational time as well as sample size needed for accurate estimation. Analysis of the NF-κB single—cell signaling responses to TNF-α reveals that NF-κB signaling dynamics improves discrimination of high concentrations of TNF-α with a relatively modest impact on discrimination of low concentrations. Provided R-package allows the approach to be used by computational biologists with only elementary knowledge of information theory.

In light of single-cell, live-imaging experiments understanding of how cells transmit information about identity and quantity of stimuli is incomplete. When exposed to the same stimulus individual cells exhibit substantial cell-to-cell heterogeneity. Besides, stimuli have been shown to regulate temporal profiles of signaling effectors. Therefore, it is, for instance, not entirely clear whether single-cell responses are binary or contain more information about the quantity of stimuli. The above questions resulted in a considerable interest to study cellular signaling within the framework of information theory. Unfortunately, the utilization of the information-theoretic perspective is handicapped in part by the lack of suitable methods that account for multivariate signaling data. Here, we propose a novel algorithm that breaks a considerable computational barrier by allowing the effective information-theoretic analysis of highly-dimensional single-cell measurements. Our approach is computationally efficient, robust and straightforward to use. Moreover, we provide a simple R-package implementation.

A Simple and Flexible Computational Framework for Inferring Sources of Heterogeneity from Single-Cell Dynamics

Single-cell time-lapse data provide the means for disentangling sources of cell-to-cell and intra-cellular variability, a key step for understanding heterogeneity in cell populations. However, single-cell analysis with dynamic models is a challenging open problem: current inference methods address only single-gene expression or neglect parameter correlations. We report on a simple, flexible, and scalable method for estimating cell-specific and population-average parameters of non-linear mixed-effects models of cellular networks, demonstrating its accuracy with a published model and dataset. We also propose sensitivity analysis for identifying which biological sub-processes quantitatively and dynamically contribute to cell-to-cell variability. Our application to endocytosis in yeast demonstrates that dynamic models of realistic size can be developed for the analysis of single-cell data and that shifting the focus from single reactions or parameters to nuanced and time-dependent contributions of sub-processes helps biological interpretation. Generality and simplicity of the approach will facilitate customized extensions for analyzing single-cell dynamics.

Emergent memory in cell signaling: Persistent adaptive dynamics in cascades can arise from the diversity of relaxation time-scales

The mitogen-activated protein kinase (MAPK) signaling cascade, an evolutionarily conserved motif present in all eukaryotic cells, is involved in coordinating crucial cellular functions. While the asymptotic dynamical behavior of the pathway stimulated by a time-invariant signal is relatively well-understood, we show using a computational model that it exhibits a rich repertoire of transient adaptive responses to changes in stimuli. When the signal is switched on, the response is characterized by long-lived modulations in frequency as well as amplitude. On withdrawing the stimulus, the activity decays over long timescales, exhibiting reverberations characterized by repeated spiking in the activated MAPK concentration. The long-term persistence of such post-stimulus activity suggests that the cascade retains memory of the signal for a significant duration following its removal. The molecular mechanism underlying the reverberatory activity is related to the existence of distinct relaxation rates for the different cascade components. This results in the imbalance of fluxes between different layers of the cascade, with the reuse of activated kinases as enzymes when they are released from sequestration in complexes. The persistent adaptive response, indicative of a cellular “short-term” memory, suggests that this ubiquitous signaling pathway plays an even more central role in information processing by eukaryotic cells.

Transition state characteristics during cell differentiation

Models describing the process of stem-cell differentiation are plentiful, and may offer insights into the underlying mechanisms and experimentally observed behaviour. Waddington's epigenetic landscape has been providing a conceptual framework for differentiation processes since its inception. It also allows, however, for detailed mathematical and quantitative analyses, as the landscape can, at least in principle, be related to mathematical models of dynamical systems. Here we focus on a set of dynamical systems features that are intimately linked to cell differentiation, by considering exemplar dynamical models that capture important aspects of stem cell differentiation dynamics. These models allow us to map the paths that cells take through gene expression space as they move from one fate to another, e.g. from a stem-cell to a more specialized cell type. Our analysis highlights the role of the transition state (TS) that separates distinct cell fates, and how the nature of the TS changes as the underlying landscape changes-change that can be induced by e.g. cellular signaling. We demonstrate that models for stem cell differentiation may be interpreted in terms of either a static or transitory landscape. For the static case the TS represents a particular transcriptional profile that all cells approach during differentiation. Alternatively, the TS may refer to the commonly observed period of heterogeneity as cells undergo stochastic transitions.

In Vivo Decoding Mechanisms of the Temporal Patterns of Blood Insulin by the Insulin-AKT Pathway in the Liver

Cells respond to various extracellular stimuli through a limited number of signaling pathways. One strategy to process such stimuli is to code the information into the temporal patterns of molecules. Although we showed that insulin selectively regulated molecules depending on its temporal patterns using Fao cells, the in vivo mechanism remains unknown. Here, we show how the insulin-AKT pathway processes the information encoded into the temporal patterns of blood insulin. We performed hyperinsulinemic-euglycemic clamp experiments and found that, in the liver, all temporal patterns of insulin are encoded into the insulin receptor, and downstream molecules selectively decode them through AKT. S6K selectively decodes the additional secretion information. G6Pase interprets the basal secretion information through FoxO1, while GSK3β decodes all secretion pattern information. Mathematical modeling revealed the mechanism via differences in network structures and from sensitivity and time constants. Given that almost all hormones exhibit distinct temporal patterns, temporal coding may be a general principle of system homeostasis by hormones.

Mean-Independent Noise Control of Cell Fates via Intermediate States

Stochasticity affects accurate signal detection and robust generation of correct cell fates. Although many known regulatory mechanisms may reduce fluctuations in signals, most simultaneously influence their mean dynamics, leading to unfaithful cell fates. Through analysis and computation, we demonstrate that a reversible signaling mechanism acting through intermediate states can reduce noise while maintaining the mean. This mean-independent noise control (MINC) mechanism is investigated in the context of an intracellular binding protein that regulates retinoic acid (RA) signaling during zebrafish hindbrain development. By comparing our models with experimental data, we find that the MINC mechanism allows for sharp boundaries of gene expression without sacrificing boundary accuracy. In addition, this MINC mechanism can modulate noise to levels that we show are beneficial to spatial patterning through noise-induced cell fate switching. These results reveal a design principle that may be important for noise regulation in many systems that control cell fate determination.

•

Mean-independent noise control allows noise attenuation without affecting the mean

•

Intermediate states enable such control through proportional coupling

•

This controls spatial gene expression noise without shifting boundary locations

•

Specific noise levels are required for successful downstream boundary sharpening

Computational Methods for Estimating Molecular System from Membrane Potential Recordings in Nerve Growth Cone

Biological cells express intracellular biomolecular information to the extracellular environment as various physical responses. We show a novel computational approach to estimate intracellular biomolecular pathways from growth cone electrophysiological responses. Previously, it was shown that cGMP signaling regulates membrane potential (MP) shifts that control the growth cone turning direction during neuronal development. We present here an integrated deterministic mathematical model and Bayesian reversed-engineering framework that enables estimation of the molecular signaling pathway from electrical recordings and considers both the system uncertainty and cell-to-cell variability. Our computational method selects the most plausible molecular pathway from multiple candidates while satisfying model simplicity and considering all possible parameter ranges. The model quantitatively reproduces MP shifts depending on cGMP levels and MP variability potential in different experimental conditions. Lastly, our model predicts that chloride channel inhibition by cGMP-dependent protein kinase (PKG) is essential in the core system for regulation of the MP shifts.

Gonadotropin-releasing hormone signaling: An information theoretic approach

Gonadotropin-releasing hormone (GnRH) is a peptide hormone that mediates central control of reproduction, acting via G-protein coupled receptors that are primarily G_q coupled and mediate GnRH effects on the synthesis and secretion of luteinizing hormone and follicle-stimulating hormone. A great deal is known about the GnRH receptor signaling network but GnRH is secreted in short pulses and much less is known about how gonadotropes decode this pulsatile signal. Similarly, single cell measures reveal considerable cell-cell heterogeneity in responses to GnRH but the impact of this variability on signaling is largely unknown. Ordinary differential equation-based mathematical models have been used to explore the decoding of pulse dynamics and information theory-derived statistical measures are increasingly used to address the influence of cell-cell variability on the amount of information transferred by signaling pathways. Here, we describe both approaches for GnRH signaling, with emphasis on novel insights gained from the information theoretic approach and on the fundamental question of why GnRH is secreted in pulses.

Epigenetic regulation of cell fate reprogramming in aging and disease: A predictive computational model

Understanding the control of epigenetic regulation is key to explain and modify the aging process. Because histone-modifying enzymes are sensitive to shifts in availability of cofactors (e.g. metabolites), cellular epigenetic states may be tied to changing conditions associated with cofactor variability. The aim of this study is to analyse the relationships between cofactor fluctuations, epigenetic landscapes, and cell state transitions. Using Approximate Bayesian Computation, we generate an ensemble of epigenetic regulation (ER) systems whose heterogeneity reflects variability in cofactor pools used by histone modifiers. The heterogeneity of epigenetic metabolites, which operates as regulator of the kinetic parameters promoting/preventing histone modifications, stochastically drives phenotypic variability. The ensemble of ER configurations reveals the occurrence of distinct epi-states within the ensemble. Whereas resilient states maintain large epigenetic barriers refractory to reprogramming cellular identity, plastic states lower these barriers, and increase the sensitivity to reprogramming. Moreover, fine-tuning of cofactor levels redirects plastic epigenetic states to re-enter epigenetic resilience, and vice versa. Our ensemble model agrees with a model of metabolism-responsive loss of epigenetic resilience as a cellular aging mechanism. Our findings support the notion that cellular aging, and its reversal, might result from stochastic translation of metabolic inputs into resilient/plastic cell states via ER systems.

Modeling of cytometry data in logarithmic space: When is a bimodal distribution not bimodal?

Recent efforts in systems immunology lead researchers to build quantitative models of cell activation and differentiation. One goal is to account for the distributions of proteins from single-cell measurements by flow cytometry or mass cytometry as readout of biological regulation. In that context, large cell-to-cell variability is often observed in biological quantities. We show here that these readouts, viewed in logarithmic scale may result in two easily-distinguishable modes, while the underlying distribution (in linear scale) is unimodal. We introduce a simple mathematical test to highlight this mismatch. We then dissect the flow of influence of cell-to-cell variability proposing a graphical model which motivates higher-dimensional analysis of the data. Finally we show how acquiring additional biological information can be used to reduce uncertainty introduced by cell-to-cell variability, helping to clarify whether the data is uni- or bimodal. This communication has cautionary implications for manual and automatic gating strategies, as well as clustering and modeling of single-cell measurements. © 2018 International Society for Advancement of Cytometry.

[Mathematical modeling: an essential tool for the study of therapeutic targeting in solid tumors]

Recent progress in biology has made the study of the medical treatment of cancer more effective, but it has also revealed the large complexity of carcinogenesis and cell signaling. For many types of cancer, several therapeutic targets are known and in some cases drugs against these targets exist. Unfortunately, the target proteins often work in networks, resulting in functional adaptation and the development of resilience/resistance to medical treatment. The use of mathematical modeling makes it possible to carry out system-level analyses for improved study of therapeutic targeting in solid tumours. We present the main types of mathematical models used in cancer research and we provide examples illustrating the relevance of these approaches in molecular oncobiology.

System identification of signaling dependent gene expression with different time-scale data

Cells decode information of signaling activation at a scale of tens of minutes by downstream gene expression with a scale of hours to days, leading to cell fate decisions such as cell differentiation. However, no system identification method with such different time scales exists. Here we used compressed sensing technology and developed a system identification method using data of different time scales by recovering signals of missing time points. We measured phosphorylation of ERK and CREB, immediate early gene expression products, and mRNAs of decoder genes for neurite elongation in PC12 cell differentiation and performed system identification, revealing the input-output relationships between signaling and gene expression with sensitivity such as graded or switch-like response and with time delay and gain, representing signal transfer efficiency. We predicted and validated the identified system using pharmacological perturbation. Thus, we provide a versatile method for system identification using data with different time scales.

Quantitative phosphoproteomics reveals involvement of multiple signaling pathways in early phagocytosis by the retinal pigmented epithelium

One of the major biological functions of the retinal pigmented epithelium (RPE) is the clearance of shed photoreceptor outer segments (POS) through a multistep process resembling phagocytosis. RPE phagocytosis helps maintain the viability of photoreceptors that otherwise could succumb to the high metabolic flux and photo-oxidative stress associated with visual processing. The regulatory mechanisms underlying phagocytosis in the RPE are not fully understood, although dysfunction of this process contributes to the pathogenesis of multiple human retinal degenerative disorders, including age-related macular degeneration. Here, we present an integrated transcriptomic, proteomic, and phosphoproteomic analysis of phagocytosing RPE cells, utilizing three different experimental models: the human-derived RPE-like cell line ARPE-19, cultured murine primary RPE cells, and RPE samples from live mice. Our combined results indicated that early stages of phagocytosis in the RPE are mainly characterized by pronounced changes in the protein phosphorylation level. Global phosphoprotein enrichment analysis revealed involvement of PI3K/Akt, mechanistic target of rapamycin (mTOR), and MEK/ERK pathways in the regulation of RPE phagocytosis, confirmed by immunoblot analyses and in vitro phagocytosis assays. Most strikingly, phagocytosis of POS by cultured RPE cells was almost completely blocked by pharmacological inhibition of phosphorylation of Akt. Our findings, along with those of previous studies, indicate that these phosphorylation events allow the RPE to integrate multiple signals instigated by shed POS at different stages of the phagocytic process.

Control mechanisms for stochastic biochemical systems via computation of reachable sets

Controlling the behaviour of cells by rationally guiding molecular processes is an overarching aim of much of synthetic biology. Molecular processes, however, are notoriously noisy and frequently nonlinear. We present an approach to studying the impact of control measures on motifs of molecular interactions that addresses the problems faced in many biological systems: stochasticity, parameter uncertainty and nonlinearity. We show that our reachability analysis formalism can describe the potential behaviour of biological (naturally evolved as well as engineered) systems, and provides a set of bounds on their dynamics at the level of population statistics: for example, we can obtain the possible ranges of means and variances of mRNA and protein expression levels, even in the presence of uncertainty about model parameters.

How to deal with parameters for whole-cell modelling

Dynamical systems describing whole cells are on the verge of becoming a reality. But as models of reality, they are only useful if we have realistic parameters for the molecular reaction rates and cell physiological processes. There is currently no suitable framework to reliably estimate hundreds, let alone thousands, of reaction rate parameters. Here, we map out the relative weaknesses and promises of different approaches aimed at redressing this issue. While suitable procedures for estimation or inference of the whole (vast) set of parameters will, in all likelihood, remain elusive, some hope can be drawn from the fact that much of the cellular behaviour may be explained in terms of smaller sets of parameters. Identifying such parameter sets and assessing their behaviour is now becoming possible even for very large systems of equations, and we expect such methods to become central tools in the development and analysis of whole-cell models.

Inferring extrinsic noise from single-cell gene expression data using approximate Bayesian computation

Background

Gene expression is known to be an intrinsically stochastic process which can involve single-digit numbers of mRNA molecules in a cell at any given time. The modelling of such processes calls for the use of exact stochastic simulation methods, most notably the Gillespie algorithm. However, this stochasticity, also termed “intrinsic noise”, does not account for all the variability between genetically identical cells growing in a homogeneous environment.

Despite substantial experimental efforts, determining appropriate model parameters continues to be a challenge. Methods based on approximate Bayesian computation can be used to obtain posterior parameter distributions given the observed data. However, such inference procedures require large numbers of simulations of the model and exact stochastic simulation is computationally costly.

In this work we focus on the specific case of trying to infer model parameters describing reaction rates and extrinsic noise on the basis of measurements of molecule numbers in individual cells at a given time point.

Results

To make the problem computationally tractable we develop an exact, model-specific, stochastic simulation algorithm for the commonly used two-state model of gene expression. This algorithm relies on certain assumptions and favourable properties of the model to forgo the simulation of the whole temporal trajectory of protein numbers in the system, instead returning only the number of protein and mRNA molecules present in the system at a specified time point. The computational gain is proportional to the number of protein molecules created in the system and becomes significant for systems involving hundreds or thousands of protein molecules.

Conclusions

We employ this simulation algorithm with approximate Bayesian computation to jointly infer the model’s rate and noise parameters from published gene expression data. Our analysis indicates that for most genes the extrinsic contributions to noise will be small to moderate but certainly are non-negligible.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0324-x) contains supplementary material, which is available to authorized users.