Mathematical modeling of the formation of apoptosome in intrinsic pathway of apoptosis.
Ontology highlight
ABSTRACT:
SUBMITTER: Ryu S
PROVIDER: S-EPMC2671592 | biostudies-literature | 2008 Jun
REPOSITORIES: biostudies-literature
Publications
Mathematical modeling of the formation of apoptosome in intrinsic pathway of apoptosis.
Ryu Seongho S Lin Shih-Chieh SC Ugel Nadia N Antoniotti Marco M Mishra Bud B
Systems and synthetic biology 20080601 1-2
Caspase-9 is the protease that mediates the intrinsic pathway of apoptosis, a type of cell death. Activation of caspase-9 is a multi-step process that requires dATP or ATP and involves at least two proteins, cytochrome c and Apaf-1. In this study, we mathematically model caspase-9 activation by using a system of ordinary differential equations (an ODE model) generated by a systems biology tool Simpathica-a simulation and reasoning system, developed to study biological pathways. A rudimentary ver ...[more]
Similar Datasets
Project description:BackgroundApoptosis is a cell suicide mechanism that enables multicellular organisms to maintain homeostasis and to eliminate individual cells that threaten the organism's survival. Dependent on the type of stimulus, apoptosis can be propagated by extrinsic pathway or intrinsic pathway. The comprehensive understanding of the molecular mechanism of apoptotic signaling allows for development of mathematical models, aiming to elucidate dynamical and systems properties of apoptotic signaling networks. There have been extensive efforts in modeling deterministic apoptosis network accounting for average behavior of a population of cells. Cellular networks, however, are inherently stochastic and significant cell-to-cell variability in apoptosis response has been observed at single cell level.ResultsTo address the inevitable randomness in the intrinsic apoptosis mechanism, we develop a theoretical and computational modeling framework of intrinsic apoptosis pathway at single-cell level, accounting for both deterministic and stochastic behavior. Our deterministic model, adapted from the well-accepted Fussenegger model, shows that an additional positive feedback between the executioner caspase and the initiator caspase plays a fundamental role in yielding the desired property of bistability. We then examine the impact of intrinsic fluctuations of biochemical reactions, viewed as intrinsic noise, and natural variation of protein concentrations, viewed as extrinsic noise, on behavior of the intrinsic apoptosis network. Histograms of the steady-state output at varying input levels show that the intrinsic noise could elicit a wider region of bistability over that of the deterministic model. However, the system stochasticity due to intrinsic fluctuations, such as the noise of steady-state response and the randomness of response delay, shows that the intrinsic noise in general is insufficient to produce significant cell-to-cell variations at physiologically relevant level of molecular numbers. Furthermore, the extrinsic noise represented by random variations of two key apoptotic proteins, namely Cytochrome C and inhibitor of apoptosis proteins (IAP), is modeled separately or in combination with intrinsic noise. The resultant stochasticity in the timing of intrinsic apoptosis response shows that the fluctuating protein variations can induce cell-to-cell stochastic variability at a quantitative level agreeing with experiments. Finally, simulations illustrate that the mean abundance of fluctuating IAP protein is positively correlated with the degree of cellular stochasticity of the intrinsic apoptosis pathway.ConclusionsOur theoretical and computational study shows that the pronounced non-genetic heterogeneity in intrinsic apoptosis responses among individual cells plausibly arises from extrinsic rather than intrinsic origin of fluctuations. In addition, it predicts that the IAP protein could serve as a potential therapeutic target for suppression of the cell-to-cell variation in the intrinsic apoptosis responsiveness.
| S-EPMC3733900 | biostudies-literature
Project description:Introduction: Mathematical model can be used to model complex biological processes, and have shown potential in describing apoptosis in chondrocytes. Method: In order to investigate the regulatory mechanisms of TNF signaling pathway in regulating chondrocyte apoptosis, a fractional-order differential equation model is proposed to describe the dynamic behavior and mutual interaction of apoptosis-related genes under the activation of TNF signaling pathway. Compared with the traditional molecular biology techniques, the proposed mathematical modeling has advantages to providing a more comprehensive understanding of the regulatory mechanisms of TNF signaling pathway in chondrocyte apoptosis. Result: In this paper, differentially expressed genes induced by IL-1β in human chondrocyte apoptosis are screened using high-throughput sequencing. It is found that they were significantly enriched in the TNF signaling pathway. Therefore, a mathematical model of the TNF signaling pathway is built. Using real-time PCR experiments, mRNA data is measured and used to identify the model parameters, as well as the correlation coefficient. Finally, the sensitivity of the model parameters is discussed by using numerical simulation methods, which can be used to predict the effects of different interventions and explore the optimal intervention strategies for regulating chondrocyte apoptosis. Discussion: Therefore, fractional-order differential equation modeling plays an important role in understanding the regulatory mechanisms of TNF signaling pathway in chondrocyte apoptosis and its potential clinical applications.
| S-EPMC10652750 | biostudies-literature
Project description:Trypanosome RNA polymerase II transcription is polycistronic, individual mRNAs being excised by trans splicing and polyadenylation. In this study, we refined the previously published mathematical model for bloodstream form parasites and extended it to the procyclic form. We used the model, together with known mRNA half-lives, to predict the abundances of individual mRNAs, assuming rapid, unregulated mRNA processing; then we compared the results with measured mRNA abundances. Remarkably, the abundances of most mRNAs in procyclic forms are predicted quite well by the model, being largely explained by variations in mRNA decay rates and length. In bloodstream forms substantially more mRNAs are less abundant than predicted. We list mRNAs that are likely to show particularly slow or inefficient processing, either in both forms or with developmental regulation. We also measured ribosome occupancies of all mRNAs in trypanosomes grown in the same conditions as were used to measure mRNA turnover. In procyclic forms there was a weak positive correlation between ribosome density and mRNA half-life, suggesting cross-talk between translation and mRNA decay; ribosome density was related to the proportion of the mRNA on polysomes, indicating control of translation initiation. Ribosomal protein mRNAs in procyclics appeared to be exceptionally rapidly processed but poorly translated. Through this study, we conclude that lLevels of mRNAs in procyclic form trypanosomes are determined mainly by length and mRNA decay, with some control of precursor processing. In bloodstream forms variations in nuclear events play a larger role in transcriptome regulation, suggesting acquisition of new control mechanisms during adaptation to mammalian parasitism. Ribosome profiling and mRNA libraries were constructed in triplicate from in vitro PCF and in duplicate from in vitro T. brucei Lister427, to understand global differntial gene transcription.
2016-03-16 | E-GEOD-79208 | biostudies-arrayexpress
Project description:Trypanosome RNA polymerase II transcription is polycistronic, individual mRNAs being excised by trans splicing and polyadenylation. In this study, we refined the previously published mathematical model for bloodstream form parasites and extended it to the procyclic form. We used the model, together with known mRNA half-lives, to predict the abundances of individual mRNAs, assuming rapid, unregulated mRNA processing; then we compared the results with measured mRNA abundances. Remarkably, the abundances of most mRNAs in procyclic forms are predicted quite well by the model, being largely explained by variations in mRNA decay rates and length. In bloodstream forms substantially more mRNAs are less abundant than predicted. We list mRNAs that are likely to show particularly slow or inefficient processing, either in both forms or with developmental regulation. We also measured ribosome occupancies of all mRNAs in trypanosomes grown in the same conditions as were used to measure mRNA turnover. In procyclic forms there was a weak positive correlation between ribosome density and mRNA half-life, suggesting cross-talk between translation and mRNA decay; ribosome density was related to the proportion of the mRNA on polysomes, indicating control of translation initiation. Ribosomal protein mRNAs in procyclics appeared to be exceptionally rapidly processed but poorly translated. Through this study, we conclude that lLevels of mRNAs in procyclic form trypanosomes are determined mainly by length and mRNA decay, with some control of precursor processing. In bloodstream forms variations in nuclear events play a larger role in transcriptome regulation, suggesting acquisition of new control mechanisms during adaptation to mammalian parasitism.
2016-03-16 | GSE79208 | GEO
Project description:The protease caspase-9 is activated on the apoptosome, a multiprotein signal transduction platform that assembles in response to mitochondria-dependent apoptosis initiation. Despite extensive molecular research, the assembly of the holo-apoptosome and the process of caspase-9 activation remain incompletely understood. Here, we therefore integrated quantitative data on the molecular interactions and proteolytic processes during apoptosome formation and apoptosis execution and conducted mathematical simulations to investigate the resulting biochemical signaling, quantitatively and kinetically. Interestingly, when implementing the homodimerization of procaspase-9 as a prerequisite for activation, the calculated kinetics of apoptosis execution and the efficacy of caspase-3 activation failed to replicate experimental data. In contrast, assuming a scenario in which procaspase-9 is activated allosterically upon binding to the apoptosome backbone, the mathematical simulations quantitatively and kinetically reproduced all experimental data. These data included a XIAP threshold concentration at which apoptosis execution is suppressed in HeLa cervical cancer cells, half-times of procaspase-9 processing, as well as the molecular timer function of the apoptosome. Our study therefore provides novel mechanistic insight into apoptosome-dependent apoptosis execution and suggests that caspase-9 is activated allosterically by binding to the apoptosome backbone. Our findings challenge the currently prevailing dogma that all initiator procaspases require homodimerization for activation.
| S-EPMC4176236 | biostudies-literature
Project description:The mechanistic target of rapamycin (mTOR) is a central regulatory pathway that integrates a variety of environmental cues to control cellular growth and homeostasis by intricate molecular feedbacks. In spite of extensive knowledge about its components, the molecular understanding of how these function together in space and time remains poor and there is a need for Systems Biology approaches to perform systematic analyses. In this work, we review the recent progress how the combined efforts of mathematical models and quantitative experiments shed new light on our understanding of the mTOR signaling pathway. In particular, we discuss the modeling concepts applied in mTOR signaling, the role of multiple feedbacks and the crosstalk mechanisms of mTOR with other signaling pathways. We also discuss the contribution of principles from information and network theory that have been successfully applied in dissecting design principles of the mTOR signaling network. We finally propose to classify the mTOR models in terms of the time scale and network complexity, and outline the importance of the classification toward the development of highly comprehensive and predictive models. WIREs Syst Biol Med 2017, 9:e1379. doi: 10.1002/wsbm.1379 For further resources related to this article, please visit the WIREs website.
| S-EPMC5573916 | biostudies-literature
Project description:Understanding the dynamics of dormant cells in microbial biofilms, in which the bacteria are embedded in extracellular matrix, is important for developing successful antibiotic therapies against pathogenic bacteria. Although some of the molecular mechanisms leading to bacterial persistence have been speculated in planktonic bacterial cell, how dormant cells emerge in the biofilms of pathogenic bacteria such as Pseudomonas aeruginosa remains unclear. The present study proposes four hypotheses of dormant cell formation; stochastic process, nutrient-dependent, oxygen-dependent, and time-dependent processes. These hypotheses were implemented into a three-dimensional individual-based model of biofilm formation. Numerical simulations of the different mechanisms yielded qualitatively different spatiotemporal distributions of dormant cells in the growing biofilm. Based on these simulation results, we discuss what kinds of experimental studies are effective for discriminating dormant cell formation mechanisms in biofilms.
| S-EPMC4446547 | biostudies-literature
Project description:Early studies of the free-living nematode C. elegans informed us how BCL-2-regulated apoptosis in humans is regulated. However, subsequent studies showed C. elegans apoptosis has several unique features compared with human apoptosis. To date, there has been no detailed analysis of apoptosis regulators in nematodes other than C. elegans. Here, we discovered BCL-2 orthologues in 89 free-living and parasitic nematode taxa representing four evolutionary clades (I, III, IV and V). Unlike in C. elegans, 15 species possess multiple (two to five) BCL-2-like proteins, and some do not have any recognisable BCL-2 sequences. Functional studies provided no evidence that BAX/BAK proteins have evolved in nematodes, and structural studies of a BCL-2 protein from the basal clade I revealed it lacks a functionally important feature of the C. elegans orthologue. Clade I CED-4/APAF-1 proteins also possess WD40-repeat sequences associated with apoptosome assembly, not present in C. elegans, or other nematode taxa studied.
| S-EPMC7456325 | biostudies-literature
Project description:K-Ras functions as a critical node in the mitogen-activated protein kinase (MAPK) pathway that regulates key cellular functions including proliferation, differentiation, and apoptosis. Following growth factor receptor activation K-Ras.GTP forms nanoclusters on the plasma membrane through interaction with the scaffold protein galectin-3. The generation of nanoclusters is essential for high fidelity signal transduction via the MAPK pathway. To explore the mechanisms underlying K-Ras.GTP nanocluster formation, we developed a mathematical model of K-Ras-galectin-3 interactions. We designed a computational method to calculate protein collision rates based on experimentally determined protein diffusion rates and diffusion mechanisms and used a genetic algorithm to search the values of key model parameters. The optimal estimated model parameters were validated using experimental data. The resulting model accurately replicates critical features of K-Ras nanoclustering, including a fixed ratio of clustered K-Ras.GTP to monomeric K-Ras.GTP that is independent of the concentration of K-Ras.GTP. The model reproduces experimental results showing that the cytosolic level of galectin-3 determines the magnitude of the K-Ras.GTP clustered fraction and illustrates that nanoclustering is regulated by key nonequilibrium processes. Our kinetic model identifies a potential biophysical mechanism for K-Ras nanoclustering and suggests general principles that may be relevant for other plasma-membrane-localized proteins.
| S-EPMC2905126 | biostudies-literature
Project description:Synaptic activity promotes resistance to diverse apoptotic insults, the mechanism behind which is incompletely understood. We show here that a coordinated downregulation of core components of the intrinsic apoptosis pathway by neuronal activity forms a key part of the underlying mechanism. Activity-dependent protection against apoptotic insults is associated with inhibition of cytochrome c release in most but not all neurons, indicative of anti-apoptotic signaling both upstream and downstream of this step. We find that enhanced firing activity suppresses expression of the proapoptotic BH3-only member gene Puma in a NMDA receptor-dependent, p53-independent manner. Puma expression is sufficient to induce cytochrome c loss and neuronal apoptosis. Puma deficiency protects neurons against apoptosis and also occludes the protective effect of synaptic activity, while blockade of physiological NMDA receptor activity in the developing mouse brain induces neuronal apoptosis that is preceded by upregulation of Puma. However, enhanced activity can also confer resistance to Puma-induced apoptosis, acting downstream of cytochrome c release. This mechanism is mediated by transcriptional suppression of apoptosome components Apaf-1 and procaspase-9, and limiting caspase-9 activity, since overexpression of procaspase-9 accelerates the rate of apoptosis in active neurons back to control levels. Synaptic activity does not exert further significant anti-apoptotic effects downstream of caspase-9 activation, since an inducible form of caspase-9 overrides the protective effect of synaptic activity, despite activity-induced transcriptional suppression of caspase-3. Thus, suppression of apoptotic gene expression may synergize with other activity-dependent events such as enhancement of antioxidant defenses to promote neuronal survival.
| S-EPMC2834927 | biostudies-literature