Cell‐specific responses to the cytokine TGFβ are determined by variability in protein levels

The cytokine TGFβ provides important information during embryonic development, adult tissue homeostasis, and regeneration. Alterations in the cellular response to TGFβ are involved in severe human diseases. To understand how cells encode the extracellular input and transmit its information to elicit appropriate responses, we acquired quantitative time‐resolved measurements of pathway activation at the single‐cell level. We established dynamic time warping to quantitatively compare signaling dynamics of thousands of individual cells and described heterogeneous single‐cell responses by mathematical modeling. Our combined experimental and theoretical study revealed that the response to a given dose of TGFβ is determined cell specifically by the levels of defined signaling proteins. This heterogeneity in signaling protein expression leads to decomposition of cells into classes with qualitatively distinct signaling dynamics and phenotypic outcome. Negative feedback regulators promote heterogeneous signaling, as a SMAD7 knock‐out specifically affected the signal duration in a subpopulation of cells. Taken together, we propose a quantitative framework that allows predicting and testing sources of cellular signaling heterogeneity.     

The “mechanostat” principle in cell differentiation. The osteochondroprogenitor paradigm

The “mechanostat” principle may be depicted as an oscillating signal of a signaling molecule, in which the amplitude, frequency, cumulative level, delay, and duration of the curve encode the information for concrete cellular responses and biological activities. When the oscillating signal is kept sustained (present delay), cell exit may be performed, whereas when the oscillating signal remains robust, cell proliferation may take place. B-catenin–Wnt signaling pathway has a key role in the differentiation of osteochondroprogenitor cells. Sustained downregulation of the β-catenin–Wnt pathway forces osteochondroprogenitors to a chondrogenic fate instead of an osteoblastic one. Other signaling, for example, bone morphogenetic protein and Notch signaling pathways interact with the Wnt pathway. The crosstalk between biochemical and mechanical stimuli produces the final information that leads to the final cell fate decisions, through the “mechanostat” principle.     

Modeling reveals that dynamic regulation of c-FLIP levels determines cell-to-cell distribution of CD95-mediated apoptosis

The expression levels of caspase-8 inhibitory c-FLIP proteins play an important role in regulating death receptor-mediated apoptosis, as their concentration at the moment when the death-inducing signaling complex (DISC) is formed determines the outcome of the DISC signal. Experimental studies have shown that c-FLIP proteins are subject to dynamic turnover and that their stability and expression levels can be rapidly altered. Even though the influence of c-FLIP on the apoptotic behavior of a single cell has been captured in mathematical simulation studies, the effect of c-FLIP turnover and stability has not been investigated. In this study, a mathematical model of apoptosis was developed to analyze how the dynamic turnover and stability of the c-FLIP isoforms regulate apoptotic signaling for both individual cells and cell populations. Intercellular parameter and concentration distributions were used to describe the behavior of cell populations. Monte-Carlo simulations of cell populations showed that c-FLIP turnover is a key determinant of death receptor responses. The fact that the developed model simulates the state of whole cell populations makes it possible to validate it by comparison with empirical data. The proposed modeling approach can be used to further determine limiting factors in the DISC signaling process.     

Brat promotes stem cell differentiation via control of a bistable switch that restricts BMP signaling

Drosophila ovarian germline stem cells (GSCs) are maintained by Dpp signaling and the Pumilio (Pum) and Nanos (Nos) translational repressors. Upon division, Dpp signaling is extinguished, and Nos is downregulated in one daughter cell, causing it to switch to a differentiating cystoblast (CB). However, downstream effectors of Pum-Nos remain unknown, and how CBs lose their responsiveness to Dpp is unclear. Here, we identify Brain Tumor (Brat) as a potent differentiation factor and target of Pum-Nos regulation. Brat is excluded from GSCs by Pum-Nos but functions with Pum in CBs to translationally repress distinct targets, including the Mad and dMyc mRNAs. Regulation of both targets simultaneously lowers cellular responsiveness to Dpp signaling, forcing the cell to become refractory to the self-renewal signal. Mathematical modeling elucidates bistability of cell fate in the Brat-mediated system, revealing how autoregulation of GSC number can arise from Brat coupling extracellular Dpp regulation to intracellular interpretation.     

