Mathematical models of cytotoxic effects in endpoint tumor cell line assays: critical assessment of the application of a single parametric value as a standard criterion to quantify the dose-response effects and new unexplored proposal formats.

The development of convenient tools for describing and quantifying the effects of standard and novel therapeutic agents is essential for the research community, to perform more precise evaluations. Although mathematical models and quantification criteria have been exchanged in the last decade between different fields of study, there are relevant methodologies that lack proper mathematical descriptions and standard criteria to quantify their responses. Therefore, part of the relevant information that can be drawn from the experimental results obtained and the quantification of its statistical reliability are lost. Despite its relevance, there is not a standard form for the in vitro endpoint tumor cell lines' assays (TCLA) that enables the evaluation of the cytotoxic dose-response effects of anti-tumor drugs. The analysis of all the specific problems associated with the diverse nature of the available TCLA used is unfeasible. However, since most TCLA share the main objectives and similar operative requirements, we have chosen the sulforhodamine B (SRB) colorimetric assay for cytotoxicity screening of tumor cell lines as an experimental case study. In this work, the common biological and practical non-linear dose-response mathematical models are tested against experimental data and, following several statistical analyses, the model based on the Weibull distribution was confirmed as the convenient approximation to test the cytotoxic effectiveness of anti-tumor compounds. Then, the advantages and disadvantages of all the different parametric criteria derived from the model, which enable the quantification of the dose-response drug-effects, are extensively discussed. Therefore, model and standard criteria for easily performing the comparisons between different compounds are established. The advantages include a simple application, provision of parametric estimations that characterize the response as standard criteria, economization of experimental effort and enabling rigorous comparisons among the effects of different compounds and experimental approaches. In all experimental data fitted, the calculated parameters were always statistically significant, the equations proved to be consistent and the correlation coefficient of determination was, in most of the cases, higher than 0.98.

[1]  C. Winsor,et al.  The Gompertz Curve as a Growth Curve. , 1932, Proceedings of the National Academy of Sciences of the United States of America.

[2]  Mohammad Fallahi-Sichani,et al.  Metrics other than potency reveal systematic variation in responses to cancer drugs. , 2013, Nature chemical biology.

[3]  J. Gibbs Mechanism-based target identification and drug discovery in cancer research. , 2000, Science.

[4]  M. Murado,et al.  The notion of hormesis and the dose-response theory: a unified approach. , 2007, Journal of theoretical biology.

[5]  Brian A. Nosek,et al.  An open investigation of the reproducibility of cancer biology research , 2014, eLife.

[6]  M. Prieto,et al.  Hydrolysis optimization of mannan, curdlan and cell walls from Endomyces fibuliger grown in mussel processing wastewaters , 2011 .

[7]  W. Weibull A statistical theory of the strength of materials , 1939 .

[8]  P. Rosin The Laws Governing the Fineness of Powdered Coal , 1933 .

[9]  M. Prieto,et al.  A critical point: the problems associated with the variety of criteria to quantify the antioxidant capacity. , 2014, Journal of agricultural and food chemistry.

[10]  A. Tsoularis,et al.  Analysis of logistic growth models. , 2002, Mathematical biosciences.

[11]  P. Verhulst Recherches mathématiques sur la loi d’accroissement de la population , 2022, Nouveaux mémoires de l'Académie royale des sciences et belles-lettres de Bruxelles.

[12]  M. Forster,et al.  Key Concepts in Model Selection: Performance and Generalizability. , 2000, Journal of mathematical psychology.

[13]  R. Goody,et al.  The original Michaelis constant: translation of the 1913 Michaelis-Menten paper. , 2011, Biochemistry.

[14]  David L. Wilson,et al.  The analysis of survival (mortality) data: Fitting Gompertz, Weibull, and logistic functions , 1994, Mechanisms of Ageing and Development.

[15]  P. Sorger,et al.  Growth rate inhibition metrics correct for confounders in measuring sensitivity to cancer drugs , 2016, Nature Methods.

[16]  D. Rodbard,et al.  Simultaneous analysis of families of sigmoidal curves: application to bioassay, radioligand assay, and physiological dose-response curves. , 1978, The American journal of physiology.

