A Minimal Model of Signaling Network Elucidates Cell-to-Cell Stochastic Variability in Apoptosis

Signaling networks are designed to sense an environmental stimulus and adapt to it. We propose and study a minimal model of signaling network that can sense and respond to external stimuli of varying strength in an adaptive manner. The structure of this minimal network is derived based on some simple assumptions on its differential response to external stimuli. We employ stochastic differential equations and probability distributions obtained from stochastic simulations to characterize differential signaling response in our minimal network model. We show that the proposed minimal signaling network displays two distinct types of response as the strength of the stimulus is decreased. The signaling network has a deterministic part that undergoes rapid activation by a strong stimulus in which case cell-to-cell fluctuations can be ignored. As the strength of the stimulus decreases, the stochastic part of the network begins dominating the signaling response where slow activation is observed with characteristic large cell-to-cell stochastic variability. Interestingly, this proposed stochastic signaling network can capture some of the essential signaling behaviors of a complex apoptotic cell death signaling network that has been studied through experiments and large-scale computer simulations. Thus we claim that the proposed signaling network is an appropriate minimal model of apoptosis signaling. Elucidating the fundamental design principles of complex cellular signaling pathways such as apoptosis signaling remains a challenging task. We demonstrate how our proposed minimal model can help elucidate the effect of a specific apoptotic inhibitor Bcl-2 on apoptotic signaling in a cell-type independent manner. We also discuss the implications of our study in elucidating the adaptive strategy of cell death signaling pathways.Comment: 9 pages, 6 figure

Essential versus accessory aspects of cell death: recommendations of the NCCD 2015

Cells exposed to extreme physicochemical or mechanical stimuli die in an uncontrollable manner, as a result of their immediate structural breakdown. Such an unavoidable variant of cellular demise is generally referred to as ‘accidental cell death’ (ACD). In most settings, however, cell death is initiated by a genetically encoded apparatus, correlating with the fact that its course can be altered by pharmacologic or genetic interventions. ‘Regulated cell death’ (RCD) can occur as part of physiologic programs or can be activated once adaptive responses to perturbations of the extracellular or intracellular microenvironment fail. The biochemical phenomena that accompany RCD may be harnessed to classify it into a few subtypes, which often (but not always) exhibit stereotyped morphologic features. Nonetheless, efficiently inhibiting the processes that are commonly thought to cause RCD, such as the activation of executioner caspases in the course of apoptosis, does not exert true cytoprotective effects in the mammalian system, but simply alters the kinetics of cellular demise as it shifts its morphologic and biochemical correlates. Conversely, bona fide cytoprotection can be achieved by inhibiting the transduction of lethal signals in the early phases of the process, when adaptive responses are still operational. Thus, the mechanisms that truly execute RCD may be less understood, less inhibitable and perhaps more homogeneous than previously thought. Here, the Nomenclature Committee on Cell Death formulates a set of recommendations to help scientists and researchers to discriminate between essential and accessory aspects of cell death

Tuning hardness in calcite by incorporation of amino acids

Structural biominerals are inorganic/organic composites that exhibit remarkable mechanical properties. However, the structure–property relationships of even the simplest building unit—mineral single crystals containing embedded macromolecules—remain poorly understood. Here, by means of a model biomineral made from calcite single crystals containing glycine (0–7 mol%) or aspartic acid (0–4 mol%), we elucidate the origin of the superior hardness of biogenic calcite. We analysed lattice distortions in these model crystals by using X-ray diffraction and molecular dynamics simulations, and by means of solid-state nuclear magnetic resonance show that the amino acids are incorporated as individual molecules. We also demonstrate that nanoindentation hardness increased with amino acid content, reaching values equivalent to their biogenic counterparts. A dislocation pinning model reveals that the enhanced hardness is determined by the force required to cut covalent bonds in the molecules

Ultrahigh specificity in a network of computationally designed protein-interaction pairs

