Investigating clinical pharmacokinetics of brivaracetam by using a pharmacokinetic modeling approach

The development of technology and the processing speed of computing machines have facilitated the evaluation of advanced pharmacokinetic (PK) models, making modeling processes simple and faster. The present model aims to analyze the PK of brivaracetam (BRV) in healthy and diseased populations. A comprehensive literature review was conducted to incorporate the BRV plasma concentration data and its input parameters into PK-Sim software, leading to the creation of intravenous (IV) and oral models for both populations. The developed physiologically based pharmacokinetic (PBPK) model of BRV was then assessed using the visual predictive checks, mean observed/predicted ratios (Robs/pre), and average fold error for PK parameters including the maximum systemic concentration (Cmax), the area under the curve at time 0 to t (AUC0-∞), and drug clearance (CL). The PBPK model of BRV demonstrated that mean Robs/pre ratios of the PK parameters remained within the acceptable limits when assessed against a twofold error margin. Furthermore, model predictions were carried out to assess how AUC0-∞ is affected following the administration of BRV in individuals with varying degrees of liver cirrhosis, ranging from different child-pugh (CP) scores like A, B, and C. Moreover, dose adjustments were recommended by considering the variations in Cmax and CL in various kidney disease stages (mild to severe).

Physiologically-based pharmacokinetic modelling guided dose evaluations of nirmatrelvir/ritonavir in renal impairment for the management of COVID-19

We aimed to address factors contributing to the pharmacokinetic changes of nirmatrelvir/ritonavir in renal impaired (RI) patients and recommend dosing adjustment via a physiologically-based pharmacokinetic (PBPK) modelling approach. A PBPK model of nirmatrelvir/ritonavir was developed via Simcyp® Simulator. Sensitivity analysis of the influence of hepatic CYP3A4 intrinsic clearance and abundance, as well as hepatic non-CYP3A4 metabolism (other human liver microsomes [HLM] CLint) was performed to evaluate the effects of RI on oral clearance of nirmatrelvir. Other HLM CLint, the most sensitive parameter, was adjusted, and the simulated plasma concentration profiles of nirmatrelvir in severe RI subjects were within the therapeutic index of 292-10 000 ng/mL for dosing regimens of loading doses of 300/100 mg followed by 150/100 mg or 75/100 mg twice daily of nirmatrelvir/ritonavir. Considering that nirmatrelvir is available as a 150 mg tablet, we recommend 300/100 mg followed by 150/100 mg twice daily as the dosing regimen to be investigated in severe RI.

The Effect of Inflammatory Bowel Disease and Irritable Bowel Syndrome on Pravastatin Oral Bioavailability: In vivo and in silico evaluation using bottom-up wbPBPK modeling

The common disorders irritable bowel syndrome (IBS) and inflammatory bowel disease (IBD) can modify the drugs' pharmacokinetics via their induced pathophysiological changes. This work aimed to investigate the impact of these two diseases on pravastatin oral bioavailability. Rat models for IBS and IBD were used to experimentally test the effects of IBS and IBD on pravastatin pharmacokinetics. Then, the observations made in rats were extrapolated to humans using a mechanistic whole-body physiologically-based pharmacokinetic (wbPBPK) model. The rat in vivo studies done herein showed that IBS and IBD decreased serum albumin (> 11% for both), decreased PRV binding in plasma, and increased pravastatin absolute oral bioavailability (0.17 and 0.53 compared to 0.01) which increased plasma, muscle, and liver exposure. However, the wbPBPK model predicted muscle concentration was much lower than the pravastatin toxicity thresholds for myotoxicity and rhabdomyolysis. Overall, IBS and IBD can significantly increase pravastatin oral bioavailability which can be due to a combination of increased pravastatin intestinal permeability and decreased pravastatin gastric degradation resulting in higher exposure. This is the first study in the literature investigating the effects of IBS and IBD on pravastatin pharmacokinetics. The high interpatient variability in pravastatin concentrations as induced by IBD and IBS can be reduced by oral administration of pravastatin using enteric-coated tablets. Such disease (IBS and IBD)-drug interaction can have more drastic consequences for narrow therapeutic index drugs prone to gastric degradation, especially for drugs with low intestinal permeability.

Pan-genome mediated therapeutic target mining in Kingella kingae and inhibition assessment using traditional Chinese medicinal compounds: an informatics approach

Kingella kingae causes bacteremia, endocarditis, osteomyelitis, septic arthritis, meningitis, spondylodiscitis, and lower respiratory tract infections in pediatric patients. Usually it demonstrates disease after inflammation of mouth, lips or infections of the upper respiratory tract. To date, therapeutic targets in this bacterium remain unexplored. We have utilized a battery of bioinformatics tools to mine these targets in this study. Core genes were initially inferred from 55 genomes of K. kingae and 39 therapeutic targets were mined using an in-house pipeline. We selected aroG product (KDPG aldolase) involved in chorismate pathway, for inhibition analysis of this bacterium using lead-like metabolites from traditional Chinese medicinal plants. Pharmacophore generation was done using control ZINC36444158 (1,16-bis[(dihydroxyphosphinyl)oxy]hexadecane), followed by molecular docking of top hits from a library of 36,000 compounds. Top prioritized compounds were ZINC95914016, ZINC33833283 and ZINC95914219. ADME profiling and simulation of compound dosing (100 mg tablet) was done to infer compartmental pharmacokinetics in a population of 300 individuals in fasting state. PkCSM based toxicity analysis revealed the compounds ZINC95914016 and ZINC95914219 as safe and with almost similar bioavailability. However, ZINC95914016 takes less time to reach maximum concentration in the plasma and shows several optimal parameters compared to other leads. In light of obtained data, we recommend this compound for further testing and induction in experimental drug design pipeline.Communicated by Ramaswamy H. Sarma.

A Comprehensive Physiologically Based Pharmacokinetic Model of Nadolol in Adults with Renal Disease and Pediatrics with Supraventricular Tachycardia

Nadolol is a long-acting non-selective β-adrenergic antagonist that helps treat angina and hypertension. The current study aimed to develop and validate the physiologically based pharmacokinetic model (PBPK) of nadolol in healthy adults, renal-compromised, and pediatric populations. A comprehensive PBPK model was established by utilizing a PK-Sim simulator. After establishing and validating the model in healthy adults, pathophysiological changes i.e., blood flow, hematocrit, and GFR that occur in renal failure were incorporated in the developed model, and the drug exposure was assessed through Box plots. The pediatric model was also developed and evaluated by considering the renal maturation process. The validation of the models was carried out by visual predictive checks, calculating predicted to observed (Rpre/obs) and the average fold error (AFE) of PK parameters i.e., the area under the concentration-time curve (AUC0-t), the maximum concentration in plasma (Cmax), and CL (clearance). The presented PBPK model successfully simulates the nadolol PK in healthy adults, renal-impaired, and pediatric populations, as the Rpre/obs values of all PK parameters fall within the acceptable range. The established PBPK model can be useful in nadolol dose optimization in patients with renal failure and children with supraventricular tachycardia.

Chronic kidney disease and physiologically based pharmacokinetic modeling: a critical review of existing models

INTRODUCTION

Physiologically based pharmacokinetic (PBPK) modeling is a paradigm shift in this era for determining the exposure of drugs in pediatrics, geriatrics, and patients with chronic diseases where clinical trials are difficult to conduct.

AREAS COVERED

This review has collated data regarding published PBPK models on chronic kidney disease (CKD), including the drug and system-specific input model parameters and model evaluation criteria. Four databases were used from 13th June 2023 to 10th July 2023 for identifying the relevant studies that met the inclusion/exclusion criteria. Alterations in plasma protein (albumin/alpha-1 acid glycoprotein), gastric emptying time, hematocrit, small intestinal transit time, the abundance of cytochrome (CYP) 450 enzymes, glomerular filtration rate, and physicochemical parameters for different drugs were explicitly elaborated from earlier reported studies. Moreover, model evaluation depicted that models in CKD for most of the included drugs were within the allowed two-fold error range.

EXPERT OPINION

This review will provide insights for researchers on applying PBPK models in managing patients with different levels of CKD to prevent undesirable side effects and increase the effectiveness of drug therapy.

Physiologically based pharmacokinetic modeling of candesartan to predict the exposure in hepatic and renal impairment and elderly populations

Background

Candesartan cilexetil is a widely used angiotensin II receptor blocker with minimal adverse effects and high tolerability for the treatment of hypertension. Candesartan is administered orally as the prodrug candesartan cilexetil, which is wholly and swiftly converted to the active metabolite candesartan by carboxylesterase during absorption in the intestinal tract. In populations with renal or hepatic impairment, candesartan's pharmacokinetic (PK) behavior may be altered, necessitating dosage adjustments.

Objectives

This study was conducted to examine how the physiologically based PK (PBPK) model characterizes the PKs of candesartan in adult and geriatric populations and to predict the PKs of candesartan in elderly populations with renal and hepatic impairment.

Design

After developing PBPK models using the reported physicochemical properties of candesartan and clinical data, these models were validated using data from clinical investigations involving various dose ranges.

Methods

Comparing predicted and observed blood concentration data and PK parameters was used to assess the fit performance of the models.

Results

Doses should be reduced to approximately 94% of Chinese healthy adults for the Chinese healthy elderly population; approximately 92%, 68%, and 64% of that of the Chinese healthy adult dose in elderly populations with mild, moderate, and severe renal impairment, respectively; and approximately 72%, 71%, and 52% of that of the Chinese healthy adult dose in elderly populations with Child-Pugh-A, Child-Pugh-B, and Child-Pugh-C hepatic impairment, respectively.

