Significant impact of time-of-day variation on metformin pharmacokinetics.

AIMS/HYPOTHESIS: The objective was to investigate if metformin pharmacokinetics is modulated by time-of-day in humans using empirical and mechanistic pharmacokinetic modelling techniques on a large clinical dataset. This study also aimed to generate and test hypotheses on the underlying mechanisms, including evidence for chronotype-dependent interindividual differences in metformin plasma and efficacy-related tissue concentrations. METHODS: A large clinical dataset consisting of individual metformin plasma and urine measurements was analysed using a newly developed empirical pharmacokinetic model. Causes of daily variation of metformin pharmacokinetics and interindividual variability were further investigated by a literature-informed mechanistic modelling analysis. RESULTS: A significant effect of time-of-day on metformin pharmacokinetics was found. Daily rhythms of gastrointestinal, hepatic and renal processes are described in the literature, possibly affecting drug pharmacokinetics. Observed metformin plasma levels were best described by a combination of a rhythm in GFR, renal plasma flow (RPF) and organic cation transporter (OCT) 2 activity. Furthermore, the large interindividual differences in measured metformin concentrations were best explained by individual chronotypes affecting metformin clearance, with impact on plasma and tissue concentrations that may have implications for metformin efficacy. CONCLUSIONS/INTERPRETATION: Metformin's pharmacology significantly depends on time-of-day in humans, determined with the help of empirical and mechanistic pharmacokinetic modelling, and rhythmic GFR, RPF and OCT2 were found to govern intraday variation. Interindividual variation was found to be partly dependent on individual chronotype, suggesting diurnal preference as an interesting, but so-far underappreciated, topic with regard to future personalised chronomodulated therapy in people with type 2 diabetes.

Physiologically Based Pharmacokinetic (PBPK) Modeling of Clopidogrel and Its Four Relevant Metabolites for CYP2B6, CYP2C8, CYP2C19, and CYP3A4 Drug-Drug-Gene Interaction Predictions.

The antiplatelet agent clopidogrel is listed by the FDA as a strong clinical index inhibitor of cytochrome P450 (CYP) 2C8 and weak clinical inhibitor of CYP2B6. Moreover, clopidogrel is a substrate of-among others-CYP2C19 and CYP3A4. This work presents the development of a whole-body physiologically based pharmacokinetic (PBPK) model of clopidogrel including the relevant metabolites, clopidogrel carboxylic acid, clopidogrel acyl glucuronide, 2-oxo-clopidogrel, and the active thiol metabolite, with subsequent application for drug-gene interaction (DGI) and drug-drug interaction (DDI) predictions. Model building was performed in PK-Sim® using 66 plasma concentration-time profiles of clopidogrel and its metabolites. The comprehensive parent-metabolite model covers biotransformation via carboxylesterase (CES) 1, CES2, CYP2C19, CYP3A4, and uridine 5'-diphospho-glucuronosyltransferase 2B7. Moreover, CYP2C19 was incorporated for normal, intermediate, and poor metabolizer phenotypes. Good predictive performance of the model was demonstrated for the DGI involving CYP2C19, with 17/19 predicted DGI AUClast and 19/19 predicted DGI Cmax ratios within 2-fold of their observed values. Furthermore, DDIs involving bupropion, omeprazole, montelukast, pioglitazone, repaglinide, and rifampicin showed 13/13 predicted DDI AUClast and 13/13 predicted DDI Cmax ratios within 2-fold of their observed ratios. After publication, the model will be made publicly accessible in the Open Systems Pharmacology repository.

Renal Transporter-Mediated Drug-Biomarker Interactions of the Endogenous Substrates Creatinine and N1 -Methylnicotinamide: A PBPK Modeling Approach.

Endogenous biomarkers for transporter-mediated drug-drug interaction (DDI) predictions represent a promising approach to facilitate and improve conventional DDI investigations in clinical studies. This approach requires high sensitivity and specificity of biomarkers for the targets of interest (e.g., transport proteins), as well as rigorous characterization of their kinetics, which can be accomplished utilizing physiologically-based pharmacokinetic (PBPK) modeling. Therefore, the objective of this study was to develop PBPK models of the endogenous organic cation transporter (OCT)2 and multidrug and toxin extrusion protein (MATE)1 substrates creatinine and N1 -methylnicotinamide (NMN). Additionally, this study aimed to predict kinetic changes of the biomarkers during administration of the OCT2 and MATE1 perpetrator drugs trimethoprim, pyrimethamine, and cimetidine. Whole-body PBPK models of creatinine and NMN were developed utilizing studies investigating creatinine or NMN exogenous administration and endogenous synthesis. The newly developed models accurately describe and predict observed plasma concentration-time profiles and urinary excretion of both biomarkers. Subsequently, models were coupled to the previously built and evaluated perpetrator models of trimethoprim, pyrimethamine, and cimetidine for interaction predictions. Increased creatinine plasma concentrations and decreased urinary excretion during the drug-biomarker interactions with trimethoprim, pyrimethamine, and cimetidine were well-described. An additional inhibition of NMN synthesis by trimethoprim and pyrimethamine was hypothesized, improving NMN plasma and urine interaction predictions. To summarize, whole-body PBPK models of creatinine and NMN were built and evaluated to better assess creatinine and NMN kinetics while uncovering knowledge gaps for future research. The models can support investigations of renal transporter-mediated DDIs during drug development.

