A Computational Systems Biology Software Platform for Multiscale Modeling and Simulation: Integrating Whole-Body Physiology, Disease Biology, and Molecular Reaction Networks

Today, in silico studies and trial simulations already complement experimental approaches in pharmaceutical R&D and have become indispensable tools for decision making and communication with regulatory agencies. While biology is multiscale by nature, project work, and software tools usually focus on isolated aspects of drug action, such as pharmacokinetics at the organism scale or pharmacodynamic interaction on the molecular level. We present a modeling and simulation software platform consisting of PK-Sim® and MoBi® capable of building and simulating models that integrate across biological scales. A prototypical multiscale model for the progression of a pancreatic tumor and its response to pharmacotherapy is constructed and virtual patients are treated with a prodrug activated by hepatic metabolization. Tumor growth is driven by signal transduction leading to cell cycle transition and proliferation. Free tumor concentrations of the active metabolite inhibit Raf kinase in the signaling cascade and thereby cell cycle progression. In a virtual clinical study, the individual therapeutic outcome of the chemotherapeutic intervention is simulated for a large population with heterogeneous genomic background. Thereby, the platform allows efficient model building and integration of biological knowledge and prior data from all biological scales. Experimental in vitro model systems can be linked with observations in animal experiments and clinical trials. The interplay between patients, diseases, and drugs and topics with high clinical relevance such as the role of pharmacogenomics, drug–drug, or drug–metabolite interactions can be addressed using this mechanistic, insight driven multiscale modeling approach.

[1]  L. Kuepfer Towards whole-body systems physiology , 2010, Molecular systems biology.

[2]  W. Schmitt,et al.  A Physiologic Model for Simulating Gastrointestinal Flow and Drug Absorption in Rats , 2003, Pharmaceutical Research.

[3]  James P. Freyer,et al.  Tumor growthin vivo and as multicellular spheroids compared by mathematical models , 1994, Bulletin of mathematical biology.

[4]  Ivan Nestorov,et al.  Whole-body physiologically based pharmacokinetic models , 2007, Expert opinion on drug metabolism & toxicology.

[5]  Benveniste,et al.  Cytochrome P450 , 1993, Handbook of Experimental Pharmacology.

[6]  J Kirchheiner,et al.  The clinical role of genetic polymorphisms in drug-metabolizing enzymes , 2008, The Pharmacogenomics Journal.

[7]  M. Rowland,et al.  Physiologically based pharmacokinetic modelling 2: predicting the tissue distribution of acids, very weak bases, neutrals and zwitterions. , 2006, Journal of pharmaceutical sciences.

[8]  H M Sauro,et al.  SCAMP: a general-purpose simulator and metabolic control analysis program , 1993, Comput. Appl. Biosci..

[9]  Dominique Tytgat,et al.  Physiologically based pharmacokinetics (PBPK) , 2009, Drug metabolism reviews.

[10]  E Schöll,et al.  Classifying the expansion kinetics and critical surface dynamics of growing cell populations. , 2007, Physical review letters.

[11]  Norah Lynn Henry,et al.  Cytochrome P450 2D6 activity predicts discontinuation of tamoxifen therapy in breast cancer patients , 2009, The Pharmacogenomics Journal.

[12]  Stefan Willmann,et al.  Theoretical Biology and Medical Modelling Open Access Dynamically Simulating the Interaction of Midazolam and the Cyp3a4 Inhibitor Itraconazole Using Individual Coupled Whole-body Physiologically-based Pharmacokinetic (wb-pbpk) Models , 2022 .

[13]  Jeffrey W. Smith,et al.  Stochastic Gene Expression in a Single Cell , .

[14]  H. Greenspan On the growth and stability of cell cultures and solid tumors. , 1976, Journal of theoretical biology.

[15]  Reinhart Heinrich,et al.  A linear steady-state treatment of enzymatic chains. A mathematical model of glycolysis of human erythrocytes. , 1974, European journal of biochemistry.

[16]  P. Sorger,et al.  Non-genetic origins of cell-to-cell variability in TRAIL-induced apoptosis , 2009, Nature.

[17]  Zhike Zi,et al.  SBML-SAT: a systems biology markup language (SBML) based sensitivity analysis tool , 2008, BMC Bioinformatics.

[18]  Stefan Willmann,et al.  Development and Validation of a Physiology-based Model for the Prediction of Oral Absorption in Monkeys , 2007, Pharmaceutical Research.

[19]  Patrick Poulin,et al.  Prediction of pharmacokinetics prior to in vivo studies. II. Generic physiologically based pharmacokinetic models of drug disposition. , 2002, Journal of pharmaceutical sciences.

[20]  Walter Schmitt,et al.  Application of physiology-based pharmacokinetic and pharmacodynamic modeling to individualized target-controlled propofol infusions , 2006, Advances in therapy.

[21]  Reinhart Heinrich,et al.  A linear steady-state treatment of enzymatic chains. General properties, control and effector strength. , 1974, European journal of biochemistry.

[22]  Malcolm Rowland,et al.  Mechanistic Approaches to Volume of Distribution Predictions: Understanding the Processes , 2007, Pharmaceutical Research.

[23]  Kannan Krishnan,et al.  Physiologically-based pharmacokinetic and toxicokinetic models in cancer risk assessment. , 2005, Journal of environmental science and health. Part C, Environmental carcinogenesis & ecotoxicology reviews.

[24]  E. Kunkel Systems biology in drug discovery , 2004, Nature Biotechnology.

[25]  M Schwab,et al.  Comprehensive analysis of the genetic factors determining expression and function of hepatic CYP2D6. , 2001, Pharmacogenetics.

[26]  Walter Schmitt,et al.  PK-Sim®: a physiologically based pharmacokinetic ‘whole-body’ model , 2003 .