Distant activation of Notch signaling induces stem cell niche assembly

Here we show that multiple modes of Notch signaling activation specify the complexity of spatial cellular interactions necessary for stem cell niche assembly. In particular, we studied the formation of the germline stem cell niche in Drosophila ovaries, which is a two-step process whereby terminal filaments are formed first. Then, terminal filaments signal to the adjacent cap cell precursors, resulting in Notch signaling activation, which is necessary for the lifelong acquisition of stem cell niche cell fate. The genetic data suggest that in order to initiate the process of stem cell niche assembly, Notch signaling is activated among non-equipotent cells via distant induction, where germline Delta is delivered to somatic cells located several diameters away via cellular projections generated by primordial germ cells. At the same time, to ensure the robustness of niche formation, terminal filament cell fate can also be induced by somatic Delta via cis- or trans-inhibition. This exemplifies a double security mechanism that guarantees that the germline stem cell niche is formed, since it is indispensable for the adjacent germline precursor cells to acquire and maintain stemness necessary for successful reproduction. These findings contribute to our understanding of the formation of stem cell niches in their natural environment, which is important for stem cell biology and regenerative medicine.     

Geometric control of myogenic cell fate

This work combines expertise in stem cell biology and bioengineering to define the system for geometric control of proliferation and differentiation of myogenic progenitor cells. We have created an artificial niche of myogenic progenitor cells, namely, modified extracellular matrix (ECM) substrates with spatially embedded growth or differentiation factors (GF, DF) that predictably direct muscle cell fate in a geometric pattern. Embedded GF and DF signal progenitor cells from specifically defined areas on the ECM successfully competed against culture media for myogenic cell fate determination at a clearly defined boundary. Differentiation of myoblasts into myotubes is induced in growth-promoting medium, myotube formation is delayed in differentiation-promoting medium, and myogenic cells, at different stages of proliferation and differentiation, can be induced to coexist adjacently in identical culture media. This method can be used to identify molecular interactions between cells in different stages of myogenic differentiation, which are likely to be important determinants of tissue repair. The designed ECM niches can be further developed into a vehicle for transplantation of myogenic progenitor cells maintaining their regenerative potential. Additionally, this work may also serve as a general model to engineer synthetic cellular niches to harness the regenerative potential of organ stem cells.     

Some mechanisms of lymphoid cell differentiation

Lymphoid cells differentiation is a well coordinated and highly integrated process, which is defined by the interaction of genes and cell signaling pathways of the maturing cell with the microenvironment factors. In the course of such interaction in the maturing cell, new ways of signal transduction appear, and on the basis of such new ways new gene networks are formed, which start to play rapidly, even during the course of the present or following cell specialization stage, the key role in the further cell fate.     

Theoretical aspects and modelling of cellular decision making,cell killing and information-processing in photodynamic therapy of cancer

Background

The aim of this report is to provide a mathematical model of the mechanism for making binary fate decisions about cell death or survival, during and after Photodynamic Therapy (PDT) treatment, and to supply the logical design for this decision mechanism as an application of rate distortion theory to the biochemical processing of information by the physical system of a cell.

Methods

Based on system biology models of the molecular interactions involved in the PDT processes previously established, and regarding a cellular decision-making system as a noisy communication channel, we use rate distortion theory to design a time dependent Blahut-Arimoto algorithm where the input is a stimulus vector composed of the time dependent concentrations of three PDT related cell death signaling molecules and the output is a cell fate decision. The molecular concentrations are determined by a group of rate equations. The basic steps are: initialize the probability of the cell fate decision, compute the conditional probability distribution that minimizes the mutual information between input and output, compute the cell probability of cell fate decision that minimizes the mutual information and repeat the last two steps until the probabilities converge. Advance to the next discrete time point and repeat the process.

Results