Conclusion

The results suggest that the PBPK model of candesartan can be utilized to optimize dosage regimens for special populations.

The Trends and Future Prospective of In Silico Models from the Viewpoint of ADME Evaluation in Drug Discovery

Drug discovery and development are aimed at identifying new chemical molecular entities (NCEs) with desirable pharmacokinetic profiles for high therapeutic efficacy. The plasma concentrations of NCEs are a biomarker of their efficacy and are governed by pharmacokinetic processes such as absorption, distribution, metabolism, and excretion (ADME). Poor ADME properties of NCEs are a major cause of attrition in drug development. ADME screening is used to identify and optimize lead compounds in the drug discovery process. Computational models predicting ADME properties have been developed with evolving model-building technologies from a simplified relationship between ADME endpoints and physicochemical properties to machine learning, including support vector machines, random forests, and convolution neural networks. Recently, in the field of in silico ADME research, there has been a shift toward evaluating the in vivo parameters or plasma concentrations of NCEs instead of using predictive results to guide chemical structure design. Another research hotspot is the establishment of a computational prediction platform to strengthen academic drug discovery. Bioinformatics projects have produced a series of in silico ADME models using free software and open-access databases. In this review, we introduce prediction models for various ADME parameters and discuss the currently available academic drug discovery platforms.

Physiologically Based Pharmacokinetic Model To Predict Metoprolol Disposition in Healthy and Disease Populations

The evolution in the development of drugs has increased the popularity of physiologically based pharmacokinetic (PBPK) models. This study seeks to assess the PK of metoprolol in populations with healthy, chronic kidney disease (CKD), and acute myocardial infarction (AMI) conditions by developing and evaluating PBPK models. An extensive literature review for identifying and selecting plasma concentration vs time profile data and other drug-related parameters was undergone for their integration into the PK-Sim program followed by the development of intravenous, oral, and diseased models. The developed PBPK model of metoprolol was then evaluated using the visual predictive checks, mean observed/predicted ratios (Robs/pre), and average fold error for all PK parameters, i.e., the area under the curve (AUC), maximal plasma concentration, and clearance. The model evaluation depicted that none of the PK parameters were out of the allowed range (2-fold error) in the case of the mean Robs/pre ratios. The model anticipations were executed to determine the influence of diseases on unbound and total AUC after the application of metoprolol in healthy, moderate, and severe CKD. The dosage reductions were also suggested based on differences in unbound and total AUC in different stages of CKD. The developed PBPK models have successfully elaborated the PK changes of metoprolol occurring in healthy individuals and those with renal and heart diseases (CKD & AMI), which may be fruitful for dose optimization among diseased patients.

A physiologically based pharmacokinetic model of cefepime to predict its pharmacokinetics in healthy, pediatric and disease populations

The physiologically based pharmacokinetic modeling (PBPK) approach can predict drug pharmacokinetics (PK) by combining changes in blood flow and pathophysiological alterations for developing drug-disease models. Cefepime hydrochloride is a parenteral cephalosporin that is used to treat pneumonia, sepsis, and febrile neutropenia, among other things. The current study sought to identify the factors that impact cefepime pharmacokinetics (PK) following dosing in healthy, diseased (CKD and obese), and pediatric populations. For model construction and simulation, the modeling tool PK-SIM was utilized. Estimating cefepime PK following intravenous (IV) application in healthy subjects served as the primary step in the model-building procedure. The prediction of cefepime PK in chronic kidney disease (CKD) and obese populations were performed after the integration of the relevant pathophysiological changes. Visual predictive checks and a comparison of the observed and predicted values of the PK parameters were used to verify the developed model. The results of the PK parameters were consistent with the reported clinical data in healthy subjects. The developed PBPK model successfully predicted cefepime PK as observed from the ratio of the observed and predicted PK parameters as they were within a two-fold error range.

Physiologically-based pharmacokinetic modeling for single and multiple dosing regimens of ceftriaxone in healthy and chronic kidney disease populations: a tool for model-informed precision dosing

Introduction: Ceftriaxone is one of commonly prescribed beta-lactam antibiotics with several label and off-label clinical indications. A high fraction of administered dose of ceftriaxone is excreted renally in an unchanged form, and it may accumulate significantly in patients with impaired renal functions, which may lead to toxicity. Methods: In this study, we employed a physiologically-based pharmacokinetic (PBPK) modeling, as a tool for precision dosing, to predict the biological exposure of ceftriaxone in a virtually-constructed healthy and chronic kidney disease patient populations, with subsequent dosing optimizations. We started developing the model by integrating the physicochemical properties of the drug with biological system information in a PBPK software platform. A PBPK model in an adult healthy population was developed and evaluated visually and numerically with respect to experimental pharmacokinetic data. The model performance was evaluated based on the fold error criteria of the predicted and reported values for different pharmacokinetic parameters. Then, the model was applied to predict drug exposure in CKD patient populations with various degrees of severity. Results: The developed PBPK model was able to precisely describe the pharmacokinetic behavior of ceftriaxone in adult healthy population and in mild, moderate, and severe CKD patient populations. Decreasing the dose by approximately 25% in mild and 50% in moderate to severe renal disease provided a comparable exposure to the healthy population. Based on the simulation of multiple dosing regimens in severe CKD population, it has been found that accumulation of 2 g every 24 h is lower than the accumulation of 1 g every 12 h dosing regimen. Discussion: In this study, the observed concentration time profiles and pharmacokinetic parameters for ceftriaxone were successfully reproduced by the developed PBPK model and it has been shown that PBPK modeling can be used as a tool for precision dosing to suggest treatment regimens in population with renal impairment.

Physiologically Based Pharmacokinetic Modelling to Predict Imatinib Exposures in Cancer Patients with Renal Dysfunction: A Case Study

Imatinib is mainly metabolised by CYP3A4 and CYP2C8 and is extensively bound to α-acid glycoprotein (AAG). A physiologically based pharmacokinetic (PBPK) model for imatinib describing the CYP3A4-mediated autoinhibition during multiple dosing in gastrointestinal stromal tumor patients with normal renal function was previously reported. After performing additional verification, the PBPK model was applied to predict the exposure of imatinib after multiple dosing in cancer patients with varying degrees of renal impairment. In agreement with the clinical data, there was a positive correlation between AAG levels and imatinib exposure. A notable finding was that for recovery of the observed data in cancer patients with moderate RI (CrCL 20 to 39 mL/min), reductions of hepatic CYP3A4 and CYP2C8 abundances, which reflect the effects of RI, had to be included in the simulations. This was not the case for mild RI (CrCL 40 to 50 mL/min). The results support the finding of the clinical study, which demonstrated that both AAG levels and the degree of renal impairment are key components that contribute to the interpatient variability associated with imatinib exposure. As indicated in the 2020 FDA draft RI guidance, PBPK modelling could be used to support an expanded inclusion of patients with RI in clinical studies.

Predicting Hydroxychloroquine Clearance in Healthy and Diseased Populations Using a Physiologically Based Pharmacokinetic Approach

Hydroxychloroquine (HCQ), a congener of chloroquine, is widely used in prophylaxis and the treatment of malaria, and also as a cure for rheumatoid arthritis, systemic lupus erythematosus, and various other diseases. Physiologically based pharmacokinetic modeling (PBPK) has attracted great interest in the past few years in predicting drug pharmacokinetics (PK). This study focuses on predicting the PK of HCQ in the healthy population and extrapolating it to the diseased populations, i.e., liver cirrhosis and chronic kidney disease (CKD), utilizing a systematically built whole-body PBPK model. The time vs. concentration profiles and drug-related parameters were obtained from the literature after a laborious search and in turn were integrated into PK-Sim software for designing healthy intravenous, oral, and diseased models. The model's evaluation was performed using observed-to-predicted ratios (Robs/Rpre) and visual predictive checks within a 2-fold error range. The healthy model was then extrapolated to liver cirrhosis and CKD populations after incorporating various disease-specific pathophysiological changes. Box-whisker plots showed an increase in AUC_0-t in liver cirrhosis, whereas a decrease in AUC_0-t was seen in the CKD population. These model predictions may assist clinicians in adjusting the administered HCQ doses in patients with different degrees of hepatic and renal impairment.

The role of roxadustat in chronic kidney disease patients complicated with anemia

The incidence of chronic kidney disease (CKD) is increasing worldwide and the current prevalence rate is 13.4%. There are > 120 million CKD patients in China and this number is expected to increase. One of the main abnormalities in patients with CKD and kidney impairment is decreased synthesis of erythropoietin (EPO), which causes anemia and affects iron metabolism. The probability of developing is higher in anemia patients with CKD than in the general population, and the incidence increases as kidney function decreases. Deficient EPO production by the kidney is the most important cause of renal anemia. Notably, anemia in patients with CKD has multiple causes, such as bleeding caused by platelet dysfunction, iron deficiency due to digestive and absorption disorders of the gastrointestinal tract, and shorter red blood cell life. Anemia is also a leading cause of hospitalization in patients with CKD. A new oral medication to treat renal anemia, the hypoxia-inducible factor prolyl hydroxylase inhibitor called roxadustat (FG-4592), regulates iron metabolism and promotes erythropoiesis. This drug has a therapeutic effect on patients with CKD. Roxadustat showed advantages over EPO in clinical experiments. This review summarizes the mechanisms of action, clinical applications, effectiveness, and safety of roxadustat.

