The primary goal of drug developers is to establish efficient and effective therapeutic protocols. Multifactorial pathologies, including dynamical diseases and complex disorders, can be difficult to treat, given the high degree of inter- and intra-patient variability and nonlinear physiological relationships. Quantitative approaches combining mechanistic disease modeling and computational strategies are increasingly leveraged to rationalize pre-clinical and clinical studies and to establish effective treatment strategies. The development of clinical trials has led to new computational methods that allow for large clinical data sets to be combined with pharmacokinetic and pharmacodynamic models of diseases. Here, we discuss recent progress using in silico clinical trials to explore treatments for a variety of complex diseases, ultimately demonstrating the immense utility of quantitative methods in drug development and medicine.
REFERENCES
1
R. M.
May
, “Dynamical diseases
,” Nature
272
, 673
–674
(1978
). 2
L.
Glass
, M. C.
Mackey
, and L.
Glass
, “Oscillation and chaos in physiological control systems
,” Science
197
, 287
–289
(1977
). 3
P.
Sorger
et al., “Quantitative and systems pharmacology in the post-genomic era: new approaches to discovering drugs and understanding therapeutic mechanisms
,” in An NIH White Paper by the QSP Workshop Group—October 2011. An NIH White Paper by the QSP Workshop Group
(National Institutes of Health
, Bethesda, MD
, 2011
).4
J. E.
Sager
, J.
Yu
, I.
Ragueneau-Majlessi
, and N.
Isoherranen
, “Physiologically based pharmacokinetic (PBPK) modeling and simulation approaches: A systematic review of published models, applications, and model verification
,” Drug Metab. Dispos.
43
, 1823
–1837
(2015
). 5
A. N.
Edginton
, F.-P.
Theil
, W.
Schmitt
, and S.
Willmann
, “Whole body physiologically-based pharmacokinetic models: Their use in clinical drug development
,” Expert Opin. Drug Metab. Toxicol.
4
, 1143
–1152
(2008
). 6
I. P.
Androulakis
, “Quantitative systems pharmacology: A framework for context
,” Curr. Pharmacol. Rep.
2
, 152
–160
(2016
). 7
M.
Craig
, “Towards quantitative systems pharmacology models of chemotherapy-induced neutropenia
,” CPT Pharmacometrics Syst. Pharmacol.
6
, 293
–304
(2017
). 8
J.
Zhang
, J. J.
Cunningham
, J. S.
Brown
, and R. A.
Gatenby
, “Integrating evolutionary dynamics into treatment of metastatic castrate-resistant prostate cancer
,” Nat. Commun.
8
, 1816
(2017
). 9
A.
Beuter
, L.
Glass
, M. C.
Mackey
, and M. S.
Titcombe
, Nonlinear Dynamics in Physiology and Medicine
(Springer
, 2003
).10
L.
Glass
and M. C.
Mackey
, “Pathological conditions resulting from instabilities in physiological control systems
,” Ann. N. Y. Acad. Sci.
316
, 214
–235
(1979
).11
C.
Foley
and M. C.
Mackey
, “Dynamic hematological disease: A review
,” J. Math. Biol.
58
, 285
–322
(2009
). 12
L.
Glass
and M. C.
Mackey
, “Pathological conditions resulting from instabilities in physiological control systems
,” Ann. N. Y. Acad. Sci.
316
, 214
–235
(1979
). 13
M.
Craig
et al, “Neutrophil dynamics during concurrent chemotherapy and G-CSF administration: Mathematical modelling guides dose optimisation to minimise neutropenia
,” J. Theor. Biol.
385
, 77
–89
(2015
). 14
G.
Brooks
et al, “Neutrophil dynamics after chemotherapy and G-CSF: The role of pharmacokinetics in shaping the response
,” J. Theor. Biol.
315
, 97
–109
(2012
). 15
I.
Badralexi
and A. A
Halanay
, “Complex model for blood cells’ evolution in chronic myelogenous leukemia
,” in 2015 20th International Conference on Control Systems and Computer Science
(IEEE, 2015
), pp. 611
–617
.16
A. V.
Hill
, “The possible effects of the aggregation of the molecules of haemoglobin on its dissociation curves
,” J. Physiol.
40
(suppl
), i
–vii
(1910
).17
M.
Santillán
, “On the use of the Hill functions in mathematical models of gene regulatory networks
,” Math. Model. Nat. Phenom.
3
, 85
–97
(2008
). 18
J.
Sebaugh
, “Guidelines for accurate EC50/IC50 estimation
,” Pharm. Stat.
10
, 128
–134
(2011
). 19
J.
Macdougall
, “Analysis of dose–response studies—E max model
,” in Dose Finding in Drug Development
(Springer, 2006
), pp. 127
–145
.20
F.
Vuletić
, V.
Zajec
, L. B.
Vuletić
, and S.
Seiwerth
, “Intratumoral heterogeneity
,” Diagn. Pathol.
4
, 2018DO
(2018
). 21
J.
Rosenfeld
and M. J.
Strong
, “Challenges in the understanding and treatment of amyotrophic lateral sclerosis/motor neuron disease
,” Neurotherapeutics
12
, 317
–325
(2015
). 22
F.
Boem
, I.
Malagrinò
, and M.
Bertolaso
, “In silico clinical trials: A possible response to complexity in pharmacology
,” in Uncertainty in Pharmacology
(Springer
, 2020
), pp. 135
–152
.23
B.
Ribba
et al, “A model of vascular tumour growth in mice combining longitudinal tumour size data with histological biomarkers
,” Eur. J. Cancer
47
, 479
–490
(2011
). 24
P.
Girard
, “Clinical trial simulation: A tool for understanding study failures and preventing them
,” Basic Clin. Pharmacol. Toxicol.
96
, 228
–234
(2005
). 25
M.
Jamei
, G. L.
Dickinson
, and A.
Rostami-Hodjegan
, “A framework for assessing inter-individual variability in pharmacokinetics using virtual human populations and integrating general knowledge of physical chemistry, biology, anatomy, physiology and genetics: A tale of ‘bottom-up’ vs ‘top-down’ recognition
,” Drug Metab. Pharmacokinet.
24
, 53
–75
(2009
). 26
Y.
Wang
et al, “Leveraging prior quantitative knowledge to guide drug development decisions and regulatory science recommendations: Impact of FDA pharmacometrics during 2004-2006
,” J. Clin. Pharmacol.
48
, 146
–156
(2008
). 27
T.
Cassidy
and M.
Craig
, “Determinants of combination GM-CSF immunotherapy and oncolytic virotherapy success identified through in silico treatment personalization
,” PLOS Comput. Biol.
15
, e1007495
(2019
). 28
V. M.
Pérez-García
et al, “Computational design of improved standardized chemotherapy protocols for grade II oligodendrogliomas
,” PLoS Comput. Biol.
15
, 1
–17
(2019
). 29
S.
Barish
, M. F.
Ochs
, E. D.
Sontag
, and J. L.
Gevertz
, “Evaluating optimal therapy robustness by virtual expansion of a sample population, with a case study in cancer immunotherapy
,” Proc. Natl. Acad. Sci. U.S.A.
114
, E6277
–E6286
(2017
). 30
P. A. J.
Valitalo
et al, “Novel model-based dosing guidelines for gentamicin and tobramycin in preterm and term neonates
,” J. Antimicrob. Chemother.
70
, 2074
–2077
(2015
). 31
D. J.
Klinke
, “Integrating epidemiological data into a mechanistic model of type 2 diabetes: Validating the prevalence of virtual patients
,” Ann. Biomed. Eng.
36
, 321
–334
(2008
). 32
A.
Kansal
and J.
Trimmer
, “Application of predictive biosimulation within pharmaceutical clinical development: Examples of significance for translational medicine and clinical trial design
,” IEE Proc. Syst. Biol.
152
, 214
–220
(2005
). 33
Z.
Barter
et al, “Determination of a quantitative relationship between hepatic CYP3A5*1/*3 and CYP3A4 expression for use in the prediction of metabolic clearance in virtual populations
,” Biopharm. Drug Dispos.
31
, 516
–532
(2010
). 34
B. J.
Schmidt
, F. P.
Casey
, T.
Paterson
, and J. R.
Chan
, “Alternate virtual populations elucidate the type I interferon signature predictive of the response to rituximab in rheumatoid arthritis
,” BMC Bioinformatics
14
, 221
(2013
). 35
D.
Brown
et al, “Trauma in silico: Individual-specific mathematical models and virtual clinical populations
,” Sci. Transl. Med.
285
, 285ra61
(2015
).36
D.
Teutonico
et al, “Generating virtual patients by multivariate and discrete re-sampling techniques
,” Pharm. Res.
32
, 3228
–3237
(2015
). 37
M.
Viceconti
, A.
Henney
, and E.
Morley-Fletcher
, “In silico clinical trials: How computer simulation will transform the biomedical industry
,” Int. J. Clin. Trials
3
, 37
–46
(2016
). 38
R. N.
Upton
and D. R.
Mould
, “Basic concepts in population modeling, simulation, and model-based drug development: Part 3-introduction to pharmacodynamic modeling methods
,” CPT Pharmacometrics Syst. Pharmacol.
3
, e88
(2014
). 39
I. P.
Androulakis
, “Systems engineering meets quantitative systems pharmacology: From low-level targets to engaging the host defenses
,” WIREs Syst. Biol. Med.
7
, 101
–112
(2015
). 40
O.
Milberg
et al, “A QSP model for predicting clinical responses to monotherapy, combination and sequential therapy following CTLA-4, PD-1, and PD-L1 checkpoint blockade
,” Sci. Rep.
9
, 11286
(2019
). 41
H.
Wang
et al, “In silico simulation of a clinical trial with anti-CTLA-4 and anti-PD-L1 immunotherapies in metastatic breast cancer using a systems pharmacology model
,” R. Soc. Open Sci.
6
, 190366
(2020
). 42
R. J.
Allen
, T. R.
Rieger
, and C. J.
Musante
, “Efficient generation and selection of virtual populations in quantitative systems pharmacology models
,” CPT Pharmacometrics Syst. Pharmacol.
5
, 140
–146
(2016
). 43
M. C.
Mackey
et al., “The timing of cyclic cytotoxic chemotherapy can worsen neutropenia and neutrophilia
,” Br. J. Clin. Pharmacol.
(published online 2020
). 44
I.
Dagogo-Jack
and A. T.
Shaw
, “Tumour heterogeneity and resistance to cancer therapies
,” Nat. Rev. Clin. Oncol.
15
, 81
–94
(2018
). 45
M.
Craig
et al, “Cooperative adaptation to therapy (CAT) confers resistance in heterogeneous non-small cell lung cancer
,” PLOS Comput. Biol.
15
, 1
–19
(2019
). 46
I.
Bozic
et al, “Accumulation of driver and passenger mutations during tumor progression
,” Proc. Natl. Acad. Sci. U.S.A.
107
, 18545
–18550
(2010
). 47
I.
Tirosh
et al, “Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq
,” Science
352
, 189
–196
(2016
). 48
E.
Kozłowska
, T.
Vallius
, J.
Hynninen
, S.
Hietanen
, and A.
Färkkilä
, “Virtual clinical trials identify effective combination therapies in ovarian cancer
,” Sci. Rep.
9
, 1
–9
(2019
). 49
E.
Kozłowska et al
., “Mathematical modeling predicts response to chemotherapy and drug combinations in ovarian cancer
,” Cancer Res.
78
, 4036
–4044
(2018
). 50
K.
Han
et al, “Simulations to predict clinical trial outcome of bevacizumab plus chemotherapy vs. chemotherapy alone in patients with first-line gastric cancer and elevated plasma VEGF-A
,” CPT Pharmacometrics Syst. Pharmacol.
5
, 352
–358
(2016
). 51
Y. Y.
Lee
and P.
Van Tassel
, “Intracranial oligodendrogliomas: Imaging findings in 35 untreated cases
,” Am. J. Roentgenol.
152
, 361
–369
(1989
). 52
D. N.
Louis
et al, “The 2016 World Health Organization classification of tumors of the central nervous system: A summary
,” Acta Neuropathol.
131
, 803
–820
(2016
). 53
C.
Conill
, E.
Verger
, and M.
Salamero
, “Performance status assessment in cancer patients
,” Cancer
65
, 1864
–1866
(1990
). 54
T. C.
Ryken
, I.
Parney
, J.
Buatti
, S. N.
Kalkanis
, and J. J.
Olson
, “The role of radiotherapy in the management of patients with diffuse low grade glioma: A systematic review and evidence-based clinical practice guideline
,” J. Neuro-Oncol.
125
, 551
–583
(2015
). 55
R.
Stupp
et al, “Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma
,” N. Engl. J. Med.
352
, 987
–996
(2005
). 56
B. J.
Fisher
et al, “Phase 2 study of temozolomide-based chemoradiation therapy for high-risk low-grade gliomas: Preliminary results of radiation therapy oncology group 0424
,” Int. J. Radiat. Oncol. Biol. Phys.
91
, 497
–504
(2015
). 57
S.
Kesari
et al, “Phase II study of protracted daily temozolomide for low-grade gliomas in adults
,” Clin. Cancer Res.
15
, 330
–337
(2009
). 58
G.
Kaloshi
et al, “Temozolomide for low-grade gliomas: Predictive impact of 1p/19q loss on response and outcome
,” Neurology
68
, 1831
–1836
(2007
). 59
L. B.
Kennedy
and A. K.
Salama
, “A review of cancer immunotherapy toxicity
,” CA. Cancer J. Clin.
70
, 86
–104
(2020
). 60
R. H. I. I.
Andtbacka
et al, “Talimogene laherparepvec improves durable response rate in patients with advanced melanoma
,” J. Clin. Oncol.
33
, 2780
–2788
(2015
). 61
S.
Adams
et al, “Current landscape of immunotherapy in breast cancer: A review
,” JAMA Oncol.
5
, 1205
–1214
(2019
). 62
J.
Bell
and G.
McFadden
, “Viruses for tumor therapy
,” Cell Host Microbe
15
, 260
–265
(2014
). 63
L. G.
dePillis
, H.
Savage
, and A. E.
Radunskaya
, “Mathematical model of colorectal cancer with monoclonal antibody treatments
,” Br. J. Med. Med. Res.
4
, 3101
–3131
(2013
). 64
A. L.
Jenner
, F.
Frascoli
, C.-O.
Yun
, and P. S.
Kim
, “Optimising hydrogel release profiles for viro-immunotherapy using oncolytic adenovirus expressing IL-12 and GM-CSF with immature dendritic cells
,” Appl. Sci.
10
, 2872
(2020
). 65
M.
Jafarnejad
et al, “A computational model of neoadjuvant PD-1 inhibition in non-small cell lung cancer
,” AAPS J.
21
, 79
(2019
). 66
G. P.
Adams
and L. M.
Weiner
, “Monoclonal antibody therapy of cancer
,” Nat. Biotechnol.
23
, 1147
–1157
(2005
). 67
L.
de Pillis
et al, “Mathematical model creation for cancer chemo-immunotherapy
,” Comput. Math. Methods Med.
10
, 165
–184
(2009
). 68
D. M.
Pardoll
, “The blockade of immune checkpoints in cancer immunotherapy
,” Nat. Rev. Cancer
12
, 252
–264
(2012
). 69
M.
Craig
, A. L.
Jenner
, B.
Namgung
, L. P.
Lee
, and A.
Goldman
, “Engineering in medicine to address the challenge of cancer drug resistance: From micro- and nanotechnologies to computational and mathematical modeling
,” Chem. Rev.
(published online 2020
). 70
T.
McGranahan
, K. E.
Therkelsen
, S.
Ahmad
, and S.
Nagpal
, “Current state of immunotherapy for treatment of glioblastoma
,” Curr. Treat. Options Oncol.
20
, 24
(2019
). 71
P.
Darvin
, S. M.
Toor
, V.
Sasidharan Nair
, and E.
Elkord
, “Immune checkpoint inhibitors: Recent progress and potential biomarkers
,” Exp. Mol. Med.
50
, 165
(2018
). 72
L. A.
Torre
, R. L.
Siegel
, and A.
Jemal
, “Lung cancer statistics
,” in Lung Cancer and Personalized Medicine
(Springer
, 2016
), pp. 1
–19
.73
E.
Sharon
, H.
Streicher
, P.
Goncalves
, and H. X.
Chen
, “Immune checkpoint inhibitors in clinical trials
,” Chin. J. Cancer
33
, 434
–444
(2014
). 74
P.
Bulutay
et al, “The importance of histological patterns on PD-L1 staining heterogeneity: Should we use pattern-based approach for selecting tumor samples for PD-L1 testing in lung adenocarcinomas?
,” Turk. J. Med. Sci.
(in press
) (2020
). 75
Y.
Cao
and W. J.
Jusko
, “Survey of monoclonal antibody disposition in man utilizing a minimal physiologically-based pharmacokinetic model
,” J. Pharmacokinet. Pharmacodyn.
41
, 571
–580
(2014
). 76
C. A.
Santa-Maria
et al, “A pilot study of durvalumab and tremelimumab and immunogenomic dynamics in metastatic breast cancer
,” Oncotarget
9
, 18985
(2018
). 77
J. J.
Havel
, D.
Chowell
, and T. A.
Chan
, “The evolving landscape of biomarkers for checkpoint inhibitor immunotherapy
,” Nat. Rev. Cancer
19
, 133
–150
(2019
). 78
H.
Ma
et al, “Combination therapy with T cell engager and PD-L1 blockade enhances the antitumor potency of T cells as predicted by a QSP model
,” J. Immunother. Cancer
8
, e001141
(2020
). 79
M.
Danhof
, E. C.
De Lange
, O. E.
Della Pasqua
, B. A.
Ploeger
, and R. A.
Voskuyl
, “Mechanism-based pharmacokinetic-pharmacodynamic (PK-PD) modeling in translational drug research
,” Trends Pharmacol. Sci.
29
, 186
–191
(2008
). 80
C.
Oelkrug
and J. M.
Ramage
, “Enhancement of T cell recruitment and infiltration into tumours
,” Clin. Exp. Immunol.
178
, 1
–8
(2014
). 81
T. M.
Pearl
, J. M.
Markert
, K. A.
Cassady
, and M. G.
Ghonime
, “Oncolytic virus-based cytokine expression to improve immune activity in brain and solid tumors
,” Mol. Ther. Oncolytics
13
, 14
–21
(2019
). 82
H.
Fukuhara
, Y.
Ino
, and T.
Todo
, “Oncolytic virus therapy: A new era of cancer treatment at Dawn
,” Cancer Sci.
107
, 1373
(2016
). 83
V. D.
Lima
et al, “Expanded access to highly active antiretroviral therapy: A potentially powerful strategy to curb the growth of the HIV epidemic
,” J. Infect. Dis.
198
, 59
–67
(2008
). 84
R.
Loddenkemper
and B.
Hauer
, “Drug-resistant tuberculosis: A worldwide epidemic poses a new challenge
,” Dtsch. Arztebl. Int.
107
, 10
–19
(2010
). 85
K. M.
Stadeli
and D. D.
Richman
, “Rates of emergence of HIV drug resistance in resource-limited settings: A systematic review
,” Antiretrovir. Ther.
18
, 115
–123
(2013
). 86
M.
Mustafa
et al, “Herpes simplex virus infections, pathophysiology and management
,” IOSR J. Dent. Med. Sci.
15
, 85
–91
(2016
). 87
J. T.
Schiffer
et al, “Mathematical modeling of herpes simplex virus-2 suppression with pritelivir predicts trial outcomes
,” Sci. Transl. Med.
8
, 324ra15
(2016
). 88
C. M.
Byrne
, S.
Gantt
, and D.
Coombs
, “Effects of spatiotemporal HSV-2 lesion dynamics and antiviral treatment on the risk of HIV-1 acquisition
,” PLoS Comput. Biol.
14
, e1006129
(2018
). 89
R. K. R.
Rajoli
et al, “In silico dose prediction for long-acting rilpivirine and cabotegravir administration to children and adolescents
,” Clin. Pharmacokinet.
57
, 255
–266
(2018
). 90
H.
Kishimoto
and M.
Maehara
, “Compliance and persistence with daily, weekly, and monthly bisphosphonates for osteoporosis in Japan: Analysis of data from the CISA
,” Arch. Osteoporos.
10
, 231
(2015
). 91
O.
Barrière
, J.
Li
, and F.
Nekka
, “A Bayesian approach for the estimation of patient compliance based on the last sampling information
,” J. Pharmacokinet. Pharmacodyn.
38
, 333
–351
(2011
). 92
P.
Assawasuwannakit
, R.
Braund
, and S.
Duffull
, “Quantification of the forgiveness of drugs to imperfect adherence
,” CPT Pharmacometrics Syst. Pharmacol.
4
, 204
–211
(2015
). 93
L.
Rong
, Z.
Feng
, and A. S.
Perelson
, “Emergence of HIV-1 drug resistance during antiretroviral treatment
,” Bull. Math. Biol.
69
, 2027
–2060
(2007
). 94
A. N.
Edginton
, W.
Schmitt
, and S.
Willmann
, “Development and evaluation of a generic physiologically based pharmacokinetic model for children
,” Clin. Pharmacokinet.
45
, 1013
–1034
(2006
). 95
A. R.
Kirtane
et al, “Development of an oral once-weekly drug delivery system for HIV antiretroviral therapy
,” Nat. Commun.
9
, 2
(2018
). 96
R.
Dantzer
, J. C.
O’Connor
, G. G.
Freund
, R. W.
Johnson
, and K. W.
Kelley
, “From inflammation to sickness and depression: When the immune system subjugates the brain
,” Nat. Rev. Neurosci.
9
, 46
–56
(2008
). 97
B.
Vishwanath
and J. L.
Moran
, “SIRS, sepsis and multiorgan failure
,” in Mechanisms of Vascular Disease: A Reference Book for Vascular Specialists
(University of Adelaide Press
, 2011
), pp. 315
–330
, 98
P. S.
Pérez
et al, “Extracellular vesicles and chronic inflammation during HIV infection
,” J. Extracell. Vesicles
8
, 1687275
(2019
). 99
G.
An
, “In silico experiments of existing and hypothetical cytokine-directed clinical trials using agent-based modeling
,” Crit. Care Med.
32
, 2050
–2060
(2004
).100
J.
Metzcar
, Y.
Wang
, R.
Heiland
, and P.
Macklin
, “A review of cell-based computational modeling in cancer biology
,” JCO Clin. Cancer Informatics
2
, 1
–13
(2019
). 101
A. L.
Jenner
, F.
Frascoli
, A. C. F.
Coster
, and P. S.
Kim
, “Enhancing oncolytic virotherapy: Observations from a voronoi cell-based model
,” J. Theor. Biol.
485
, 110052
(2020
). 102
A. J.
Lamparello
et al, “A conceptual time window-based model for the early stratification of trauma patients
,” J. Intern. Med.
286
, 2
–15
(2019
). 103
G.
Clermont
et al, “In silico design of clinical trials: A method coming of age
,” Crit. Care Med.
32
, 2061
(2004
). 104
G.
Arango Duque
and A.
& Descoteaux
, “Macrophage cytokines: Involvement in immunity and infectious diseases
,” Front. Immunol.
5
, 491
(2014
). 105
K.
Reinhart
et al, “Randomized, placebo-controlled trial of the anti-tumor necrosis factor antibody fragment afelimomab in hyperinflammatory response during severe sepsis: The RAMSES study
,” Crit. Care Med.
29
, 765
–769
(2001
). 106
E.
Panacek
, J.
Marshall
, C. J.
Supreme
, S.
Fischkoff
, W.
Barchuk
, and L. B. T.-C.
Teoh
, “Neutralization of TNF by a monoclonal antibody improves survival and reduces organ dysfunction in human sepsis: Results of the MONARCS trial
,” Chest
118
, 88S
(2000
).© 2020 Author(s).
2020
Author(s)
You do not currently have access to this content.