Novel models for the prediction of drug-gene interactions.

INTRODUCTION: Adverse drug reactions (ADRs) are among the leading causes of death, and frequently associated with drug-gene interactions (DGIs). In addition to pharmacogenomic programs for implementation of genetic preemptive testing into clinical practice, mathematical modeling can help to understand, quantify and predict the effects of DGIs in vivo. Moreover, modeling can contribute to optimize prospective clinical drug trial activities and to reduce DGI-related ADRs. AREAS COVERED: Approaches and challenges of mechanistical DGI implementation and model parameterization are discussed for population pharmacokinetic and physiologically based pharmacokinetic models. The broad spectrum of published DGI models and their applications is presented, focusing on the investigation of DGI effects on pharmacology and model-based dose adaptations. EXPERT OPINION: Mathematical modeling provides an opportunity to investigate complex DGI scenarios and can facilitate the development process of safe and efficient personalized dosing regimens. However, reliable DGI model input data from in vivo and in vitro measurements are crucial. For this, collaboration among pharmacometricians, laboratory scientists and clinicians is important to provide homogeneous datasets and unambiguous model parameters. For a broad adaptation of validated DGI models in clinical practice, interdisciplinary cooperation should be promoted and qualification toolchains must be established.

A Physiologically-Based Pharmacokinetic Model of Trimethoprim for MATE1, OCT1, OCT2, and CYP2C8 Drug-Drug-Gene Interaction Predictions.

Trimethoprim is a frequently-prescribed antibiotic and therefore likely to be co-administered with other medications, but it is also a potent inhibitor of multidrug and toxin extrusion protein (MATE) and a weak inhibitor of cytochrome P450 (CYP) 2C8. The aim of this work was to develop a physiologically-based pharmacokinetic (PBPK) model of trimethoprim to investigate and predict its drug-drug interactions (DDIs). The model was developed in PK-Sim®, using a large number of clinical studies (66 plasma concentration-time profiles with 36 corresponding fractions excreted in urine) to describe the trimethoprim pharmacokinetics over the entire published dosing range (40 to 960 mg). The key features of the model include intestinal efflux via P-glycoprotein (P-gp), metabolism by CYP3A4, an unspecific hepatic clearance process, and a renal clearance consisting of glomerular filtration and tubular secretion. The DDI performance of this new model was demonstrated by prediction of DDIs and drug-drug-gene interactions (DDGIs) of trimethoprim with metformin, repaglinide, pioglitazone, and rifampicin, with all predicted DDI and DDGI AUClast and Cmax ratios within 1.5-fold of the clinically-observed values. The model will be freely available in the Open Systems Pharmacology model repository, to support DDI studies during drug development.

A Mechanistic, Enantioselective, Physiologically Based Pharmacokinetic Model of Verapamil and Norverapamil, Built and Evaluated for Drug-Drug Interaction Studies.

The calcium channel blocker and antiarrhythmic agent verapamil is recommended by the FDA for drug-drug interaction (DDI) studies as a moderate clinical CYP3A4 index inhibitor and as a clinical Pgp inhibitor. The purpose of the presented work was to develop a mechanistic whole-body physiologically based pharmacokinetic (PBPK) model to investigate and predict DDIs with verapamil. The model was established in PK-Sim®, using 45 clinical studies (dosing range 0.1-250 mg), including literature as well as unpublished Boehringer Ingelheim data. The verapamil R- and S-enantiomers and their main metabolites R- and S-norverapamil are represented in the model. The processes implemented to describe the pharmacokinetics of verapamil and norverapamil include enantioselective plasma protein binding, enantioselective metabolism by CYP3A4, non-stereospecific Pgp transport, and passive glomerular filtration. To describe the auto-inhibitory and DDI potential, mechanism-based inactivation of CYP3A4 and non-competitive inhibition of Pgp by the verapamil and norverapamil enantiomers were incorporated based on in vitro literature. The resulting DDI performance was demonstrated by prediction of DDIs with midazolam, digoxin, rifampicin, and cimetidine, with 21/22 predicted DDI AUC ratios or Ctrough ratios within 1.5-fold of the observed values. The thoroughly built and qualified model will be freely available in the Open Systems Pharmacology model repository to support model-informed drug discovery and development.

