Physiologically Based Pharmacokinetic Modeling to Predict the Impact of Liver Cirrhosis on Glucuronidation via UGT1A4 and UGT2B7/2B4-A Case Study with Midazolam.

Hepatic impairment, due to liver cirrhosis, decreases the activity of cytochrome P450 enzymes (CYPs). The use of physiologically based pharmacokinetic (PBPK) modeling to predict this effect for CYP substrates has been well-established, but the effect of cirrhosis on uridine-glucuronosyltransferase (UGT) activities is less studied and few PBPK models have been reported. UGT enzymes are involved in primary N-glucuronidation of midazolam and glucuronidation of 1'-OH-midazolam following CYP3A hydroxylation. In this study, Simcyp was used to establish PBPK models for midazolam, its primary metabolites midazolam-N-glucuronide (UGT1A4) and 1'-OH midazolam (CYP3A4/3A5), and the secondary metabolite 1'-OH-midazolam-O-glucuronide (UGT2B7/2B4), allowing to simulate the impact of liver cirrhosis on the primary and secondary glucuronidation of midazolam. The model was verified in noncirrhotic subjects before extrapolation to cirrhotic patients of Child-Pugh (CP) classes A, B, and C. Our model successfully predicted the exposures of midazolam and its metabolites in noncirrhotic and cirrhotic patients, with 86% of observed plasma concentrations within 5th-95th percentiles of predictions and observed geometrical mean of area under the plasma concentration curve between 0 hours to infinity and maximal plasma concentration within 0.7- to 1.43-fold of predictions. The simulated metabolic ratio defined as the ratio of the glucuronide metabolite AUC over the parent compound AUC (AUCglucuronide/AUCparent, metabolic ratio [MR]), was calculated for midazolam-N-glucuronide to midazolam (indicative of UGT1A4 activity) and decreased by 40% (CP A), 48% (CP B), and 75% (CP C). For 1'-OH-midazolam-O-glucuronide to 1'-OH-midazolam, the MR (indicative of UGT2B7/2B4 activity) dropped by 35% (CP A), 51% (CP B), and 64% (CP C). These predicted MRs were corroborated by the observed data. This work thus increases confidence in Simcyp predictions of the effect of liver cirrhosis on the pharmacokinetics of UGT1A4 and UGT2B7/UGT2B4 substrates. SIGNIFICANCE STATEMENT: This article presents a physiologically based pharmacokinetic model for midazolam and its metabolites and verifies the accurate simulation of pharmacokinetic profiles when using the Simcyp hepatic impairment population models. Exposure changes of midazolam-N-glucuronide and 1'-OH-midazolam-O-glucuronide reflect the impact of decreases in UGT1A4 and UGT2B7/2B4 glucuronidation activity in cirrhotic patients. The approach used in this study may be extended to verify the modeling of other uridine glucuronosyltransferase enzymes affected by liver cirrhosis.

Quantitative Cytochrome P450 3A4 Induction Risk Assessment Using Human Hepatocytes Complemented with Pregnane X Receptor-Activating Profiles.

Reliable in vitro to in vivo translation of cytochrome P450 (CYP) 3A4 induction potential is essential to support risk mitigation for compounds during pharmaceutical discovery and development. In this study, a linear correlation of CYP3A4 mRNA induction potential in human hepatocytes with the respective pregnane-X receptor (PXR) activation in a reporter gene assay using DPX2 cells was successfully demonstrated for 13 clinically used drugs. Based on this correlation, using rifampicin as a positive control, the magnitude of CYP3A4 mRNA induction for 71 internal compounds at several concentrations up to 10 µM (n = 90) was predicted within 2-fold error for 64% of cases with only a few false positives (19%). Furthermore, the in vivo area under the curve reduction of probe CYP substrates was reasonably predicted for eight marketed drugs (carbamazepine, dexamethasone, enzalutamide, nevirapine, phenobarbital, phenytoin, rifampicin, and rufinamide) using the static net effect model using both the PXR activation and CYP3A4 mRNA induction data. The liver exit concentrations were used for the model in place of the inlet concentrations to avoid false positive predictions and the concentration achieving twofold induction (F2) was used to compensate for the lack of full induction kinetics due to cytotoxicity and solubility limitations in vitro. These findings can complement the currently available induction risk mitigation strategy and potentially influence the drug interaction modeling work conducted at clinical stages. SIGNIFICANCE STATEMENT: The established correlation of CYP3A4 mRNA in human hepatocytes to PXR activation provides a clear cut-off to identify a compound showing an in vitro induction risk, complementing current regulatory guidance. Also, the demonstrated in vitro-in vivo translation of induction data strongly supports a clinical development program although limitations remain for drug candidates showing complex disposition pathways, such as involvement of auto-inhibition/induction, active transport and high protein binding.

A novel application of hydrophilic interaction liquid chromatography for the identification of compounds with intramolecular hydrogen bonds.

The incorporation of intramolecular hydrogen bonds (IMHB) into small molecules constitutes an interesting optimization strategy to afford potential drug candidates with enhanced solubility as well as permeability and consequently improved bioavailability (if metabolic stability is high). Common methods to assess IMHB rely on spectroscopic or diffraction techniques, which, however, have limited throughput when screening for hit compounds in early phases of drug discovery. Inspired by literature findings using supercritical fluid chromatography (SFC) as an indirect method for IMHB identification in a screening context, we aimed at developing a secondary chromatographic methodology taking advantage of commonly used HPLC-MS instrumentation. In this work, we explored hydrophilic interaction liquid chromatography (HILIC) and developed a method for discriminating compounds based on their hydrogen bonding features. By quantifying retention of different matched molecular pairs (MMP) and using information about their low energy conformations from quantum-mechanical calculations, we defined a hydrogen bonding-driven adsorption (kads) chromatographic parameter to assess a compound's propensity to forming IMHB. In addition to the MMP analysis, we found that the kads parameter allows for the differentiation of analytes forming IMHB regardless of the comparison with control compounds.