Bridging Time Scales in Cellular Decision Making with a Stochastic Bistable Switch

Cellular transformations which involve a significant phenotypical change of the cell's state use bistable biochemical switches as underlying decision systems. In this work, we aim at linking cellular decisions taking place on a time scale of years to decades with the biochemical dynamics in signal transduction and gene regulation, occuring on a time scale of minutes to hours. We show that a stochastic bistable switch forms a viable biochemical mechanism to implement decision processes on long time scales. As a case study, the mechanism is applied to model the initiation of follicle growth in mammalian ovaries, where the physiological time scale of follicle pool depletion is on the order of the organism's lifespan. We construct a simple mathematical model for this process based on experimental evidence for the involved genetic mechanisms. Despite the underlying stochasticity, the proposed mechanism turns out to yield reliable behavior in large populations of cells subject to the considered decision process. Our model explains how the physiological time constant may emerge from the intrinsic stochasticity of the underlying gene regulatory network. Apart from ovarian follicles, the proposed mechanism may also be of relevance for other physiological systems where cells take binary decisions over a long time scale.Comment: 14 pages, 4 figure

Systems analysis of apoptosis protein expression allows the case-specific prediction of cell death responsiveness of melanoma cells.

Many cancer entities and their associated cell line models are highly heterogeneous in their responsiveness to apoptosis inducers and, despite a detailed understanding of the underlying signaling networks, cell death susceptibility currently cannot be predicted reliably from protein expression profiles. Here, we demonstrate that an integration of quantitative apoptosis protein expression data with pathway knowledge can predict the cell death responsiveness of melanoma cell lines. By a total of 612 measurements, we determined the absolute expression (nM) of 17 core apoptosis regulators in a panel of 11 melanoma cell lines, and enriched these data with systems-level information on apoptosis pathway topology. By applying multivariate statistical analysis and multi-dimensional pattern recognition algorithms, the responsiveness of individual cell lines to tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) or dacarbazine (DTIC) could be predicted with very high accuracy (91 and 82% correct predictions), and the most effective treatment option for individual cell lines could be pre-determined in silico. In contrast, cell death responsiveness was poorly predicted when not taking knowledge on protein-protein interactions into account (55 and 36% correct predictions). We also generated mathematical predictions on whether anti-apoptotic Bcl-2 family members or x-linked inhibitor of apoptosis protein (XIAP) can be targeted to enhance TRAIL responsiveness in individual cell lines. Subsequent experiments, making use of pharmacological Bcl-2/Bcl-xL inhibition or siRNA-based XIAP depletion, confirmed the accuracy of these predictions. We therefore demonstrate that cell death responsiveness to TRAIL or DTIC can be predicted reliably in a large number of melanoma cell lines when investigating expression patterns of apoptosis regulators in the context of their network-level interplay. The capacity to predict responsiveness at the cellular level may contribute to personalizing anti-cancer treatments in the future

Mechanistic Oral Absorption Modeling and Simulation for Formulation Development and Bioequivalence Evaluation: Report of an FDA Public Workshop

Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/138394/1/psp412204.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/138394/2/psp412204_am.pd

Why Are Outcomes Different for Registry Patients Enrolled Prospectively and Retrospectively? Insights from the Global Anticoagulant Registry in the FIELD-Atrial Fibrillation (GARFIELD-AF).

Background: Retrospective and prospective observational studies are designed to reflect real-world evidence on clinical practice, but can yield conflicting results. The GARFIELD-AF Registry includes both methods of enrolment and allows analysis of differences in patient characteristics and outcomes that may result. Methods and Results: Patients with atrial fibrillation (AF) and ≥1 risk factor for stroke at diagnosis of AF were recruited either retrospectively (n = 5069) or prospectively (n = 5501) from 19 countries and then followed prospectively. The retrospectively enrolled cohort comprised patients with established AF (for a least 6, and up to 24 months before enrolment), who were identified retrospectively (and baseline and partial follow-up data were collected from the emedical records) and then followed prospectively between 0-18 months (such that the total time of follow-up was 24 months; data collection Dec-2009 and Oct-2010). In the prospectively enrolled cohort, patients with newly diagnosed AF (≤6 weeks after diagnosis) were recruited between Mar-2010 and Oct-2011 and were followed for 24 months after enrolment. Differences between the cohorts were observed in clinical characteristics, including type of AF, stroke prevention strategies, and event rates. More patients in the retrospectively identified cohort received vitamin K antagonists (62.1% vs. 53.2%) and fewer received non-vitamin K oral anticoagulants (1.8% vs . 4.2%). All-cause mortality rates per 100 person-years during the prospective follow-up (starting the first study visit up to 1 year) were significantly lower in the retrospective than prospectively identified cohort (3.04 [95% CI 2.51 to 3.67] vs . 4.05 [95% CI 3.53 to 4.63]; p = 0.016). Conclusions: Interpretations of data from registries that aim to evaluate the characteristics and outcomes of patients with AF must take account of differences in registry design and the impact of recall bias and survivorship bias that is incurred with retrospective enrolment. Clinical Trial Registration: - URL: http://www.clinicaltrials.gov . Unique identifier for GARFIELD-AF (NCT01090362)