A Comprehensive Whole-Body Physiologically Based Pharmacokinetic Drug-Drug-Gene Interaction Model of Metformin and Cimetidine in Healthy Adults and Renally Impaired Individuals.

BACKGROUND: Metformin is a widely prescribed antidiabetic BCS Class III drug (low permeability) that depends on active transport for its absorption and disposition. It is recommended by the US Food and Drug Administration as a clinical substrate of organic cation transporter 2/multidrug and toxin extrusion protein for drug-drug interaction studies. Cimetidine is a potent organic cation transporter 2/multidrug and toxin extrusion protein inhibitor. OBJECTIVE: The objective of this study was to provide mechanistic whole-body physiologically based pharmacokinetic models of metformin and cimetidine, built and evaluated to describe the metformin-SLC22A2 808G>T drug-gene interaction, the cimetidine-metformin drug-drug interaction, and the impact of renal impairment on metformin exposure. METHODS: Physiologically based pharmacokinetic models were developed in PK-Sim® (version 8.0). Thirty-nine clinical studies (dosing range 0.001-2550 mg), providing metformin plasma and urine data, positron emission tomography measurements of tissue concentrations, studies in organic cation transporter 2 polymorphic volunteers, drug-drug interaction studies with cimetidine, and data from patients in different stages of chronic kidney disease, were used to develop the metformin model. Twenty-seven clinical studies (dosing range 100-800 mg), reporting cimetidine plasma and urine concentrations, were used for the cimetidine model development. RESULTS: The established physiologically based pharmacokinetic models adequately describe the available clinical data, including the investigated drug-gene interaction, drug-drug interaction, and drug-drug-gene interaction studies, as well as the metformin exposure during renal impairment. All modeled drug-drug interaction area under the curve and maximum concentration ratios are within 1.5-fold of the observed ratios. The clinical data of renally impaired patients shows the expected increase in metformin exposure with declining kidney function, but also indicates counter-regulatory mechanisms in severe renal disease; these mechanisms were implemented into the model based on findings in preclinical species. CONCLUSIONS: Whole-body physiologically based pharmacokinetic models of metformin and cimetidine were built and qualified for the prediction of metformin pharmacokinetics during drug-gene interaction, drug-drug interaction, and different stages of renal disease. The model files will be freely available in the Open Systems Pharmacology model repository. Current guidelines for metformin treatment of renally impaired patients should be reviewed to avoid overdosing in CKD3 and to allow metformin therapy of CKD4 patients.

Physiologically Based Pharmacokinetic Models for Prediction of Complex CYP2C8 and OATP1B1 (SLCO1B1) Drug-Drug-Gene Interactions: A Modeling Network of Gemfibrozil, Repaglinide, Pioglitazone, Rifampicin, Clarithromycin and Itraconazole.

BACKGROUND: Drug-drug interactions (DDIs) and drug-gene interactions (DGIs) pose a serious health risk that can be avoided by dose adaptation. These interactions are investigated in strictly controlled setups, quantifying the effect of one perpetrator drug or polymorphism at a time, but in real life patients frequently take more than two medications and are very heterogenous regarding their genetic background. OBJECTIVES: The first objective of this study was to provide whole-body physiologically based pharmacokinetic (PBPK) models of important cytochrome P450 (CYP) 2C8 perpetrator and victim drugs, built and evaluated for DDI and DGI studies. The second objective was to apply these models to describe complex interactions with more than two interacting partners. METHODS: PBPK models of the CYP2C8 and organic-anion-transporting polypeptide (OATP) 1B1 perpetrator drug gemfibrozil (parent-metabolite model) and the CYP2C8 victim drugs repaglinide (also an OATP1B1 substrate) and pioglitazone were developed using a total of 103 clinical studies. For evaluation, these models were applied to predict 34 different DDI studies, establishing a CYP2C8 and OATP1B1 PBPK DDI modeling network. RESULTS: The newly developed models show a good performance, accurately describing plasma concentration-time profiles, area under the plasma concentration-time curve (AUC) and maximum plasma concentration (Cmax) values, DDI studies as well as DGI studies. All 34 of the modeled DDI AUC ratios (AUC during DDI/AUC control) and DDI Cmax ratios (Cmax during DDI/Cmax control) are within twofold of the observed values. CONCLUSIONS: Whole-body PBPK models of gemfibrozil, repaglinide, and pioglitazone have been built and qualified for DDI and DGI prediction. PBPK modeling is applicable to investigate complex interactions between multiple drugs and genetic polymorphisms.