Protein networks in all organisms comprise homologous interacting pairs. In these networks, some proteins are specific, interacting with one or a few binding partners, whereas others are multispecific and bind a range of targets. We describe an algorithm that starts from an interacting pair and designs dozens of new pairs with diverse backbone conformations at the binding site as well as new binding orientations and sequences. Applied to a high-affinity bacterial pair, the algorithm results in 18 new ones, with cognate affinities from pico- to micromolar. Three pairs exhibit 3-5 orders of magnitude switch in specificity relative to the wild type, whereas others are multispecific, collectively forming a protein-interaction network. Crystallographic analysis confirms design accuracy, including in new backbones and polar interactions. Preorganized polar interaction networks are responsible for high specificity, thus defining design principles that can be applied to program synthetic cellular interaction networks of desired affinity and specificity

Relaxation oscillations and hierarchy of feedbacks in MAPK signaling

We formulated a computational model for a MAPK signaling cascade downstream of the EGF receptor to investigate how interlinked positive and negative feedback loops process EGF signals into ERK pulses of constant amplitude but dose-dependent duration and frequency. A positive feedback loop involving RAS and SOS, which leads to bistability and allows for switch-like responses to inputs, is nested within a negative feedback loop that encompasses RAS and RAF, MEK, and ERK that inhibits SOS via phosphorylation. This negative feedback, operating on a longer time scale, changes switch-like behavior into oscillations having a period of 1 hour or longer. Two auxiliary negative feedback loops, from ERK to MEK and RAF, placed downstream of the positive feedback, shape the temporal ERK activity profile but are dispensable for oscillations. Thus, the positive feedback introduces a hierarchy among negative feedback loops, such that the effect of a negative feedback depends on its position with respect to the positive feedback loop. Furthermore, a combination of the fast positive feedback involving slow-diffusing membrane components with slower negative feedbacks involving faster diffusing cytoplasmic components leads to local excitation/global inhibition dynamics, which allows the MAPK cascade to transmit paracrine EGF signals into spatially non-uniform ERK activity pulses.Peer reviewe

Modeling the TNFα-Induced Apoptosis Pathway in Hepatocytes

The proinflammatory cytokine TNFα fails to provoke cell death in isolated hepatocytes but has been implicated in hepatocyte apoptosis during liver diseases associated with chronic inflammation. Recently, we showed that TNFα is able to sensitize primary murine hepatocytes cultured on collagen to Fas ligand-induced apoptosis and presented a mathematical model of the sensitizing effect. Here, we analyze how TNFα induces apoptosis in combination with the transcriptional inhibitor actinomycin D (ActD). Accumulation of reactive oxygen species (ROS) in response to TNFR activation turns out to be critical for sustained activation of JNK which then triggers mitochondrial pathway-dependent apoptosis. In addition, the amount of JNK is strongly upregulated in a ROS-dependent way. In contrast to TNFα plus cycloheximide no cFLIP degradation is observed suggesting a different apoptosis pathway in which the Itch-mediated cFLIP degradation and predominantly caspase-8 activation is not involved. Time-resolved data of the respective pro- and antiapoptotic factors are obtained and subjected to mathematical modeling. On the basis of these data we developed a mathematical model which reproduces the complex interplay regulating the phosphorylation status of JNK and generation of ROS. This model was fully integrated with our model of TNFα/Fas ligand sensitizing as well as with a published NF-κB-model. The resulting comprehensive model delivers insight in the dynamical interplay between the TNFα and FasL pathways, NF-κB and ROS and gives an example for successful model integration

Elevated cerebrospinal fluid pressure in patients with Alzheimer's disease

