
-
Previous Article
Flow optimization in vascular networks
- MBE Home
- This Issue
- Next Article
Modeling and simulation for toxicity assessment
1. | Department of Mathematics and Statistics, Grant MacEwan University, Edmonton, Alberta, T5P2P7, Canada |
2. | Department of Mathematical and statistical Sciences, University of Alberta, Edmonton, Alberta, T6G2G1, Canada |
3. | Alberta Health, Edmonton, Alberta, T5J1S6, Canada |
4. | Department of Laboratory Medicine and Pathology, University of Alberta, Edmonton, Alberta, T6G2B7, Canada |
5. | Alberta Centre for Toxicology, University of Calgary, Calgary, Alberta, T2N4N1, Canada |
6. | ACEA Biosciences Inc, San Diego, California, 92121, USA |
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.
References:
[1] |
C. Biernacki, G. Celeux and G. Govaert,
Choosing starting values for the EM algorithm for getting the highest likelihood in multivariate gaussian mixture models, Comput. Statist. Data Anal., 41 (2003), 561-575.
doi: 10.1016/S0167-9473(02)00163-9. |
[2] |
F. Cannavó,
Sensitivity analysis for volcanic source modeling quality assessment and model selection, Computers & Geosciences,, 44 (2012), 52-59.
|
[3] |
Z. Ghahramani and S. Roweis, Learning nonlinear dynamical systems using an EM algorithm, in Advances in Neural Information Processing Systems (eds. M. Kearns, S. Solla and C. D. A.), MIT Press, 1999,599-605. |
[4] |
T. Hallam, C. Clark and G. Jordan,
Effects of toxicants on populations: A qualitative approach Ⅱ. First order kinetics, J. Math. Biology, 18 (1983), 25-37.
|
[5] |
J. He and K. Wang,
The survival analysis for a population in a polluted environment, Nonlinear Analysis: Real World Applications, 10 (2009), 1555-1571.
doi: 10.1016/j.nonrwa.2008.01.027. |
[6] |
B. Huang and J. Xing,
Dynamic modeling and prediction of cytotoxicity on microelectronic cell sensor array, The Canadian Journal of Chemical Engineering, 84 (2006), 393-405.
|
[7] |
Q. Huang, L. Parshotam, H. Wang, C. Bampfylde and M. Lewis,
A model for the impact of contaminants on fish population dynamics, Journal of Theoretical Biology, 334 (2013), 71-79.
doi: 10.1016/j.jtbi.2013.05.018. |
[8] |
F. Ibrahim, B. Huang, J. Xing and S. Gabos,
Early determination of toxicant concentration in water supply using MHE, Water Research, 44 (2010), 3252-3260.
|
[9] |
A.M. Jarrett, Y. Liu, N. Cogan and M.Y. Hussaini,
Global sensitivity analysis used to interpret biological experimental results, Journal of Mathematical Biology, 71 (2015), 151-170.
doi: 10.1007/s00285-014-0818-3. |
[10] |
J. Jiao, W. Long and L. Chen,
A single stage-structured population model with mature individuals in a polluted environment and pulse input of environmental toxin, Nonlinear Analysis: Real World Applications, 10 (2009), 3073-3081.
doi: 10.1016/j.nonrwa.2008.10.007. |
[11] |
S. Julier, J. Uhlmann and H. Durrant-White,
A new method for the nonlinear transformation of means and covariances in filters and estimators, IEEE Trans. Aut. Control, 45 (2000), 477-482.
doi: 10.1109/9.847726. |
[12] |
S. Julier, J. Uhlmann and H. Durrant-Whyte, A new approach for filtering nonlinear systems, in American Control Conference, Seattle, Washington, 1995,1628–1632.
doi: 10.1109/ACC.1995.529783. |
[13] |
A. Kiparissides, S. Kucherenko, A. Mantalaris and E.N. Pistikopoulos,
Global sensitivity analysis challenges in biological systems modeling, Industrial & Engineering Chemistry Research, 48 (2009), 7168-7180.
doi: 10.1021/ie900139x. |
[14] |
K. Kothawad, A. Pathan and M. Logad,
Evaluation of in vitro anti-cancer activity of fruit lagenaria siceraria against MCF7, HOP62 and DU145 cell line, Int. J. Pharm. & Technol, 4 (2012), 3909-4392.
|
[15] |
M. Liu and K. Wang,
Survival analysis of stochastic single-species population models in polluted environments, Ecological Modeling, 220 (2009), 1347-1357.
|
[16] |
M. Liu, K. Wang and X. Liu,
Long term behaviors of stochastic single-species growth models in a polluted environment, Applied Mathematical Modelling, 35 (2011), 752-762.
doi: 10.1016/j.apm.2010.07.031. |
[17] |
X. Meng and D. Van Dyk,
The EM algorithm -an old folk-song to a fast new tune, J.R. Statist. Soc.B, 59 (1997), 511-567.
doi: 10.1111/1467-9868.00082. |
[18] |
R. Neal and G. Hinton, A view of the EM algorithm that justifies incremental, sparse, an other variants, in Learning in Graphical Models (ed. M. Jordan), 89 (1998), 355-368.
doi: 10.1007/978-94-011-5014-9_12. |
[19] |
T. Pan, B. Huang, W. Zhang, S. Gabos, D. Huang and V. Devendran,
Cytotoxicity assessment based on the AUC50 using multi-concentration time-dependent cellular response curves, Anal. Chim. Acta, 764 (2013), 44-52.
|
[20] |
T. Pan, S. Khare, F. Ackah, B. Huang, W. Zhang, S. Gabos, C. Jin and M. Stampfl,
In vitro cytotoxicity assessment based on KC50 with real-time cell analyzer (RTCA) assay, Comp. Biol. Chem., 47 (2013), 113-120.
|
[21] |
L. Perko, Differential Equations and Dynamical Systems, Springer, New York, 2001.
doi: 10.1007/978-1-4613-0003-8. |
[22] |
R. Shumway and D. Stoffer,
An approach to time series smoothing and forecasting using the EM algorithm, J. Time Ser. Anal., 3 (1982), 253-264.
doi: 10.1111/j.1467-9892.1982.tb00349.x. |
[23] |
I.M. Sobol,
Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates, Mathematics and Computers in Simulation, 55 (2001), 271-280.
doi: 10.1016/S0378-4754(00)00270-6. |
[24] |
H. Thieme, Mathematics in Population Biology, Princeton Series in theoretical and Computational Biology., 2003 |
[25] |
E. A. Wan, R. Van der Merwe and A. T. Nelson, Dual estimation and the unscented transformation, in Advances in Neural Information Processing Systems (ed. M. I. J. et al.), MIT Press, 2000. |
[26] |
C. Wu,
On the convergence properties of the EM algorithm, The Annals of Statistics, 11 (1983), 95-103.
doi: 10.1214/aos/1176346060. |
[27] |
Z. Xi, S. Khare, A. Cheung, B. Huang, T. Pan, W. Zhang, F. Ibrahim, C. Jin and S. Gabos,
Mode of action classification of chemicals using multi-concentration time-dependent cellular response profiles, Comp. Biol. Chem., 49 (2014), 23-35.
doi: 10.1016/j.compbiolchem.2013.12.004. |
[28] |
J. Xing, L. Zhu, S. Gabos and L. Xie,
Microelectronic cell sensor assay for detection of cytotoxicity and prediction of acute toxicity, Toxicology in Vitro, 20 (2006), 995-1004.
doi: 10.1016/j.tiv.2005.12.008. |
[29] |
M. Zhang, D. Aguilera, C. Das, H. Vasquez, P. Zage, V. Gopalakrishnan and J. Wolff,
Measuring cytotoxicity: A new perspective on LC50, Anticancer Res., 27 (2007), 35-38.
|
show all references
References:
[1] |
C. Biernacki, G. Celeux and G. Govaert,
Choosing starting values for the EM algorithm for getting the highest likelihood in multivariate gaussian mixture models, Comput. Statist. Data Anal., 41 (2003), 561-575.
doi: 10.1016/S0167-9473(02)00163-9. |
[2] |
F. Cannavó,
Sensitivity analysis for volcanic source modeling quality assessment and model selection, Computers & Geosciences,, 44 (2012), 52-59.
|
[3] |
Z. Ghahramani and S. Roweis, Learning nonlinear dynamical systems using an EM algorithm, in Advances in Neural Information Processing Systems (eds. M. Kearns, S. Solla and C. D. A.), MIT Press, 1999,599-605. |
[4] |
T. Hallam, C. Clark and G. Jordan,
Effects of toxicants on populations: A qualitative approach Ⅱ. First order kinetics, J. Math. Biology, 18 (1983), 25-37.
|
[5] |
J. He and K. Wang,
The survival analysis for a population in a polluted environment, Nonlinear Analysis: Real World Applications, 10 (2009), 1555-1571.
doi: 10.1016/j.nonrwa.2008.01.027. |
[6] |
B. Huang and J. Xing,
Dynamic modeling and prediction of cytotoxicity on microelectronic cell sensor array, The Canadian Journal of Chemical Engineering, 84 (2006), 393-405.
|
[7] |
Q. Huang, L. Parshotam, H. Wang, C. Bampfylde and M. Lewis,
A model for the impact of contaminants on fish population dynamics, Journal of Theoretical Biology, 334 (2013), 71-79.
doi: 10.1016/j.jtbi.2013.05.018. |
[8] |
F. Ibrahim, B. Huang, J. Xing and S. Gabos,
Early determination of toxicant concentration in water supply using MHE, Water Research, 44 (2010), 3252-3260.
|
[9] |
A.M. Jarrett, Y. Liu, N. Cogan and M.Y. Hussaini,
Global sensitivity analysis used to interpret biological experimental results, Journal of Mathematical Biology, 71 (2015), 151-170.
doi: 10.1007/s00285-014-0818-3. |
[10] |
J. Jiao, W. Long and L. Chen,
A single stage-structured population model with mature individuals in a polluted environment and pulse input of environmental toxin, Nonlinear Analysis: Real World Applications, 10 (2009), 3073-3081.
doi: 10.1016/j.nonrwa.2008.10.007. |
[11] |
S. Julier, J. Uhlmann and H. Durrant-White,
A new method for the nonlinear transformation of means and covariances in filters and estimators, IEEE Trans. Aut. Control, 45 (2000), 477-482.
doi: 10.1109/9.847726. |
[12] |
S. Julier, J. Uhlmann and H. Durrant-Whyte, A new approach for filtering nonlinear systems, in American Control Conference, Seattle, Washington, 1995,1628–1632.
doi: 10.1109/ACC.1995.529783. |
[13] |
A. Kiparissides, S. Kucherenko, A. Mantalaris and E.N. Pistikopoulos,
Global sensitivity analysis challenges in biological systems modeling, Industrial & Engineering Chemistry Research, 48 (2009), 7168-7180.
doi: 10.1021/ie900139x. |
[14] |
K. Kothawad, A. Pathan and M. Logad,
Evaluation of in vitro anti-cancer activity of fruit lagenaria siceraria against MCF7, HOP62 and DU145 cell line, Int. J. Pharm. & Technol, 4 (2012), 3909-4392.
|
[15] |
M. Liu and K. Wang,
Survival analysis of stochastic single-species population models in polluted environments, Ecological Modeling, 220 (2009), 1347-1357.
|
[16] |
M. Liu, K. Wang and X. Liu,
Long term behaviors of stochastic single-species growth models in a polluted environment, Applied Mathematical Modelling, 35 (2011), 752-762.
doi: 10.1016/j.apm.2010.07.031. |
[17] |
X. Meng and D. Van Dyk,
The EM algorithm -an old folk-song to a fast new tune, J.R. Statist. Soc.B, 59 (1997), 511-567.
doi: 10.1111/1467-9868.00082. |
[18] |
R. Neal and G. Hinton, A view of the EM algorithm that justifies incremental, sparse, an other variants, in Learning in Graphical Models (ed. M. Jordan), 89 (1998), 355-368.
doi: 10.1007/978-94-011-5014-9_12. |
[19] |
T. Pan, B. Huang, W. Zhang, S. Gabos, D. Huang and V. Devendran,
Cytotoxicity assessment based on the AUC50 using multi-concentration time-dependent cellular response curves, Anal. Chim. Acta, 764 (2013), 44-52.
|
[20] |
T. Pan, S. Khare, F. Ackah, B. Huang, W. Zhang, S. Gabos, C. Jin and M. Stampfl,
In vitro cytotoxicity assessment based on KC50 with real-time cell analyzer (RTCA) assay, Comp. Biol. Chem., 47 (2013), 113-120.
|
[21] |
L. Perko, Differential Equations and Dynamical Systems, Springer, New York, 2001.
doi: 10.1007/978-1-4613-0003-8. |
[22] |
R. Shumway and D. Stoffer,
An approach to time series smoothing and forecasting using the EM algorithm, J. Time Ser. Anal., 3 (1982), 253-264.
doi: 10.1111/j.1467-9892.1982.tb00349.x. |
[23] |
I.M. Sobol,
Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates, Mathematics and Computers in Simulation, 55 (2001), 271-280.
doi: 10.1016/S0378-4754(00)00270-6. |
[24] |
H. Thieme, Mathematics in Population Biology, Princeton Series in theoretical and Computational Biology., 2003 |
[25] |
E. A. Wan, R. Van der Merwe and A. T. Nelson, Dual estimation and the unscented transformation, in Advances in Neural Information Processing Systems (ed. M. I. J. et al.), MIT Press, 2000. |
[26] |
C. Wu,
On the convergence properties of the EM algorithm, The Annals of Statistics, 11 (1983), 95-103.
doi: 10.1214/aos/1176346060. |
[27] |
Z. Xi, S. Khare, A. Cheung, B. Huang, T. Pan, W. Zhang, F. Ibrahim, C. Jin and S. Gabos,
Mode of action classification of chemicals using multi-concentration time-dependent cellular response profiles, Comp. Biol. Chem., 49 (2014), 23-35.
doi: 10.1016/j.compbiolchem.2013.12.004. |
[28] |
J. Xing, L. Zhu, S. Gabos and L. Xie,
Microelectronic cell sensor assay for detection of cytotoxicity and prediction of acute toxicity, Toxicology in Vitro, 20 (2006), 995-1004.
doi: 10.1016/j.tiv.2005.12.008. |
[29] |
M. Zhang, D. Aguilera, C. Das, H. Vasquez, P. Zage, V. Gopalakrishnan and J. Wolff,
Measuring cytotoxicity: A new perspective on LC50, Anticancer Res., 27 (2007), 35-38.
|














