1
|
He Y, Battista NA, Waldrop LD. Mixed uncertainty analysis on pumping by peristaltic hearts using Dempster-Shafer theory. J Math Biol 2024; 89:13. [PMID: 38879850 DOI: 10.1007/s00285-024-02116-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2024] [Revised: 05/22/2024] [Accepted: 06/04/2024] [Indexed: 06/28/2024]
Abstract
In this paper, we introduce the numerical strategy for mixed uncertainty propagation based on probability and Dempster-Shafer theories, and apply it to the computational model of peristalsis in a heart-pumping system. Specifically, the stochastic uncertainty in the system is represented with random variables while epistemic uncertainty is represented using non-probabilistic uncertain variables with belief functions. The mixed uncertainty is propagated through the system, resulting in the uncertainty in the chosen quantities of interest (QoI, such as flow volume, cost of transport and work). With the introduced numerical method, the uncertainty in the statistics of QoIs will be represented using belief functions. With three representative probability distributions consistent with the belief structure, global sensitivity analysis has also been implemented to identify important uncertain factors and the results have been compared between different peristalsis models. To reduce the computational cost, physics constrained generalized polynomial chaos method is adopted to construct cheaper surrogates as approximations for the full simulation.
Collapse
Affiliation(s)
- Yanyan He
- Department of Mathematics, and of Computer Science and Engineering, University of North Texas, 1155 Union Circle, Denton, TX, 76203, USA.
| | - Nicholas A Battista
- Department of Mathematics and Statistics, The College of New Jersey, Pennington Rd, Ewing Township, NJ, 08618, USA
| | - Lindsay D Waldrop
- Schmid College of Science and Technology, Chapman University, One University Drive, Orange, CA, 92866, USA
| |
Collapse
|
2
|
Aggarwal M, Lewis O, Jarrett A, Hussaini MY, Cogan NG. A Model of Gastric Mucosal pH Regulation: Extending Sensitivity Analysis Using Sobol' Indices to Understand Higher Moments. Bull Math Biol 2024; 86:77. [PMID: 38775877 DOI: 10.1007/s11538-024-01308-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2023] [Accepted: 04/25/2024] [Indexed: 06/25/2024]
Abstract
Several recent theoretical studies have indicated that a relatively simple secretion control mechanism in the epithelial cells lining the stomach may be responsible for maintaining a neutral (healthy) pH adjacent to the stomach wall, even in the face of enormous electrodiffusive acid transport from the interior of the stomach. Subsequent work used Sobol' Indices (SIs) to quantify the degree to which this secretion mechanism is "self-regulating" i.e. the degree to which the wall pH is held neutral as mathematical parameters vary. However, questions remain regarding the nature of the control that specific parameters exert over the maintenance of a healthy stomach wall pH. Studying the sensitivity of higher moments of the statistical distribution of a model output can provide useful information, for example, how one parameter may skew the distribution towards or away from a physiologically advantageous regime. In this work, we prove a relationship between SIs and the higher moments and show how it can potentially reduce the cost of computing sensitivity of said moments. We define γ -indices to quantify sensitivity of variance, skewness, and kurtosis to the choice of value of a parameter, and we propose an efficient strategy that uses both SIs and γ -indices for a more comprehensive sensitivity analysis. Our analysis uncovers a control parameter which governs the "tightness of control" that the secretion mechanism exerts on wall pH. Finally, we discuss how uncertainty in this parameter can be reduced using expert information about higher moments, and speculate about the physiological advantage conferred by this control mechanism.
Collapse
Affiliation(s)
- Manu Aggarwal
- Department of Mathematics, Florida State University, Tallahassee, USA.
| | - Owen Lewis
- Department of Mathematics & Statistics, University of New Mexico, Albuquerque, New Mexico, USA
| | - Angela Jarrett
- Institute for Computational Engineering and Sciences, University of Texas at Austin, Austin, Texas, USA
| | - M Y Hussaini
- Department of Mathematics, Florida State University, Tallahassee, USA
| | - N G Cogan
- Department of Mathematics, Florida State University, Tallahassee, USA
| |
Collapse
|
3
|
Gasior KI, Cogan NG. Untangling the Molecular Interactions Underlying Intracellular Phase Separation Using Combined Global Sensitivity Analyses. Bull Math Biol 2024; 86:60. [PMID: 38641666 DOI: 10.1007/s11538-024-01288-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2023] [Accepted: 03/27/2024] [Indexed: 04/21/2024]
Abstract
Liquid-liquid phase separation is an intracellular mechanism by which molecules, usually proteins and RNAs, interact and then rapidly demix from the surrounding matrix to form membrane-less compartments necessary for cellular function. Occurring in both the cytoplasm and the nucleus, properties of the resulting droplets depend on a variety of characteristics specific to the molecules involved, such as valency, density, and diffusion within the crowded environment. Capturing these complexities in a biologically relevant model is difficult. To understand the nuanced dynamics between proteins and RNAs as they interact and form droplets, as well as the impact of these interactions on the resulting droplet properties, we turn to sensitivity analysis. In this work, we examine a previously published mathematical model of two RNA species competing for the same protein-binding partner. We use the combined analyses of Morris Method and Sobol' sensitivity analysis to understand the impact of nine molecular parameters, subjected to three different initial conditions, on two observable LLPS outputs: the time of phase separation and the composition of the droplet field. Morris Method is a screening method capable of highlighting the most important parameters impacting a given output, while the variance-based Sobol' analysis can quantify both the importance of a given parameter, as well as the other model parameters it interacts with, to produce the observed phenomena. Combining these two techniques allows Morris Method to identify the most important dynamics and circumvent the large computational expense associated with Sobol', which then provides more nuanced information about parameter relationships. Together, the results of these combined methodologies highlight the complicated protein-RNA relationships underlying both the time of phase separation and the composition of the droplet field. Sobol' sensitivity analysis reveals that observed spatial and temporal dynamics are due, at least in part, to high-level interactions between multiple (3+) parameters. Ultimately, this work discourages using a single measurement to extrapolate the value of any single rate or parameter value, while simultaneously establishing a framework in which to analyze and assess the impact of these small-scale molecular interactions on large-scale droplet properties.
Collapse
Affiliation(s)
- Kelsey I Gasior
- Department of Applied and Computational Mathematics and Statistics, University of Notre Dame, Notre Dame, IN, USA.
| | - Nicholas G Cogan
- Department of Mathematics, Florida State University, Tallahassee, FL, USA
| |
Collapse
|
4
|
Wanika L, Evans ND, Johnson M, Tomkinson H, Chappell MJ. In vitro PK/PD modeling of tyrosine kinase inhibitors in non-small cell lung cancer cell lines. Clin Transl Sci 2024; 17:e13714. [PMID: 38477045 PMCID: PMC10933606 DOI: 10.1111/cts.13714] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2023] [Revised: 11/22/2023] [Accepted: 12/14/2023] [Indexed: 03/14/2024] Open
Abstract
Tyrosine kinase inhibitors (TKIs) are routinely prescribed for the treatment of non-small cell lung cancer (NSCLC). As with all medications, patients can experience adverse events due to TKIs. Unfortunately, the relationship between many TKIs and the occurrence of certain adverse events remains unclear. There are limited in vivo studies which focus on TKIs and their effects on different regulation pathways. Many in vitro studies, however, that investigate the effects of TKIs observe additional changes, such as changes in gene activations or protein expressions. These studies could potentially help to gain greater understanding of the mechanisms for TKI induced adverse events. However, in order to utilize these pathways in a pharmacokinetic/pharmacodynamic (PK/PD) framework, an in vitro PK/PD model needs to be developed, in order to characterize the effects of TKIs in NSCLC cell lines. Through the use of ordinary differential equations, cell viability data and nonlinear mixed effects modeling, an in vitro TKI PK/PD model was developed with estimated PK and PD parameter values for the TKIs alectinib, crizotinib, erlotinib, and gefitinib. The relative standard errors for the population parameters are all less than 25%. The inclusion of random effects enabled the model to predict individual parameter values which provided a closer fit to the observed response. It is hoped that this model can be extended to include in vitro data of certain pathways that may potentially be linked with adverse events and provide a better understanding of TKI-induced adverse events.
Collapse
Affiliation(s)
- Linda Wanika
- School of EngineeringUniversity of WarwickCoventryUK
| | - Neil D. Evans
- School of EngineeringUniversity of WarwickCoventryUK
| | | | | | | |
Collapse
|
5
|
Lazarus A, Dalton D, Husmeier D, Gao H. Sensitivity analysis and inverse uncertainty quantification for the left ventricular passive mechanics. Biomech Model Mechanobiol 2022; 21:953-982. [PMID: 35377030 PMCID: PMC9132878 DOI: 10.1007/s10237-022-01571-8] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2021] [Accepted: 02/28/2022] [Indexed: 01/08/2023]
Abstract
Personalized computational cardiac models are considered to be a unique and powerful tool in modern cardiology, integrating the knowledge of physiology, pathology and fundamental laws of mechanics in one framework. They have the potential to improve risk prediction in cardiac patients and assist in the development of new treatments. However, in order to use these models for clinical decision support, it is important that both the impact of model parameter perturbations on the predicted quantities of interest as well as the uncertainty of parameter estimation are properly quantified, where the first task is a priori in nature (meaning independent of any specific clinical data), while the second task is carried out a posteriori (meaning after specific clinical data have been obtained). The present study addresses these challenges for a widely used constitutive law of passive myocardium (the Holzapfel-Ogden model), using global sensitivity analysis (SA) to address the first challenge, and inverse-uncertainty quantification (I-UQ) for the second challenge. The SA is carried out on a range of different input parameters to a left ventricle (LV) model, making use of computationally efficient Gaussian process (GP) surrogate models in place of the numerical forward simulator. The results of the SA are then used to inform a low-order reparametrization of the constitutive law for passive myocardium under consideration. The quality of this parameterization in the context of an inverse problem having observed noisy experimental data is then quantified with an I-UQ study, which again makes use of GP surrogate models. The I-UQ is carried out in a Bayesian manner using Markov Chain Monte Carlo, which allows for full uncertainty quantification of the material parameter estimates. Our study reveals insights into the relation between SA and I-UQ, elucidates the dependence of parameter sensitivity and estimation uncertainty on external factors, like LV cavity pressure, and sheds new light on cardio-mechanic model formulation, with particular focus on the Holzapfel-Ogden myocardial model.
Collapse
Affiliation(s)
- Alan Lazarus
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - David Dalton
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - Dirk Husmeier
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - Hao Gao
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| |
Collapse
|
6
|
Johnson KE, Howard GR, Morgan D, Brenner EA, Gardner AL, Durrett RE, Mo W, Al’Khafaji A, Sontag ED, Jarrett AM, Yankeelov TE, Brock A. Integrating transcriptomics and bulk time course data into a mathematical framework to describe and predict therapeutic resistance in cancer. Phys Biol 2020; 18:016001. [PMID: 33215611 PMCID: PMC8156495 DOI: 10.1088/1478-3975/abb09c] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
A significant challenge in the field of biomedicine is the development of methods to integrate the multitude of dispersed data sets into comprehensive frameworks to be used to generate optimal clinical decisions. Recent technological advances in single cell analysis allow for high-dimensional molecular characterization of cells and populations, but to date, few mathematical models have attempted to integrate measurements from the single cell scale with other types of longitudinal data. Here, we present a framework that actionizes static outputs from a machine learning model and leverages these as measurements of state variables in a dynamic model of treatment response. We apply this framework to breast cancer cells to integrate single cell transcriptomic data with longitudinal bulk cell population (bulk time course) data. We demonstrate that the explicit inclusion of the phenotypic composition estimate, derived from single cell RNA-sequencing data (scRNA-seq), improves accuracy in the prediction of new treatments with a concordance correlation coefficient (CCC) of 0.92 compared to a prediction accuracy of CCC = 0.64 when fitting on longitudinal bulk cell population data alone. To our knowledge, this is the first work that explicitly integrates single cell clonally-resolved transcriptome datasets with bulk time-course data to jointly calibrate a mathematical model of drug resistance dynamics. We anticipate this approach to be a first step that demonstrates the feasibility of incorporating multiple data types into mathematical models to develop optimized treatment regimens from data.
Collapse
Affiliation(s)
- Kaitlyn E Johnson
- Department of Biomedical Engineering, The University of
Texas at Austin, Austin, TX, 78712, United States of America
| | - Grant R Howard
- Department of Biomedical Engineering, The University of
Texas at Austin, Austin, TX, 78712, United States of America
| | - Daylin Morgan
- Department of Biomedical Engineering, The University of
Texas at Austin, Austin, TX, 78712, United States of America
| | - Eric A Brenner
- Department of Biomedical Engineering, The University of
Texas at Austin, Austin, TX, 78712, United States of America
- Institute for Cellular and Molecular Biology, The
University of Texas at Austin, Austin, TX, 78712, United States of America
| | - Andrea L Gardner
- Department of Biomedical Engineering, The University of
Texas at Austin, Austin, TX, 78712, United States of America
| | - Russell E Durrett
- Department of Biomedical Engineering, The University of
Texas at Austin, Austin, TX, 78712, United States of America
- Institute for Cellular and Molecular Biology, The
University of Texas at Austin, Austin, TX, 78712, United States of America
| | - William Mo
- Department of Biomedical Engineering, The University of
Texas at Austin, Austin, TX, 78712, United States of America
| | - Aziz Al’Khafaji
- Department of Biomedical Engineering, The University of
Texas at Austin, Austin, TX, 78712, United States of America
- Institute for Cellular and Molecular Biology, The
University of Texas at Austin, Austin, TX, 78712, United States of America
| | - Eduardo D Sontag
- Department of Electrical and Computer Engineering,
Northeastern University, Boston, MA, 02115, United States of America
- Department of Bioengineering, Northeastern University,
Boston, MA, 02115, United States of America
- Laboratory of Systems Pharmacology, Program in Therapeutics
Science, Harvard Medical School, Boston, MA, 02115, United States of America
| | - Angela M Jarrett
- Livestrong Cancer Institutes, Dell Medical School, The
University of Texas at Austin, Austin, TX, 78712, United States of America
- Oden Institute for Computational Engineering and Sciences,
The University of Texas at Austin
| | - Thomas E Yankeelov
- Department of Biomedical Engineering, The University of
Texas at Austin, Austin, TX, 78712, United States of America
- Livestrong Cancer Institutes, Dell Medical School, The
University of Texas at Austin, Austin, TX, 78712, United States of America
- Oden Institute for Computational Engineering and Sciences,
The University of Texas at Austin
- Department of Diagnostic Medicine, The University of Texas
at Austin, Austin, TX, 78712, United States of America
- Department of Oncology, The University of Texas at Austin,
Austin, TX, 78712, United States of America
- Department of Imaging Physics, The MD Anderson Cancer
Center Houston, TX, 77030, United States of America
| | - Amy Brock
- Department of Biomedical Engineering, The University of
Texas at Austin, Austin, TX, 78712, United States of America
- Institute for Cellular and Molecular Biology, The
University of Texas at Austin, Austin, TX, 78712, United States of America
- Livestrong Cancer Institutes, Dell Medical School, The
University of Texas at Austin, Austin, TX, 78712, United States of America
| |
Collapse
|
7
|
Hormuth DA, Jarrett AM, Lima EA, McKenna MT, Fuentes DT, Yankeelov TE. Mechanism-Based Modeling of Tumor Growth and Treatment Response Constrained by Multiparametric Imaging Data. JCO Clin Cancer Inform 2019; 3:1-10. [PMID: 30807209 PMCID: PMC6535803 DOI: 10.1200/cci.18.00055] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/05/2018] [Indexed: 12/19/2022] Open
Abstract
Multiparametric imaging is a critical tool in the noninvasive study and assessment of cancer. Imaging methods have evolved over the past several decades to provide quantitative measures of tumor and healthy tissue characteristics related to, for example, cell number, blood volume fraction, blood flow, hypoxia, and metabolism. Mechanistic models of tumor growth also have matured to a point where the incorporation of patient-specific measures could provide clinically relevant predictions of tumor growth and response. In this review, we identify and discuss approaches that use multiparametric imaging data, including diffusion-weighted magnetic resonance imaging, dynamic contrast-enhanced magnetic resonance imaging, diffusion tensor imaging, contrast-enhanced computed tomography, [18F]fluorodeoxyglucose positron emission tomography, and [18F]fluoromisonidazole positron emission tomography to initialize and calibrate mechanistic models of tumor growth and response. We focus the discussion on brain and breast cancers; however, we also identify three emerging areas of application in kidney, pancreatic, and lung cancers. We conclude with a discussion of the future directions for incorporating multiparametric imaging data and mechanistic modeling into clinical decision making for patients with cancer.
Collapse
|
8
|
Jarrett AM, Lima EABF, Hormuth DA, McKenna MT, Feng X, Ekrut DA, Resende ACM, Brock A, Yankeelov TE. Mathematical models of tumor cell proliferation: A review of the literature. Expert Rev Anticancer Ther 2018; 18:1271-1286. [PMID: 30252552 PMCID: PMC6295418 DOI: 10.1080/14737140.2018.1527689] [Citation(s) in RCA: 51] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
INTRODUCTION A defining hallmark of cancer is aberrant cell proliferation. Efforts to understand the generative properties of cancer cells span all biological scales: from genetic deviations and alterations of metabolic pathways to physical stresses due to overcrowding, as well as the effects of therapeutics and the immune system. While these factors have long been studied in the laboratory, mathematical and computational techniques are being increasingly applied to help understand and forecast tumor growth and treatment response. Advantages of mathematical modeling of proliferation include the ability to simulate and predict the spatiotemporal development of tumors across multiple experimental scales. Central to proliferation modeling is the incorporation of available biological data and validation with experimental data. Areas covered: We present an overview of past and current mathematical strategies directed at understanding tumor cell proliferation. We identify areas for mathematical development as motivated by available experimental and clinical evidence, with a particular emphasis on emerging, non-invasive imaging technologies. Expert commentary: The data required to legitimize mathematical models are often difficult or (currently) impossible to obtain. We suggest areas for further investigation to establish mathematical models that more effectively utilize available data to make informed predictions on tumor cell proliferation.
Collapse
Affiliation(s)
- Angela M Jarrett
- a Institute for Computational Engineering and Sciences , The University of Texas at Austin , Austin , USA
- b Livestrong Cancer Institutes , The University of Texas at Austin , Austin , USA
| | - Ernesto A B F Lima
- a Institute for Computational Engineering and Sciences , The University of Texas at Austin , Austin , USA
| | - David A Hormuth
- a Institute for Computational Engineering and Sciences , The University of Texas at Austin , Austin , USA
- b Livestrong Cancer Institutes , The University of Texas at Austin , Austin , USA
| | - Matthew T McKenna
- c Department of Biomedical Engineering , Vanderbilt University , Nashville , USA
| | - Xinzeng Feng
- a Institute for Computational Engineering and Sciences , The University of Texas at Austin , Austin , USA
| | - David A Ekrut
- a Institute for Computational Engineering and Sciences , The University of Texas at Austin , Austin , USA
| | - Anna Claudia M Resende
- a Institute for Computational Engineering and Sciences , The University of Texas at Austin , Austin , USA
- d Department of Computational Modeling , National Laboratory for Scientific Computing , Petrópolis , Brazil
| | - Amy Brock
- b Livestrong Cancer Institutes , The University of Texas at Austin , Austin , USA
- e Department of Biomedical Engineering , The University of Texas at Austin , Austin , USA
| | - Thomas E Yankeelov
- a Institute for Computational Engineering and Sciences , The University of Texas at Austin , Austin , USA
- b Livestrong Cancer Institutes , The University of Texas at Austin , Austin , USA
- e Department of Biomedical Engineering , The University of Texas at Austin , Austin , USA
- f Department of Diagnostic Medicine , The University of Texas at Austin , Austin , USA
- g Department of Oncology , The University of Texas at Austin , Austin , USA
| |
Collapse
|
9
|
Kostrytsia A, Papirio S, Mattei MR, Frunzo L, Lens PNL, Esposito G. Sensitivity analysis for an elemental sulfur-based two-step denitrification model. WATER SCIENCE AND TECHNOLOGY : A JOURNAL OF THE INTERNATIONAL ASSOCIATION ON WATER POLLUTION RESEARCH 2018; 78:1296-1303. [PMID: 30388086 DOI: 10.2166/wst.2018.398] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
A local sensitivity analysis was performed for a chemically synthesized elemental sulfur (S0)-based two-step denitrification model, accounting for nitrite (NO2 -) accumulation, biomass growth and S0 hydrolysis. The sensitivity analysis was aimed at verifying the model stability, understanding the model structure and individuating the model parameters to be further optimized. The mass specific area of the sulfur particles (a*) and hydrolysis kinetic constant (k1) were identified as the dominant parameters on the model outputs, i.e. nitrate (NO3 -), NO2 - and sulfate (SO4 2-) concentrations, confirming that the microbially catalyzed S0 hydrolysis is the rate-limiting step during S0-driven denitrification. Additionally, the maximum growth rates of the denitrifying biomass on NO3 - and NO2 - were detected as the most sensitive kinetic parameters.
Collapse
Affiliation(s)
- A Kostrytsia
- Department of Civil and Mechanical Engineering, University of Cassino and Southern Lazio, via Di Biasio 43, 03043 Cassino, FR, Italy E-mail:
| | - S Papirio
- Department of Civil, Architectural and Environmental Engineering, University of Naples Federico II, via Claudio 21, 80125 Naples, Italy
| | - M R Mattei
- Department of Mathematics and Applications Renato Caccioppoli, University of Naples Federico II, via Cintia, Monte S. Angelo, 1-80126 Naples, Italy
| | - L Frunzo
- Department of Mathematics and Applications Renato Caccioppoli, University of Naples Federico II, via Cintia, Monte S. Angelo, 1-80126 Naples, Italy
| | - P N L Lens
- UNESCO-IHE, Institute for Water Education, P.O. Box 3015, 2601 DA Delft, The Netherlands
| | - G Esposito
- Department of Civil and Mechanical Engineering, University of Cassino and Southern Lazio, via Di Biasio 43, 03043 Cassino, FR, Italy E-mail:
| |
Collapse
|
10
|
Combining Two Methods of Global Sensitivity Analysis to Investigate MRSA Nasal Carriage Model. Bull Math Biol 2017; 79:2258-2272. [PMID: 28752384 DOI: 10.1007/s11538-017-0329-7] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2017] [Accepted: 07/20/2017] [Indexed: 10/19/2022]
Abstract
We apply two different sensitivity techniques to a model of bacterial colonization of the anterior nares to better understand the dynamics of Staphylococcus aureus nasal carriage. Specifically, we use partial rank correlation coefficients to investigate sensitivity as a function of time and identify a reduced model with fewer than half of the parameters of the full model. The reduced model is used for the calculation of Sobol' indices to identify interacting parameters by their additional effects indices. Additionally, we found that the model captures an interesting characteristic of the biological phenomenon related to the initial population size of the infection; only two parameters had any significant additional effects, and these parameters have biological evidence suggesting they are connected but not yet completely understood. Sensitivity is often applied to elucidate model robustness, but we show that combining sensitivity measures can lead to synergistic insight into both model and biological structures.
Collapse
|
11
|
Anton C, Deng J, Wong YS, Zhang Y, Zhang W, Gabos S, Huang DY, Jin C. Modeling and simulation for toxicity assessment. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2017; 14:581-606. [PMID: 28092954 DOI: 10.3934/mbe.2017034] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
The effect of various toxicants on growth/death and morphology of human cells is investigated using the xCELLigence Real-Time Cell Analysis High Troughput in vitro assay. The cell index is measured as a proxy for the number of cells, and for each test substance in each cell line, time-dependent concentration response curves (TCRCs) are generated. In this paper we propose a mathematical model to study the effect of toxicants with various initial concentrations on the cell index. This model is based on the logistic equation and linear kinetics. We consider a three dimensional system of differential equations with variables corresponding to the cell index, the intracellular concentration of toxicant, and the extracellular concentration of toxicant. To efficiently estimate the model's parameters, we design an Expectation Maximization algorithm. The model is validated by showing that it accurately represents the information provided by the TCRCs recorded after the experiments. Using stability analysis and numerical simulations, we determine the lowest concentration of toxin that can kill the cells. This information can be used to better design experimental studies for cytotoxicity profiling assessment.
Collapse
Affiliation(s)
- Cristina Anton
- Department of Mathematics and Statistics, Grant MacEwan University, Edmonton, Alberta, T5P2P7, Canada .
| | | | | | | | | | | | | | | |
Collapse
|
12
|
Sensitivity Analysis of a Pharmacokinetic Model of Vaginal Anti-HIV Microbicide Drug Delivery. J Pharm Sci 2016; 105:1772-1778. [DOI: 10.1016/j.xphs.2016.02.015] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2015] [Revised: 01/30/2016] [Accepted: 02/02/2016] [Indexed: 12/14/2022]
|
13
|
McCarthy GD, Drewell RA, Dresch JM. Global sensitivity analysis of a dynamic model for gene expression in Drosophila embryos. PeerJ 2015; 3:e1022. [PMID: 26157608 PMCID: PMC4476099 DOI: 10.7717/peerj.1022] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2015] [Accepted: 05/25/2015] [Indexed: 11/20/2022] Open
Abstract
It is well known that gene regulation is a tightly controlled process in early organismal development. However, the roles of key processes involved in this regulation, such as transcription and translation, are less well understood, and mathematical modeling approaches in this field are still in their infancy. In recent studies, biologists have taken precise measurements of protein and mRNA abundance to determine the relative contributions of key factors involved in regulating protein levels in mammalian cells. We now approach this question from a mathematical modeling perspective. In this study, we use a simple dynamic mathematical model that incorporates terms representing transcription, translation, mRNA and protein decay, and diffusion in an early Drosophila embryo. We perform global sensitivity analyses on this model using various different initial conditions and spatial and temporal outputs. Our results indicate that transcription and translation are often the key parameters to determine protein abundance. This observation is in close agreement with the experimental results from mammalian cells for various initial conditions at particular time points, suggesting that a simple dynamic model can capture the qualitative behavior of a gene. Additionally, we find that parameter sensitivites are temporally dynamic, illustrating the importance of conducting a thorough global sensitivity analysis across multiple time points when analyzing mathematical models of gene regulation.
Collapse
Affiliation(s)
| | | | - Jacqueline M Dresch
- Department of Mathematics, Amherst College , Amherst, MA , USA ; Department of Mathematics and Computer Science, Clark University , Worcester, MA , USA
| |
Collapse
|