[27]  Hiroaki Kitano,et al.  The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models , 2003, Bioinform..

[28]  R. Heinrich,et al.  The Regulation of Cellular Systems , 1996, Springer US.

[29]  Walter Schmitt,et al.  General approach for the calculation of tissue to plasma partition coefficients. , 2008, Toxicology in vitro : an international journal published in association with BIBRA.

[30]  Jörg Lippert,et al.  From physicochemistry to absorption and distribution: predictive mechanistic modelling and computational tools , 2005, Expert opinion on drug metabolism & toxicology.

[31]  D. Fell,et al.  Differential feedback regulation of the MAPK cascade underlies the quantitative differences in EGF and NGF signalling in PC12 cells , 2000, FEBS letters.

[32]  T. G. Coleman,et al.  Systems biology and integrative physiological modelling , 2011, The Journal of physiology.

[33]  Stefan Willmann,et al.  Whole-body physiologically based pharmacokinetic population modelling of oral drug administration: inter-individual variability of cimetidine absorption. , 2009, The Journal of pharmacy and pharmacology.

[34]  F. Theil,et al.  Prediction of adipose tissue: plasma partition coefficients for structurally unrelated drugs. , 2001, Journal of pharmaceutical sciences.

[35]  I. Sobol,et al.  Global sensitivity indices for nonlinear mathematical models. Review , 2005 .

[36]  Walter Schmitt,et al.  Whole body physiologically-based pharmacokinetic models: their use in clinical drug development , 2008, Expert opinion on drug metabolism & toxicology.

[37]  Malcolm Rowland,et al.  Tissue distribution of basic drugs: accounting for enantiomeric, compound and regional differences amongst beta-blocking drugs in rat. , 2005, Journal of pharmaceutical sciences.

[38]  W. Kalow,et al.  Pharmacogenetics and pharmacogenomics: origin, status, and the hope for personalized medicine , 2006, The Pharmacogenomics Journal.

[39]  H. Kacser,et al.  The control of flux. , 1995, Biochemical Society transactions.

[40]  Walter Schmitt,et al.  Development of a Physiology-Based Whole-Body Population Model for Assessing the Influence of Individual Variability on the Pharmacokinetics of Drugs , 2007, Journal of Pharmacokinetics and Pharmacodynamics.

[41]  L. Berezhkovskiy,et al.  Volume of distribution at steady state for a linear pharmacokinetic system with peripheral elimination. , 2004, Journal of pharmaceutical sciences.

[42]  Walter Schmitt,et al.  A Mechanistic Approach for the Scaling of Clearance in Children , 2006, Clinical pharmacokinetics.

[43]  Stefan Willmann,et al.  Physiology-Based Simulations of a Pathological Condition , 2008, Clinical pharmacokinetics.

[44]  Edda Klipp,et al.  Systems Biology , 1994 .

[45]  Walter Schmitt,et al.  A physiological model for the estimation of the fraction dose absorbed in humans. , 2004, Journal of medicinal chemistry.

[46]  Patrick Poulin,et al.  Prediction of pharmacokinetics prior to in vivo studies. 1. Mechanism-based prediction of volume of distribution. , 2002, Journal of pharmaceutical sciences.

[47]  Walter Schmitt,et al.  Development and Evaluation of a Generic Physiologically Based Pharmacokinetic Model for Children , 2006, Clinical pharmacokinetics.

[48]  A. Brú,et al.  The universal dynamics of tumor growth. , 2003, Biophysical journal.

[49]  A. Hindmarsh,et al.  CVODE, a stiff/nonstiff ODE solver in C , 1996 .

[50]  Henrik Madsen,et al.  Non-linear mixed-effects pharmacokinetic/pharmacodynamic modelling in NLME using differential equations , 2004, Comput. Methods Programs Biomed..

[51]  Jukka Corander,et al.  CYP2D6 worldwide genetic variation shows high frequency of altered activity variants and no continental structure , 2007, Pharmacogenetics and genomics.

[52]  Susanne Briest,et al.  Tamoxifen metabolism and its effect on endocrine treatment of breast cancer. , 2009, Clinical advances in hematology & oncology : H&O.

[53]  C. Rao,et al.  Control, exploitation and tolerance of intracellular noise , 2002, Nature.

[54]  S. Willmann,et al.  Risk to the Breast‐Fed Neonate From Codeine Treatment to the Mother: A Quantitative Mechanistic Modeling Study , 2009, Clinical pharmacology and therapeutics.

[55]  L. Shampine Solving ODEs and DDEs with residual control , 2005 .

[56]  D L S McElwain,et al.  A history of the study of solid tumour growth: The contribution of mathematical modelling , 2004, Bulletin of mathematical biology.

[57]  Uri Alon,et al.  Dynamics and variability of ERK2 response to EGF in individual living cells. , 2009, Molecular cell.

[58]  Edward C Stites,et al.  Network Analysis of Oncogenic Ras Activation in Cancer , 2007, Science.

[59]  R. Milo,et al.  Variability and memory of protein levels in human cells , 2006, Nature.

[60]  I. Sobol Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates , 2001 .

[61]  M. Rowland,et al.  Physiologically based pharmacokinetic modeling 1: predicting the tissue distribution of moderate-to-strong bases. , 2005, Journal of pharmaceutical sciences.

[62]  George Loizou,et al.  Development of good modelling practice for physiologically based pharmacokinetic models for use in risk assessment: the first steps. , 2008, Regulatory toxicology and pharmacology : RTP.

[63]  E Schöll,et al.  Comparing the growth kinetics of cell populations in two and three dimensions. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[64]  Thomas S Deisboeck,et al.  Personalizing medicine: a systems biology perspective , 2009, Molecular systems biology.