[17]  Ana Maria Carvalho,et al.  Chemical composition of wild and commercial Achillea millefolium L. and bioactivity of the methanolic extract, infusion and decoction. , 2013, Food chemistry.

[18]  Tyler Sa,et al.  Dynamics of normal growth. , 1965 .

[19]  M. Delignette-Muller,et al.  Estimating the bacterial lag time: which model, which precision? , 2004, International journal of food microbiology.

[20]  Joshua A. Bittker,et al.  Correlating chemical sensitivity and basal gene expression reveals mechanism of action , 2015, Nature chemical biology.

[21]  A. Hill The Combinations of Haemoglobin with Oxygen and with Carbon Monoxide. I. , 1913, The Biochemical journal.

[22]  R. D. Berger Comparison of the Gompertz and Logistic Equations to Describe Plant Disease Progress , 1981 .

[23]  J. Lorenzo,et al.  Evaluation of non-linear equations to model different animal growths with mono and bisigmoid profiles. , 2012, Journal of theoretical biology.

[24]  P. Polese,et al.  SOLVERSTAT: a new utility for multipurpose analysis. An application to the investigation of dioxygenated Co(II) complex formation in dimethylsulfoxide solution. , 2003, Talanta.

[25]  A. Anderson,et al.  Evolution of cell motility in an individual-based model of tumour growth. , 2009, Journal of theoretical biology.

[26]  Smirnov Ip The organizing function of medical technology in public health , 1972 .

[27]  José Antonio Vázquez,et al.  Dose-response relationships: an overview, a generative model and its application to the verification of descriptive models , 2002 .

[28]  F. Rombouts,et al.  Modeling of the Bacterial Growth Curve , 1990, Applied and environmental microbiology.

[29]  Marcello Farina,et al.  Modular model of TNFα cytotoxicity , 2011, Bioinform..

[30]  Francisco Taboada,et al.  Comparative study of four sigmoid models of pressure-volume curve in acute lung injury , 2007, Biomedical engineering online.

[31]  Martin A. Nowak,et al.  A spatial model predicts that dispersal and cell turnover limit intratumour heterogeneity , 2015, Nature.

[32]  Chih-Ling Tsai,et al.  Regression model selection—a residual likelihood approach , 2002 .

[33]  Benjamin Gompertz,et al.  XXIV. On the nature of the function expressive of the law of human mortality, and on a new mode of determining the value of life contingencies. In a letter to Francis Baily, Esq. F. R. S. &c , 1825, Philosophical Transactions of the Royal Society of London.

[34]  S. Keller,et al.  Nonlinear least-squares data fitting in Excel spreadsheets , 2010, Nature Protocols.

[35]  M. Prieto,et al.  Comparison of several mathematical models for describing the joint effect of temperature and ph on glucanex activity , 2012, Biotechnology progress.

[36]  M. Prieto,et al.  A Time-Dose Model to Quantify the Antioxidant Responses of the Oxidative Hemolysis Inhibition Assay (OxHLIA) and Its Extension to Evaluate Other Hemolytic Effectors , 2014, BioMed research international.

[37]  Krister Wennerberg,et al.  Quantitative scoring of differential drug sensitivity for individually optimized anticancer therapies , 2014, Scientific Reports.

[38]  F. J. Richards A Flexible Growth Function for Empirical Use , 1959 .

[39]  Miguel Ángel Prieto Lage,et al.  Dose-response analysis in the joint action of two effectors. A new approach to simulation, identification and modelling of some basic interactions. , 2013 .

[40]  D. Banerjee,et al.  Cytotoxicity and Cell Growth Assays , 2006 .

[41]  Ana Maria Carvalho,et al.  Bioactivity and chemical characterization in hydrophilic and lipophilic compounds of Chenopodium ambrosioides L , 2013 .

[42]  R. Rust,et al.  Model selection criteria: an investigation of relative accuracy, posterior probabilities, and combinations of criteria , 1995 .

[43]  Dejian Huang,et al.  The chemistry behind antioxidant capacity assays. , 2005, Journal of agricultural and food chemistry.

