Physiologically based pharmacokinetic modeling to predict the impact of liver cirrhosis on glucuronidation via UGT1A4 & 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 non-cirrhotic 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 non-cirrhotic and cirrhotic patients, with 86% of observed plasma concentrations within 5 th -95 th percentiles of predictions and observed geometrical mean of AUC inf and C max within 0.7-1.43-fold of predictions. The simulated metabolic ratio (AUC glucuronide /AUC parent , 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 paper presents a PBPK 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 UGT enzymes affected by liver cirrhosis.