The investigation of the complex population-drug-drug interaction between ritonavir-boosted lopinavir and chloroquine or ivermectin using physiologically-based pharmacokinetic modeling

OBJECTIVES

Therapy failure caused by complex population-drug-drug (PDDI) interactions including CYP3A4 can be predicted using mechanistic physiologically-based pharmacokinetic (PBPK) modeling. A synergy between ritonavir-boosted lopinavir (LPVr), ivermectin, and chloroquine was suggested to improve COVID-19 treatment. This work aimed to study the PDDI of the two CYP3A4 substrates (ivermectin and chloroquine) with LPVr in mild-to-moderate COVID-19 adults, geriatrics, and pregnancy populations.

METHODS

The PDDI of LPVr with ivermectin or chloroquine was investigated. Pearson's correlations between plasma, saliva, and lung interstitial fluid (ISF) levels were evaluated. Target site (lung epithelial lining fluid [ELF]) levels of ivermectin and chloroquine were estimated.

RESULTS

Upon LPVr coadministration, while the chloroquine plasma levels were reduced by 30, 40, and 20%, the ivermectin plasma levels were increased by a minimum of 425, 234, and 453% in adults, geriatrics, and pregnancy populations, respectively. The established correlation equations can be useful in therapeutic drug monitoring (TDM) and dosing regimen optimization.

CONCLUSIONS

Neither chloroquine nor ivermectin reached therapeutic ELF levels in the presence of LPVr despite reaching toxic ivermectin plasma levels. PBPK modeling, guided with TDM in saliva, can be advantageous to evaluate the probability of reaching therapeutic ELF levels in the presence of PDDI, especially in home-treated patients.

Renin-angiotensin-aldosterone pathway modulators in chronic kidney disease: A comparative review

Chronic kidney disease presents a health challenge that has a complex underlying pathophysiology, both acquired and inherited. The pharmacotherapeutic treatment options available today lower the progression of the disease and improve the quality of life but cannot completely cure it. This poses a challenge to the healthcare provider to choose, from the available options, the best way to manage the disease as per the presentation of the patient. As of now, the recommended first line of treatment to control the blood pressure in chronic kidney disease is the administration of renin-angiotensin-aldosterone system modulators. These are represented mainly by the direct renin inhibitor, angiotensin-converting enzyme inhibitors, and angiotensin II receptor blockers. These modulators are varied in their structure and mechanisms of action, hence showing varying treatment outcomes. The choice of administration of these modulators is determined by the presentation and the co-morbidities of the patient, the availability and affordability of the treatment option, and the expertise of the healthcare provider. A direct head-to-head comparison between these significant renin-angiotensin-aldosterone system modulators is lacking, which can benefit healthcare providers and researchers. In this review, a comparison has been drawn between the direct renin inhibitor (aliskiren), angiotensin-converting enzyme inhibitors, and angiotensin II receptor blockers. This can be of significance for healthcare providers and researchers to find the particular loci of interest, either in structure or mechanism, and to intervene as per the case presentation to obtain the best possible treatment option.

Evaluation of Renal Impairment Influence on Metabolic Drug Clearance using a Modelling Approach

BACKGROUND AND OBJECTIVE

Chronic kidney disease (CKD) may alter drug renal elimination but is also known for interacting with hepatic metabolism via multiple uremic components. However, few global models, considering the five major cytochromes, have been published, and none specifically address the decrease in cytochrome P450 (CYP450) activity. The aim of our study was to estimate the possibility of quantifying residual cytochrome activity as a function of filtration rate, according to the data available in the literature.

METHODS

For each drug in the DDI-predictor database, we collected available pharmacokinetic data comparing drug exposition in the healthy patient and in various stages of CKD, before building a model capable of predicting the variation of exposure according to the degree of renal damage. We followed an In vivo Mechanistic Static Model (IMSM) approach, previously validated for predicting change in liver clearance. We estimated the remaining fraction parameters at glomerular filtration rate (GFR) = 0 and the alpha value of GFR to 50% impairment for the 5 major cytochromes using a non-linear constrained regression using Matlab software.

RESULTS

Thirty-one compounds had usable pharmacokinetic data, with 51 AUC ratios between healthy and renal impaired patients. The remaining CYP3A4 activity was estimated to be 0.4 when CYP2D6, 2C9, 2C19 and 1A2 activity was estimated to be 0.43; 1; 0.73 and 0.7, respectively. The alpha value was estimated to be at 6.62; 25; 9.8; 1.38 and 11.04 for each cytochrome. In comparison with published data, all estimates but one were correctly predicted in the range of 0.5-2.

CONCLUSION

Our approach was able to describe the impact of CKD on metabolic elimination. Modelling this process makes it possible to anticipate changes in clearance and drug exposure in CKD patients, with the advantage of greater simplicity than approaches based on physiologically-based pharmacokinetic modelling. However, a precise estimation of the impact of renal failure is not possible with an IMSM approach due to the large variability of the published data, and thus should rely on specific pharmacokinetic modelling for narrow therapeutic margin drugs.

Development and Evaluation of a Physiologically Based Pharmacokinetic Model of Labetalol in Healthy and Diseased Populations

Labetalol is a drug that exhibits both alpha and beta-adrenergic receptor-blocking properties. The American Heart Association/American Stroke Association (AHA/ASA) has recommended labetalol as an initial treatment option for the management of severe hypertension. The physiologically based pharmacokinetic (PBPK) model is an in silico approach to determining the pharmacokinetics (PK) of a drug by incorporating blood flow and tissue composition of the organs. This study was conducted to evaluate the primary reasons for the difference in PK after intravenous (IV) and oral administration in healthy and diseased (renal and hepatic) populations. A comprehensive literature search was done using two databases, PubMed and Google Scholar. Various PK parameters were screened for the development of the PBPK model utilizing a population-based PK-Sim simulator. Simulations were performed after creating building blocks firstly in healthy individuals and then in diseased patients after IV and oral administration. The disposition of labetalol after IV and oral administration occurring in patients with the hepatic and renal disease was predicted. The model was evaluated by calculating the Robs/pred ratio and average fold error (AFE), which was in the two-fold error range. Moreover, Box-whisker plots were made to compare the overall concentration of the drug in the body at various stages of disease severity. The presented model provides useful quantitative estimates of drug dosing in patients fighting against severe chronic diseases.

Determining the Effects of Chronic Kidney Disease on Organic Anion Transporter1/3 Activity Through Physiologically Based Pharmacokinetic Modeling

BACKGROUND AND OBJECTIVE

The renal excretion of drugs via organic anion transporters 1 and 3 (OAT1/3) is significantly decreased in patients with renal impairment. This study uses physiologically based pharmacokinetic models to quantify the reduction in OAT1/3-mediated secretion of drugs throughout varying stages of chronic kidney disease.

METHODS

Physiologically based pharmacokinetic models were constructed for four OAT1/3 substrates in healthy individuals: acyclovir, meropenem, furosemide, and ciprofloxacin. Observed data from drug-drug interaction studies with probenecid, a potent OAT1/3 inhibitor, were used to parameterize the contribution of OAT1/3 to the renal elimination of each drug. The models were then translated to patients with chronic kidney disease by accounting for changes in glomerular filtration rate, kidney volume, renal blood flow, plasma protein binding, and hematocrit. Additionally, a relationship was derived between the estimated glomerular filtration rate and the reduction in OAT1/3-mediated secretion of drugs based on the renal extraction ratios of ƿ-aminohippuric acid in patients with varying degrees of renal impairment. The relationship was evaluated in silico by evaluating the predictive performance of each final model in describing the pharmacokinetics (PK) of drugs across stages of chronic kidney disease.

RESULTS

OAT1/3-mediated renal excretion of drugs was found to be decreased by 27-49%, 50-68%, and 70-96% in stage 3, stage 4, and stage 5 of chronic kidney disease, respectively. In support of the parameterization, physiologically based pharmacokinetic models of four OAT1/3 substrates were able to adequately characterize the PK in patients with different degrees of renal impairment. Total exposure after intravenous administration was predicted within a 1.5-fold error and 85% of the observed data points fell within a 1.5-fold prediction error. The models modestly under-predicted plasma concentrations in patients with end-stage renal disease undergoing intermittent hemodialysis. However, results should be interpreted with caution because of the limited number of molecules analyzed and the sparse sampling in observed chronic kidney disease pharmacokinetic studies.

CONCLUSIONS

A quantitative understanding of the reduction in OAT1/3-mediated excretion of drugs in differing stages of renal impairment will contribute to better predictive accuracy for physiologically based pharmacokinetic models in drug development, assisting with clinical trial planning and potentially sparing this population from unnecessary toxic exposures.

Coproporphyrin I as an Endogenous Biomarker to Detect Reduced OATP1B Activity and Shift in Elimination Route in Chronic Kidney Disease