[44]  E. Tjørve Shapes and functions of species–area curves: a review of possible models , 2003 .

[45]  Carlos F. Lopez,et al.  An unbiased metric of antiproliferative drug effect in vitro , 2016, Nature Methods.

[46]  W. Weibull A Statistical Distribution Function of Wide Applicability , 1951 .

[47]  R. C. Whiting,et al.  When is simple good enough: a comparison of the Gompertz, Baranyi, and three-phase linear models for fitting bacterial growth curves , 1997 .

[48]  M. Prieto,et al.  Preparation of marine silage of swordfish, ray and shark visceral waste by lactic acid bacteria , 2011 .

[49]  Omer Dushek,et al.  Phenotypic models of T cell activation , 2014, Nature Reviews Immunology.

[50]  Emmanuel Barillot,et al.  Mathematical Modelling of Molecular Pathways Enabling Tumour Cell Invasion and Migration , 2015, PLoS Comput. Biol..

[51]  R. Shoemaker The NCI60 human tumour cell line anticancer drug screen , 2006, Nature Reviews Cancer.

[52]  Kanyawim Kirtikara,et al.  Sulforhodamine B colorimetric assay for cytotoxicity screening , 2006, Nature Protocols.

[53]  J. Black,et al.  Operational models of pharmacological agonism , 1983, Proceedings of the Royal Society of London. Series B. Biological Sciences.

[54]  M. Prieto,et al.  Oversimplification and Overstandardization in Biological Methods: Sperm Bioassays in Ecotoxicology as a Case of Study and a Proposal for Their Reformulation , 2014, TheScientificWorldJournal.

[55]  M. Prieto,et al.  NOEC and LOEC as merely concessive expedients: two unambiguous alternatives and some criteria to maximize the efficiency of dose-response experimental designs. , 2013, The Science of the total environment.

[56]  P. Clemons,et al.  Target identification and mechanism of action in chemical biology and drug discovery. , 2013, Nature chemical biology.

[57]  P. Zage,et al.  Measuring cytotoxicity: a new perspective on LC50. , 2007, Anticancer research.

[58]  P. Campbell,et al.  Enzyme changes in neonatal skeletal muscle: effect of thyroid deficiency. , 1978, The American journal of physiology.

[59]  M. Maurin,et al.  REVIEW ARTICLE doi: 10.1111/j.1472-8206.2008.00633.x The Hill equation: a review of its capabilities in pharmacological modelling , 2008 .

[60]  M. Fréchet Sur la loi de probabilité de l'écart maximum , 1928 .

[61]  Joshua M. Stuart,et al.  Subtype and pathway specific responses to anticancer compounds in breast cancer , 2011, Proceedings of the National Academy of Sciences.

[62]  D. Rivers,et al.  Model Selection Tests for Nonlinear Dynamic Models , 2002 .

[63]  Britta Basse,et al.  Modelling cell population growth with applications to cancer therapy in human tumour cell lines. , 2004, Progress in biophysics and molecular biology.

[64]  I. Puri,et al.  Mathematical model for chemotherapeutic drug efficacy in arresting tumour growth based on the cancer stem cell hypothesis , 2007, Cell proliferation.

[65]  Mustafa Özilgen,et al.  Kinetic Model of Lipid Oxidation in Foods , 1990 .

[66]  Benjamin Haibe-Kains,et al.  Inconsistency in large pharmacogenomic studies , 2013, Nature.

[67]  Philip Hall,et al.  Integrated Stochastic Model of DNA Damage Repair by Non-homologous End Joining and p53/p21- Mediated Early Senescence Signalling , 2015, PLoS Comput. Biol..

[68]  Adam A. Margolin,et al.  The Cancer Cell Line Encyclopedia enables predictive modeling of anticancer drug sensitivity , 2012, Nature.

[69]  G. Peters,et al.  Comparison of the sulforhodamine B protein and tetrazolium (MTT) assays for in vitro chemosensitivity testing. , 1991, European journal of cancer.

[70]  Erik Sahai,et al.  Differing modes of tumour cell invasion have distinct requirements for Rho/ROCK signalling and extracellular proteolysis , 2003, Nature Cell Biology.