PBPK modeling of irbesartan: incorporation of hepatic uptake

Physiologically based pharmacokinetic (PBPK) modeling to predict the pharmacokinetics of irbesartan in different CYP2C9 genotypes

Irbesartan, a potent and selective angiotensin II type-1 (AT1) receptor blocker (ARB), is one of the representative medications for the treatment of hypertension. Cytochrome P450 (CYP) 2C9 is primarily involved in the oxidation of irbesartan. CYP2C9 is highly polymorphic, and genetic polymorphism of this enzyme is the leading cause of significant alterations in the pharmacokinetics of irbesartan. This study aimed to establish the physiologically based pharmacokinetic (PBPK) model to predict the pharmacokinetics of irbesartan in different CYP2C9 genotypes. The irbesartan PBPK model was established using the PK-Sim® software. Our previously reported pharmacogenomic data for irbesartan was leveraged in the development of the PBPK model and collected clinical pharmacokinetic data for irbesartan was used for the validation of the model. Physicochemical and ADME properties of irbesartan were obtained from previously reported data, predicted by the modeling software, or optimized to fit the observed plasma concentration-time profiles. Model evaluation was performed by comparing the predicted plasma concentration-time profiles and pharmacokinetic parameters to the observed results. Predicted plasma concentration-time profiles were visually similar to observed profiles. Predicted AUCinf in CYP2C9*1/*3 and CYP2C9*1/*13 genotypes were increased by 1.54- and 1.62-fold compared to CYP2C9*1/*1 genotype, respectively. All fold error values for AUC and Cmax in non-genotyped and CYP2C9 genotyped models were within the two-fold error criterion. We properly established the PBPK model of irbesartan in different CYP2C9 genotypes. It can be used to predict the pharmacokinetics of irbesartan for personalized pharmacotherapy in individuals of various races, ages, and CYP2C9 genotypes.

A Genetically Encoded Magnetic Resonance Imaging Reporter Enables Sensitive Detection and Tracking of Spontaneous Metastases in Deep Tissues

Metastasis is the leading cause of cancer-related death. However, it remains a poorly understood aspect of cancer biology, and most preclinical cancer studies do not examine metastasis, focusing solely on the primary tumor. One major factor contributing to this paradox is a gap in available tools for accurate spatiotemporal measurements of metastatic spread in vivo. Here, our objective was to develop an imaging reporter system that offers sensitive three-dimensional (3D) detection of cancer cells at high resolutions in live mice. An organic anion-transporting polypeptide 1b3 (oatp1b3) was used as an MRI reporter gene, and its sensitivity was systematically optimized for in vivo tracking of viable cancer cells in a spontaneous metastasis model. Metastases with oatp1b3-MRI could be observed at the single lymph node level and tracked over time as cancer cells spread to multiple lymph nodes and different organ systems in individual animals. While initial single lesions were successfully imaged in parallel via bioluminescence, later metastases were largely obscured by light scatter from the initial node. Importantly, MRI could detect micrometastases in lung tissue comprised on the order of 1,000 cancer cells. In summary, oatp1b3-MRI enables longitudinal tracking of cancer cells with combined high resolution and high sensitivity that provides 3D spatial information and the surrounding anatomical context.

SIGNIFICANCE

An MRI reporter gene system optimized for tracking metastasis in deep tissues at high resolutions and able to detect spontaneous micrometastases in lungs of mice provides a useful tool for metastasis research.

Integrated multi-biomarker responses of juvenile rainbow trout (Oncorhynchus mykiss) to an environmentally relevant pharmaceutical mixture

Pharmaceuticals are emerging pollutants of concern for aquatic ecosystems where they are occurring in complex mixtures. In the present study, the chronic toxicity of an environmentally relevant pharmaceutical mixture on juvenile rainbow trout (Oncorhynchus mykiss) was investigated. Five pharmaceuticals (paracetamol, carbamazepine, diclofenac, naproxen and irbesartan) were selected based on their detection frequency and concentration levels in the Meuse river (Belgium). Fish were exposed for 42 days to three different concentrations of the mixture, the median one detected in the Meuse river, 10-times and 100-times this concentration. Effects on the nervous, immune, antioxidant, and detoxification systems were evaluated throughout the exposure period and their response standardized using the Integrated Biomarker Response (IBRv2) index. IBRv2 scores increased over time in the fish exposed to the highest concentration. After 42 days, fish exposed to the highest concentration displayed significantly higher levels in lysozyme activity (p < 0.01). The mixture also caused significant changes in brain serotonin turnover (p < 0.05). In short, our results indicate that the subchronic waterborne exposure to a pharmaceutical mixture commonly occurring in freshwater ecosystems may affect the neuroendocrine and immune systems of juvenile rainbow trout.

Incorporating renal excretion via the OCT2 transporter in physiologically based kinetic modelling to predict in vivo kinetics of mepiquat in rat

The present study aimed at incorporating active renal excretion via the organic cation transporter 2 (OCT2) into a generic rat physiologically based kinetic (PBK) model using an in vitro human renal proximal tubular epithelial cell line (SA7K) and mepiquat chloride (MQ) as the model compound. The Vmax (10.5 pmol/min/mg protein) and Km (20.6 μM) of OCT2 transport of MQ were determined by concentration-dependent uptake in SA7K cells using doxepin as inhibitor. PBK model predictions incorporating these values in the PBK model were 6.7-8.4-fold different from the reported in vivo data on the blood concentration of MQ in rat. Applying an overall scaling factor that also corrects for potential differences in OCT2 activity in the SA7K cells and in vivo kidney cortex and species differences resulted in adequate predictions for in vivo kinetics of MQ in rat (2.3-3.2-fold). The results indicate that using SA7K cells to define PBK parameters for active renal OCT2 mediated excretion with adequate scaling enables incorporation of renal excretion via the OCT2 transporter in PBK modelling to predict in vivo kinetics of mepiquat in rat. This study demonstrates a proof-of-principle on how to include active renal excretion into generic PBK models.

Novel testing strategy for prediction of rat biliary excretion of intravenously administered estradiol-17β glucuronide

The aim of the present study was to develop a generic rat physiologically based kinetic (PBK) model that includes a novel testing strategy where active biliary excretion is incorporated using estradiol-17β glucuronide (E217βG) as the model substance. A major challenge was the definition of the scaling factor for the in vitro to in vivo conversion of the PBK-model parameter Vmax. In vitro values for the Vmax and Km for transport of E217βG were found in the literature in four different studies based on experiments with primary rat hepatocytes. The required scaling factor was defined based on fitting the PBK model-based predicted values to reported experimental data on E217βG blood levels and cumulative biliary E217βG excretion. This resulted in a scaling factor of 129 mg protein/g liver. With this scaling factor the PBK model predicted the in vivo data for blood and cumulative biliary E217βG levels with on average of less than 1.8-fold deviation. The study provides a proof of principle on how biliary excretion can be included in a generic PBK model using primary hepatocytes to define the kinetic parameters that describe the biliary excretion.

Cardiovascular Risk Management and Hepatitis C: Combining Drugs

Direct-acting antivirals (DAAs) are known victims (substrate) and perpetrators (cause) of drug–drug interactions (DDIs). These DAAs are used for the treatment of hepatitis C virus (HCV) infections and are highly effective drugs. Drugs used for cardiovascular risk management are frequently used by HCV-infected patients, whom also are treated with DAAs. Therefore, the aim of this review was to describe DDIs between cardiovascular drugs (CVDs) and DAAs. An extensive literature search was performed containing search terms for the marketed DAAs and CVDs (β-blocking agents, ACE inhibitors, angiotensin II antagonists, renin inhibitors, diuretics, calcium channel blockers, statins/ezetimibe, fibrates, platelet aggregation inhibitors, vitamin K antagonists, heparins, direct Xa inhibitors, nitrates, amiodarone, and digoxin). In particular, the drug labels from the European Medicines Agency and the US Food and Drug Administration were used. A main finding of this review is that CVDs are mostly victims of DDIs with DAAs. Therefore, when possible, monitoring of pharmacodynamics is recommended when coadministering these drugs with DAAs. Nevertheless, it is sometimes better to discontinue a drug on a temporary basis (statins, ezetimide). The DAAs are victims of DDIs in combination with bisoprolol, carvedilol, labetalol, verapamil, and gemfibrozil. Despite there are many DDIs predicted in this review, most of these DDIs can be managed by monitoring the efficacy and toxicity of the victim drug or by switching to another CVD/DAA.

Understanding the Oral Absorption of Irbesartan Using Biorelevant Dissolution Testing and PBPK Modeling

Poorly soluble weak bases form a significant proportion of the drugs available in the market thereby making it imperative to understand their absorption behavior. This work aims to mechanistically understand the oral absorption behavior for a weakly basic drug, Irbesartan (IRB), by investigating its pH dependent solubility, supersaturation, and precipitation behavior. Simulations performed using the equilibrium solubility could not accurately predict oral absorption. A multi-compartmental biorelevant dissolution testing model was used to evaluate dissolution in the stomach and duodenal compartment and mimic oral drug administration. This model exhibited sustained intestinal supersaturation (2-4-fold) even upon varying flow rates (4 mL/min, 7 mL/min, and mono-exponential transfer) from gastric to intestinal compartment. Simulation of oral absorption using GastroPlus™ and dissolution data collectively predicted plasma exposure with higher accuracy (% prediction error values within ± 15%), thereby indicating that multi-compartment dissolution testing enabled an improved prediction for oral pharmacokinetics of Irbesartan. Additionally, precipitates obtained in the intestinal compartment were characterized to determine the factors underlying intestinal supersaturation of Irbesartan. The solid form of these precipitates was amorphous with considerable particle size reduction. This indicated that following gastric transit, precipitate formation in the amorphous form coupled with an approximately 10 times particle size reduction could be potential factors leading to the generation and sustenance of intestinal drug supersaturation.

Positron Emission Tomography Imaging of [11 C]Rosuvastatin Hepatic Concentrations and Hepatobiliary Transport in Humans in the Absence and Presence of Cyclosporin A

Using positron emission tomography imaging, we determined the hepatic concentrations and hepatobiliary transport of [¹¹ C]rosuvastatin (RSV; i.v. injection) in the absence (n = 6) and presence (n = 4 of 6) of cyclosporin A (CsA; i.v. infusion) following a therapeutic dose of unlabeled RSV (5 mg, p.o.) in healthy human volunteers. The sinusoidal uptake, sinusoidal efflux, and biliary efflux clearance (CL; mL/minute) of [¹¹ C]RSV, estimated through compartment modeling were 1,205.6 ± 384.8, 16.2 ± 11.2, and 5.1 ± 1.8, respectively (n = 6). CsA (blood concentration: 2.77 ± 0.24 μM), an organic-anion-transporting polypeptide, Na⁺ -taurocholate cotransporting polypeptide, and breast cancer resistance protein inhibitor increased [¹¹ C]RSV systemic blood exposure (45%; P < 0.05), reduced its biliary efflux CL (52%; P < 0.05) and hepatic uptake (25%; P > 0.05) but did not affect its distribution into the kidneys. CsA increased plasma concentrations of coproporphyrin I and III and total bilirubin by 297 ± 69%, 384 ± 102%, and 81 ± 39%, respectively (P < 0.05). These data can be used in the future to verify predictions of hepatic concentrations and hepatobiliary transport of RSV.

Physiologically-based pharmacokinetic modelling of infant exposure to efavirenz through breastfeeding

Background: Very little is known about the level of infant exposure to many drugs commonly used during breastfeeding. The aim of this study was to develop a physiologically-based pharmacokinetic (PBPK) model for predicting infant exposure to maternal efavirenz through breastmilk. Methods: A breastfeeding PBPK model combining whole-body maternal and infant sub-models was constructed from drug-specific and system parameters affecting drug disposition using mathematical descriptions. The model was validated against published data on the pharmacokinetics of efavirenz in nursing mother-infant pairs. Further simulations were conducted to assess exposure in the context of the 400 mg reduced dose of efavirenz as well as best- and worse-case scenarios. Results: The model adequately described efavirenz pharmacokinetics, with over 80% of observed data points (203 matched breast milk and plasma pairs) within the predictive interval. All parameters were within 2-fold difference of clinical data. Median (range) predicted versus observed breast milk AUC0-24, Cmax and Cmin at the standard 600 mg dose were 75.0 (18.5-324) versus 68.5 (26.3-257) µg.hr/mL, 4.56 (1.17-16.0) versus 5.39 (1.43-18.4) µg/mL, and 2.11 (0.38-12.3) versus 1.68 (0.316-9.57) µg/mL, respectively. Predicted plasma AUC0-24, Cmax and Cmin at 400 mg reduced dose were similar to clinical data from non-breastfeeding adults. Model-predicted infant plasma concentrations were similar to clinical data, 0.15 (0.026–0.78) μg/mL at the 400 mg maternal dose in pooled analysis, approximately 25% lower than simulated exposure at 600 mg. The maximum exposure index was observed in the youngest infants, 5.9% (2.2-20) at 400 mg and 8.7% (3.2-29) at 600 mg. Thirteen and 36% of 10 days-1 month old infants were predicted to have exposure index above the 10% recommended threshold at 400 mg and 600 mg maternal dose, respectively. Conclusions: This application of PBPK modelling opens up opportunities for expanding our understanding of infant exposure to maternal drugs through breastfeeding.

Prediction and Quantification of Hepatic Transporter-Mediated Uptake of Pitavastatin Utilizing a Combination of the Relative Activity Factor Approach and Mechanistic Modeling

Quantification of the fraction transported (f_t) by a particular transporter will facilitate more robust estimations of transporter interactions. Using pitavastatin as a model uptake transporter substrate, we investigated the utility of the relative activity factor (RAF) approach and mechanistic modeling to estimate f_t in hepatocytes. The transporters evaluated were organic anion-transporting polypeptides OATP1B1 and OATP1B3 and sodium-taurocholate cotransporting polypeptide. Transporter-expressing human embryonic kidney 293 cells and human hepatocytes were used for determining RAF values, which were then incorporated into the mechanistic model to simulate hepatocyte uptake of pitavastatin over time. There was excellent agreement between simulated and observed hepatocyte uptake of pitavastatin, indicating the suitability of this approach for translation of uptake from individual transporter-expressing cells to more holistic in vitro models. Subsequently, f_t values were determined. The largest contributor to hepatocyte uptake of pitavastatin was OATP1B1, which correlates with what is known about the in vivo disposition of pitavastatin. The f_t values were then used for evaluating in vitro-in vivo correlations of hepatic uptake inhibition with OATP inhibitors rifampicin and cyclosporine. Predictions were compared with previously reported plasma exposure changes of pitavastatin with these inhibitors. Although hepatic uptake inhibition of pitavastatin was 2-3-fold underpredicted, incorporation of scaling factors (SFs) into RAF values significantly improved the predictive ability. We propose that calibration of hepatocytes with standard transporter substrates and inhibitors would allow for determination of system-specific SFs, which could subsequently be used for refining predictions of clinical DDI potential for new chemical entities that undergo active hepatic uptake.

Successful Prediction of In Vivo Hepatobiliary Clearances and Hepatic Concentrations of Rosuvastatin Using Sandwich-Cultured Rat Hepatocytes, Transporter-Expressing Cell Lines, and Quantitative Proteomics

We determined whether in vivo transporter-mediated hepatobiliary clearance (CL) and hepatic concentrations of rosuvastatin (RSV) in the rat could be predicted by transport activity in sandwich-cultured rat hepatocytes (SCRHs) and/or transporter-expressing cell lines scaled by differences in transporter protein expression between SCRHs, cell lines, and rat liver. The predicted hepatobiliary CLs and hepatic concentrations of RSV were compared with our previously published positron emission tomography imaging data. Sinusoidal uptake CL ([Formula: see text]) and efflux (canalicular and sinusoidal) CLs of [³H]-RSV in SCRHs were evaluated in the presence and absence of Ca²⁺ and in the absence and presence of 1 mM unlabeled RSV (to estimate passive diffusion CL). [Formula: see text] of RSV into cells expressing organic anion transporting polypeptide (Oatp) 1a1, 1a4, and 1b2 was also determined. Protein expression of Oatps in SCRHs and Oatp-expressing cells was quantified by liquid chromatography tandem mass spectrometry. SCRHs well predicted the in vivo RSV sinusoidal and canalicular efflux CLs but significantly underestimated in vivo [Formula: see text]. Oatp expression in SCRHs was significantly lower than that in the rat liver. [Formula: see text], based on RSV [Formula: see text] into Oatp-expressing cells (active transport) plus passive diffusion CL in SCRHs, scaled by the difference in protein expression in Oatp cells versus SCRH versus rat liver, was within 2-fold of that observed in SCRHs or in vivo. In vivo hepatic RSV concentrations were well predicted by Oatp-expressing cells after correcting [Formula: see text] for Oatp protein expression. This is the first demonstration of the successful prediction of in vivo hepatobiliary CLs and hepatic concentrations of RSV using transporter-expressing cells and SCRHs.

Impaired N-linked glycosylation of uptake and efflux transporters in human non-alcoholic fatty liver disease

BACKGROUND & AIMS

N-linked glycosylation of proteins is critical for proper protein folding and trafficking to the plasma membrane. Drug transporters are one class of proteins that have reduced function when glycosylation is impaired. N-linked glycosylation of plasma proteins has also been investigated as a biomarker for several liver diseases, including non-alcoholic fatty liver disease (NAFLD). The purpose of this study was to assess the transcriptomic expression of genes involved in protein processing and glycosylation, and to determine the glycosylation status of key drug transporters during human NAFLD progression.

METHODS

Human liver samples diagnosed as healthy, steatosis, and non-alcoholic steatohepatitis (NASH) were analysed for gene expression of glycosylation-related genes and for protein glycosylation using immunoblot.

RESULTS

Genes involved in protein processing in the ER and biosynthesis of N-glycans were significantly enriched for down-regulation in NAFLD progression. Included in the down regulated N-glycan biosynthesis category were genes involved in the oligosaccharyltransferase complex, N-glycan quality control, N-glycan precursor biosynthesis, N-glycan trimming to the core, and N-glycan extension from the core. N-glycan degradation genes were unaltered in the progression to NASH. Immunoblot analysis of the uptake transporters organic anion transporting polypeptide-1B1 (OATP1B1), OATP1B3, OATP2B1, and Sodium/Taurocholate Co-transporting Polypeptide (NTCP) and the efflux transporter multidrug resistance-associated protein 2 (MRP2) demonstrated a significant loss of glycosylation following the progression to NASH.

CONCLUSIONS

These data suggest that the loss of glycosylation of key uptake and efflux transporters in humans NASH may influence transporter function and contribute to altered drug disposition observed in NASH.

Application of the MechPeff model to predict passive effective intestinal permeability in the different regions of the rodent small intestine and colon

A major component of physiologically based pharmacokinetic (PBPK) models is the prediction of the rate and extent of absorption of orally dosed drugs for which knowledge of effective passive intestinal permeability (P_eff ) is essential. Single-pass intestinal perfusion (SPIP) studies are used to establish effective permeability in vivo but are difficult to perform in rodents, while mechanistic models to predict drug P_eff in rat and mouse have not been published. This work evaluates the predictive performance of the 'MechPeff' model to predict P_eff in the rodent intestine based upon knowledge of regional gut physiology and drug-specific physicochemical parameters. The 'MechPeff' model, built-in to the Simcyp Rat and Mouse Simulators, predicts transcellular, paracellular and mucus layer permeabilities and combines these to give the overall P_eff . The jejunal and/or ileal P_eff was predicted for 12 (4) acidic, 13 (12) basic, 10 (8) neutral and 2 (0) ampholytic drugs in the rat (mouse), spanning a wide range of MW and logP_o:w , and compared with experimental P_eff obtained using SPIP. A key input is the intrinsic transcellular permeability (P_trans,0 ) which can be derived from modelling of appropriate in vitro permeability experiments or predicted from physicochemical properties. The P_eff predictions were reasonably good when experimentally derived P_trans,0 was used; from 42 P_eff,rat values, 24 (57%) were within 3-fold, and of 19 P_eff,mouse values, 12 (63%) were within 3-fold, of observed P_eff . Considering the lack of alternative models to predict P_eff in preclinical species, and the minimal drug-specific inputs required, this model provides a valuable tool within drug discovery and development programmes. Copyright © 2017 John Wiley & Sons, Ltd.