Based on the model from communication theory described in this work, and assuming that the activation of the death signal processing occurs when any of the molecular stimulants increases higher than a predefined threshold (50% of the maximum concentrations), for 1800s of treatment, the cell undergoes necrosis within the first 30 minutes with probability range 90.0%-99.99% and in the case of repair/survival, it goes through apoptosis within 3-4 hours with probability range 90.00%-99.00%. Although, there is no experimental validation of the model at this moment, it reproduces some patterns of survival ratios of predicted experimental data.

Conclusions

Analytical modeling based on cell death signaling molecules has been shown to be an independent and useful tool for prediction of cell surviving response to PDT. The model can be adjusted to provide important insights for cellular response to other treatments such as hyperthermia, and diseases such as neurodegeneration.

     

Mechanotransduction from the ECM to the genome: Are the pieces now in place?

A multitude of biochemical signaling processes have been characterized that affect gene expression and cellular activity. However, living cells often need to integrate biochemical signals with mechanical information from their microenvironment as they respond. In fact, the signals received by shape alone can dictate cell fate. This mechanotrasduction of information is powerful, eliciting proliferation, differentiation, or apoptosis in a manner dependent upon the extent of physical deformation. The cells internal "prestressed" structure and its "hardwired" interaction with the extra-cellular matrix (ECM) appear to confer this ability to filter biochemical signals and decide between divergent cell functions influenced by the nature of signals from the mechanical environment. In some instances mechanical signaling through the tissue microenvironment has been shown to be dominant over genomic defects, imparting a normal phenotype on cells that otherwise have transforming genetic lesions. This mechanical control of phenotype is postulated to have a central role in embryogenesis, tissue physiology as well as the pathology of a wide variety of diseases, including cancer. We will briefly review studies showing physical continuity between the external cellular microenvironment and the interior of the cell nucleus. Newly characterized structures, termed nuclear envelope lamina spanning complexes (NELSC), and their interactions will be described as part of a model for mechanical transduction of extracellular cues from the ECM to the genome.     

Dickkopf1 regulates fate decision and drives breast cancer stem cells to differentiation: an experimentally supported mathematical model

Background

Modulation of cellular signaling pathways can change the replication/differentiation balance in cancer stem cells (CSCs), thus affecting tumor growth and recurrence. Analysis of a simple, experimentally verified, mathematical model suggests that this balance is maintained by quorum sensing (QS).

Methodology/Principal Findings

To explore the mechanism by which putative QS cellular signals in mammary stem cells (SCs) may regulate SC fate decisions, we developed a multi-scale mathematical model, integrating extra-cellular and intra-cellular signal transduction within the mammary tissue dynamics. Preliminary model analysis of the single cell dynamics indicated that Dickkopf1 (Dkk1), a protein known to negatively regulate the Wnt pathway, can serve as anti-proliferation and pro-maturation signal to the cell. Simulations of the multi-scale tissue model suggested that Dkk1 may be a QS factor, regulating SC density on the level of the whole tissue: relatively low levels of exogenously applied Dkk1 have little effect on SC numbers, whereas high levels drive SCs into differentiation. To verify these model predictions, we treated the MCF-7 cell line and primary breast cancer (BC) cells from 3 patient samples with different concentrations and dosing regimens of Dkk1, and evaluated subsequent formation of mammospheres (MS) and the mammary SC marker CD44⁺CD24^lo. As predicted by the model, low concentrations of Dkk1 had no effect on primary BC cells, or even increased MS formation among MCF-7 cells, whereas high Dkk1 concentrations decreased MS formation among both primary BC cells and MCF-7 cells.

Conclusions/Significance

Our study suggests that Dkk1 treatment may be more robust than other methods for eliminating CSCs, as it challenges a general cellular homeostasis mechanism, namely, fate decision by QS. The study also suggests that low dose Dkk1 administration may be counterproductive; we showed experimentally that in some cases it can stimulate CSC proliferation, although this needs validating in vivo.     

Predictive modeling of signaling crosstalk during C. elegans vulval development

