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

Bistability in Apoptosis by Receptor Clustering

Apoptosis is a highly regulated cell death mechanism involved in many physiological processes. A key component of extrinsically activated apoptosis is the death receptor Fas, which, on binding to its cognate ligand FasL, oligomerize to form the death-inducing signaling complex. Motivated by recent experimental data, we propose a mathematical model of death ligand-receptor dynamics where FasL acts as a clustering agent for Fas, which form locally stable signaling platforms through proximity-induced receptor interactions. Significantly, the model exhibits hysteresis, providing an upstream mechanism for bistability and robustness. At low receptor concentrations, the bistability is contingent on the trimerism of FasL. Moreover, irreversible bistability, representing a committed cell death decision, emerges at high concentrations, which may be achieved through receptor pre-association or localization onto membrane lipid rafts. Thus, our model provides a novel theory for these observed biological phenomena within the unified context of bistability. Importantly, as Fas interactions initiate the extrinsic apoptotic pathway, our model also suggests a mechanism by which cells may function as bistable life/death switches independently of any such dynamics in their downstream components. Our results highlight the role of death receptors in deciding cell fate and add to the signal processing capabilities attributed to receptor clustering.Comment: Accepted by PLoS Comput Bio

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

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

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

XIAP impairs Smac release from the mitochondria during apoptosis

X-linked inhibitor of apoptosis protein (XIAP) is a potent inhibitor of caspases 3, 7 and 9, and mitochondrial Smac (second mitochondria-derived activator of caspase) release during apoptosis inhibits the activity of XIAP. In this study we show that cytosolic XIAP also feeds back to mitochondria to impair Smac release. We constructed a fluorescent XIAP-fusion protein by labelling NH2- and COOH-termini with Cerulean fluorescent protein (C-XIAP-C). Immunoprecipitation confirmed that C-XIAP-C retained the ability to interact with Smac and impaired extrinsically and intrinsically activated apoptosis in response to tumour necrosis factor-related apoptosis-inducing ligand/cycloheximide and staurosporine. In C-XIAP-C-expressing cells, cytochrome c release from mitochondria proceeded normally, whereas Smac release was significantly prolonged and incomplete. In addition, physiological expression of native XIAP prolonged or limited Smac release in HCT-116 colon cancer cells and primary mouse cortical neurons. The Smac-binding capacity of XIAP, but not caspase inhibition, was central for mitochondrial Smac retention, as evidenced in experiments using XIAP mutants that cannot bind to Smac or effector caspases. Similarly, the release of a Smac mutant that cannot bind to XIAP was not impaired by C-XIAP-C expression. Full Smac release could however be provoked by rapid cytosolic C-XIAP-C depletion upon digitonin-induced plasma membrane permeabilization. Our findings suggest that although mitochondria may already contain pores sufficient for cytochrome c release, elevated amounts of XIAP can selectively impair and limit the release of Smac

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

VASCo: computation and visualization of annotated protein surface contacts

<p>Abstract</p> <p>Background</p> <p>Structural data from crystallographic analyses contain a vast amount of information on protein-protein contacts. Knowledge on protein-protein interactions is essential for understanding many processes in living cells. The methods to investigate these interactions range from genetics to biophysics, crystallography, bioinformatics and computer modeling. Also crystal contact information can be useful to understand biologically relevant protein oligomerisation as they rely in principle on the same physico-chemical interaction forces. Visualization of crystal and biological contact data including different surface properties can help to analyse protein-protein interactions.</p> <p>Results</p> <p>VASCo is a program package for the calculation of protein surface properties and the visualization of annotated surfaces. Special emphasis is laid on protein-protein interactions, which are calculated based on surface point distances. The same approach is used to compare surfaces of two aligned molecules. Molecular properties such as electrostatic potential or hydrophobicity are mapped onto these surface points. Molecular surfaces and the corresponding properties are calculated using well established programs integrated into the package, as well as using custom developed programs. The modular package can easily be extended to include new properties for annotation. The output of the program is most conveniently displayed in PyMOL using a custom-made plug-in.</p> <p>Conclusion</p> <p>VASCo supplements other available protein contact visualisation tools and provides additional information on biological interactions as well as on crystal contacts. The tool provides a unique feature to compare surfaces of two aligned molecules based on point distances and thereby facilitates the visualization and analysis of surface differences.</p