@article{KnoechelKloftHuisinga2018, author = {Kn{\"o}chel, Jane and Kloft, Charlotte and Huisinga, Wilhelm}, title = {Understanding and reducing complex systems pharmacology models based on a novel input-response index}, series = {Journal of pharmacokinetics and pharmacodynamics}, volume = {45}, journal = {Journal of pharmacokinetics and pharmacodynamics}, number = {1}, publisher = {Springer Science + Business Media B.V.}, address = {New York}, issn = {1567-567X}, doi = {10.1007/s10928-017-9561-x}, pages = {139 -- 157}, year = {2018}, abstract = {A growing understanding of complex processes in biology has led to large-scale mechanistic models of pharmacologically relevant processes. These models are increasingly used to study the response of the system to a given input or stimulus, e.g., after drug administration. Understanding the input-response relationship, however, is often a challenging task due to the complexity of the interactions between its constituents as well as the size of the models. An approach that quantifies the importance of the different constituents for a given input-output relationship and allows to reduce the dynamics to its essential features is therefore highly desirable. In this article, we present a novel state- and time-dependent quantity called the input-response index that quantifies the importance of state variables for a given input-response relationship at a particular time. It is based on the concept of time-bounded controllability and observability, and defined with respect to a reference dynamics. In application to the brown snake venom-fibrinogen (Fg) network, the input-response indices give insight into the coordinated action of specific coagulation factors and about those factors that contribute only little to the response. We demonstrate how the indices can be used to reduce large-scale models in a two-step procedure: (i) elimination of states whose dynamics have only minor impact on the input-response relationship, and (ii) proper lumping of the remaining (lower order) model. In application to the brown snake venom-fibrinogen network, this resulted in a reduction from 62 to 8 state variables in the first step, and a further reduction to 5 state variables in the second step. We further illustrate that the sequence, in which a recursive algorithm eliminates and/or lumps state variables, has an impact on the final reduced model. The input-response indices are particularly suited to determine an informed sequence, since they are based on the dynamics of the original system. In summary, the novel measure of importance provides a powerful tool for analysing the complex dynamics of large-scale systems and a means for very efficient model order reduction of nonlinear systems.}, language = {en} } @article{DemarisWidigsonIlvemarketal.2022, author = {D{\´e}maris, Alix and Widigson, Ella S. K. and Ilvemark, Johan F. K. F. and Steenholdt, Casper and Seidelin, Jakob B. and Huisinga, Wilhelm and Michelet, Robin and Aulin, Linda B. S. and Kloft, Charlotte}, title = {Ulcerative colitis and acute severe ulcerative colitis patients are overlooked in infliximab population pharmacokinetic models}, series = {Pharmaceutics / Molecular Diversity Preservation International}, volume = {14}, journal = {Pharmaceutics / Molecular Diversity Preservation International}, number = {10}, publisher = {MDPI}, address = {Basel}, issn = {1999-4923}, doi = {10.3390/pharmaceutics14102095}, pages = {32}, year = {2022}, abstract = {Ulcerative colitis (UC) is part of the inflammatory bowels diseases, and moderate to severe UC patients can be treated with anti-tumour necrosis alpha monoclonal antibodies, including infliximab (IFX). Even though treatment of UC patients by IFX has been in place for over a decade, many gaps in modelling of IFX PK in this population remain. This is even more true for acute severe UC (ASUC) patients for which early prediction of IFX pharmacokinetic (PK) could highly improve treatment outcome. Thus, this review aims to compile and analyse published population PK models of IFX in UC and ASUC patients, and to assess the current knowledge on disease activity impact on IFX PK. For this, a semi-systematic literature search was conducted, from which 26 publications including a population PK model analysis of UC patients receiving IFX therapy were selected. Amongst those, only four developed a model specifically for UC patients, and only three populations included severe UC patients. Investigations of disease activity impact on PK were reported in only 4 of the 14 models selected. In addition, the lack of reported model codes and assessment of predictive performance make the use of published models in a clinical setting challenging. Thus, more comprehensive investigation of PK in UC and ASUC is needed as well as more adequate reports on developed models and their evaluation in order to apply them in a clinical setting.}, language = {en} } @article{WichaHuisingaKloft2017, author = {Wicha, Sebastian G. and Huisinga, Wilhelm and Kloft, Charlotte}, title = {Translational pharmacometric evaluation of typical antibiotic broad-spectrum combination therapies against staphylococcus aureus exploiting in vitro information}, series = {CPT: pharmacometrics \& systems pharmacology}, volume = {6}, journal = {CPT: pharmacometrics \& systems pharmacology}, publisher = {Wiley}, address = {Hoboken}, issn = {2163-8306}, doi = {10.1002/psp4.12197}, pages = {512 -- 522}, year = {2017}, abstract = {Broad-spectrum antibiotic combination therapy is frequently applied due to increasing resistance development of infective pathogens. The objective of the present study was to evaluate two common empiric broad-spectrum combination therapies consisting of either linezolid (LZD) or vancomycin (VAN) combined with meropenem (MER) against Staphylococcus aureus (S. aureus) as the most frequent causative pathogen of severe infections. A semimechanistic pharmacokinetic-pharmacodynamic (PK-PD) model mimicking a simplified bacterial life-cycle of S. aureus was developed upon time-kill curve data to describe the effects of LZD, VAN, and MER alone and in dual combinations. The PK-PD model was successfully (i) evaluated with external data from two clinical S. aureus isolates and further drug combinations and (ii) challenged to predict common clinical PK-PD indices and breakpoints. Finally, clinical trial simulations were performed that revealed that the combination of VAN-MER might be favorable over LZD-MER due to an unfavorable antagonistic interaction between LZD and MER.}, language = {en} } @article{MuellerSchoellGroenlandScherfClaveletal.2020, author = {Mueller-Schoell, Anna and Groenland, Stefanie L. and Scherf-Clavel, Oliver and van Dyk, Madele and Huisinga, Wilhelm and Michelet, Robin and Jaehde, Ulrich and Steeghs, Neeltje and Huitema, Alwin D. R. and Kloft, Charlotte}, title = {Therapeutic drug monitoring of oral targeted antineoplastic drugs}, series = {European journal of clinical pharmacology}, volume = {77}, journal = {European journal of clinical pharmacology}, number = {4}, publisher = {Springer}, address = {Heidelberg}, issn = {0031-6970}, doi = {10.1007/s00228-020-03014-8}, pages = {441 -- 464}, year = {2020}, abstract = {Purpose This review provides an overview of the current challenges in oral targeted antineoplastic drug (OAD) dosing and outlines the unexploited value of therapeutic drug monitoring (TDM). Factors influencing the pharmacokinetic exposure in OAD therapy are depicted together with an overview of different TDM approaches. Finally, current evidence for TDM for all approved OADs is reviewed. Methods A comprehensive literature search (covering literature published until April 2020), including primary and secondary scientific literature on pharmacokinetics and dose individualisation strategies for OADs, together with US FDA Clinical Pharmacology and Biopharmaceutics Reviews and the Committee for Medicinal Products for Human Use European Public Assessment Reports was conducted. Results OADs are highly potent drugs, which have substantially changed treatment options for cancer patients. Nevertheless, high pharmacokinetic variability and low treatment adherence are risk factors for treatment failure. TDM is a powerful tool to individualise drug dosing, ensure drug concentrations within the therapeutic window and increase treatment success rates. After reviewing the literature for 71 approved OADs, we show that exposure-response and/or exposure-toxicity relationships have been established for the majority. Moreover, TDM has been proven to be feasible for individualised dosing of abiraterone, everolimus, imatinib, pazopanib, sunitinib and tamoxifen in prospective studies. There is a lack of experience in how to best implement TDM as part of clinical routine in OAD cancer therapy. Conclusion Sub-therapeutic concentrations and severe adverse events are current challenges in OAD treatment, which can both be addressed by the application of TDM-guided dosing, ensuring concentrations within the therapeutic window.}, language = {en} } @article{MoldenhawerMorenoSchindleretal.2022, author = {Moldenhawer, Ted and Moreno, Eduardo and Schindler, Daniel and Flemming, Sven and Holschneider, Matthias and Huisinga, Wilhelm and Alonso, Sergio and Beta, Carsten}, title = {Spontaneous transitions between amoeboid and keratocyte-like modes of migration}, series = {Frontiers in Cell and Developmental Biology}, volume = {10}, journal = {Frontiers in Cell and Developmental Biology}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-634X}, doi = {10.3389/fcell.2022.898351}, pages = {13}, year = {2022}, abstract = {The motility of adherent eukaryotic cells is driven by the dynamics of the actin cytoskeleton. Despite the common force-generating actin machinery, different cell types often show diverse modes of locomotion that differ in their shape dynamics, speed, and persistence of motion. Recently, experiments in Dictyostelium discoideum have revealed that different motility modes can be induced in this model organism, depending on genetic modifications, developmental conditions, and synthetic changes of intracellular signaling. Here, we report experimental evidence that in a mutated D. discoideum cell line with increased Ras activity, switches between two distinct migratory modes, the amoeboid and fan-shaped type of locomotion, can even spontaneously occur within the same cell. We observed and characterized repeated and reversible switchings between the two modes of locomotion, suggesting that they are distinct behavioral traits that coexist within the same cell. We adapted an established phenomenological motility model that combines a reaction-diffusion system for the intracellular dynamics with a dynamic phase field to account for our experimental findings.}, language = {en} } @article{HenrichJoergerKraffetal.2017, author = {Henrich, Andrea and Joerger, Markus and Kraff, Stefanie and Jaehde, Ulrich and Huisinga, Wilhelm and Kloft, Charlotte and Parra-Guillen, Zinnia Patricia}, title = {Semimechanistic Bone Marrow Exhaustion Pharmacokinetic/Pharmacodynamic Model for Chemotherapy-Induced Cumulative Neutropenia}, series = {Journal of Pharmacology and Experimental Therapeutics}, volume = {362}, journal = {Journal of Pharmacology and Experimental Therapeutics}, number = {2}, publisher = {American Society for Pharmacology and Experimental Therapeutics}, address = {Bethesda}, issn = {0022-3565}, doi = {10.1124/jpet.117.240309}, pages = {347 -- 358}, year = {2017}, abstract = {Paclitaxel is a commonly used cytotoxic anticancer drug with potentially life-threatening toxicity at therapeutic doses and high interindividual pharmacokinetic variability. Thus, drug and effect monitoring is indicated to control dose-limiting neutropenia. Joerger et al. (2016) developed a dose individualization algorithm based on a pharmacokinetic (PK)/pharmacodynamic (PD) model describing paclitaxel and neutrophil concentrations. Furthermore, the algorithm was prospectively compared in a clinical trial against standard dosing (Central European Society for Anticancer Drug Research Study of Paclitaxel Therapeutic Drug Monitoring; 365 patients, 720 cycles) but did not substantially improve neutropenia. This might be caused by misspecifications in the PK/PD model underlying the algorithm, especially without consideration of the observed cumulative pattern of neutropenia or the platinum-based combination therapy, both impacting neutropenia. This work aimed to externally evaluate the original PK/PD model for potential misspecifications and to refine the PK/PD model while considering the cumulative neutropenia pattern and the combination therapy. An underprediction was observed for the PK (658 samples), the PK parameters, and these parameters were re-estimated using the original estimates as prior information. Neutrophil concentrations (3274 samples) were over-predicted by the PK/PD model, especially for later treatment cycles when the cumulative pattern aggravated neutropenia. Three different modeling approaches (two from the literature and one newly developed) were investigated. The newly developed model, which implemented the bone marrow hypothesis semiphysiologically, was superior. This model further included an additive effect for toxicity of carboplatin combination therapy. Overall, a physiologically plausible PK/PD model was developed that can be used for dose adaptation simulations and prospective studies to further improve paclitaxel/ carboplatin combination therapy.}, language = {en} } @article{EhmannZollerMinichmayretal.2017, author = {Ehmann, Lisa and Zoller, Michael and Minichmayr, Iris K. and Scharf, Christina and Maier, Barbara and Schmitt, Maximilian V. and Hartung, Niklas and Huisinga, Wilhelm and Vogeser, Michael and Frey, Lorenz and Zander, Johannes and Kloft, Charlotte}, title = {Role of renal function in risk assessment of target non-attainment after standard dosing of meropenem in critically ill patients}, series = {Critical care}, volume = {21}, journal = {Critical care}, publisher = {BioMed Central}, address = {London}, issn = {1466-609X}, doi = {10.1186/s13054-017-1829-4}, pages = {14}, year = {2017}, abstract = {Background: Severe bacterial infections remain a major challenge in intensive care units because of their high prevalence and mortality. Adequate antibiotic exposure has been associated with clinical success in critically ill patients. The objective of this study was to investigate the target attainment of standard meropenem dosing in a heterogeneous critically ill population, to quantify the impact of the full renal function spectrum on meropenem exposure and target attainment, and ultimately to translate the findings into a tool for practical application. Methods: A prospective observational single-centre study was performed with critically ill patients with severe infections receiving standard dosing of meropenem. Serial blood samples were drawn over 4 study days to determine meropenem serum concentrations. Renal function was assessed by creatinine clearance according to the Cockcroft and Gault equation (CLCRCG). Variability in meropenem serum concentrations was quantified at the middle and end of each monitored dosing interval. The attainment of two pharmacokinetic/pharmacodynamic targets (100\% T->MIC, 50\% T->4xMIC) was evaluated for minimum inhibitory concentration (MIC) values of 2 mg/L and 8 mg/L and standard meropenem dosing (1000 mg, 30-minute infusion, every 8 h). Furthermore, we assessed the impact of CLCRCG on meropenem concentrations and target attainment and developed a tool for risk assessment of target non-attainment. Results: Large inter-and intra-patient variability in meropenem concentrations was observed in the critically ill population (n = 48). Attainment of the target 100\% T->MIC was merely 48.4\% and 20.6\%, given MIC values of 2 mg/L and 8 mg/L, respectively, and similar for the target 50\% T->4xMIC. A hyperbolic relationship between CLCRCG (25-255 ml/minute) and meropenem serum concentrations at the end of the dosing interval (C-8h) was derived. For infections with pathogens of MIC 2 mg/L, mild renal impairment up to augmented renal function was identified as a risk factor for target non-attainment (for MIC 8 mg/L, additionally, moderate renal impairment). Conclusions: The investigated standard meropenem dosing regimen appeared to result in insufficient meropenem exposure in a considerable fraction of critically ill patients. An easy-and free-to-use tool (the MeroRisk Calculator) for assessing the risk of target non-attainment for a given renal function and MIC value was developed.}, language = {en} } @article{GrisicEserHuisingaetal.2020, author = {Grisic, Ana-Marija and Eser, Alexander and Huisinga, Wilhelm and Reinisch, Walter and Kloft, Charlotte}, title = {Quantitative relationship between infliximab exposure and inhibition of C-reactive protein synthesis to support inflammatory bowel disease management}, series = {British journal of clinical pharmacology}, volume = {87}, journal = {British journal of clinical pharmacology}, number = {5}, publisher = {Wiley}, address = {Hoboken}, issn = {0306-5251}, doi = {10.1111/bcp.14648}, pages = {2374 -- 2384}, year = {2020}, abstract = {Aim Quantitative and kinetic insights into the drug exposure-disease response relationship might enhance our knowledge on loss of response and support more effective monitoring of inflammatory activity by biomarkers in patients with inflammatory bowel disease (IBD) treated with infliximab (IFX). This study aimed to derive recommendations for dose adjustment and treatment optimisation based on mechanistic characterisation of the relationship between IFX serum concentration and C-reactive protein (CRP) concentration.
Methods Data from an investigator-initiated trial included 121 patients with IBD during IFX maintenance treatment. Serum concentrations of IFX, antidrug antibodies (ADA), CRP, and disease-related covariates were determined at the mid-term and end of a dosing interval. Data were analysed using a pharmacometric nonlinear mixed-effects modelling approach. An IFX exposure-CRP model was generated and applied to evaluate dosing regimens to achieve CRP remission.
Results The generated quantitative model showed that IFX has the potential to inhibit up to 72\% (9\% relative standard error [RSE]) of CRP synthesis in a patient. IFX concentration leading to 90\% of the maximum CRP synthesis inhibition was 18.4 mu g/mL (43\% RSE). Presence of ADA was the most influential factor on IFX exposure. With standard dosing strategy, >= 55\% of ADA+ patients experienced CRP nonremission. Shortening the dosing interval and co-therapy with immunomodulators were found to be the most beneficial strategies to maintain CRP remission.
Conclusions With the generated model we could for the first time establish a robust relationship between IFX exposure and CRP synthesis inhibition, which could be utilised for treatment optimisation in IBD patients.}, language = {en} } @article{WeisseMiddletonHuisinga2010, author = {Weiße, Andrea Y. and Middleton, Richard H. and Huisinga, Wilhelm}, title = {Quantifying uncertainty, variability and likelihood for ordinary differential equation models}, issn = {1752-0509}, doi = {10.1186/1752-0509-4-144}, year = {2010}, abstract = {Background: In many applications, ordinary differential equation (ODE) models are subject to uncertainty or variability in initial conditions and parameters. Both, uncertainty and variability can be quantified in terms of a probability density function on the state and parameter space. Results: The partial differential equation that describes the evolution of this probability density function has a form that is particularly amenable to application of the well- known method of characteristics. The value of the density at some point in time is directly accessible by the solution of the original ODE extended by a single extra dimension (for the value of the density). This leads to simple methods for studying uncertainty, variability and likelihood, with significant advantages over more traditional Monte Carlo and related approaches especially when studying regions with low probability. Conclusions: While such approaches based on the method of characteristics are common practice in other disciplines, their advantages for the study of biological systems have so far remained unrecognized. Several examples illustrate performance and accuracy of the approach and its limitations.}, language = {en} } @article{MakaravaMenzThevesetal.2014, author = {Makarava, Natallia and Menz, Stephan and Theves, Matthias and Huisinga, Wilhelm and Beta, Carsten and Holschneider, Matthias}, title = {Quantifying the degree of persistence in random amoeboid motion based on the Hurst exponent of fractional Brownian motion}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {90}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {4}, publisher = {American Physical Society}, address = {College Park}, issn = {1539-3755}, doi = {10.1103/PhysRevE.90.042703}, pages = {6}, year = {2014}, abstract = {Amoebae explore their environment in a random way, unless external cues like, e. g., nutrients, bias their motion. Even in the absence of cues, however, experimental cell tracks show some degree of persistence. In this paper, we analyzed individual cell tracks in the framework of a linear mixed effects model, where each track is modeled by a fractional Brownian motion, i.e., a Gaussian process exhibiting a long-term correlation structure superposed on a linear trend. The degree of persistence was quantified by the Hurst exponent of fractional Brownian motion. Our analysis of experimental cell tracks of the amoeba Dictyostelium discoideum showed a persistent movement for the majority of tracks. Employing a sliding window approach, we estimated the variations of the Hurst exponent over time, which allowed us to identify points in time, where the correlation structure was distorted ("outliers"). Coarse graining of track data via down-sampling allowed us to identify the dependence of persistence on the spatial scale. While one would expect the (mode of the) Hurst exponent to be constant on different temporal scales due to the self-similarity property of fractional Brownian motion, we observed a trend towards stronger persistence for the down-sampled cell tracks indicating stronger persistence on larger time scales.}, language = {en} } @article{NassarHohmannMicheletetal.2022, author = {Nassar, Yomna M. and Hohmann, Nicolas and Michelet, Robin and Gottwalt, Katharina and Meid, Andreas D. and Burhenne, J{\"u}rgen and Huisinga, Wilhelm and Haefeli, Walter E. and Mikus, Gerd and Kloft, Charlotte}, title = {Quantification of the Time Course of CYP3A Inhibition, Activation, and Induction Using a Population Pharmacokinetic Model of Microdosed Midazolam Continuous Infusion}, series = {Clinical Pharmacokinetics}, volume = {61}, journal = {Clinical Pharmacokinetics}, number = {11}, publisher = {Springer}, address = {Northcote}, issn = {0312-5963}, doi = {10.1007/s40262-022-01175-6}, pages = {1595 -- 1607}, year = {2022}, abstract = {Background Cytochrome P450 (CYP) 3A contributes to the metabolism of many approved drugs. CYP3A perpetrator drugs can profoundly alter the exposure of CYP3A substrates. However, effects of such drug-drug interactions are usually reported as maximum effects rather than studied as time-dependent processes. Identification of the time course of CYP3A modulation can provide insight into when significant changes to CYP3A activity occurs, help better design drug-drug interaction studies, and manage drug-drug interactions in clinical practice. Objective We aimed to quantify the time course and extent of the in vivo modulation of different CYP3A perpetrator drugs on hepatic CYP3A activity and distinguish different modulatory mechanisms by their time of onset, using pharmacologically inactive intravenous microgram doses of the CYP3A-specific substrate midazolam, as a marker of CYP3A activity. Methods Twenty-four healthy individuals received an intravenous midazolam bolus followed by a continuous infusion for 10 or 36 h. Individuals were randomized into four arms: within each arm, two individuals served as a placebo control and, 2 h after start of the midazolam infusion, four individuals received the CYP3A perpetrator drug: voriconazole (inhibitor, orally or intravenously), rifampicin (inducer, orally), or efavirenz (activator, orally). After midazolam bolus administration, blood samples were taken every hour (rifampicin arm) or every 15 min (remaining study arms) until the end of midazolam infusion. A total of 1858 concentrations were equally divided between midazolam and its metabolite, 1'-hydroxymidazolam. A nonlinear mixed-effects population pharmacokinetic model of both compounds was developed using NONMEM (R). CYP3A activity modulation was quantified over time, as the relative change of midazolam clearance encountered by the perpetrator drug, compared to the corresponding clearance value in the placebo arm. Results Time course of CYP3A modulation and magnitude of maximum effect were identified for each perpetrator drug. While efavirenz CYP3A activation was relatively fast and short, reaching a maximum after approximately 2-3 h, the induction effect of rifampicin could only be observed after 22 h, with a maximum after approximately 28-30 h followed by a steep drop to almost baseline within 1-2 h. In contrast, the inhibitory impact of both oral and intravenous voriconazole was prolonged with a steady inhibition of CYP3A activity followed by a gradual increase in the inhibitory effect until the end of sampling at 8 h. Relative maximum clearance changes were +59.1\%, +46.7\%, -70.6\%, and -61.1\% for efavirenz, rifampicin, oral voriconazole, and intravenous voriconazole, respectively. Conclusions We could distinguish between different mechanisms of CYP3A modulation by the time of onset. Identification of the time at which clearance significantly changes, per perpetrator drug, can guide the design of an optimal sampling schedule for future drug-drug interaction studies. The impact of a short-term combination of different perpetrator drugs on the paradigm CYP3A substrate midazolam was characterized and can define combination intervals in which no relevant interaction is to be expected.}, language = {en} } @article{KrippendorffOyarzunHuisinga2012, author = {Krippendorff, Ben-Fillippo and Oyarz{\´u}n, Diego A. and Huisinga, Wilhelm}, title = {Predicting the F(ab)-mediated effect of monoclonal antibodies in vivo by combining cell-level kinetic and pharmacokinetic modelling}, series = {Journal of pharmacokinetics and pharmacodynamics}, volume = {39}, journal = {Journal of pharmacokinetics and pharmacodynamics}, number = {2}, publisher = {Springer}, address = {New York}, issn = {1567-567X}, doi = {10.1007/s10928-012-9243-7}, pages = {125 -- 139}, year = {2012}, abstract = {Cell-level kinetic models for therapeutically relevant processes increasingly benefit the early stages of drug development. Later stages of the drug development processes, however, rely on pharmacokinetic compartment models while cell-level dynamics are typically neglected. We here present a systematic approach to integrate cell-level kinetic models and pharmacokinetic compartment models. Incorporating target dynamics into pharmacokinetic models is especially useful for the development of therapeutic antibodies because their effect and pharmacokinetics are inherently interdependent. The approach is illustrated by analysing the F(ab)-mediated inhibitory effect of therapeutic antibodies targeting the epidermal growth factor receptor. We build a multi-level model for anti-EGFR antibodies by combining a systems biology model with in vitro determined parameters and a pharmacokinetic model based on in vivo pharmacokinetic data. Using this model, we investigated in silico the impact of biochemical properties of anti-EGFR antibodies on their F(ab)-mediated inhibitory effect. The multi-level model suggests that the F(ab)-mediated inhibitory effect saturates with increasing drug-receptor affinity, thereby limiting the impact of increasing antibody affinity on improving the effect. This indicates that observed differences in the therapeutic effects of high affinity antibodies in the market and in clinical development may result mainly from Fc-mediated indirect mechanisms such as antibody-dependent cell cytotoxicity.}, language = {en} } @article{MelinParraGuillenHartungetal.2018, author = {Melin, Johanna and Parra-Guillen, Zinnia Patricia and Hartung, Niklas and Huisinga, Wilhelm and Ross, Richard J. and Whitaker, Martin J. and Kloft, Charlotte}, title = {Predicting Cortisol Exposure from Paediatric Hydrocortisone Formulation Using a Semi-Mechanistic Pharmacokinetic Model Established in Healthy Adults}, series = {Clinical Pharmacokinetics}, volume = {57}, journal = {Clinical Pharmacokinetics}, number = {4}, publisher = {Springer}, address = {Northcote}, issn = {0312-5963}, doi = {10.1007/s40262-017-0575-8}, pages = {515 -- 527}, year = {2018}, abstract = {Background and objective Optimisation of hydrocortisone replacement therapy in children is challenging as there is currently no licensed formulation and dose in Europe for children under 6 years of age. In addition, hydrocortisone has non-linear pharmacokinetics caused by saturable plasma protein binding. A paediatric hydrocortisone formulation, Infacort (R) oral hydrocortisone granules with taste masking, has therefore been developed. The objective of this study was to establish a population pharmacokinetic model based on studies in healthy adult volunteers to predict hydrocortisone exposure in paediatric patients with adrenal insufficiency. Methods Cortisol and binding protein concentrations were evaluated in the absence and presence of dexamethasone in healthy volunteers (n = 30). Dexamethasone was used to suppress endogenous cortisol concentrations prior to and after single doses of 0.5, 2, 5 and 10 mg of Infacort (R) or 20 mg of Infacort (R)/hydrocortisone tablet/hydrocortisone intravenously. A plasma protein binding model was established using unbound and total cortisol concentrations, and sequentially integrated into the pharmacokinetic model. Results Both specific (non-linear) and non-specific (linear) protein binding were included in the cortisol binding model. A two-compartment disposition model with saturable absorption and constant endogenous cortisol baseline (Baseline (cort),15.5 nmol/L) described the data accurately. The predicted cortisol exposure for a given dose varied considerably within a small body weight range in individuals weighing < 20 kg. Conclusions Our semi-mechanistic population pharmacokinetic model for hydrocortisone captures the complex pharmacokinetics of hydrocortisone in a simplified but comprehensive framework. The predicted cortisol exposure indicated the importance of defining an accurate hydrocortisone dose to mimic physiological concentrations for neonates and infants weighing < 20 kg.}, language = {en} } @article{HetheyHartungWangorschetal.2021, author = {Hethey, Christoph Philipp and Hartung, Niklas and Wangorsch, Gaby and Weisser, Karin and Huisinga, Wilhelm}, title = {Physiology-based toxicokinetic modelling of aluminium in rat and man}, series = {Archives of toxicology : official journal of EUROTOX}, volume = {95}, journal = {Archives of toxicology : official journal of EUROTOX}, number = {9}, publisher = {Springer}, address = {Berlin ; Heidelberg}, issn = {0340-5761}, doi = {10.1007/s00204-021-03107-y}, pages = {2977 -- 3000}, year = {2021}, abstract = {A sufficient quantitative understanding of aluminium (Al) toxicokinetics (TK) in man is still lacking, although highly desirable for risk assessment of Al exposure. Baseline exposure and the risk of contamination severely limit the feasibility of TK studies administering the naturally occurring isotope Al-27, both in animals and man. These limitations are absent in studies with Al-26 as a tracer, but tissue data are limited to animal studies. A TK model capable of inter-species translation to make valid predictions of Al levels in humans-especially in toxicological relevant tissues like bone and brain-is urgently needed. Here, we present: (i) a curated dataset which comprises all eligible studies with single doses of Al-26 tracer administered as citrate or chloride salts orally and/or intravenously to rats and humans, including ultra-long-term kinetic profiles for plasma, blood, liver, spleen, muscle, bone, brain, kidney, and urine up to 150 weeks; and (ii) the development of a physiology-based (PB) model for Al TK after intravenous and oral administration of aqueous Al citrate and Al chloride solutions in rats and humans. Based on the comprehensive curated Al-26 dataset, we estimated substance-dependent parameters within a non-linear mixed-effect modelling context. The model fitted the heterogeneous Al-26 data very well and was successfully validated against datasets in rats and humans. The presented PBTK model for Al, based on the most extensive and diverse dataset of Al exposure to date, constitutes a major advancement in the field, thereby paving the way towards a more quantitative risk assessment in humans.}, language = {en} } @article{KluweMicheletMuellerSchoelletal.2020, author = {Kluwe, Franziska and Michelet, Robin and M{\"u}ller-Sch{\"o}ll, Anna and Maier, Corinna and Klopp-Schulze, Lena and van Dyk, Madele and Mikus, Gerd and Huisinga, Wilhelm and Kloft, Charlotte}, title = {Perspectives on model-informed precision dosing in the digital health era}, series = {Clinical pharmacology \& therapeutics}, volume = {109}, journal = {Clinical pharmacology \& therapeutics}, number = {1}, publisher = {Wiley}, address = {Hoboken}, issn = {0009-9236}, doi = {10.1002/cpt.2049}, pages = {29 -- 36}, year = {2020}, language = {en} } @article{HartungWahlRastogietal.2021, author = {Hartung, Niklas and Wahl, Martin and Rastogi, Abhishake and Huisinga, Wilhelm}, title = {Nonparametric goodness-of-fit testing for parametric covariate models in pharmacometric analyses}, series = {CPT: pharmacometrics \& systems pharmacology}, volume = {10}, journal = {CPT: pharmacometrics \& systems pharmacology}, number = {6}, publisher = {Nature Publ. Group}, address = {London}, issn = {2163-8306}, doi = {10.1002/psp4.12614}, pages = {564 -- 576}, year = {2021}, abstract = {The characterization of covariate effects on model parameters is a crucial step during pharmacokinetic/pharmacodynamic analyses. Although covariate selection criteria have been studied extensively, the choice of the functional relationship between covariates and parameters, however, has received much less attention. Often, a simple particular class of covariate-to-parameter relationships (linear, exponential, etc.) is chosen ad hoc or based on domain knowledge, and a statistical evaluation is limited to the comparison of a small number of such classes. Goodness-of-fit testing against a nonparametric alternative provides a more rigorous approach to covariate model evaluation, but no such test has been proposed so far. In this manuscript, we derive and evaluate nonparametric goodness-of-fit tests for parametric covariate models, the null hypothesis, against a kernelized Tikhonov regularized alternative, transferring concepts from statistical learning to the pharmacological setting. The approach is evaluated in a simulation study on the estimation of the age-dependent maturation effect on the clearance of a monoclonal antibody. Scenarios of varying data sparsity and residual error are considered. The goodness-of-fit test correctly identified misspecified parametric models with high power for relevant scenarios. The case study provides proof-of-concept of the feasibility of the proposed approach, which is envisioned to be beneficial for applications that lack well-founded covariate models.}, language = {en} } @article{FrontonPilariHuisinga2014, author = {Fronton, Ludivine and Pilari, Sabine and Huisinga, Wilhelm}, title = {Monoclonal antibody disposition: a simplified PBPK model and its implications for the derivation and interpretation of classical compartment models}, series = {Journal of pharmacokinetics and pharmacodynamics}, volume = {41}, journal = {Journal of pharmacokinetics and pharmacodynamics}, number = {2}, publisher = {Springer}, address = {New York}, issn = {1567-567X}, doi = {10.1007/s10928-014-9349-1}, pages = {87 -- 107}, year = {2014}, abstract = {The structure, interpretation and parameterization of classical compartment models as well as physiologically-based pharmacokinetic (PBPK) models for monoclonal antibody (mAb) disposition are very diverse, with no apparent consensus. In addition, there is a remarkable discrepancy between the simplicity of experimental plasma and tissue profiles and the complexity of published PBPK models. We present a simplified PBPK model based on an extravasation rate-limited tissue model with elimination potentially occurring from various tissues and plasma. Based on model reduction (lumping), we derive several classical compartment model structures that are consistent with the simplified PBPK model and experimental data. We show that a common interpretation of classical two-compartment models for mAb disposition-identifying the central compartment with the total plasma volume and the peripheral compartment with the interstitial space (or part of it)-is not consistent with current knowledge. Results are illustrated for the monoclonal antibodies 7E3 and T84.66 in mice.}, language = {en} } @article{MicheletBindelliniMelinetal.2023, author = {Michelet, Robin and Bindellini, Davide and Melin, Johanna and Neumann, Uta and Blankenstein, Oliver and Huisinga, Wilhelm and Johnson, Trevor N. and Whitaker, Martin J. and Ross, Richard and Kloft, Charlotte}, title = {Insights in the maturational processes influencing hydrocortisone pharmacokinetics in congenital adrenal hyperplasia patients using a middle-out approach}, series = {Frontiers in Pharmacology}, volume = {13}, journal = {Frontiers in Pharmacology}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {1663-9812}, doi = {10.3389/fphar.2022.1090554}, pages = {14}, year = {2023}, abstract = {Introduction: Hydrocortisone is the standard of care in cortisol replacement therapy for congenital adrenal hyperplasia patients. Challenges in mimicking cortisol circadian rhythm and dosing individualization can be overcome by the support of mathematical modelling. Previously, a non-linear mixed-effects (NLME) model was developed based on clinical hydrocortisone pharmacokinetic (PK) pediatric and adult data. Additionally, a physiologically-based pharmacokinetic (PBPK) model was developed for adults and a pediatric model was obtained using maturation functions for relevant processes. In this work, a middle-out approach was applied. The aim was to investigate whether PBPK-derived maturation functions could provide a better description of hydrocortisone PK inter-individual variability when implemented in the NLME framework, with the goal of providing better individual predictions towards precision dosing at the patient level. Methods: Hydrocortisone PK data from 24 adrenal insufficiency pediatric patients and 30 adult healthy volunteers were used for NLME model development, while the PBPK model and maturation functions of clearance and cortisol binding globulin (CBG) were developed based on previous studies published in the literature. Results: Clearance (CL) estimates from both approaches were similar for children older than 1 year (CL/F increasing from around 150 L/h to 500 L/h), while CBG concentrations differed across the whole age range (CBG(NLME) stable around 0.5 mu M vs. steady increase from 0.35 to 0.8 mu M for CBG (PBPK)). PBPK-derived maturation functions were subsequently included in the NLME model. After inclusion of the maturation functions, none, a part of, or all parameters were re-estimated. However, the inclusion of CL and/or CBG maturation functions in the NLME model did not result in improved model performance for the CL maturation function (\& UDelta;OFV > -15.36) and the re-estimation of parameters using the CBG maturation function most often led to unstable models or individual CL prediction bias. Discussion: Three explanations for the observed discrepancies could be postulated, i) non-considered maturation of processes such as absorption or first-pass effect, ii) lack of patients between 1 and 12 months, iii) lack of correction of PBPK CL maturation functions derived from urinary concentration ratio data for the renal function relative to adults. These should be investigated in the future to determine how NLME and PBPK methods can work towards deriving insights into pediatric hydrocortisone PK.}, language = {en} } @article{FuhrmannKloftHuisinga2017, author = {Fuhrmann, Saskia and Kloft, Charlotte and Huisinga, Wilhelm}, title = {Impact of altered endogenous IgG on unspecific mAb clearance}, series = {Journal of pharmacokinetics and pharmacodynamics}, volume = {44}, journal = {Journal of pharmacokinetics and pharmacodynamics}, publisher = {Springer}, address = {New York}, issn = {1567-567X}, doi = {10.1007/s10928-017-9524-2}, pages = {351 -- 374}, year = {2017}, abstract = {Immunodeficient mice are crucial models to evaluate the efficacy of monoclonal antibodies (mAbs). When studying mAb pharmacokinetics (PK), protection from elimination by binding to the neonatal Fc receptor (FcRn) is known to be a major process influencing the unspecific clearance of endogenous and therapeutic IgG. The concentration of endogenous IgG in immunodeficient mice, however is reduced, and this effect on the FcRn protection mechanism and subsequently on unspecific mAb clearance is unknown, yet of great importance for the interpretation of mAb PK data. We used a PBPK modelling approach to elucidate the influence of altered endogenous IgG concentrations on unspecific mAb clearance. To this end, we used PK data in immunodeficient mice, i.e. nude and severe combined immunodeficiency mice. To avoid impact of target-mediated clearance processes, we focussed on mAbs without affinity to a target antigen in these mice. In addition, intravenous immunoglobulin (IVIG) data of immunocompetent mice was used to study the impact of increased total IgG concentrations on unspecific therapeutic antibody clearance. The unspecific clearance is linear, whenever therapeutic IgG concentrations, i.e. mAb and IVIG concentrations are lower than FcRn; it can be non-linear if therapeutic IgG concentrations are larger than FcRn and endogenous IgG concentrations (e.g., under IVIG therapy). Unspecific mAb clearance of immunodeficient mice is effectively linear (under mAb doses as typically used in human). Studying the impact of reduced endogenous IgG concentrations on unspecific mAb clearance is of great relevance for the extrapolation to clinical species, e.g., when predicting mAb PK in immunosuppressed cancer patients.}, language = {en} } @article{MenzLatorreSchuetteetal.2012, author = {Menz, Stephan and Latorre, Juan C. and Sch{\"u}tte, Christof and Huisinga, Wilhelm}, title = {Hybrid stochastic-deterministic solution of the chemical master equation}, series = {Multiscale modeling \& simulation : a SIAM interdisciplinary journal}, volume = {10}, journal = {Multiscale modeling \& simulation : a SIAM interdisciplinary journal}, number = {4}, publisher = {Society for Industrial and Applied Mathematics}, address = {Philadelphia}, issn = {1540-3459}, doi = {10.1137/110825716}, pages = {1232 -- 1262}, year = {2012}, abstract = {The chemical master equation (CME) is the fundamental evolution equation of the stochastic description of biochemical reaction kinetics. In most applications it is impossible to solve the CME directly due to its high dimensionality. Instead, indirect approaches based on realizations of the underlying Markov jump process are used, such as the stochastic simulation algorithm (SSA). In the SSA, however, every reaction event has to be resolved explicitly such that it becomes numerically inefficient when the system's dynamics include fast reaction processes or species with high population levels. In many hybrid approaches, such fast reactions are approximated as continuous processes or replaced by quasi-stationary distributions in either a stochastic or a deterministic context. Current hybrid approaches, however, almost exclusively rely on the computation of ensembles of stochastic realizations. We present a novel hybrid stochastic-deterministic approach to solve the CME directly. Our starting point is a partitioning of the molecular species into discrete and continuous species that induces a partitioning of the reactions into discrete-stochastic and continuous-deterministic processes. The approach is based on a WKB (Wentzel-Kramers-Brillouin) ansatz for the conditional probability distribution function (PDF) of the continuous species (given a discrete state) in combination with Laplace's method of integral approximation. The resulting hybrid stochastic-deterministic evolution equations comprise a CME with averaged propensities for the PDF of the discrete species that is coupled to an evolution equation of the related expected levels of the continuous species for each discrete state. In contrast to indirect hybrid methods, the impact of the evolution of discrete species on the dynamics of the continuous species has to be taken into account explicitly. The proposed approach is efficient whenever the number of discrete molecular species is small. We illustrate the performance of the new hybrid stochastic-deterministic approach in an application to model systems of biological interest.}, language = {en} }