BACKGROUND: Abnormalities in cerebrospinal fluid (CSF) production and turnover, seen in normal pressure hydrocephalus (NPH) and in Alzheimer's disease (AD), may be an important cause of amyloid retention in the brain and may relate the two diseases. There is a high incidence of AD pathology in patients being shunted for NPH, the AD-NPH syndrome. We now report elevated CSF pressure (CSFP), consistent with very early hydrocephalus, in a subset of AD patients enrolled in a clinical trial of chronic low-flow CSF drainage. Our objective was to determine the frequency of elevated CSFP in subjects meeting National Institutes of Neurological and Communicative Diseases and Stroke – Alzheimer's Disease and Related Disorders Association (NINCDS-ADRDA) criteria for AD, excluding those with signs of concomitant NPH. METHODS: AD subjects by NINCDS-ADRDA criteria (n = 222), were screened by history, neurological examination, and radiographic imaging to exclude those with clinical or radiographic signs of NPH. As part of this exclusion process, opening CSFP was measured supine under general anesthesia during device implantation surgery at a controlled pCO(2 )of 40 Torr (40 mmHg). RESULTS: Of the 222 AD subjects 181 had pressure measurements recorded. Seven subjects (3.9%) enrolled in the study had CSFP of 220 mmH(2)0 or greater, mean 249 ± 20 mmH(2)0 which was significantly higher than 103 ± 47 mmH(2)O for the AD-only group. AD-NPH patients were significantly younger and significantly less demented on the Mattis Dementia Rating Scale (MDRS). CONCLUSION: Of the AD subjects who were carefully screened to exclude those with clinical NPH, 4% had elevated CSFP. These subjects were presumed to have the AD-NPH syndrome and were withdrawn from the remainder of the study

Bcl-2 inhibits apoptosis by increasing the time-to-death and intrinsic cell-to-cell variations in the mitochondrial pathway of cell death

BH3 mimetics have been proposed as new anticancer therapeutics. They target anti-apoptotic Bcl-2 proteins, up-regulation of which has been implicated in the resistance of many cancer cells, particularly leukemia and lymphoma cells, to apoptosis. Using probabilistic computational modeling of the mitochondrial pathway of apoptosis, verified by single-cell experimental observations, we develop a model of Bcl-2 inhibition of apoptosis. Our results clarify how Bcl-2 imparts its anti-apoptotic role by increasing the time-to-death and cell-to-cell variability. We also show that although the commitment to death is highly impacted by differences in protein levels at the time of stimulation, inherent stochastic fluctuations in apoptotic signaling are sufficient to induce cell-to-cell variability and to allow single cells to escape death. This study suggests that intrinsic cell-to-cell stochastic variability in apoptotic signaling is sufficient to cause fractional killing of cancer cells after exposure to BH3 mimetics. This is an unanticipated facet of cancer chemoresistance.Comment: 11 pages, In pres

Stochastic Competition between Mechanistically Independent Slippage and Death Pathways Determines Cell Fate during Mitotic Arrest

Variability in cell-to-cell behavior within clonal populations can be attributed to the inherent stochasticity of biochemical reactions. Most single-cell studies have examined variation in behavior due to randomness in gene transcription. Here we investigate the mechanism of cell fate choice and the origin of cell-to-cell variation during mitotic arrest, when transcription is silenced. Prolonged mitotic arrest is commonly observed in cells treated with anti-mitotic drugs. Cell fate during mitotic arrest is determined by two alternative pathways, one promoting cell death, the other promoting cyclin B1 degradation, which leads to mitotic slippage and survival. It has been unclear whether these pathways are mechanistically coupled or independent. In this study we experimentally uncoupled these two pathways using zVAD-fmk to block cell death or Cdc20 knockdown to block slippage. We then used time-lapse imaging to score the kinetics of single cells adopting the remaining fate. We also used kinetic simulation to test whether the behaviors of death versus slippage in cell populations where both pathways are active can be quantitatively recapitulated by a model that assumes stochastic competition between the pathways. Our data are well fit by a model where the two pathways are mechanistically independent, and cell fate is determined by a stochastic kinetic competition between them that results in cell-to-cell variation