Coproporphyrin I (CPI) is an endogenous biomarker of organic anion transporting polypeptide 1B transporter (OATP1B). CPI plasma baseline was reported to increase with severity of chronic kidney disease (CKD). Further, ratio of CPI area under the plasma concentration-time curve (AUCR) in the presence/absence of OATP1B inhibitor rifampin was higher in patients with CKD compared with healthy participants, in contrast to pitavastatin (a clinical OATP1B probe). This study investigated mechanism(s) contributing to altered CPI baseline in patients with CKD by extending a previously developed physiologically-based pharmacokinetic (PBPK) model to this patient population. CKD-related covariates were evaluated in a stepwise manner on CPI fraction unbound in plasma (f_u,p ), OATP1B-mediated hepatic uptake clearance (CL_active ), renal clearance (CL_R ), and endogenous synthesis (k_syn ). The CPI model successfully recovered increased baseline and rifampin-mediated AUCR in patients with CKD by accounting for the following disease-related changes: 13% increase in f_u,p , 29% and 39% decrease in CL_active in mild and moderate to severe CKD, respectively, decrease in CL_R proportional to decline in glomerular filtration rate, and 27% decrease in k_syn in severe CKD. Almost complete decline in CPI renal elimination in severe CKD increased its fraction transported by OATP1B, rationalizing differences in the CPI-rifampin interaction observed between healthy participants and patients with CKD. In conclusion, mechanistic modeling performed here supports CKD-related decrease in OATP1B function to inform prospective PBPK modeling of OATP1B-mediated drug-drug interaction in these patients. Monitoring of CPI allows detection of CKD-drug interaction risk for OATP1B drugs with combined hepatic and renal elimination which may be underestimated by extrapolating the interaction risk based on pitavastatin data in healthy participants.

Enhancing Atorvastatin In Vivo Oral Bioavailability in the Presence of Inflammatory Bowel Disease and Irritable Bowel Syndrome Using Supercritical Fluid Technology Guided by wbPBPK Modeling in Rat and Human

Inflammatory bowel disease (IBD) and irritable bowel syndrome (IBS) are common disorders that can change the body's physiology and drugs pharmacokinetics. Solid dispersion (SD) preparation using supercritical fluid technology (SFT) has many advantages. Our study aimed to explore the effect of IBS and IBD on atorvastatin (ATV) pharmacokinetics, enhance ATV oral bioavailability (BCS II drug) using SFT, and analyze drug-disease-formulation interaction using a whole-body physiologically based pharmacokinetic (wbPBPK) model in rat and human. A novel ATV formulation was prepared using SFT and characterized in vitro and in vivo in healthy, IBS, and IBD rats. The resulting ATV plasma levels were analyzed using a combination of conventional and wbPBPK approaches. The novel formulation increased ATV solubility by 20-fold and resulted in a zero-order release of up to 95%. Both IBS and IBD increased ATV exposure after oral and intravenous administration by more than 30%. The novel SFT formulation increased ATV bioavailability by 28, 14, and 18% in control, IBD, and IBD rat groups and resulted in more consistent exposure as compared to raw ATV solution. Higher improvements in ATV bioavailability of more than 2-fold upon receiving the novel SFT formulation were predicted by the human wbPBPK model as compared to receiving the conventional tablets. Finally, the established wbPBPK model could describe ATV ADME in the presence of IBS and IBD after oral administration of raw ATV and using the novel SFT formula and can help scale the optimized ATV dosing regimens in the presence of IBS and IBD from rats to humans.

The Development of a PBPK Model for Atomoxetine Using Levels in Plasma, Saliva and Brain Extracellular Fluid in Patients with Normal and Deteriorated Kidney Function

BACKGROUND

Atomoxetine is a treatment for attention-deficit hyperactivity disorder. It inhibits Norepinephrine Transporters (NET) in the brain. Renal impairment can reduce hepatic CYP2D6 activity and atomoxetine elimination which may increase its body exposure. Atomoxetine can be secreted in saliva.

OBJECTIVE

The objective of this work was to test the hypothesis that atomoxetine saliva levels (sATX) can be used to predict ATX brain Extracellular Fluid (bECF) levels and their pharmacological effects in healthy subjects and those with End-Stage Renal Disease (ESRD).

METHODS

The pharmacokinetics of atomoxetine after intravenous administration to rats with chemically induced acute and chronic renal impairments were investigated. A physiologically-based pharmacokinetic (PBPK) model was built and verified in rats using previously published measured atomoxetine levels in plasma and brain tissue. The rat PBPK model was then scaled to humans and verified using published measured atomoxetine levels in plasma, saliva, and bECF.

RESULTS

The rat PBPK model predicted the observed reduced atomoxetine clearance due to renal impairment in rats. The PBPK model predicted atomoxetine exposure in human plasma, sATX and bECF. Additionally, it predicted that ATX bECF levels needed to inhibit NET are achieved at 80 mg dose. In ESRD patients, the developed PBPK model predicted that the previously reported 65% increase in plasma exposure in these patients can be associated with a 63% increase in bECF. The PBPK simulations showed that there is a significant correlation between sATX and bECF in human.

CONCLUSION

Saliva levels can be used to predict atomoxetine pharmacological response.

Effect of Chronic Kidney Disease on the Renal Secretion via Organic Anion Transporters 1/3: Implications for Physiologically-Based Pharmacokinetic Modeling and Dose Adjustment

There is growing evidence that active tubular secretory clearance (CL_s) may not decline proportionally with the glomerular filtration rate (GFR) in chronic kidney disease (CKD), leading to the overestimation of renal clearance (CL_r) when using solely GFR to approximate disease effect on renal elimination. The clinical pharmacokinetic data of 33 renally secreted OAT1/3 substrates were collated to investigate the impact of mild, moderate, and severe CKD on CL_r, tubular secretion and protein binding (f_u,p). The f_u,p of the collated substrates ranged from 0.0026 to 1.0 in healthy populations; observed CKD‐related increase in the f_u,p (up to 2.7‐fold) of 8 highly bound substrates (f_u,p ≤ 0.2) was accounted for in the analysis. Use of prediction equation based on disease‐related changes in albumin resulted in underprediction of the CKD‐related increase in f_u,p of highly bound substrates, highlighting the necessity to measure protein binding in severe CKD. The critical analysis of clinical data for 33 OAT1/3 probes established that decrease in OAT1/3 activity proportional to the changes in GFR was insufficient to recapitulate effects of severe CKD on unbound tubular secretion clearance. OAT1/3‐mediated CL_s was estimated to decline by an additional 50% relative to the GFR decline in severe CKD, whereas change in active secretion in mild and moderate CKD was proportional to GFR. Consideration of this additional 50% decline in OAT1/3‐mediated CL_s is recommended for physiologically‐based pharmacokinetic models and dose adjustment of OAT1/3 substrates in severe CKD, especially for substrates with high contribution of the active secretion to CL_r.

Physiologically-Based Pharmacokinetic Modeling Approaches for Patients with SARS-CoV-2 Infection: A Case Study with Imatinib

Severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) infection, which causes coronavirus disease 2019 (COVID‐19), manifests as mild respiratory symptoms to severe respiratory failure and is associated with inflammation and other physiological changes. Of note, substantial increases in plasma concentrations of α₁‐acid‐glycoprotein and interleukin‐6 have been observed among patients admitted to the hospital with advanced SARS‐CoV‐2 infection. A physiologically based pharmacokinetic (PBPK) approach is a useful tool to evaluate and predict disease‐related changes on drug pharmacokinetics. A PBPK model of imatinib has previously been developed and verified in healthy people and patients with cancer. In this study, the PBPK model of imatinib was successfully extrapolated to patients with SARS‐CoV‐2 infection by accounting for disease‐related changes in plasma α₁‐acid‐glycoprotein concentrations and the potential drug interaction between imatinib and dexamethasone. The model demonstrated a good predictive performance in describing total and unbound imatinib concentrations in patients with SARS‐CoV‐2 infection. PBPK simulations highlight that an equivalent dose of imatinib may lead to substantially higher total drug concentrations in patients with SARS‐CoV‐2 infection compared to that in patients with cancer, while the unbound concentrations remain comparable between the 2 patient populations. This supports the notion that unbound trough concentration is a better exposure metric for dose adjustment of imatinib in patients with SARS‐CoV‐2 infection, compared to the corresponding total drug concentration. Potential strategies for refinement and generalization of the PBPK modeling approach in the patient population with SARS‐CoV‐2 are also provided in this article, which could be used to guide study design and inform dose adjustment in the future.

A physiologically based pharmacokinetic model of clopidogrel in populations of European and Japanese ancestry: An evaluation of CYP2C19 activity

Treatment response to clopidogrel is associated with CYP2C19 activity through the formation of the active H4 metabolite. The aims of this study were to develop a physiologically based pharmacokinetic (PBPK) model of clopidogrel and its metabolites for populations of European ancestry, to predict the pharmacokinetics in the Japanese population by CYP2C19 phenotype, and to investigate the effect of clinical and demographic factors. A PBPK model was developed and verified to describe the two metabolic pathways of clopidogrel (H4 metabolite, acyl glucuronide metabolite) for a population of European ancestry using plasma data from published studies. Subsequently, model predictions in the Japanese population were evaluated. The effects of CYP2C19 activity, fluvoxamine coadministration (CYP2C19 inhibitor), and population‐specific factors (age, sex, BMI, body weight, cancer, hepatic, and renal dysfunction) on the pharmacokinetics of clopidogrel and its metabolites were then characterized. The predicted/observed ratios for clopidogrel and metabolite exposure parameters were acceptable (twofold acceptance criteria). For all CYP2C19 phenotypes, steady‐state AUC_0‐τ of the H4 metabolite was lower for the Japanese (e.g., EM, 7.69 [6.26–9.45] ng·h/ml; geometric mean [95% CI]) than European (EM, 24.8 [20.4–30.1] ng·h/ml, p < .001) population. In addition to CYP2C19‐poor metabolizer phenotype, fluvoxamine coadministration, hepatic, and renal dysfunction were found to reduce H4 metabolite but not acyl glucuronide metabolite concentrations. This is the first PBPK model describing the two major metabolic pathways of clopidogrel, which can be applied to populations of European and Japanese ancestry by CYP2C19 phenotype. The differences between the two populations appear to be determined primarily by the effect of varying CYP2C19 liver activity.

