1
|
Balti A, Zugaj D, Fenneteau F, Tremblay PO, Nekka F. Dynamical systems analysis as an additional tool to inform treatment outcomes: The case study of a quantitative systems pharmacology model of immuno-oncology. CHAOS (WOODBURY, N.Y.) 2021; 31:023124. [PMID: 33653032 DOI: 10.1063/5.0022238] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/18/2020] [Accepted: 01/22/2021] [Indexed: 06/12/2023]
Abstract
Quantitative systems pharmacology (QSP) proved to be a powerful tool to elucidate the underlying pathophysiological complexity that is intensified by the biological variability and overlapped by the level of sophistication of drug dosing regimens. Therapies combining immunotherapy with more traditional therapeutic approaches, including chemotherapy and radiation, are increasingly being used. These combinations are purposed to amplify the immune response against the tumor cells and modulate the suppressive tumor microenvironment (TME). In order to get the best performance from these combinatorial approaches and derive rational regimen strategies, a better understanding of the interaction of the tumor with the host immune system is needed. The objective of the current work is to provide new insights into the dynamics of immune-mediated TME and immune-oncology treatment. As a case study, we will use a recent QSP model by Kosinsky et al. [J. Immunother. Cancer 6, 17 (2018)] that aimed to reproduce the dynamics of interaction between tumor and immune system upon administration of radiation therapy and immunotherapy. Adopting a dynamical systems approach, we here investigate the qualitative behavior of the representative components of this QSP model around its key parameters. The ability of T cells to infiltrate tumor tissue, originally identified as responsible for individual therapeutic inter-variability [Y. Kosinsky et al., J. Immunother. Cancer 6, 17 (2018)], is shown here to be a saddle-node bifurcation point for which the dynamical system oscillates between two states: tumor-free or maximum tumor volume. By performing a bifurcation analysis of the physiological system, we identified equilibrium points and assessed their nature. We then used the traditional concept of basin of attraction to assess the performance of therapy. We showed that considering the therapy as input to the dynamical system translates into the changes of the trajectory shapes of the solutions when approaching equilibrium points and thus providing information on the issue of therapy.
Collapse
Affiliation(s)
- Aymen Balti
- Faculty of pharmacy, Université de Montréal, Montreal, Quebec H3C 3J7, Canada
| | - Didier Zugaj
- Syneos Health, Clinical Pharmacology, Quebec, Quebec G1P 0A2, Canada
| | | | | | - Fahima Nekka
- Faculty of pharmacy, Université de Montréal, Montreal, Quebec H3C 3J7, Canada
| |
Collapse
|
2
|
Fornari C, Pin C, Yates JWT, Mettetal JT, Collins TA. Importance of Stability Analysis When Using Nonlinear Semimechanistic Models to Describe Drug-Induced Hematotoxicity. CPT-PHARMACOMETRICS & SYSTEMS PHARMACOLOGY 2020; 9:498-508. [PMID: 32453487 PMCID: PMC7499189 DOI: 10.1002/psp4.12514] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/23/2020] [Accepted: 03/31/2020] [Indexed: 12/27/2022]
Abstract
Stability analysis, often overlooked in pharmacometrics, is essential to explore dynamical systems. The model developed by Friberg et al.1 to describe drug‐induced hematotoxicity is widely used to support decisions across drug development, and parameter values are often identified from observed blood counts. We use stability analysis to study the parametric dependence of stable and unstable solutions of several Friberg‐type models and highlight the risks associated with system instability in the context of nonlinear mixed effects modeling. We emphasize the consequences of unstable solutions on prediction performance by demonstrating nonbiological system behaviors in a real case study of drug‐induced thrombocytopenia. Ultimately, we provide simple criteria for identifying parameters associated with stable solutions of Friberg‐type models. For instance, in the original Friberg model, we find that stability depends only on the parameter that governs the feedback from peripheral cells to progenitors and provide the exact range of values that results in stable solutions.
Collapse
Affiliation(s)
- Chiara Fornari
- Clinical Pharmacology and Safety Sciences, BioPharmaceuticals R&D, AstraZeneca, Cambridge, UK
| | - Carmen Pin
- Clinical Pharmacology and Safety Sciences, BioPharmaceuticals R&D, AstraZeneca, Cambridge, UK
| | - James W T Yates
- Drug Metabolism and Pharmacokinetic, Oncology R&D, AstraZeneca, Cambridge, UK
| | | | - Teresa A Collins
- Clinical Pharmacology and Safety Sciences, BioPharmaceuticals R&D, AstraZeneca, Cambridge, UK
| |
Collapse
|
3
|
Pillai N, Schwartz SL, Ho T, Dokoumetzidis A, Bies R, Freedman I. Estimating parameters of nonlinear dynamic systems in pharmacology using chaos synchronization and grid search. J Pharmacokinet Pharmacodyn 2019; 46:193-210. [PMID: 30929120 PMCID: PMC6491657 DOI: 10.1007/s10928-019-09629-4] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2018] [Accepted: 03/23/2019] [Indexed: 12/04/2022]
Abstract
Bridging fundamental approaches to model optimization for pharmacometricians, systems pharmacologists and statisticians is a critical issue. These fields rely primarily on Maximum Likelihood and Extended Least Squares metrics with iterative estimation of parameters. Our research combines adaptive chaos synchronization and grid search to estimate physiological and pharmacological systems with emergent properties by exploring deterministic methods that are more appropriate and have potentially superior performance than classical numerical approaches, which minimize the sum of squares or maximize the likelihood. We illustrate these issues with an established model of cortisol in human with nonlinear dynamics. The model describes cortisol kinetics over time, including its chaotic oscillations, by a delay differential equation. We demonstrate that chaos synchronization helps to avoid the tendency of the gradient-based optimization algorithms to end up in a local minimum. The subsequent analysis illustrates that the hybrid adaptive chaos synchronization for estimation of linear parameters with coarse-to-fine grid search for optimal values of non-linear parameters can be applied iteratively to accurately estimate parameters and effectively track trajectories for a wide class of noisy chaotic systems.
Collapse
Affiliation(s)
- Nikhil Pillai
- Computational and Data-Enabled Science, University at Buffalo, Buffalo, USA
| | - Sorell L Schwartz
- Department of Pharmacology and Physiology, Georgetown University Medical Center, Georgetown, USA
| | - Thang Ho
- Vertex Pharmaceuticals, Boston, USA
| | - Aris Dokoumetzidis
- Department of Pharmaceutical Technology, University of Athens, Athens, Greece
| | - Robert Bies
- Computational and Data-Enabled Science, University at Buffalo, Buffalo, USA
- Pharmaceutical Science, University at Buffalo, Buffalo, USA
| | | |
Collapse
|
4
|
Michalaki LI, Goussis DA. Asymptotic analysis of a TMDD model: when a reaction contributes to the destruction of its product. J Math Biol 2018; 77:821-855. [DOI: 10.1007/s00285-018-1234-x] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2017] [Revised: 01/28/2018] [Indexed: 02/04/2023]
|
5
|
Taneja A, Vermeulen A, Huntjens DRH, Danhof M, De Lange ECM, Proost JH. Modeling of prolactin response following dopamine D 2 receptor antagonists in rats: can it be translated to clinical dosing? Pharmacol Res Perspect 2017; 5. [PMID: 29226628 PMCID: PMC5723698 DOI: 10.1002/prp2.364] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2017] [Accepted: 09/06/2017] [Indexed: 01/23/2023] Open
Abstract
Prolactin release is a side effect of antipsychotic therapy with dopamine antagonists, observed in rats as well as humans. We examined whether two semimechanistic models could describe prolactin response in rats and subsequently be translated to predict pituitary dopamine D2 receptor occupancy and plasma prolactin concentrations in humans following administration of paliperidone or remoxipride. Data on male Wistar rats receiving single or multiple doses of risperidone, paliperidone, or remoxipride was described by two semimechanistic models, the precursor pool model and the agonist–antagonist interaction model. Using interspecies scaling approaches, human D2 receptor occupancy and plasma prolactin concentrations were predicted for a range of clinical paliperidone and remoxipride doses. The predictions were compared with corresponding observations described in literature as well as with predictions from published models developed on human data. The pool model could predict D2 receptor occupancy and prolactin response in humans following single doses of paliperidone and remoxipride. Tolerance of prolactin release was predicted following multiple doses. The interaction model underpredicted both D2 receptor occupancy and prolactin response. Prolactin elevation may be deployed as a suitable biomarker for interspecies translation and can inform the clinical safe and effective dose range of antipsychotic drugs. While the pool model was more predictive than the interaction model, it overpredicted tolerance on multiple dosing. Shortcomings of the translations reflect the need for better mechanistic models.
Collapse
Affiliation(s)
- Amit Taneja
- Division of Pharmacokinetics, Toxicology and TargetingGroningen Research Institute of PharmacyUniversity of GroningenAntonius Deusinglaan 19713 AVGroningenThe Netherlands
| | - An Vermeulen
- Division of Janssen Pharmaceutica NVClinical Pharmacology and PharmacometricsJanssen Research and DevelopmentBeerseBelgium
| | - Dymphy R. H. Huntjens
- Division of Janssen Pharmaceutica NVClinical Pharmacology and PharmacometricsJanssen Research and DevelopmentBeerseBelgium
| | - Meindert Danhof
- Department of PharmacologyLeiden Academic Center for Drug ResearchLeiden UniversityLeidenThe Netherlands
| | - Elizabeth C. M. De Lange
- Department of PharmacologyLeiden Academic Center for Drug ResearchLeiden UniversityLeidenThe Netherlands
| | - Johannes H. Proost
- Division of Pharmacokinetics, Toxicology and TargetingGroningen Research Institute of PharmacyUniversity of GroningenAntonius Deusinglaan 19713 AVGroningenThe Netherlands
| |
Collapse
|
6
|
Schulthess P, Post TM, Yates J, van der Graaf PH. Frequency-Domain Response Analysis for Quantitative Systems Pharmacology Models. CPT Pharmacometrics Syst Pharmacol 2017; 7:111-123. [PMID: 29193852 PMCID: PMC5824121 DOI: 10.1002/psp4.12266] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2017] [Revised: 10/30/2017] [Accepted: 11/02/2017] [Indexed: 01/06/2023] Open
Abstract
Drug dosing regimen can significantly impact drug effect and, thus, the success of treatments. Nevertheless, trial and error is still the most commonly used method by conventional pharmacometric approaches to optimize dosing regimen. In this tutorial, we utilize four distinct classes of quantitative systems pharmacology models to introduce frequency-domain response analysis, a method widely used in electrical and control engineering that allows the analytical optimization of drug treatment regimen from the dynamics of the model.
Collapse
Affiliation(s)
- Pascal Schulthess
- Systems Biomedicine & PharmacologyLACDR, Leiden UniversityLeidenThe Netherlands
| | - Teun M. Post
- Systems Biomedicine & PharmacologyLACDR, Leiden UniversityLeidenThe Netherlands
- Leiden Experts on Advanced Pharmacokinetics and Pharmacodynamics (LAP&P)LeidenThe Netherlands
| | - James Yates
- DMPK, Oncology, Innovative Medicines and Early Development, AstraZeneca, Chesterford Research ParkUK
| | - Piet H. van der Graaf
- Systems Biomedicine & PharmacologyLACDR, Leiden UniversityLeidenThe Netherlands
- Certara QSP, Canterbury Innovation HouseCanterburyUK
| |
Collapse
|
7
|
Peletier LA, Gabrielsson J. Impact of mathematical pharmacology on practice and theory: four case studies. J Pharmacokinet Pharmacodyn 2017; 45:3-21. [PMID: 28884259 PMCID: PMC5847232 DOI: 10.1007/s10928-017-9539-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2017] [Accepted: 08/18/2017] [Indexed: 11/25/2022]
Abstract
Drug-discovery has become a complex discipline in which the amount of knowledge about human biology, physiology, and biochemistry have increased. In order to harness this complex body of knowledge mathematics can play a critical role, and has actually already been doing so. We demonstrate through four case studies, taken from previously published data and analyses, what we can gain from mathematical/analytical techniques when nonlinear concentration-time courses have to be transformed into their equilibrium concentration-response (target or complex) relationships and new structures of drug potency have to be deciphered; when pattern recognition needs to be carried out for an unconventional response-time dataset; when what-if? predictions beyond the observational concentration-time range need to be made; or when the behaviour of a semi-mechanistic model needs to be elucidated or challenged. These four examples are typical situations when standard approaches known to the general community of pharmacokineticists prove to be inadequate.
Collapse
Affiliation(s)
| | - Johan Gabrielsson
- Division of Pharmacology and Toxicology, Department of Biomedical Sciences and Veterinary Public Health, Swedish University of Agricultural Sciences, Box 7028, 750 07 Uppsala, Sweden
| |
Collapse
|
8
|
Peletier LA, de Winter W. Impact of saturable distribution in compartmental PK models: dynamics and practical use. J Pharmacokinet Pharmacodyn 2017; 44:1-16. [PMID: 28050672 PMCID: PMC5306145 DOI: 10.1007/s10928-016-9500-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2016] [Accepted: 12/05/2016] [Indexed: 11/02/2022]
Abstract
We explore the impact of saturable distribution over the central and the peripheral compartment in pharmacokinetic models, whilst assuming that back flow into the central compartiment is linear. Using simulations and analytical methods we demonstrate characteristic tell-tale differences in plasma concentration profiles of saturable versus linear distribution models, which can serve as a guide to their practical applicability. For two extreme cases, relating to (i) the size of the peripheral compartment with respect to the central compartment and (ii) the magnitude of the back flow as related to direct elimination from the central compartment, we derive explicit approximations which make it possible to give quantitative estimates of parameters. In three appendices we give detailed explanations of how these estimates are derived. They demonstrate how singular perturbation methods can be successfully employed to gain insight in the dynamics of multi-compartment pharmacokinetic models. These appendices are also intended to serve as an introductory tutorial to these ideas.
Collapse
Affiliation(s)
- Lambertus A Peletier
- Mathematical Institute, Leiden University, PB 9512, 2300 RA, Leiden, The Netherlands.
| | - Willem de Winter
- Janssen Research & Development, Janssen Prevention Center, Archimedesweg 6, 2333 CN, Leiden, The Netherlands
| |
Collapse
|
9
|
van den Brink WJ, Wong YC, Gülave B, van der Graaf PH, de Lange ECM. Revealing the Neuroendocrine Response After Remoxipride Treatment Using Multi-Biomarker Discovery and Quantifying It by PK/PD Modeling. AAPS JOURNAL 2016; 19:274-285. [PMID: 27785749 DOI: 10.1208/s12248-016-0002-3] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/29/2016] [Accepted: 10/03/2016] [Indexed: 01/10/2023]
Abstract
To reveal unknown and potentially important mechanisms of drug action, multi-biomarker discovery approaches are increasingly used. Time-course relationships between drug action and multi-biomarker profiles, however, are typically missing, while such relationships will provide increased insight in the underlying body processes. The aim of this study was to investigate the effect of the dopamine D2 antagonist remoxipride on the neuroendocrine system. Different doses of remoxipride (0, 0.7, 5.2, or 14 mg/kg) were administered to rats by intravenous infusion. Serial brain extracellular fluid (brainECF) and plasma samples were collected and analyzed for remoxipride pharmacokinetics (PK). Plasma samples were analyzed for concentrations of the eight pituitary-related hormones as a function of time. A Mann-Whitney test was used to identify the responding hormones, which were further analyzed by pharmacokinetic/pharmacodynamic (PK/PD) modeling. A three-compartment PK model adequately described remoxipride PK in plasma and brainECF. Not only plasma PRL, but also adrenocorticotrophic hormone (ACTH) concentrations were increased, the latter especially at higher concentrations of remoxipride. Brain-derived neurotropic factor (BDNF), follicle stimulating hormone (FSH), growth hormone (GH), luteinizing hormone (LH), and thyroid stimulating hormones (TSH) did not respond to remoxipride at the tested doses, while oxytocin (OXT) measurements were below limit of quantification. Precursor pool models were linked to brainECF remoxipride PK by Emax drug effect models, which could accurately describe the PRL and ACTH responses. To conclude, this study shows how a multi-biomarker identification approach combined with PK/PD modeling can reveal and quantify a neuroendocrine multi-biomarker response for single drug action.
Collapse
Affiliation(s)
- Willem J van den Brink
- Systems Pharmacology, Division of Pharmacology, Leiden Academic Center for Drug Research, Leiden University, PO box 9502, 2300 RA, Leiden, The Netherlands
| | - Yin C Wong
- Systems Pharmacology, Division of Pharmacology, Leiden Academic Center for Drug Research, Leiden University, PO box 9502, 2300 RA, Leiden, The Netherlands
| | - Berfin Gülave
- Systems Pharmacology, Division of Pharmacology, Leiden Academic Center for Drug Research, Leiden University, PO box 9502, 2300 RA, Leiden, The Netherlands
| | - Piet H van der Graaf
- Systems Pharmacology, Division of Pharmacology, Leiden Academic Center for Drug Research, Leiden University, PO box 9502, 2300 RA, Leiden, The Netherlands.,Certara QSP, Canterbury Innovation House, Canterbury, UK
| | - Elizatbeth C M de Lange
- Systems Pharmacology, Division of Pharmacology, Leiden Academic Center for Drug Research, Leiden University, PO box 9502, 2300 RA, Leiden, The Netherlands.
| |
Collapse
|