Symbol | Definition |
cell index ≈ cell population | |
toxicant concentration inside the cell | |
toxicant concentration outside the cell | |
cell growth rate in the absence of toxicant | |
capacity volume | |
effect coefficient of toxicant on the cell's growth | |
the uptake rate of the toxicant from environment | |
the toxicant uptake rate from cells | |
the toxicant input rate to the environment | |
the losses rate of toxicant absorbed by cells |
Symbol | Definition |
cell index ≈ cell population | |
toxicant concentration inside the cell | |
toxicant concentration outside the cell | |
cell growth rate in the absence of toxicant | |
capacity volume | |
effect coefficient of toxicant on the cell's growth | |
the uptake rate of the toxicant from environment | |
the toxicant uptake rate from cells | |
the toxicant input rate to the environment | |
the losses rate of toxicant absorbed by cells |
Initialize the model parameters |
Repeat until the log likelihood has converged |
The E step |
For k=1 to N |
Run the UF filter to compute |
For k=N to 1 |
Calculate the smoothed values |
The M step |
Update the values of the parameters |
Initialize the model parameters |
Repeat until the log likelihood has converged |
The E step |
For k=1 to N |
Run the UF filter to compute |
For k=N to 1 |
Calculate the smoothed values |
The M step |
Update the values of the parameters |
Toxicant | Cluster | β | K | |||||
PF431396 | Ⅹ | 0.077 | 21.912 | 0.273 | 0.058 | 0 | 0.008 | 0.238 |
monastrol | Ⅹ | 0.074 | 18.17 | 0.209 | 0.177 | 0.204 | 0.5 | 0.016 |
ABT888 | Ⅰ | 0.083 | 17.543 | 0.079 | 0.177 | 0.205 | 0.5 | 0.005 |
HA1100 hydrochloride | Ⅰ | 0.077 | 21.913 | 0.143 | 0.0098 | 0.0786 | 0.147 | 0.351 |
Toxicant | Cluster | β | K | |||||
PF431396 | Ⅹ | 0.077 | 21.912 | 0.273 | 0.058 | 0 | 0.008 | 0.238 |
monastrol | Ⅹ | 0.074 | 18.17 | 0.209 | 0.177 | 0.204 | 0.5 | 0.016 |
ABT888 | Ⅰ | 0.083 | 17.543 | 0.079 | 0.177 | 0.205 | 0.5 | 0.005 |
HA1100 hydrochloride | Ⅰ | 0.077 | 21.913 | 0.143 | 0.0098 | 0.0786 | 0.147 | 0.351 |
[1] |
Azmy S. Ackleh, Jeremy J. Thibodeaux. Parameter estimation in a structured erythropoiesis model. Mathematical Biosciences & Engineering, 2008, 5 (4) : 601-616. doi: 10.3934/mbe.2008.5.601 |
[2] |
Alex Capaldi, Samuel Behrend, Benjamin Berman, Jason Smith, Justin Wright, Alun L. Lloyd. Parameter estimation and uncertainty quantification for an epidemic model. Mathematical Biosciences & Engineering, 2012, 9 (3) : 553-576. doi: 10.3934/mbe.2012.9.553 |
[3] |
Marcello Delitala, Tommaso Lorenzi. A mathematical model for value estimation with public information and herding. Kinetic and Related Models, 2014, 7 (1) : 29-44. doi: 10.3934/krm.2014.7.29 |
[4] |
Gianni Gilioli, Sara Pasquali, Fabrizio Ruggeri. Nonlinear functional response parameter estimation in a stochastic predator-prey model. Mathematical Biosciences & Engineering, 2012, 9 (1) : 75-96. doi: 10.3934/mbe.2012.9.75 |
[5] |
Sebastian Springer, Heikki Haario, Vladimir Shemyakin, Leonid Kalachev, Denis Shchepakin. Robust parameter estimation of chaotic systems. Inverse Problems and Imaging, 2019, 13 (6) : 1189-1212. doi: 10.3934/ipi.2019053 |
[6] |
Simon Hubmer, Andreas Neubauer, Ronny Ramlau, Henning U. Voss. On the parameter estimation problem of magnetic resonance advection imaging. Inverse Problems and Imaging, 2018, 12 (1) : 175-204. doi: 10.3934/ipi.2018007 |
[7] |
Robert Azencott, Yutheeka Gadhyan. Accurate parameter estimation for coupled stochastic dynamics. Conference Publications, 2009, 2009 (Special) : 44-53. doi: 10.3934/proc.2009.2009.44 |
[8] |
Francisco de la Hoz, Anna Doubova, Fernando Vadillo. Persistence-time estimation for some stochastic SIS epidemic models. Discrete and Continuous Dynamical Systems - B, 2015, 20 (9) : 2933-2947. doi: 10.3934/dcdsb.2015.20.2933 |
[9] |
Blaise Faugeras, Olivier Maury. An advection-diffusion-reaction size-structured fish population dynamics model combined with a statistical parameter estimation procedure: Application to the Indian Ocean skipjack tuna fishery. Mathematical Biosciences & Engineering, 2005, 2 (4) : 719-741. doi: 10.3934/mbe.2005.2.719 |
[10] |
Azmy S. Ackleh, H.T. Banks, Keng Deng, Shuhua Hu. Parameter Estimation in a Coupled System of Nonlinear Size-Structured Populations. Mathematical Biosciences & Engineering, 2005, 2 (2) : 289-315. doi: 10.3934/mbe.2005.2.289 |
[11] |
Chongyang Liu, Meijia Han, Zhaohua Gong, Kok Lay Teo. Robust parameter estimation for constrained time-delay systems with inexact measurements. Journal of Industrial and Management Optimization, 2021, 17 (1) : 317-337. doi: 10.3934/jimo.2019113 |
[12] |
Krzysztof Fujarewicz, Krzysztof Łakomiec. Parameter estimation of systems with delays via structural sensitivity analysis. Discrete and Continuous Dynamical Systems - B, 2014, 19 (8) : 2521-2533. doi: 10.3934/dcdsb.2014.19.2521 |
[13] |
Dominique Chapelle, Philippe Moireau, Patrick Le Tallec. Robust filtering for joint state-parameter estimation in distributed mechanical systems. Discrete and Continuous Dynamical Systems, 2009, 23 (1&2) : 65-84. doi: 10.3934/dcds.2009.23.65 |
[14] |
Andrea Arnold, Daniela Calvetti, Erkki Somersalo. Vectorized and parallel particle filter SMC parameter estimation for stiff ODEs. Conference Publications, 2015, 2015 (special) : 75-84. doi: 10.3934/proc.2015.0075 |
[15] |
Alessandro Corbetta, Adrian Muntean, Kiamars Vafayi. Parameter estimation of social forces in pedestrian dynamics models via a probabilistic method. Mathematical Biosciences & Engineering, 2015, 12 (2) : 337-356. doi: 10.3934/mbe.2015.12.337 |
[16] |
Ferenc Hartung. Parameter estimation by quasilinearization in differential equations with state-dependent delays. Discrete and Continuous Dynamical Systems - B, 2013, 18 (6) : 1611-1631. doi: 10.3934/dcdsb.2013.18.1611 |
[17] |
Jiangqi Wu, Linjie Wen, Jinglai Li. Resampled ensemble Kalman inversion for Bayesian parameter estimation with sequential data. Discrete and Continuous Dynamical Systems - S, 2022, 15 (4) : 837-850. doi: 10.3934/dcdss.2021045 |
[18] |
Suqi Ma. Low viral persistence of an immunological model. Mathematical Biosciences & Engineering, 2012, 9 (4) : 809-817. doi: 10.3934/mbe.2012.9.809 |
[19] |
Abraão D. C. Nascimento, Leandro C. Rêgo, Raphaela L. B. A. Nascimento. Compound truncated Poisson normal distribution: Mathematical properties and Moment estimation. Inverse Problems and Imaging, 2019, 13 (4) : 787-803. doi: 10.3934/ipi.2019036 |
[20] |
Laura Martín-Fernández, Gianni Gilioli, Ettore Lanzarone, Joaquín Míguez, Sara Pasquali, Fabrizio Ruggeri, Diego P. Ruiz. A Rao-Blackwellized particle filter for joint parameter estimation and biomass tracking in a stochastic predator-prey system. Mathematical Biosciences & Engineering, 2014, 11 (3) : 573-597. doi: 10.3934/mbe.2014.11.573 |
2018 Impact Factor: 1.313
Tools
Metrics
Other articles
by authors
[Back to Top]