Evaluation of the drug-drug interaction potential of treosulfan using a physiologically-based pharmacokinetic modelling approach.

AIMS: The aim of this work is the development of a mechanistic physiologically-based pharmacokinetic (PBPK) model using in vitro to in vivo extrapolation to conduct a drug-drug interaction (DDI) assessment of treosulfan against two cytochrome p450 (CYP) isoenzymes and P-glycoprotein (P-gp) substrates. METHODS: A PBPK model for treosulfan was developed de novo based on literature and unpublished clinical data. The PBPK DDI analysis was conducted using the U.S. Food and Drug Administration (FDA) DDI index drugs (probe substrates) midazolam, omeprazole and digoxin for CYP3A4, CYP2C19 and P-gp, respectively. Qualified and documented PBPK models of the probe substrates have been adopted from an open-source online model database. RESULTS: The PBPK model for treosulfan, based on both in vitro and in vivo data, was able to predict the plasma concentration-time profiles and exposure levels of treosulfan applied for a standard conditioning treatment. Medium and low potentials for DDI on CYP3A4 (maximum area under the concentration-time curve ratio (AUCRmax = 2.23) and CYP2C19 (AUCRmax = 1.6) were predicted, respectively, using probe substrates midazolam and omeprazole. Treosulfan was not predicted to cause a DDI on P-gp. CONCLUSION: Medicinal products with a narrow therapeutic index (eg, digoxin) that are substrates for CYP3A4, CYP2C19 or P-gp should not be given during treatment with treosulfan. However, considering the comprehensive treosulfan-based conditioning treatment schedule and the respective pharmacokinetic properties of the concomitantly used drugs (eg, half-life), the potential for interaction on all evaluated mechanisms would be low (AUCR < 1.25), if concomitantly administered drugs are dosed either 2 hours before or 8 hours after the 2-hour intravenous infusion of treosulfan.

MONALISA for stochastic simulations of Petri net models of biochemical systems.

BACKGROUND: The concept of Petri nets (PN) is widely used in systems biology and allows modeling of complex biochemical systems like metabolic systems, signal transduction pathways, and gene expression networks. In particular, PN allows the topological analysis based on structural properties, which is important and useful when quantitative (kinetic) data are incomplete or unknown. Knowing the kinetic parameters, the simulation of time evolution of such models can help to study the dynamic behavior of the underlying system. If the number of involved entities (molecules) is low, a stochastic simulation should be preferred against the classical deterministic approach of solving ordinary differential equations. The Stochastic Simulation Algorithm (SSA) is a common method for such simulations. The combination of the qualitative and semi-quantitative PN modeling and stochastic analysis techniques provides a valuable approach in the field of systems biology. RESULTS: Here, we describe the implementation of stochastic analysis in a PN environment. We extended MONALISA - an open-source software for creation, visualization and analysis of PN - by several stochastic simulation methods. The simulation module offers four simulation modes, among them the stochastic mode with constant firing rates and Gillespie's algorithm as exact and approximate versions. The simulator is operated by a user-friendly graphical interface and accepts input data such as concentrations and reaction rate constants that are common parameters in the biological context. The key features of the simulation module are visualization of simulation, interactive plotting, export of results into a text file, mathematical expressions for describing simulation parameters, and up to 500 parallel simulations of the same parameter sets. To illustrate the method we discuss a model for insulin receptor recycling as case study. CONCLUSIONS: We present a software that combines the modeling power of Petri nets with stochastic simulation of dynamic processes in a user-friendly environment supported by an intuitive graphical interface. The program offers a valuable alternative to modeling, using ordinary differential equations, especially when simulating single-cell experiments with low molecule counts. The ability to use mathematical expressions provides an additional flexibility in describing the simulation parameters. The open-source distribution allows further extensions by third-party developers. The software is cross-platform and is licensed under the Artistic License 2.0.

Predictions of Bedaquiline Central Nervous System Exposure in Patients with Tuberculosis Meningitis Using Physiologically based Pharmacokinetic Modeling.

BACKGROUND AND OBJECTIVE: The use of bedaquiline as a treatment option for drug-resistant tuberculosis meningitis (TBM) is of interest to address the increased prevalence of resistance to first-line antibiotics. To this end, we describe a whole-body physiologically based pharmacokinetic (PBPK) model for bedaquiline to predict central nervous system (CNS) exposure. METHODS: A whole-body PBPK model was developed for bedaquiline and its metabolite, M2. The model included compartments for brain and cerebrospinal fluid (CSF). Model predictions were evaluated by comparison to plasma PK time profiles following different dosing regimens and sparse CSF concentrations data from patients. Simulations were then conducted to compare CNS and lung exposures to plasma exposure at clinically relevant dosing schedules. RESULTS: The model appropriately described the observed plasma and CSF bedaquiline and M2 concentrations from patients with pulmonary tuberculosis (TB). The model predicted a high impact of tissue binding on target site drug concentrations in CNS. Predicted unbound exposures within brain interstitial exposures were comparable with unbound vascular plasma and unbound lung exposures. However, unbound brain intracellular exposures were predicted to be 7% of unbound vascular plasma and unbound lung intracellular exposures. CONCLUSIONS: The whole-body PBPK model for bedaquiline and M2 predicted unbound concentrations in brain to be significantly lower than the unbound concentrations in the lung at clinically relevant doses. Our findings suggest that bedaquiline may result in relatively inferior efficacy against drug-resistant TBM when compared with efficacy against drug-resistant pulmonary TB.

A Physiologically-Based Quantitative Systems Pharmacology Model of the Incretin Hormones GLP-1 and GIP and the DPP4 Inhibitor Sitagliptin.

Incretin hormones glucagon-like peptide-1 (GLP-1) and glucose-dependent insulinotropic polypeptide (GIP) play a major role in regulation of postprandial glucose and the development of type 2 diabetes mellitus. The incretins are rapidly metabolized, primarily by the enzyme dipeptidyl-peptidase 4 (DPP4), and the neutral endopeptidase (NEP), although the exact metabolization pathways are unknown. We developed a physiologically-based (PB) quantitative systems pharmacology model of GLP-1 and GIP and their metabolites that describes the secretion of the incretins in response to intraduodenal glucose infusions and their degradation by DPP4 and NEP. The model describes the observed data and suggests that NEP significantly contributes to the metabolization of GLP-1, and the traditional assays for the total GLP-1 and GIP forms measure yet unknown entities produced by NEP. We further extended the model with a PB pharmacokinetics/pharmacodynamics model of the DPP4 inhibitor sitagliptin that allows predictions of the effects of this medication class on incretin concentrations.

A Quantitative Systems Pharmacology Kidney Model of Diabetes Associated Renal Hyperfiltration and the Effects of SGLT Inhibitors.

The early stage of diabetes mellitus is characterized by increased glomerular filtration rate (GFR), known as hyperfiltration, which is believed to be one of the main causes leading to renal injury in diabetes. Sodium-glucose cotransporter 2 inhibitors (SGLT2i) have been shown to be able to reverse hyperfiltration in some patients. We developed a mechanistic computational model of the kidney that explains the interplay of hyperglycemia and hyperfiltration and integrates the pharmacokinetics/pharmacodynamics (PK/PD) of the SGLT2i dapagliflozin. Based on simulation results, we propose kidney growth as the necessary process for hyperfiltration progression. Further, the model indicates that renal SGLT1i could significantly improve hyperfiltration when added to SGTL2i. Integrated into a physiologically based PK/PD (PBPK/PD) Diabetes Platform, the model presents a powerful tool for aiding drug development, prediction of hyperfiltration risk, and allows the assessment of the outcomes of individualized treatments with SGLT1-inhibitors and SGLT2-inhibitors and their co-administration with insulin.