Caenorhabditis elegans vulval development provides an important paradigm for studying the process of cell fate determination and pattern formation during animal development. Although many genes controlling vulval cell fate specification have been identified, how they orchestrate themselves to generate a robust and invariant pattern of cell fates is not yet completely understood. Here, we have developed a dynamic computational model incorporating the current mechanistic understanding of gene interactions during this patterning process. A key feature of our model is the inclusion of multiple modes of crosstalk between the epidermal growth factor receptor (EGFR) and LIN-12/Notch signaling pathways, which together determine the fates of the six vulval precursor cells (VPCs). Computational analysis, using the model-checking technique, provides new biological insights into the regulatory network governing VPC fate specification and predicts novel negative feedback loops. In addition, our analysis shows that most mutations affecting vulval development lead to stable fate patterns in spite of variations in synchronicity between VPCs. Computational searches for the basis of this robustness show that a sequential activation of the EGFR-mediated inductive signaling and LIN-12 / Notch-mediated lateral signaling pathways is key to achieve a stable cell fate pattern. We demonstrate experimentally a time-delay between the activation of the inductive and lateral signaling pathways in wild-type animals and the loss of sequential signaling in mutants showing unstable fate patterns; thus, validating two key predictions provided by our modeling work. The insights gained by our modeling study further substantiate the usefulness of executing and analyzing mechanistic models to investigate complex biological behaviors.     

Dose-to-duration encoding and signaling beyond saturation in intracellular signaling networks

The cellular response elicited by an environmental cue typically varies with the strength of the stimulus. For example, in the yeast Saccharomyces cerevisiae, the concentration of mating pheromone determines whether cells undergo vegetative growth, chemotropic growth, or mating. This implies that the signaling pathways responsible for detecting the stimulus and initiating a response must transmit quantitative information about the intensity of the signal. Our previous experimental results suggest that yeast encode pheromone concentration as the duration of the transmitted signal. Here we use mathematical modeling to analyze possible biochemical mechanisms for performing this “dose-to-duration” conversion. We demonstrate that modulation of signal duration increases the range of stimulus concentrations for which dose-dependent responses are possible; this increased dynamic range produces the counterintuitive result of “signaling beyond saturation” in which dose-dependent responses are still possible after apparent saturation of the receptors. We propose a mechanism for dose-to-duration encoding in the yeast pheromone pathway that is consistent with current experimental observations. Most previous investigations of information processing by signaling pathways have focused on amplitude encoding without considering temporal aspects of signal transduction. Here we demonstrate that dose-to-duration encoding provides cells with an alternative mechanism for processing and transmitting quantitative information about their surrounding environment. The ability of signaling pathways to convert stimulus strength into signal duration results directly from the nonlinear nature of these systems and emphasizes the importance of considering the dynamic properties of signaling pathways when characterizing their behavior. Understanding how signaling pathways encode and transmit quantitative information about the external environment will not only deepen our understanding of these systems but also provide insight into how to reestablish proper function of pathways that have become dysregulated by disease.     

Coupled positive feedbacks provoke slow induction plus fast switching in apoptosis

Apoptosis is a form of a programmed cell death for multicellular organisms to remove unwanted or damaged cells. This critical choice of cellular fate is an all-or-none process, but its dynamics remains unraveled. The switch-like apoptotic decision has to be reliable, and once a pro-apoptotic fate is determined it requires fast and irreversible execution. One of the key regulators in apoptosis is caspase-3. Interestingly, activated caspase-3 quickly executes apoptosis, but it takes considerable time to activate it. Here, we have analyzed this "slow induction plus fast switching" mechanism of caspase-3 through mathematical modeling and computational simulation. First, we have shown that two positive feedbacks, composed of caspase-8 and XIAP, are essential for the "slow induction plus fast switching" behavior of caspase-3. Second, we have found that XIAP in the feedback loops primarily regulates induction time of caspase-3. In many cancer cells activation of caspase-3 is suppressed. Our results suggest that reinforcement of the positive feedback by XIAP, which relieves XIAP-mediated caspase-3 inhibition, might favor a pro-apoptotic cellular fate.