PBPK Simulation-Based Evaluation of Ganciclovir Crystalluria Risk Factors: Effect of Renal Impairment, Old Age, and Low Fluid Intake

Dosing guidance is often lacking for chronic kidney disease (CKD) due to exclusion of such patients from pivotal clinical trials. Physiologically based pharmacokinetic (PBPK) modelling supports model-informed dosing when clinical data are lacking, but application of these approaches to patients with impaired renal function is not yet at full maturity. In the current study, a ganciclovir PBPK model was developed for patients with normal renal function and extended to CKD population. CKD-related changes in tubular secretion were explored in the mechanistic kidney model and implemented either as proportional or non-proportional decline relative to GFR. Crystalluria risk was evaluated in different clinical settings (old age, severe CKD and low fluid intake) by simulating ganciclovir medullary collecting duct (MCD) concentrations. The ganciclovir PBPK model captured observed changes in systemic pharmacokinetic endpoints in mild-to-severe CKD; these trends were evident irrespective of assumed pathophysiological mechanism of altered active tubular secretion in the model. Minimal difference in simulated ganciclovir MCD concentrations was noted between young adult and geriatric populations with normal renal function and urine flow (1 mL/min), with lower concentrations predicted for severe CKD patients. High crystalluria risk was identified at reduced urine flow (0.1 mL/min) as simulated ganciclovir MCD concentrations exceeded its solubility (2.6–6 mg/mL), irrespective of underlying renal function. The analysis highlighted the importance of appropriate distribution of virtual subjects’ systems data in CKD populations. The ganciclovir PBPK model illustrates the ability of this translational tool to explore individual and combined effects of age, urine flow, and renal impairment on local drug renal exposure.

Sex-related pharmacokinetic differences with aging

PURPOSE

The proportion of women increases with advanced age, but older women are often underrepresented in clinical trials. Therefore, little is known about the combined effect of sex- and age-related physiological changes on drug pharmacokinetics.

METHODS

We compiled clinical studies, which investigated sex-related pharmacokinetic differences both in older and young women and men. The ratio women/men was calculated for various pharmacokinetic parameters across adulthood to assess sex-related differences in drug pharmacokinetics with aging. The contribution of body weight and drug characteristics to sex-related pharmacokinetic differences were explored using analysis of variance.

RESULTS

We found 67 studies reporting the pharmacokinetics for 56 drugs both in older and young women and men. Median peak concentration (C_max) (interquartile range (IQR)) and drug exposure (AUC) (IQR) were 22% (8-41%) and 20% (0-39%) higher in women compared with men whereas time to peak concentration (t_max), apparent volume of distribution (VdF) and elimination half-life (t_1/2) were not significantly different. Body weight and the drug main elimination pathway contributed to sex-related differences in C_max and AUC. Relative to men, women had a modest increase in C_max with increasing age (r = 0.19, p = 0.04). Conversely, sex-related differences in AUC remained constant with increasing age.

CONCLUSION

The pharmacokinetic differences between women and men were modest and, with the exception of C_max, remained constant with increasing age. The higher plasma concentration might be correlated to more adverse events in older women and thus, drug treatment should be started on the lower recommended dosage when appropriate particularly for drugs characterized by a narrow therapeutic index.

Physiologically based pharmacokinetic (PBPK) modelling of tamsulosin related to CYP2D6*10 allele