Global Analysis of Dynamical Decision-Making Models through Local Computation around the Hidden Saddle

Bistable dynamical switches are frequently encountered in mathematical modeling of biological systems because binary decisions are at the core of many cellular processes. Bistable switches present two stable steady-states, each of them corresponding to a distinct decision. In response to a transient signal, the system can flip back and forth between these two stable steady-states, switching between both decisions. Understanding which parameters and states affect this switch between stable states may shed light on the mechanisms underlying the decision-making process. Yet, answering such a question involves analyzing the global dynamical (i.e., transient) behavior of a nonlinear, possibly high dimensional model. In this paper, we show how a local analysis at a particular equilibrium point of bistable systems is highly relevant to understand the global properties of the switching system. The local analysis is performed at the saddle point, an often disregarded equilibrium point of bistable models but which is shown to be a key ruler of the decision-making process. Results are illustrated on three previously published models of biological switches: two models of apoptosis, the programmed cell death and one model of long-term potentiation, a phenomenon underlying synaptic plasticity

Exact model reduction of combinatorial reaction networks

Receptors and scaffold proteins usually possess a high number of distinct binding domains inducing the formation of large multiprotein signaling complexes. Due to combinatorial reasons the number of distinguishable species grows exponentially with the number of binding domains and can easily reach several millions. Even by including only a limited number of components and binding domains the resulting models are very large and hardly manageable. A novel model reduction technique allows the significant reduction and modularization of these models

Caspase-8 activity has an essential role in CD95/Fas-mediated MAPK activation

Stimulation of CD95/Fas/APO-1 results in the induction of both apoptotic and non-apoptotic signaling pathways. The processes regulating these two opposing pathways have not been thoroughly elucidated to date. In this study, using quantitative immunoblots, imaging, and mathematical modeling, we addressed the dynamics of the DED proteins of the death-inducing signaling complex (DISC), procaspase-8, and cellular FLICE inhibitory proteins (c-FLIPs) to the onset of CD95-mediated ERK1/2 and p38 mitogen-activated protein kinase (MAPK) activation. We found that CD95 DISC-induced caspase-8 activity is important for the initiation of ERK1/2 and p38 MAPK activation. The long c-FLIP isoform, c-FLIPL, and the short c-FLIP isoform, c-FLIPR, inhibited MAPK induction by blocking caspase-8 processing at the DISC. Furthermore, we built a mathematical model describing CD95 DISC-mediated MAPK activation and apoptosis. The model quantitatively defined the dynamics of DED proteins, procaspase-8, and c-FLIP, which lead to caspase-8 activation and induction of apoptotic and non-apoptotic signaling pathways. In conclusion, the combination of biochemical analysis with mathematical modeling provides evidence for an important role of caspase-8 in CD95-mediated activation of MAPKs, while c-FLIP exerts a regulatory function in this process

A comparative study of qualitative and quantitative dynamic models of biological regulatory networks

Backgroun

Exploring the Contextual Sensitivity of Factors that Determine Cell-to-Cell Variability in Receptor-Mediated Apoptosis

Stochastic fluctuations in gene expression give rise to cell-to-cell variability in protein levels which can potentially cause variability in cellular phenotype. For TRAIL (TNF-related apoptosis-inducing ligand) variability manifests itself as dramatic differences in the time between ligand exposure and the sudden activation of the effector caspases that kill cells. However, the contribution of individual proteins to phenotypic variability has not been explored in detail. In this paper we use feature-based sensitivity analysis as a means to estimate the impact of variation in key apoptosis regulators on variability in the dynamics of cell death. We use Monte Carlo sampling from measured protein concentration distributions in combination with a previously validated ordinary differential equation model of apoptosis to simulate the dynamics of receptor-mediated apoptosis. We find that variation in the concentrations of some proteins matters much more than variation in others and that precisely which proteins matter depends both on the concentrations of other proteins and on whether correlations in protein levels are taken into account. A prediction from simulation that we confirm experimentally is that variability in fate is sensitive to even small increases in the levels of Bcl-2. We also show that sensitivity to Bcl-2 levels is itself sensitive to the levels of interacting proteins. The contextual dependency is implicit in the mathematical formulation of sensitivity, but our data show that it is also important for biologically relevant parameter values. Our work provides a conceptual and practical means to study and understand the impact of cell-to-cell variability in protein expression levels on cell fate using deterministic models and sampling from parameter distributions