An automated fitting procedure and software for dose-response curves with multiphasic features

In cancer pharmacology (and many other areas), most dose-response curves are satisfactorily described by a classical Hill equation (i.e. 4 parameters logistical). Nevertheless, there are instances where the marked presence of more than one point of inflection, or the presence of combined agonist and antagonist effects, prevents straight-forward modelling of the data via a standard Hill equation. Here we propose a modified model and automated fitting procedure to describe dose-response curves with multiphasic features. The resulting general model enables interpreting each phase of the dose-response as an independent dose-dependent process. We developed an algorithm which automatically generates and ranks dose-response models with varying degrees of multiphasic features. The algorithm was implemented in new freely available Dr Fit software (sourceforge.net/projects/drfit/). We show how our approach is successful in describing dose-response curves with multiphasic features. Additionally, we analysed a large cancer cell viability screen involving 11650 dose-response curves. Based on our algorithm, we found that 28% of cases were better described by a multiphasic model than by the Hill model. We thus provide a robust approach to fit dose-response curves with various degrees of complexity, which, together with the provided software implementation, should enable a wide audience to easily process their own data.

[1]  Henggui Zhang,et al.  hERG Inhibitors with Similar Potency But Different Binding Kinetics Do Not Pose the Same Proarrhythmic Risk: Implications for Drug Safety Assessment , 2014, Journal of cardiovascular electrophysiology.

[2]  David Allard,et al.  Inhibition of Hedgehog Signaling Enhances Delivery of Chemotherapy in a Mouse Model of Pancreatic Cancer , 2009, Science.

[3]  Shu-Shen Liu,et al.  Predicting hormetic effects of ionic liquid mixtures on luciferase activity using the concentration addition model. , 2011, Environmental science & technology.

[4]  Laurence Roy,et al.  CABAS: a freely available PC program for fitting calibration curves in chromosome aberration dosimetry. , 2007, Radiation protection dosimetry.

[5]  Heino Prinz,et al.  Hill coefficients, dose–response curves and allosteric mechanisms , 2010, Journal of chemical biology.

[6]  J. Satti.  The Emerging Low-dose Therapy for Advanced Cancers , 2009, Dose-response : a publication of International Hormesis Society.

[7]  Xiaotong Shen,et al.  Adaptive Model Selection , 2002 .

[8]  Kunihiko Takahashi,et al.  Cubic Spline Regression of J-shaped Dose-Response Curves with Likelihood-based Assignments of Grouped Exposure Levels , 2013 .

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

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

[11]  Dale Bowman,et al.  A Parametric Model for Detecting Hormetic Effects in Developmental Toxicity Studies , 2004, Risk analysis : an official publication of the Society for Risk Analysis.

[12]  P. Sorger,et al.  Systems biology and combination therapy in the quest for clinical efficacy , 2006, Nature chemical biology.

[13]  Ruili Huang,et al.  A Grid Algorithm for High Throughput Fitting of Dose-Response Curve Data , 2010, Current chemical genomics.

[14]  R. Conolly,et al.  Nonmonotonic dose-response relationships: mechanistic basis, kinetic modeling, and implications for risk assessment. , 2004, Toxicological sciences : an official journal of the Society of Toxicology.

[15]  S. Davis,et al.  Pharmacological curve fitting to analyze cutaneous adrenergic responses. , 2011, Journal of applied physiology.

[16]  E. Ward.  A review and comparison of four commonly used Bayesian and maximum likelihood model selection tools , 2008 .

[17]  C. I. Bliss.  THE TOXICITY OF POISONS APPLIED JOINTLY1 , 1939 .

[18]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[19]  U. Nielsen,et al.  Preclinical activity of nanoliposomal irinotecan is governed by tumor deposition and intratumor prodrug conversion. , 2014, Cancer research.

[20]  Directed Walk Designs for Dose‐Response Problems with Competing Failure Modes , 2003, Biometrics.

[21]  R. Weinberg,et al.  The Biology of Cancer , 2006 .

[22]  Geoffrey A Head,et al.  A five-parameter logistic equation for investigating asymmetry of curvature in baroreflex studies. , 1999, American journal of physiology. Regulatory, integrative and comparative physiology.

[23]  R. Liedl,et al.  Complexity vs. Simplicity: Groundwater Model Ranking Using Information Criteria , 2014, Ground water.

[24]  E. Calabrese,et al.  Reevaluation of the Fundamental Dose–Response Relationship , 1999 .