Tamsulosin, a selective [Formula: see text]-adrenoceptor blocker, is commonly used for alleviation of lower urinary tract symptoms related to benign prostatic hyperplasia. Tamsulosin is predominantly metabolized by CYP3A4 and CYP2D6 enzymes, and several studies reported the effects of CYP2D6 genetic polymorphism on the pharmacokinetics of tamsulosin. This study aims to develop and validate the physiologically based pharmacokinetic (PBPK) model of tamsulosin in CYP2D6*wt/*wt, CYP2D6*wt/*10, and CYP2D6*10/*10 genotypes, using Simcyp® simulator. Physicochemical, and formulation properties and data for absorption, distribution, metabolism and excretion were collected from previous publications, predicted in the simulator, or optimized in different CYP2D6 genotypes. The tamsulosin PBPK model in CYP2D6*wt/*wt and CYP2D6*wt/*10 genotypes were developed based on the clinical pharmacokinetic study where a single oral dose of 0.2 mg tamsulosin was administered to 25 healthy Korean male volunteers with CYP2D6*wt/*wt and CYP2D6*wt/*10 genotypes. A previous pharmacokinetic study was used to develop the model in CYP2D6*10/*10 genotype. The developed model was validated using other clinical pharmacokinetic studies not used in development. The predicted exposures via the PBPK model in CYP2D6*wt/*10 and CYP2D6*10/*10 genotype was 1.23- and 1.76-fold higher than CYP2D6*wt/*wt genotype, respectively. The simulation profiles were visually similar to the observed profiles, and fold errors of all development and validation datasets were included within the criteria. Therefore, the tamsulosin PBPK model in different CYP2D6 genotypes with regards to CYP2D6*10 alleles was appropriately established. Our model can contribute to the implementation of personalized pharmacotherapy of patients, appropriately predicting the pharmacokinetics of tamsulosin reflecting their demographic and CYP2D6 genotype characteristics without unnecessary drug exposure.

A Physiological Approach to Pharmacokinetics in Chronic Kidney Disease

The conventional approach to approximating the pharmacokinetics of drugs in patients with chronic kidney disease (CKD) only accounts for changes in the estimated glomerular filtration rate. However, CKD is a systemic and multifaceted disease that alters many body systems. Therefore, the objective of this exercise was to develop and evaluate a whole-body mechanistic approach to predicting pharmacokinetics in patients with CKD. Physiologically based pharmacokinetic models were developed in PK-Sim v8.0 (www.open-systems-pharmacology.org) to mechanistically represent the disposition of 7 compounds in healthy human adults. The 7 compounds selected were eliminated by glomerular filtration and active tubular secretion by the organic cation transport system to varying degrees. After a literature search, the healthy adult models were adapted to patients with CKD by numerically accounting for changes in glomerular filtration rate, kidney volume, renal perfusion, hematocrit, plasma protein concentrations, and gastrointestinal transit. Literature-informed interindividual variability was applied to the physiological parameters to facilitate a population approach. Model performance in CKD was evaluated against pharmacokinetic data from 8 clinical trials in the literature. Overall, integration of the CKD parameterization enabled exposure predictions that were within 1.5-fold error across all compounds and patients with varying stages of renal impairment. Notable improvement was observed over the conventional approach to scaling exposure, which failed in all but 1 scenario in patients with advanced CKD. Further research is required to qualify its use for first-in-CKD dose selection and clinical trial planning for a wider selection of renally eliminated compounds, including those subject to anion transport.

Pooled Analysis of Gastric Emptying in Patients With Obesity: Implications for Oral Absorption Projection

PURPOSE

Gastric emptying time is one of limiting factors that determines the pharmacokinetic properties of drugs administered by mouth. Despite the high prevalence of obesity worldwide, modifications in gastric emptying time have not been systematically addressed in this set of patients. The current analysis aims to quantitatively address obesity-related changes in gastric emptying time of solids, semisolids, and liquids compared with lean individuals, highlighting the relevant pharmacokinetic implications of oral drug absorption in patients with obesity.

METHODS

We searched the Cochrane Library, PubMed, Web of Science, and Embase for all relevant articles published until November 1, 2020. Differences in gastrointestinal variables in relation to gastric emptying between obese and lean individuals were quantified by weighted mean difference (WMD) and ratio of means (RoM). Robustness of the analyses was evaluated by subgroup analysis and publication bias test.

FINDINGS

A total of 17 studies with 906 participants were included. The gastric half-emptying time of solids (WMD, -10.4 minutes; P = 0.001; RoM, 0.90; P = 0.01) and liquids (WMD, -6.14 minutes; P < 0.001; RoM, 0.83, P = 0.03) was significantly shorter in individuals with obesity compared with lean individuals. These findings were confirmed by the subgroup analyses and publication bias tests.

IMPLICATIONS

Our pooled analysis systemically quantifies the differences in gastric half-emptying time between individuals with obesity and lean individuals, facilitating better understanding and prediction of drug absorption in individuals with obesity through physiologically based pharmacokinetic approaches. Obesity is associated with a faster transit of both solids and liquids through the stomach.

Physiologically Based Pharmacokinetic Modelling to Identify Pharmacokinetic Parameters Driving Drug Exposure Changes in the Elderly

BACKGROUND

Medication use is highly prevalent with advanced age, but clinical studies are rarely conducted in the elderly, leading to limited knowledge regarding age-related pharmacokinetic changes.

OBJECTIVE

The objective of this study was to investigate which pharmacokinetic parameters determine drug exposure changes in the elderly by conducting virtual clinical trials for ten drugs (midazolam, metoprolol, lisinopril, amlodipine, rivaroxaban, repaglinide, atorvastatin, rosuvastatin, clarithromycin and rifampicin) using our physiologically based pharmacokinetic (PBPK) framework.

METHODS

PBPK models for all ten drugs were developed in young adults (20-50 years) following the best practice approach, before predicting pharmacokinetics in the elderly (≥ 65 years) without any modification of drug parameters. A descriptive relationship between age and each investigated pharmacokinetic parameter (peak concentration [C_max], time to C_max [t_max], area under the curve [AUC], clearance, volume of distribution, elimination-half-life) was derived using the final PBPK models, and verified with independent clinically observed data from 52 drugs.

RESULTS

The age-related changes in drug exposure were successfully simulated for all ten drugs. Pharmacokinetic parameters were predicted within 1.25-fold (70%), 1.5-fold (86%) and 2-fold (100%) of clinical data. AUC increased progressively by 0.9% per year throughout adulthood from the age of 20 years, which was explained by decreased clearance, while C_max, t_max and volume of distribution were not affected by human aging. Additional clinical data of 52 drugs were contained within the estimated variability of the established age-dependent correlations for each pharmacokinetic parameter.

CONCLUSION

The progressive decrease in hepatic and renal blood flow, as well as glomerular filtration, rate led to a reduced clearance driving exposure changes in the healthy elderly, independent of the drug.

Comparing Predictions of a PBPK Model for Cyclosporine With Drug Levels From Therapeutic Drug Monitoring

This study compared simulations of a physiologically based pharmacokinetic (PBPK) model implemented for cyclosporine with drug levels from therapeutic drug monitoring to evaluate the predictive performance of a PBPK model in a clinical population. Based on a literature search model parameters were determined. After calibrating the model using the pharmacokinetic profiles of healthy volunteers, 356 cyclosporine trough levels of 32 renal transplant outpatients were predicted based on their biometric parameters. Model performance was assessed by calculating absolute and relative deviations of predicted and observed trough levels. The median absolute deviation was 6 ng/ml (interquartile range: 30 to 31 ng/ml, minimum = -379 ng/ml, maximum = 139 ng/ml). 86% of predicted cyclosporine trough levels deviated less than twofold from observed values. The high intra-individual variability of observed cyclosporine levels was not fully covered by the PBPK model. Perspectively, consideration of clinical and additional patient-related factors may improve the model's performance. In summary, the current study has shown that PBPK modeling may offer valuable contributions for pharmacokinetic research in clinical drug therapy.

Design and conduct considerations for studies in patients with impaired renal function

An impaired renal function, including acute and chronic kidney disease and end‐stage renal disease, can be the result of aging, certain disease conditions, the use of some medications, or as a result of smoking. In patients with renal impairment (RI), the pharmacokinetics (PKs) of drugs or drug metabolites may change and result in increased safety risks or decreased efficacy. In order to make specific dose recommendations in the label of drugs for patients with RI, a clinical trial may have to be conducted or, when not feasible, modeling and simulations approaches, such as population PK modeling or physiologically‐based PK modelling may be applied. This tutorial aims to provide an overview of the global regulatory landscape and a practical guidance for successfully designing and conducting clinical RI trials or, alternatively, on applying modeling and simulation tools to come to a dose recommendation for patients with RI in the most efficient manner.

Physiologically-based pharmacokinetic model for alectinib, ruxolitinib, and panobinostat in the presence of cancer, renal impairment, and hepatic impairment

Renal (RIP) and hepatic (HIP) impairments are prevalent conditions in cancer patients. They can cause changes in gastric emptying time, albumin levels, hematocrit, glomerular filtration rate, hepatic functional volume, blood flow rates, and metabolic activity that can modify drug pharmacokinetics. Performing clinical studies in such populations has ethical and practical issues. Using predictive physiologically-based pharmacokinetic (PBPK) models in the evaluation of the PK of alectinib, ruxolitinib, and panobinostat exposures in the presence of cancer, RIP, and HIP can help in using optimal doses with lower toxicity in these populations. Verified PBPK models were customized under scrutiny to account for the pathophysiological changes induced in these diseases. The PBPK model-predicted plasma exposures in patients with different health conditions within average 2-fold error. The PBPK model predicted an area under the curve ratio (AUCR) of 1, and 1.8, for ruxolitinib and panobinostat, respectively, in the presence of severe RIP. On the other hand, the severe HIP was associated with AUCR of 1.4, 2.9, and 1.8 for alectinib, ruxolitinib, and panobinostat, respectively, in agreement with the observed AUCR. Moreover, the PBPK model predicted that alectinib therapeutic cerebrospinal fluid levels are achieved in patients with non-small cell lung cancer, moderate HIP, and severe HIP at 1-, 1.5-, and 1.8-fold that of healthy subjects. The customized PBPK models showed promising ethical alternatives for simulating clinical studies in patients with cancer, RIP, and HIP. More work is needed to quantify other pathophysiological changes induced by simultaneous affliction by cancer and RIP or HIP.

Estimation of Ontogeny Functions for Renal Transporters Using a Combined Population Pharmacokinetic and Physiology-Based Pharmacokinetic Approach: Application to OAT1,3

To date, information on the ontogeny of renal transporters is limited. Here, we propose to estimate the in vivo functional ontogeny of transporters using a combined population pharmacokinetic (popPK) and physiology-based pharmacokinetic (PBPK) modeling approach called popPBPK. Clavulanic acid and amoxicillin were used as probes for glomerular filtration, combined glomerular filtration, and active secretion through OAT1,3, respectively. The predictive value of the estimated OAT1,3 ontogeny function was assessed by PBPK predictions of renal clearance (CL_R) of other OAT1,3 substrates: cefazolin and piperacillin. Individual CL_Rpost-hoc values, obtained from a published popPK model on the concomitant use of clavulanic acid and amoxicillin in critically ill children between 1 month and 15 years, were used as dependent variables in the popPBPK analysis. CL_R was re-parameterized according to PBPK principles, resulting in the estimation of OAT1,3-mediated intrinsic clearance (CL_{int,OAT1,3,invivo}) and its ontogeny. CL_{int,OAT1,3,invivo} ontogeny was described by a sigmoidal function, reaching half of adult level around 7 months of age, comparable to findings based on renal transporter-specific protein expression data. PBPK-based CL_R predictions including this ontogeny function were reasonably accurate for piperacillin in a similar age range (2.5 months–15 years) as well as for cefazolin in neonates as compared to published data (%RMSPE of 21.2 and 22.8%, respectively and %PE within ±50%). Using this novel approach, we estimated an in vivo functional ontogeny profile for CL_{int,OAT1,3,invivo} that yields accurate CL_R predictions for different OAT1,3 substrates across different ages. This approach deserves further study on functional ontogeny of other transporters.

Development and evaluation of physiologically based pharmacokinetic drug-disease models for predicting captopril pharmacokinetics in chronic diseases

The advancement in the processing speeds of computing machines has facilitated the development of complex physiologically based pharmacokinetic (PBPK) models. These PBPK models can incorporate disease-specific data and could be used to predict pharmacokinetics (PK) of administered drugs in different chronic conditions. The present study aimed to develop and evaluate PBPK drug-disease models for captopril after incorporating relevant pathophysiological changes occurring in adult chronic kidney disease (CKD) and chronic heart failure (CHF) populations. The population-based PBPK simulator Simcyp was used as a modeling and simulation platform. The visual predictive checks and mean observed/predicted ratios (ratio_(Obs/pred)) of the PK parameters were used for model evaluation. The developed disease models were successful in predicting captopril PK in all three stages of CKD (mild, moderate, and severe) and CHF, as the observed and predicted PK profiles and the ratio_(obs/pred) for the PK parameters were in close agreement. The developed captopril PBPK models can assist in tailoring captopril dosages in patients with different disease severity (CKD and CHF).

Teicoplanin physiologically based pharmacokinetic modeling offers a quantitative assessment of a theoretical influence of serum albumin and renal function on its disposition

PURPOSE

Variability in teicoplanin pharmacokinetics has been explained by multiple factors such as body weight, renal function, and serum albumin level. To improve mechanistic understanding of the causes of variability, a physiologically based pharmacokinetic (PBPK) model can be used as a systematic platform. In this study, a PBPK model of teicoplanin was developed to quantitatively assess the effects of physiological changes due to disease status using virtual populations.

METHODS

Predictive performance of the models was evaluated by comparing simulated and observed concentration-time profiles of teicoplanin. Subsequently, sensitivity analyses were conducted to identify potential factors contributing to individual differences in teicoplanin PK.

RESULTS

The developed PBPK model generated concentration-time profiles that were comparable to clinical observations in healthy adults, including Caucasians and Japanese, and after single-dose and multiple-dose administration. The predicted PK parameters (i.e., C_max, AUC, clearance) were within a two-fold range of the observed data in patients with renal impairments as well as healthy adults. Changes in total and unbound teicoplanin concentrations at 72 h, after various dosing regimens (tested 4-14 mg/kg q12h for three doses as a loading dose and then 4-14 mg/kg daily as a maintenance dose), were sensitive to renal function and serum albumin concentrations.

CONCLUSION

The PBPK model of teicoplanin provides mechanistic insight into the factors altering its disposition and allows assessments of the theoretical and quantitative impact of individual changes in physiological parameters on its PK even when an actual assessment with adequate sample sizes of patients is challenging.

Physiologically based pharmacokinetic model of renally cleared antibacterial drugs in Chinese renal impairment patients

To preliminarily develop physiologically based population models for Chinese renal impairment patients and to evaluate the prediction performance of new population models by renally cleared antibacterial drugs. First, demographic data and physiological parameters of Chinese renal impairment patients were collected, and then the coefficients of the relative demographic and physiological equation were recalibrated to construct the new population models. Second, drug‐independent parameters of ceftazidime, cefodizime, vancomycin, and cefuroxime were collected and verified by Chinese healthy volunteers, Caucasian healthy volunteers, and Caucasian renal impairment population models built in Simcyp. Finally, the newly developed population models were applied to predict the plasma concentration of four antibacterial drugs in Chinese renal impairment patients. The new physiologically based pharmacokinetic (PBPK) population models can predict the main pharmacokinetic parameters, including area under the plasma concentration–time curve extrapolated to infinity (AUC_inf), renal clearance (CL_r), and peak concentration (C_max), of ceftazidime, cefodizime, vancomycin, and cefuroxime following intravenous administrations with less than twofold error in mild, moderate, and severe Chinese renal impairment patients. The accuracy and precision of the predictions were improved compared with the Chinese healthy volunteers and Caucasian renal impairment population models. The PBPK population models were preliminarily developed and the first‐step validation results of four antibacterial drugs following intravenous administration showed acceptable accuracy and precision. The population models still need more systematic validation by using more drugs and scenarios in future studies to support their applications on dosage recommendation for Chinese renal impairment patients.

Physiologically-Based Pharmacokinetic Modeling in Renal and Hepatic Impairment Populations: A Pharmaceutical Industry Perspective

The predictive performance of physiologically‐based pharmacokinetics (PBPK) models for pharmacokinetics (PK) in renal impairment (RI) and hepatic impairment (HI) populations was evaluated using clinical data from 29 compounds with 106 organ impairment study arms were collected from 19 member companies of the International Consortium for Innovation and Quality in Pharmaceutical Development. Fifty RI and 56 HI study arms with varying degrees of organ insufficiency along with control populations were evaluated. For RI, the area under the curve (AUC) ratios of RI to healthy control were predicted within twofold of the observed ratios for > 90% (N = 47/50 arms). For HI, > 70% (N = 43/56 arms) of the hepatically impaired to healthy control AUC ratios were predicted within twofold. Inaccuracies, typically overestimation of AUC ratios, occurred more in moderate and severe HI. PBPK predictions can help determine the need and timing of organ impairment study. It may be suitable for predicting the impact of RI on PK of drugs predominantly cleared by metabolism with varying contribution of renal clearance. PBPK modeling may be used to support mild impairment study waivers or clinical study design.

Application of Physiologically Based Pharmacokinetic Modeling to Predict the Effect of Renal Impairment on the Pharmacokinetics of Olanzapine and Samidorphan Given in Combination

BACKGROUND

A combination of the antipsychotic olanzapine and opioid receptor antagonist samidorphan (OLZ/SAM) is in development for the treatment of patients with schizophrenia or bipolar I disorder. The effect of severe renal impairment on the pharmacokinetics of olanzapine and samidorphan after a single oral dose of OLZ/SAM was evaluated in a clinical study. Complementary to the clinical findings, physiologically based pharmacokinetic modeling was used to assess the effects of varying degrees of renal impairment on the pharmacokinetics of olanzapine and samidorphan.

METHODS

A physiologically based pharmacokinetic model for OLZ/SAM was developed and validated by comparing model-simulated data with observed clinical data. The model was applied to predict changes in olanzapine and samidorphan pharmacokinetics after administration of OLZ/SAM in subjects with mild, moderate, and severe renal impairment relative to age-matched controls with normal renal function.

RESULTS

The model predicted 1.5- and 2.2-fold increases in olanzapine and samidorphan area under the plasma concentration-time curve (AUC), respectively, after a single dose of OLZ/SAM in subjects with severe renal impairment vs controls, which was consistent with results from the clinical study. Application of the model prediction indicated increases in steady-state olanzapine AUC of 1.2-, 1.5-, and 1.6-fold, and samidorphan AUC of 1.4-, 1.8-, and 2.2-fold, in subjects with mild, moderate, and severe renal impairment, respectively, relative to healthy controls.

CONCLUSIONS

Physiologically based pharmacokinetic modeling extended the findings from a clinical study in severe renal impairment to other untested clinical scenarios; these data could be of interest to clinicians treating patients with renal impairment.

Dose Optimization in Kidney Disease: Opportunities for PBPK Modeling and Simulation

Kidney disease affects pharmacokinetic (PK) profiles of not only renally cleared drugs but also nonrenally cleared drugs. The impact of kidney disease on drug disposition has not been fully elucidated, but describing the extent of such impact is essential for conducting dose optimization in kidney disease. Accurate evaluation of kidney function has been a clinical interest for dose optimization, and more scientists pay attention and conduct research for clarifying the role of drug transporters, metabolic enzymes, and their interplay in drug disposition as kidney disease progresses. Physiologically based pharmacokinetic (PBPK) modeling and simulation can provide valuable insights for dose optimization in kidney disease. It is a powerful tool to integrate discrete knowledge from preclinical and clinical research and mechanistically investigate system- and drug-dependent factors that may contribute to the changes in PK profiles. PBPK-based prediction of drug exposures may be used a priori to adjust dosing regimens and thereby minimize the likelihood of drug-related toxicity. With real-time clinical studies, parameter estimation may be performed with PBPK approaches that can facilitate identification of sources of interindividual variability. PBPK modeling may also facilitate biomarker research that aids dose optimization in kidney disease. U.S. Food and Drug Administration guidances related to conduction of PK studies in kidney impairment and PBPK documentation provide the foundation for facilitating model-based dose-finding research in kidney disease.

Effect of Kidney Function and Dialysis on the Pharmacokinetics and Pharmacodynamics of Roxadustat, an Oral Hypoxia-Inducible Factor Prolyl Hydroxylase Inhibitor

Background and Objectives

Roxadustat is an orally active hypoxia-inducible factor prolyl hydroxylase inhibitor for anemia in chronic kidney disease. The pharmacokinetics, metabolic profile, and pharmacodynamics of roxadustat were investigated in subjects with different degrees of kidney function.

Methods

This phase 1 open-label study enrolled subjects with normal and severely impaired kidney function, and end-stage renal disease (ESRD) on continuous ambulatory peritoneal dialysis (CAPD) or automated peritoneal dialysis (APD) or hemodialysis/hemodiafiltration (HD/HDF). All subjects received a single 100-mg dose of oral roxadustat. Within a single-sequence, two-treatment period design (P1/P2), subjects with ESRD on HD/HDF received roxadustat 2 h after (P1) and 2 h before (P2) a dialysis session. Area under the plasma concentration–time curve (AUC) from administration to infinity (AUC_inf), maximum concentration (C_max), and terminal elimination half-life (t_1/2) were assessed for roxadustat; AUC and C_max were assessed for erythropoietin.

Results

Thirty-four subjects were enrolled and received roxadustat (normal kidney function, n = 12; severely impaired kidney function, n = 9; ESRD on CAPD/APD, n = 1; ESRD on HD/HDF, n = 12). The geometric least-square mean ratio of AUC_inf was 223% and 195% in subjects with severely impaired kidney function and ESRD on HD/HDF, respectively, relative to subjects with normal kidney function; C_max and t_1/2 were comparable. The pharmacokinetic profile of roxadustat was not affected by HD/HDF. AUC_inf and t_1/2 for the metabolites of roxadustat increased in subjects with kidney impairment. The AUC and C_max of erythropoietin increased in subjects with severely impaired kidney function or ESRD on HD/HDF. Roxadustat was well tolerated.

Conclusions

Kidney function impairment increased the AUC of roxadustat and its metabolites. The C_max and t_1/2 of roxadustat were comparable among groups. Roxadustat and its metabolites were not cleared by HD/HDF.

Clinical Trials Registration Number: NCT02965040.

Electronic supplementary material

The online version of this article (10.1007/s13318-020-00658-w) contains supplementary material, which is available to authorized users.

Novel Mechanistic PBPK Model to Predict Renal Clearance in Varying Stages of CKD by Incorporating Tubular Adaptation and Dynamic Passive Reabsorption

Chronic kidney disease (CKD) has significant effects on renal clearance (CLr ) of drugs. Physiologically-based pharmacokinetic (PBPK) models have been used to predict CKD effects on transporter-mediated renal active secretion and CLr for hydrophilic nonpermeable compounds. However, no studies have shown systematic PBPK modeling of renal passive reabsorption or CLr for hydrophobic permeable drugs in CKD. The goal of this study was to expand our previously developed and verified mechanistic kidney model to develop a universal model to predict changes in CLr in CKD for permeable and nonpermeable drugs that accounts for the dramatic nonlinear effect of CKD on renal passive reabsorption of permeable drugs. The developed model incorporates physiologically-based tubular changes of reduced water reabsorption/increased tubular flow rate per remaining functional nephron in CKD. The final adaptive kidney model successfully (absolute fold error (AFE) all < 2) predicted renal passive reabsorption and CLr for 20 permeable and nonpermeable test compounds across the stages of CKD. In contrast, use of proportional glomerular filtration rate reduction approach without addressing tubular adaptation processes in CKD to predict CLr generated unacceptable CLr predictions (AFE = 2.61-7.35) for permeable compounds in severe CKD. Finally, the adaptive kidney model accurately predicted CLr of para-amino-hippuric acid and memantine, two secreted compounds, in CKD, suggesting successful integration of active secretion into the model, along with passive reabsorption. In conclusion, the developed adaptive kidney model enables mechanistic predictions of in vivo CLr through CKD progression without any empirical scaling factors and can be used for CLr predictions prior to assessment of drug disposition in renal impairment.

Physiologically Based Pharmacokinetic Modeling of Transdermal Selegiline and Its Metabolites for the Evaluation of Disposition Differences between Healthy and Special Populations

A physiologically based pharmacokinetic (PBPK) model of selegiline (SEL), and its metabolites, was developed in silico to evaluate the disposition differences between healthy and special populations. SEL is metabolized to methamphetamine (MAP) and desmethyl selegiline (DMS) by several CYP enzymes. CYP2D6 metabolizes the conversion of MAP to amphetamine (AMP), while CYP2B6 and CYP3A4 predominantly mediate the conversion of DMS to AMP. The overall prediction error in simulated PK, using the developed PBPK model, was within 0.5-1.5-fold after intravenous and transdermal dosing in healthy and elderly populations. Simulation results generated in the special populations demonstrated that a decrease in cardiac output is a potential covariate that affects the SEL exposure in renally impaired (RI) and hepatic impaired (HI) subjects. A decrease in CYP2D6 levels increased the systemic exposure of MAP. DMS exposure increased due to a reduction in the abundance of CYP2B6 and CYP3A4 in RI and HI subjects. In addition, an increase in the exposure of the primary metabolites decreased the exposure of AMP. No significant difference between the adult and adolescent populations, in terms of PK, were observed. The current PBPK model predictions indicate that subjects with HI or RI may require closer clinical monitoring to identify any untoward effects associated with the administration of transdermal SEL patch.

Impact of Disease on Plasma and Lung Exposure of Chloroquine, Hydroxychloroquine and Azithromycin: Application of PBPK Modeling

We use a mechanistic lung model to demonstrate that accumulation of chloroquine (CQ), hydroxychloroquine (HCQ), and azithromycin (AZ) in the lungs is sensitive to changes in lung pH, a parameter that can be affected in patients with coronavirus disease 2019 (COVID-19). A reduction in pH from 6.7 to 6 in the lungs, as observed in respiratory disease, led to 20-fold, 4.0-fold, and 2.7-fold increases in lung exposure of CQ, HCQ, and AZ, respectively. Simulations indicated that the relatively high concentrations of CQ and HCQ in lung tissue were sustained long after administration of the drugs had stopped. Patients with COVID-19 often present with kidney failure. Our simulations indicate that renal impairment (plus lung pH reduction) caused 30-fold, 8.0-fold, and 3.4-fold increases in lung exposures for CQ, HCQ, and AZ, respectively, with relatively small accompanying increases (20 to 30%) in systemic exposure. Although a number of different dosage regimens were assessed, the purpose of our study was not to provide recommendations for a dosing strategy, but to demonstrate the utility of a physiologically-based pharmacokinetic modeling approach to estimate lung concentrations. This, used in conjunction with robust in vitro and clinical data, can help in the assessment of COVID-19 therapeutics going forward.

The Influence of Drug Properties and Ontogeny of Transporters on Pediatric Renal Clearance through Glomerular Filtration and Active Secretion: a Simulation-Based Study

Glomerular filtration (GF) and active tubular secretion (ATS) contribute to renal drug elimination, with the latter remaining understudied across the pediatric age range. Therefore, we systematically analyzed the influence of transporter ontogeny on the relative contribution of GF and ATS to renal clearance CL_R for drugs with different properties in children. A physiology-based model for CL_R in adults was extrapolated to the pediatric population by including maturation functions for the system-specific parameters. This model was used to predict GF and ATS for hypothetical drugs with a range of drug-specific properties, including transporter-mediated intrinsic clearance (CL_int,T) values, that are substrates for renal secretion transporters with different ontogeny patterns. To assess the impact of transporter ontogeny on ATS and total CL_R, a percentage prediction difference (%PD) was calculated between the predicted CL_R in the presence and absence of transporter ontogeny. The contribution of ATS to CL_R ranges between 41 and 90% in children depending on fraction unbound and CL_int,T values. If ontogeny of renal transporters is < 0.2 of adult values, CL_R predictions are unacceptable (%PD > 50%) for the majority of drugs regardless of the pediatric age. Ignoring ontogeny patterns of secretion transporters increasing with age in children younger than 2 years results in CL_R predictions that are not systematically acceptable for all hypothetical drugs (%PD > 50% for some drugs). This analysis identified for what drug-specific properties and at what ages the contribution of ATS on total pediatric CL_R cannot be ignored. Drugs with these properties may be sensitive in vivo probes to investigate transporter ontogeny.

Effect of ageing on antiretroviral drug pharmacokinetics using clinical data combined with modelling and simulation

AIMS

The impact of ageing on antiretroviral pharmacokinetics remains uncertain, leading to missing dosing recommendations for elderly people living with human immunodeficiency virus (HIV: PLWH). The objective of this study was to investigate whether ageing leads to clinically relevant pharmacokinetic changes of antiretrovirals that would support a dose adjustment based on the age of the treated PLWH.

METHODS

Plasma concentrations for 10 first-line antiretrovirals were obtained in PLWH ≥55 years, participating in the Swiss HIV Cohort Study, and used to proof the predictive performance of our physiologically based pharmacokinetic (PBPK) model. The verified PBPK model predicted the continuous effect of ageing on HIV drug pharmacokinetics across adulthood (20-99 years). The impact of ethnicity on age-related pharmacokinetic changes between whites and other races was statistically analysed.

RESULTS

Clinically observed concentration-time profiles of all investigated antiretrovirals were generally within the 95% confidence interval of the PBPK simulations, demonstrating the predictive power of the modelling approach used. The predicted decline in drug clearance drove age-related pharmacokinetic changes of antiretrovirals, resulting in a maximal 70% [95% confidence interval: 40%, 120%] increase in antiretrovirals exposure across adulthood. Peak concentration, time to peak concentration and apparent volume of distribution were predicted to be unaltered by ageing. There was no statistically significant difference of age-related pharmacokinetic changes between studied ethnicities.

CONCLUSION

Dose adjustment for antiretrovirals based on the age of male and female PLWH is a priori not necessary in the absence of severe comorbidities considering the large safety margin of the current first-line HIV treatments.

Drugs Commonly Applied to Kidney Patients May Compromise Renal Tubular Uremic Toxins Excretion

In chronic kidney disease (CKD), the secretion of uremic toxins is compromised leading to their accumulation in blood, which contributes to uremic complications, in particular cardiovascular disease. Organic anion transporters (OATs) are involved in the tubular secretion of protein-bound uremic toxins (PBUTs). However, OATs also handle a wide range of drugs, including those used for treatment of cardiovascular complications and their interaction with PBUTs is unknown. The aim of this study was to investigate the interaction between commonly prescribed drugs in CKD and endogenous PBUTs with respect to OAT1-mediated uptake. We exposed a unique conditionally immortalized proximal tubule cell line (ciPTEC) equipped with OAT1 to a panel of selected drugs, including angiotensin-converting enzyme inhibitors (ACEIs: captopril, enalaprilate, lisinopril), angiotensin receptor blockers (ARBs: losartan and valsartan), furosemide and statins (pravastatin and simvastatin), and evaluated the drug-interactions using an OAT1-mediated fluorescein assay. We show that selected ARBs and furosemide significantly reduced fluorescein uptake, with the highest potency for ARBs. This was exaggerated in presence of some PBUTs. Selected ACEIs and statins had either no or a slight effect at supratherapeutic concentrations on OAT1-mediated fluorescein uptake. In conclusion, we demonstrate that PBUTs may compete with co-administrated drugs commonly used in CKD management for renal OAT1 mediated secretion, thus potentially compromising the residual renal function.

Modelling Age-Related Changes in the Pharmacokinetics of Risperidone and 9-Hydroxyrisperidone in Different CYP2D6 Phenotypes Using a Physiologically Based Pharmacokinetic Approach

PURPOSE

Dose-optimization strategies for risperidone are gaining in importance, especially in the elderly. Based on the genetic polymorphism of cytochrome P 450 (CYP) 2D6 genetically and age-related changes cause differences in the pharmacokinetics of risperidone and 9-hydroxyrisperidone. The goal of the study was to develop physiologically based pharmacokinetic (PBPK) models for the elderly aged 65+ years. Additionally, CYP2D6 phenotyping using metabolic ratio were applied and different pharmacokinetic parameter for different age classes predicted.

METHODS

Plasma concentrations of risperidone and 9-hydroxyrisperidone were used to phenotype 17 geriatric inpatients treated under naturalistic conditions. For this purpose, PBPK models were developed to examine age-related changes in the pharmacokinetics between CYP2D6 extensive metabolizer, intermediate metabolizer, poor metabolizer, (PM) and ultra-rapid metabolizer.

RESULTS

PBPK-based metabolic ratio was able to predict different CYP2D6 phenotypes during steady-state. One inpatient was identified as a potential PM, showing a metabolic ratio of 3.39. About 88.2% of all predicted plasma concentrations of the inpatients were within the 2-fold error range. Overall, age-related changes of the pharmacokinetics in the elderly were mainly observed in Cmax and AUC. Comparing a population of young adults with the oldest-old, Cmax of risperidone increased with 24-44% and for 9-hydroxyrisperidone with 35-37%.

CONCLUSIONS

Metabolic ratio combined with PBPK modelling can provide a powerful tool to identify potential CYP2D6 PM during therapeutic drug monitoring. Based on genetic, anatomical and physiological changes during aging, PBPK models ultimately support decision-making regarding dose-optimization strategies to ensure the best therapy for each patient over the age of 65 years.