Modeling small cell lung cancer (SCLC) biology through deterministic and stochastic mathematical models.
Ontology highlight
ABSTRACT:
SUBMITTER: Salgia R
PROVIDER: S-EPMC5995226 | biostudies-other | 2018 May
REPOSITORIES: biostudies-other
Publications
Modeling small cell lung cancer (SCLC) biology through deterministic and stochastic mathematical models.
Salgia Ravi R Mambetsariev Isa I Hewelt Blake B Achuthan Srisairam S Li Haiqing H Poroyko Valeriy V Wang Yingyu Y Sattler Martin M
Oncotarget 20180525 40
Mathematical cancer models are immensely powerful tools that are based in part on the fractal nature of biological structures, such as the geometry of the lung. Cancers of the lung provide an opportune model to develop and apply algorithms that capture changes and disease phenotypes. We reviewed mathematical models that have been developed for biological sciences and applied them in the context of small cell lung cancer (SCLC) growth, mutational heterogeneity, and mechanisms of metastasis. The u ...[more]
Similar Datasets
Project description:Hybrid deterministic-stochastic methods provide an efficient alternative to a fully stochastic treatment of models which include components with disparate levels of stochasticity. However, general-purpose hybrid solvers for spatially resolved simulations of reaction-diffusion systems are not widely available. Here we describe fundamentals of a general-purpose spatial hybrid method. The method generates realizations of a spatially inhomogeneous hybrid system by appropriately integrating capabilities of a deterministic partial differential equation solver with a popular particle-based stochastic simulator, Smoldyn. Rigorous validation of the algorithm is detailed, using a simple model of calcium 'sparks' as a testbed. The solver is then applied to a deterministic-stochastic model of spontaneous emergence of cell polarity. The approach is general enough to be implemented within biologist-friendly software frameworks such as Virtual Cell.
| S-EPMC5154471 | biostudies-literature
Project description:BackgroundStochastic and asymptotic methods are powerful tools in developing multiscale systems biology models; however, little has been done in this context to compare the efficacy of these methods. The majority of current systems biology modelling research, including that of auxin transport, uses numerical simulations to study the behaviour of large systems of deterministic ordinary differential equations, with little consideration of alternative modelling frameworks.ResultsIn this case study, we solve an auxin-transport model using analytical methods, deterministic numerical simulations and stochastic numerical simulations. Although the three approaches in general predict the same behaviour, the approaches provide different information that we use to gain distinct insights into the modelled biological system. We show in particular that the analytical approach readily provides straightforward mathematical expressions for the concentrations and transport speeds, while the stochastic simulations naturally provide information on the variability of the system.ConclusionsOur study provides a constructive comparison which highlights the advantages and disadvantages of each of the considered modelling approaches. This will prove helpful to researchers when weighing up which modelling approach to select. In addition, the paper goes some way to bridging the gap between these approaches, which in the future we hope will lead to integrative hybrid models.
| S-EPMC2873313 | biostudies-literature
Project description:Small-cell lung cancer (SCLC) accounts for 13-15% of all new lung cancer cases in the US. The tumor has a tendency to disseminate early resulting in 80-85% of patients being diagnosed with extensive disease (ES-SCLC). Chemotherapy has provided SCLC patients considerable survival benefits over the past three decades. Nonetheless, most patients relapse and rarely survive beyond 2 years. Despite consistent overall response rates of ≥50%, until recently, median survival times and 2-year survivals only ranged between 7-10 months and 10-20%, respectively. Several chemotherapy agents possess activity against SCLC, both, as single agents and in combinations but etoposide-platinum emerged as the preferred first line regimen. Upon relapse, many patients remain candidates for additional therapy. However, the sensitivity of relapsed SCLC to further therapies is markedly reduced and dependent upon the level and duration of response to the initial treatment (platinum-sensitive vs. resistant relapse). Multiple factors suggest a therapeutic role for immunotherapy in SCLC: SCLC has been associated with immune-mediated paraneoplastic processes (cerebellar degeneration, limbic encephalitis, and Lambert-Eaton syndrome) and patients presenting with these paraneoplastic syndromes have shown more favorable outcomes, suggesting an underlying immune response mechanism.Comprehensive genomic profiling of SCLC indicates that the majority lack functional p53 (90%) and Rb1 (65%). These universal genetic aberrations facilitate poor genomic stability, thus perpetuating the generation of tumor associated antigens, amenable to targeting with immunotherapy.SCLC has one of the highest mutational loads, likely a reflection of the myriad of insults inflicted by smoking-related carcinogens. The relationship between tumor mutational load and response to immune checkpoint inhibitors has been established in multiple solid tumors, including preliminary results in relapsed SCLC. In this manuscript, we review the early (some failed and discontinued, some partly successful, and still ongoing) attempts to incorporate immunotherapy (particularly vaccine based approaches) to the treatment of SCLC, and the latest attempts (mostly incorporating the use of checkpoint inhibitors), including those with favorable but preliminary results (CheckMate 032, Keynote 028 and 158), and those with more definitive positive (iMpower 133 and CASPIAN) and negative (CheckMate 331 and 451) results.
| S-EPMC7378389 | biostudies-literature
Project description:Experimental studies of the flowering of Arabidopsis thaliana have shown that a large complex gene regulatory network (GRN) is responsible for its regulation. This process has been mathematically modelled with deterministic differential equations by considering the interactions between gene activators and inhibitors (Valentim et al. in PLoS ONE 10(2):e0116973, 2015; van Mourik et al. in BMC Syst Biol 4(1):1, 2010). However, due to complexity of the model, the properties of the network and the roles of the individual genes cannot be deducted from the numerical solution the published work offers. Here, we propose simplifications of the model, based on decoupling of the original GRN to motifs, described with three and two differential equations. A stable solution of the original model is sought by linearisation of the original model which contributes to further investigation of the role of the individual genes to the flowering. Furthermore, we study the role of noise by introducing and investigating two types of stochastic elements into the model. The deterministic and stochastic nonlinear dynamic models of Arabidopsis flowering time are considered by following the deterministic delayed model introduced in Valentim et al. (2015). Steady-state regimes and stability of the deterministic original model are investigated analytically and numerically. By decoupling some concentrations, the system was reduced to emphasise the role played by the transcription factor Suppressor of Overexpression of Constants1 ([Formula: see text]) and the important floral meristem identity genes, Leafy ([Formula: see text]) and Apetala1 ([Formula: see text]). Two-dimensional motifs, based on the dynamics of [Formula: see text] and [Formula: see text], are obtained from the reduced network and parameter ranges ensuring flowering are determined. Their stability analysis shows that [Formula: see text] and [Formula: see text] are regulating each other for flowering, matching experimental findings. New sufficient conditions of mean square stability in the stochastic model are obtained using a stochastic Lyapunov approach. Our numerical simulations demonstrate that the reduced models of Arabidopsis flowering time, describing specific motifs of the GRN, can capture the essential behaviour of the full system and also introduce the conditions of flowering initiation. Additionally, they show that stochastic effects can change the behaviour of the stability region through a stability switch. This study thus contributes to a better understanding of the role of [Formula: see text] and [Formula: see text] in Arabidopsis flowering.
| S-EPMC6320361 | biostudies-other
Project description:BackgroundGillespie's stochastic simulation algorithm (SSA) for chemical reactions admits three kinds of elementary processes, namely, mass action reactions of 0th, 1st or 2nd order. All other types of reaction processes, for instance those containing non-integer kinetic orders or following other types of kinetic laws, are assumed to be convertible to one of the three elementary kinds, so that SSA can validly be applied. However, the conversion to elementary reactions is often difficult, if not impossible. Within deterministic contexts, a strategy of model reduction is often used. Such a reduction simplifies the actual system of reactions by merging or approximating intermediate steps and omitting reactants such as transient complexes. It would be valuable to adopt a similar reduction strategy to stochastic modelling. Indeed, efforts have been devoted to manipulating the chemical master equation (CME) in order to achieve a proper propensity function for a reduced stochastic system. However, manipulations of CME are almost always complicated, and successes have been limited to relative simple cases.ResultsWe propose a rather general strategy for converting a deterministic process model into a corresponding stochastic model and characterize the mathematical connections between the two. The deterministic framework is assumed to be a generalized mass action system and the stochastic analogue is in the format of the chemical master equation. The analysis identifies situations: where a direct conversion is valid; where internal noise affecting the system needs to be taken into account; and where the propensity function must be mathematically adjusted. The conversion from deterministic to stochastic models is illustrated with several representative examples, including reversible reactions with feedback controls, Michaelis-Menten enzyme kinetics, a genetic regulatory motif, and stochastic focusing.ConclusionsThe construction of a stochastic model for a biochemical network requires the utilization of information associated with an equation-based model. The conversion strategy proposed here guides a model design process that ensures a valid transition between deterministic and stochastic models.
| S-EPMC3236013 | biostudies-literature
Project description:BackgroundCyclic adenosine monophosphate (cAMP) has a key signaling role in all eukaryotic organisms. In Saccharomyces cerevisiae, it is the second messenger in the Ras/PKA pathway which regulates nutrient sensing, stress responses, growth, cell cycle progression, morphogenesis, and cell wall biosynthesis. A stochastic model of the pathway has been reported.ResultsWe have created deterministic mathematical models of the PKA module of the pathway, as well as the complete cAMP pathway. First, a simplified conceptual model was created which reproduced the dynamics of changes in cAMP levels in response to glucose addition in wild-type as well as cAMP phosphodiesterase deletion mutants. This model was used to investigate the role of the regulatory Krh proteins that had not been included previously. The Krh-containing conceptual model reproduced very well the experimental evidence supporting the role of Krh as a direct inhibitor of PKA. These results were used to develop the Complete cAMP Model. Upon simulation it illustrated several important features of the yeast cAMP pathway: Pde1p is more important than is Pde2p for controlling the cAMP levels following glucose pulses; the proportion of active PKA is not directly proportional to the cAMP level, allowing PKA to exert negative feedback; negative feedback mechanisms include activating Pde1p and deactivating Ras2 via phosphorylation of Cdc25. The Complete cAMP model is easier to simulate, and although significantly simpler than the existing stochastic one, it recreates cAMP levels and patterns of changes in cAMP levels observed experimentally in vivo in response to glucose addition in wild-type as well as representative mutant strains such as pde1Delta, pde2Delta, cyr1Delta, and others. The complete model is made available in SBML format.ConclusionWe suggest that the lower number of reactions and parameters makes these models suitable for integrating them with models of metabolism or of the cell cycle in S. cerevisiae. Similar models could be also useful for studies in the human pathogen Candida albicans as well as other less well-characterized fungal species.
| S-EPMC2719611 | biostudies-literature
Project description:BackgroundOncolytic viruses that specifically target tumor cells are promising anti-cancer therapeutic agents. The interaction between an oncolytic virus and tumor cells is amenable to mathematical modeling using adaptations of techniques employed previously for modeling other types of virus-cell interaction.ResultsA complete parametric analysis of dynamic regimes of a conceptual model of anti-tumor virus therapy is presented. The role and limitations of mass-action kinetics are discussed. A functional response, which is a function of the ratio of uninfected to infected tumor cells, is proposed to describe the spread of the virus infection in the tumor. One of the main mathematical features of ratio-dependent models is that the origin is a complicated equilibrium point whose characteristics determine the main properties of the model. It is shown that, in a certain area of parameter values, the trajectories of the model form a family of homoclinics to the origin (so-called elliptic sector). Biologically, this means that both infected and uninfected tumor cells can be eliminated with time, and complete recovery is possible as a result of the virus therapy within the framework of deterministic models.ConclusionOur model, in contrast to the previously published models of oncolytic virus-tumor interaction, exhibits all possible outcomes of oncolytic virus infection, i.e., no effect on the tumor, stabilization or reduction of the tumor load, and complete elimination of the tumor. The parameter values that result in tumor elimination, which is, obviously, the desired outcome, are compatible with some of the available experimental data.ReviewersThis article was reviewed by Mikhail Blagosklonny, David Krakauer, Erik Van Nimwegen, and Ned Wingreen.Open peer reviewReviewed by Mikhail Blagosklonny, David Krakauer, Erik Van Nimwegen, and Ned Wingreen. For the full reviews, please go to the Reviewers' comments section.
| S-EPMC1403749 | biostudies-literature
Project description:Network of signaling proteins and functional interaction between the infected cell and the leishmanial parasite, though are not well understood, may be deciphered computationally by reconstructing the immune signaling network. As we all know signaling pathways are well-known abstractions that explain the mechanisms whereby cells respond to signals, collections of pathways form networks, and interactions between pathways in a network, known as cross-talk, enables further complex signaling behaviours. In silico perturbations can help identify sensitive crosstalk points in the network which can be pharmacologically tested. In this study, we have developed a model for immune signaling cascade in leishmaniasis and based upon the interaction analysis obtained through simulation, we have developed a model network, between four signaling pathways i.e., CD14, epidermal growth factor (EGF), tumor necrotic factor (TNF) and PI3 K mediated signaling. Principal component analysis of the signaling network showed that EGF and TNF pathways can be potent pharmacological targets to curb leishmaniasis. The approach is illustrated with a proposed workable model of epidermal growth factor receptor (EGFR) that modulates the immune response. EGFR signaling represents a critical junction between inflammation related signal and potent cell regulation machinery that modulates the expression of cytokines.
| S-EPMC3824814 | biostudies-other
Project description:We study the calcium-induced vesicle release into the synaptic cleft using a deterministic algorithm and MCell, a Monte Carlo algorithm that tracks individual molecules. We compare the average vesicle release probability obtained using both algorithms and investigate the effect of the three main sources of noise: diffusion, sensor kinetics and fluctuations from the voltage-dependent calcium channels (VDCCs). We find that the stochastic opening kinetics of the VDCCs are the main contributors to differences in the release probability. Our results show that the deterministic calculations lead to reliable results, with an error of less than 20%, when the sensor is located at least 50 nm from the VDCCs, corresponding to microdomain signaling. For smaller distances, i.e. nanodomain signaling, the error becomes larger and a stochastic algorithm is necessary.
| S-EPMC2892017 | biostudies-literature
Project description:We present deterministic and stochastic programming models for the workover rig problem, one of the most challenging problems in the oil industry. In the deterministic approach, an integer linear programming model is used to determine the rig fleet size and schedule needed to service wells while maximizing oil production and minimizing rig usage cost. The stochastic approach is an extension of the deterministic method and relies on a two-stage stochastic programming model to define the optimal rig fleet size considering uncertainty in the intervention time. In this approach, different scenario-generation methods are compared. Several experiments were performed using instances based on real-world problems. The results suggest that the proposed methodology can be used to solve large instances and produces quality solutions in computationally reasonable times.
| S-EPMC6156096 | biostudies-literature