[25]  A. Reynolds,et al.  Potential Relevance of Bell-Shaped and U-Shaped Dose-Responses for the Therapeutic Targeting of Angiogenesis in Cancer , 2010, Dose-response : a publication of International Hormesis Society.

[26]  Hua Liang,et al.  A Strategy to Model Nonmonotonic Dose-Response Curve and Estimate IC50 , 2013, PloS one.

[27]  Christian Ritz,et al.  Toward a unified approach to dose–response modeling in ecotoxicology , 2010, Environmental toxicology and chemistry.

[28]  Peter Dalgaard,et al.  R Development Core Team (2010): R: A language and environment for statistical computing , 2010 .

[29]  W. I L L I A,et al.  A General Approach to Modeling Biphasic Relationships , 2008 .

[30]  Sudhindra R Gadagkar,et al.  Computational tools for fitting the Hill equation to dose-response curves. , 2015, Journal of pharmacological and toxicological methods.

[31]  D. Svendsgaard,et al.  U-shaped dose-response curves: their occurrence and implications for risk assessment. , 1990, Journal of toxicology and environmental health.

[32]  Shu-Shen Liu,et al.  Modeling non-monotonic dose-response relationships: model evaluation and hormetic quantities exploration. , 2013, Ecotoxicology and environmental safety.

[33]  J. Timbrell.  Principles of biochemical toxicology , 1991 .

[34]  J. Ramsay.  Monotone Regression Splines in Action , 1988 .

[35]  H. Akaike.  A new look at the statistical model identification , 1974 .

[36]  J. Hostettler.  Introduction to the , 1983 .

[37]  C. Stevens,et al.  INTRODUCTION TO PHARMACOLOGY , 1965, The Ulster Medical Journal.

[38]  J. Giraldo,et al.  Assessing the (a)symmetry of concentration-effect curves: empirical versus mechanistic models. , 2002, Pharmacology & therapeutics.

[39]  P. Bain,et al.  Cytotoxicity of binary mixtures of human pharmaceuticals in a fish cell line: approaches for non-monotonic concentration-response relationships. , 2014, Chemosphere.

[40]  John W. A. Findlay,et al.  Appropriate calibration curve fitting in ligand binding assays , 2007, The AAPS Journal.

[41]  A. Juška.  A minimal model of non-hyperbolic enzyme and receptor kinetics. , 2003, Biochemical and biophysical research communications.

[42]  G. Hecht,et al.  Myosin regulation of NKCC1: effects on cAMP-mediated Cl- secretion in intestinal epithelia. , 1999, American journal of physiology. Cell physiology.

[43]  Jonathan D Parsons,et al.  A high-throughput method for fitting dose-response curves using Microsoft Excel. , 2007, Analytical biochemistry.

[44]  D. Zheleva,et al.  Paclitaxel and CYC3, an aurora kinase A inhibitor, synergise in pancreatic cancer cells but not bone marrow precursor cells , 2012, British Journal of Cancer.

[45]  Shu-Shen Liu,et al.  Support vector regression and least squares support vector regression for hormetic dose-response curves fitting. , 2010, Chemosphere.

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

[47]  G. Schwarz.  Estimating the Dimension of a Model , 1978 .

[48]  Christian Ritz,et al.  Improved empirical models describing hormesis , 2005, Environmental toxicology and chemistry.

[49]  Statistical Dose-Response Models with Hormetic Effects , 2001 .

[50]  S. Ramaswamy,et al.  Systematic identification of genomic markers of drug sensitivity in cancer cells , 2012, Nature.

[51]  J. Giraldo,et al.  On the fitting of binding data when receptor dimerization is suspected , 2008, British journal of pharmacology.

[52]  Roger D. Cousens,et al.  An equation to describe dose responses where there is stimulation of growth at low doses , 1989 .

[53]  R. Prentice,et al.  A generalization of the probit and logit methods for dose response curves. , 1976, Biometrics.

[54]  Nicole Genco,et al.  Evaluation of an asymmetry parameter for curve-fitting in single-chemical and mixture toxicity assessment. , 2012, Toxicology.

[55]  Davis Jm,et al.  U‐Shaped dose‐response curves: Their occurrence and implications for risk assessment , 1990 .

[56]  W. Greco,et al.  The search for synergy: a critical review from a response surface perspective. , 1995, Pharmacological reviews.

[57]  Piero Quatto,et al.  Flexible meta-regression functions for modeling aggregate dose-response data, with an application to alcohol and mortality. , 2004, American journal of epidemiology.