Visual Overview
Abstract
Tolvaptan, a vasopressin V2-receptor antagonist, has demonstrated efficacy in slowing kidney function decline in patients with autosomal dominant polycystic kidney disease (ADPKD). In the pivotal clinical trial, the incidence of elevated liver enzymes was higher in patients receiving tolvaptan compared with placebo. Adjudication by a panel of expert hepatologists concluded a causal link of tolvaptan to liver injury in patients with ADPKD. An ex situ isolated perfused liver (IPL) study of tolvaptan disposition was undertaken in a rodent model of ADPKD, the polycystic kidney (PCK) rat (n = 5), and compared with wild-type (WT) Sprague-Dawley rats (n = 6). Livers were perfused with tolvaptan, followed by a tolvaptan-free washout phase. Total recovery (mean ± S.D. percentage of dose; PCK vs. WT) of tolvaptan and two metabolites, DM-4103 and DM-4107, quantified by liquid chromatography–tandem mass spectroscopy, was 58.14% ± 24.72% vs. 43.40% ± 18.11% in liver, 20.10% ± 9.15% vs. 21.17% ± 12.51% in outflow perfusate, and 0.08% ± 0.01% vs. 0.39% ± 0.32% in bile. DM-4103 recovery (mean ± S.D. percentage of dose) was decreased in PCK vs. WT bile (<0.01% ± <0.01% vs. 0.02% ± 0.01%; P = 0.0037), and DM-4107 recovery was increased in PCK vs. WT outflow perfusate (1.60% ± 0.57% vs. 0.43% ± 0.29%; P = 0.0017). A pharmacokinetic compartmental model assuming first-order processes was developed to describe the rate vs. time profiles of tolvaptan and DM-4103 + DM-4107 in rat IPLs. The model-derived estimate of tolvaptan’s biliary clearance was significantly decreased in PCK compared with WT IPLs. The model predicted greater hepatocellular concentrations of tolvaptan and DM-4103 + DM-4107 in PCK compared with WT IPLs. Increased hepatocellular exposure to tolvaptan and metabolites may contribute to the hepatotoxicity in patients with ADPKD treated with tolvaptan.
Introduction
Polycystic kidney disease (PKD) is a class of ciliopathies arising from primary cilium defects in epithelial cells of various organs, resulting in the development of fluid-filled cysts (Igarashi and Somlo, 2007). Autosomal dominant polycystic kidney disease (ADPKD) is the most common type of PKD and hereditary form of kidney disease (Torres et al., 2007), and a significant health concern that leads to kidney failure in approximately one-half of affected individuals (Grantham, 2008). ADPKD is primarily characterized by the progressive development of renal cysts and declining kidney function, but extrarenal cysts are common in this multiorgan disorder. The liver is affected in more than 80% of patients with ADPKD, leading to polycystic liver disease (Grantham, 2008). Despite complications arising from hepatic cysts such as cyst rupture, infection, and obstruction of bile ducts, liver function is thought to be preserved in polycystic liver disease (Cnossen and Drenth, 2014).
ADPKD arises from mutations in the PKD1 or PKD2 gene, encoding polycystin 1 and 2, respectively, which together form a receptor-channel complex (Kaimori and Germino, 2008). The illness has similarities to the less prevalent autosomal recessive PKD, which is caused by mutations in PKHD1, encoding the membrane-associated receptor-like protein fibrocystin/polyductin (Kim et al., 2008). The polycystic kidney (PCK) rat, a rodent model of human PKD that spontaneously arose in a colony of Sprague-Dawley rats, harbors a mutation in the Pkhd1 gene, similar to autosomal recessive PKD (Katsuyama et al., 2000; Masyuk et al., 2004). Since the natural history and the hepatorenal abnormalities in PCK disease and ADPKD are similar, this rodent model has been used to study ADPKD (Lager et al., 2001; Sabbatini et al., 2014).
Tolvaptan (TVP), an oral, selective vasopressin V2-receptor antagonist, has demonstrated efficacy in slowing kidney function decline in patients with ADPKD. A concerning finding in the pivotal clinical trial Tolvaptan Efficacy and Safety in Management of Autosomal Dominant Polycystic Kidney Disease and Its Outcomes (TEMPO) 3:4 was the higher incidence of elevated alanine aminotransferase (ALT) in patients receiving TVP (4.9%) compared with placebo (1.2%), with ALT > 2.5× the upper limit of normal (ULN) at any trial visit (Torres et al., 2012). Adjudication by a panel of hepatologists concluded a causal link between TVP and liver injury in patients with ADPKD, and in TEMPO 3:4 and an open-label extension trial, TEMPO 4:4, three patients met Hy’s law criteria (Watkins et al., 2015). Similar percentages (5.6% vs. 1.2%, TVP vs. placebo, respectively) in ALT elevations (>3× the ULN) but no additional Hy’s law cases were observed in a second pivotal trial, Replicating Evidence of Preserved Renal Function: an Investigation of Tolvaptan Safety and Efficacy in ADPKD with monthly monitoring of liver enzymes (Torres et al., 2017).
TVP is primarily metabolized by cytochrome P450 3A4 (CYP3A4) (Shoaf et al., 2011); CYP3A5, CYP2D6, uridine 5’-diphospho-glucuronosyltransferase 2B7 (UGT2B7), and UGT2B17 also contribute to TVP metabolite formation (Mazzarino et al., 2017). Numerous TVP metabolites have been identified, but most studies focus on the two major metabolites, DM-4103 and DM-4107 (Jiang et al., 2016; Lu et al., 2016). In humans, the main circulating metabolite is DM-4103, while a substantial portion of TVP is excreted in urine and feces as DM-4107 (Sorbera et al., 2002). Interactions between TVP, DM-4103, and DM-4107 and various human hepatic transporters have been investigated. DM-4103 is a relatively strong inhibitor of Na+-taurocholate cotransporting polypeptide [SLC10A1; half-maximal inhibitory concentration (IC50 value) = 16.3 μM], bile salt export pump (ABCB11; IC50 = 4.15 μM), and multidrug resistance-associated protein 4 (ABCC4; IC50 = 4.26 μM), while TVP is a moderate inhibitor of bile salt export pump (IC50 = 31.6 μM), and DM-4107 is a moderate inhibitor of multidrug resistance-associated protein 4 (IC50 = 37.9 μM) (Slizgi et al., 2016). Furthermore, TVP is a substrate and competitive inhibitor of P-glycoprotein (MDR1/ABCB1) (Shoaf et al., 2011). Currently, it remains unclear whether TVP, DM-4103, DM-4107, or other metabolites are responsible for the elevated liver enzymes observed in the ADPKD clinical trials.
The isolated perfused liver (IPL) is a physiologically relevant ex situ model used to assess hepatobiliary drug clearance (Brouwer and Thurman, 1996), and evaluate the hepatic transport and metabolic properties of a compound (Booth et al., 1998; Xiong et al., 2000; Chandra et al., 2005; Miranda et al., 2008; Pfeifer et al., 2013). To our knowledge, hepatobiliary function in PCK rats has not been evaluated previously using the IPL model. In the present study, the disposition of TVP and the major metabolites, DM-4103 and DM-4107, was investigated in IPLs from PCK and wild-type (WT) Sprague-Dawley rats to elucidate disease-mediated alterations in hepatobiliary function.
Materials and Methods
Chemicals.
All chemicals were purchased from Sigma-Aldrich (St. Louis, MO), unless otherwise stated. Livers were perfused with freshly prepared, protein-free, modified Krebs-Henseleit buffer supplemented with 0.373 g/l calcium chloride dihydrate and 2.1 g/l sodium bicarbonate (Fisher Scientific, Fair Lawn, NJ); 7.735 mg/l taurocholic acid sodium salt hydrate was added (15 µM) to maintain bile flow. The pH of the buffer was adjusted to 7.4 with sodium hydroxide or hydrochloric acid. TVP, DM-4103, DM-4107, and OPC-41100, which was used as the internal standard (IS), were provided by Otsuka Pharmaceutical Co., Ltd. (Tokushima, Japan).
Animals.
Male, 20- to 21-week-old, WT and PCK (PCK/CrljCrl-Pkhd1pck/Crl) Sprague-Dawley rats (Charles River Laboratories, Wilmington, MA) were used for the IPL studies. Cyst development in the PCK rat is substantial at this age, evidenced by an increase of approximately 85% in hepatic cyst volume from 12 to 20 weeks, and by significantly increased liver and kidney weights in 20-week-old PCK rats compared with age-matched WT rats (Brock et al., 2018). Rats were housed in an alternating 12-hr light/dark cycle, allowed water and food ad libitum, and acclimated for a minimum of 1 week prior to liver perfusions. All animal procedures were carried out in accordance with the Guide for the Care and Use of Laboratory Animals as adopted and promulgated by the US National Institutes of Health, and were approved by the Institutional Animal Care and Use Committee (University of North Carolina, Chapel Hill, NC). All surgeries were performed under full anesthesia with ketamine:xylazine (60:12 mg/kg, i.p. injection).
Tolvaptan Disposition in Rat Isolated Perfused Livers.
Following cannulation of the bile duct and portal vein, the livers were perfused ex situ in a single-pass manner with continuously oxygenated, albumin-free Krebs-Henseleit buffer (35 ml/min) as described previously (Brouwer and Thurman, 1996). Albumin was excluded so that more TVP, which is highly protein bound, could be taken up by the hepatocytes, resulting in more quantifiable concentrations of TVP and metabolites in outflow perfusate and bile. After a 15-min equilibration period with blank perfusate, livers were perfused for 30 min with TVP-containing perfusate (∼7 μM; the exact concentration was measured for each experiment), followed by TVP-free perfusion for 30 min. Bile was collected continuously over 5-min intervals for 60 min. Outflow perfusate was collected directly from the perfusate flowing out of the vena cava every 5 min for 60 min. Blank bile was collected during the surgery immediately after bile duct cannulation. After each perfusion, the liver was blotted and weighed. All samples were initially snap frozen, and subsequently stored at −80°C. Liver viability was assessed by monitoring portal perfusion pressure (<15 cm H2O), gross morphology, bile flow (>2 μl/min) (Zhang et al., 2016), and lactate dehydrogenase (LDH) levels in the outflow perfusate measured by a Pierce LDH Cytotoxicity Assay Kit (Thermo Fisher Scientific, Waltham, MA) per the manufacturer’s instructions. As a positive control for the LDH assay, one WT rat liver was perfused with 1% Triton X-100, which resulted in immediate liver toxicity based on gross morphology and major LDH release in outflow perfusate. Using dilutions and a standard curve for the outflow perfusate sample with the most LDH release (defined as 100% cytotoxicity), the percentage of cytotoxicity in the other perfused livers was evaluated; LDH release from these livers was negligible (<1% cytotoxicity, data not shown).
Bioanalysis.
TVP, DM-4103, and DM-4107 were quantified by liquid chromatography/electrospray ionization (ESI)–tandem mass spectroscopy. After thawing, liver tissue was mixed with saline and homogenized with a Shake Master NEO (BioMedical Science, Tokyo, Japan) prior to the pretreatment. Each sample was mixed with the IS and methanol in an ice-water bath. After centrifugation at 20,400g for 5 min at 10°C or below, the supernatant was analyzed by liquid chromatography/ESI–tandem mass spectroscopy. The liquid chromatography–tandem mass spectroscopy system consisted of a Prominence UFLC system (Shimadzu, Kyoto, Japan) and an API4000 triple quadrupole mass spectrometer equipped with an ESI source (AB SCIEX, Framingham, MA) and a switching valve (Valco Instruments, Houston, TX). Analyte separation was achieved on a reversed-phase high-performance liquid chromatography column XBridge C18 (3.5 μm, 2.1 mm i.d. × 50 mm; Waters, Milford, MA) by the binary gradient method with mobile phases of water-formic acid (1000:1, v/v) and acetonitrile-formic acid (1000:1, v/v) at a flow rate of 0.35 ml/min. The mass spectrometer was operated in ESI positive ion mode using multiple reaction monitoring transitions at mass-to-charge ratio (m/z) 449.2/252 for TVP, m/z 479.1/252 for DM-4103, m/z 481.1/252 for DM-4107, and m/z 463.2/266 for IS. Quantification of each analyte was achieved by comparison of the analyte/IS peak area ratios to those of a calibration curve (Analyst 1.5.1; AB SCIEX). Blank perfusate, blank bile, and blank liver homogenate served as matrices for the corresponding samples. TVP and metabolites were quantified in outflow perfusate, bile, and terminal liver samples from the IPL studies; these data were used to assess mass balance and conduct pharmacokinetic analyses. The lower limits of quantitation were 5nM, 50 nM, and 0.6 nmol/g tissue for the outflow perfusate, bile, and liver samples, respectively. Beal’s M5 method (Beal, 2001; Jusko, 2012) was used to manage individual values below the lower limit of quantitation.
Pharmacokinetic Modeling.
Modeling and simulation were used to evaluate the disposition of TVP and the two major metabolites in rat IPLs. A compartmental model assuming first-order processes was developed to describe the rate vs. time profiles using WinNonlin version 8.1 (Phoenix, St. Louis, MO), and using a mixed (combined additive and proportional) residual error model. The model scheme depicting the single-pass IPL system (Fig. 1) consisted of an extrahepatic and hepatic compartment, which was divided into five subcompartments to mimic liver physiology (Watanabe et al., 2009), and a bile compartment. The model was fit simultaneously to the average biliary excretion and outflow perfusate-time data for TVP and DM-4103 + DM-4107, for the WT IPLs and separately for the PCK IPLs. The individual outflow perfusate appearance and biliary excretion rates are presented in Supplemental Table 1. To describe the apparent delay in the biliary detection of TVP and DM-4103 + DM-4107 in PCK IPLs (likely due to the cyst burden in the PCK bile ducts), the bile compartment consisted of five transit compartments. The differential equations used to describe the model scheme in Fig. 1 are as follows, where the initial conditions in all compartments were assumed to be zero:
WT and PCK extracellular compartment 1 (TVP)(1)
WT and PCK extracellular compartments 2–5 (TVP)(2)
WT and PCK intracellular compartments 1–5 (TVP)(3)
WT and PCK concentration data in outflow perfusate expressed as dX/dt (TVP)(4)
WT concentration data in bile expressed as dX/dt (TVP)(5)
PCK transit compartment 1 following biliary excretion (TVP)(6)
PCK transit compartments 2–5 following biliary excretion (TVP)(7)
PCK concentration data in bile expressed as dX/dt (TVP)(8)
WT and PCK extracellular compartment 1 (DM-4103 + DM-4107)(9)
WT and PCK extracellular compartments 2–5 (DM-4103 + DM-4107)(10)
WT and PCK intracellular compartments 1–5 (DM-4103 + DM-4107)(11)
WT and PCK concentration data in outflow perfusate expressed as dX/dt (DM-4103 + DM-4107)(12)
WT concentration data in bile expressed as dX/dt (DM-4103 + DM-4107)(13)
PCK transit compartment 1 following biliary excretion (DM-4103 + DM-4107)(14)
PCK transit compartments 2–5 following biliary excretion (DM-4103 + DM-4107)(15)
PCK concentration data in bile expressed as dX/dt (DM-4103 + DM-4107)(16)
In eqs. 1–16, the variables and parameters are defined as in Fig. 1, with further explanation as follows: EC denotes extracellular; IC denotes intracellular; CL03,07 denotes the formation clearance of DM-4103 + DM-4107 from TVP; CLBile denotes biliary clearance of TVP; CLBile,03,07 denotes biliary clearance of DM-4103 + DM-4107; CLBL denotes basolateral clearance of TVP; CLBL,03,07 denotes basolateral clearance of DM-4103 + DM-4107; CLMet denotes formation clearance of other metabolites from TVP; CLUP denotes uptake clearance of TVP; KLag,Bile denotes the transit rate constant for the biliary excretion of TVP; KLag,Bile,03,07 denotes the transit rate constant for the biliary excretion of DM-4103 + DM-4107; CEC,n represents the extracellular concentration in the corresponding sinusoidal compartment, calculated as ; and VEC is the extracellular volume of the liver, which was assumed to be in equilibrium with the sinusoidal space and estimated as 20% of the total liver mass (Watanabe et al., 2009; Hobbs et al., 2012). Total liver mass was converted to total liver volume (VL) using the average density of rat liver (1.084 g/ml) (Pretlow and Pretlow, 1987). Here, Q represents the perfusate flow rate of 35 ml/min, which was assumed to be constant throughout the sinusoidal space of the liver; Cin represents the concentration of TVP in the inflow perfusate, measured for each preparation; binding of TVP to the perfusion tubing and apparatus was negligible (data not shown); and CIC,n and C03,07,IC,n represent the total intracellular concentrations of TVP and DM-4103+DM-4107, respectively, in the corresponding hepatocellular compartment. The total liver concentration was calculated as the sum of the mass in intra- and extracellular liver subcompartments 1–5 divided by VL. The total intracellular liver concentration was calculated as the sum of the mass in hepatocellular subcompartments 1–5, divided by the intracellular volume (VIC), calculated as . The perfusate flow rate multiplied by the concentration in the final extracellular liver compartment () was fit to the observed appearance rate in outflow perfusate, while the biliary excretion rate ( and for WT and PCK rats, respectively) was fit to the observed biliary excretion rate data. The hepatocellular concentrations of TVP and DM-4103 + DM-4107 were simulated using the final parameter estimates based on the extensive sampling of outflow perfusate and bile.
Data Analysis.
The majority of the data sets passed the Shapiro-Wilk normality test. All raw data are presented as mean ± S.D. from n = 5 livers in the PCK group, and n = 6 in the WT group. Statistical tests [t tests and Bonferroni-Dunn multiple comparisons tests (α = 0.05)] were performed using GraphPad Prism 7.02 for Windows (GraphPad Software, La Jolla, CA) to assess differences in: 1) concentrations of TVP, DM-4103, and DM-4107 in outflow perfusate at 30 and 60 min; 2) cumulative amounts of TVP, DM-4103, and DM-4107 in bile at the end of the perfusion; and 3) recovery (in terms of percentage of dose) of TVP, DM-4103, and DM-4107 in the outflow perfusate, bile, and liver. Final parameter estimates were considered significantly different if the 95% confidence intervals determined in Phoenix WinNonlin version 8.1 between the WT and PCK groups did not overlap. This method of judging the significance of differences is more conservative than the “standard” method (Schenker and Gentleman, 2001), which assumes consistency, asymptotic normality, and asymptotic independence of the estimates.
Results
At the end of liver perfusion, larger amounts of TVP (1.29-fold), DM-4103 (1.42-fold), and DM-4107 (1.45-fold) were recovered in PCK compared with WT rat livers, but differences were not statistically significant (Fig. 2). PCK rat livers weighed significantly more than those of WT rats (P = 0.0002; Supplemental Table 2), consistent with previous reports (Katsuyama et al., 2000; Lager et al., 2001; Shimomura et al., 2015).
Outflow Perfusate Concentrations of Tolvaptan and Metabolites.
After 30 min of TVP perfusion, when the perfusate was switched from TVP-containing to blank buffer, the TVP concentrations in outflow perfusate did not differ between PCK and WT IPLs (Fig. 3A). However, at the end of the washout phase, the outflow perfusate concentration of TVP was 2.33-fold higher in PCK compared with the WT rat livers (PCK = 522.4 ± 127.6 nM vs. WT = 224.6 ± 162.3 nM, P = 0.0089). DM-4103 concentrations in outflow perfusate trended higher for most time points in PCK compared with WT IPLs, especially toward the end of the perfusion, but these differences were not statistically significant (Fig. 3B). A similar trend was seen for DM-4107, but the concentrations at 30 min (3.94-fold increase in PCK) and at the end of the washout (3.69-fold increase in PCK) were significantly different (PCK = 62.57 ± 25.05 nM vs. WT = 15.89 ± 8.67 nM, P = 0.0020, and PCK = 93.04 ± 30.01 nM vs. WT = 25.23 ± 19.38 nM, P = 0.0014, respectively; Fig. 3C).
Biliary Excretion of Tolvaptan and Metabolites.
The overall bile flow in PCK rats was reduced ∼1.4-fold (data not shown). The cumulative amount of TVP, DM-4103, and DM-4107 recovered in the bile of PCK compared with WT rat livers tended to be lower throughout the perfusion (Fig. 4). The cumulative biliary recovery at 60 min was significantly lower for DM-4103 (P = 0.0038; Fig. 4B).
Total Recovery of Tolvaptan and Metabolites.
The total recovery of TVP, expressed as a percentage of dose, in liver tissue (Fig. 5A), outflow perfusate (Fig. 5B), and bile (Fig. 5C) at the end of the 60-min perfusions was similar between PCK and WT IPLs. DM-4107 recovery was significantly increased in the outflow perfusate of PCK (1.60% ± 0.57%) compared with WT (0.43% ± 0.29%) IPLs (P = 0.0017). DM-4103 recovery was significantly decreased in the bile of PCK (<0.01% ± <0.01%) compared with WT (0.02% ± 0.01%) IPLs (P = 0.0037). When assessing the metabolite/parent area under the curve ratios in liver tissue, outflow perfusate, and bile, only DM-4107 in outflow perfusate (PCK = 9.40% ± 2.63%; WT = 2.67% ± 1.98%) reached statistical significance (P = 0.0019), which coincides with Fig. 5B.
Pharmacokinetic Modeling.
Differential equations corresponding to that the model described in Fig. 1 were fit to the data to obtain estimates for the seven clearance parameters, and also for the two rate constants in the case of PCK rats (Table 1). The appearance rate in outflow perfusate, and the excretion rate in bile are plotted vs. time in Fig. 6. Additional model diagnostic plots are shown in Supplemental Fig. 1.
The biliary clearance value of TVP (CLBile) was significantly decreased in PCK compared with WT IPLs (2.19-fold; Table 1). Based on the confidence intervals of the remaining parameter estimates, no additional differences in clearance parameters were significant between PCK and WT IPLs. For those parameters describing clearance from a compartment of identical volume, additional comparisons were made within the same rat group (indicated by the # or † symbols in Table 1). In both PCK and WT groups, CLBL > CLBile suggesting a higher affinity of TVP for hepatic basolateral efflux transporters than for canalicular transporters. Additionally, in both groups, CLMet > CL03,07 indicating that the formation clearance of other metabolites is greater in rat IPLs than the formation clearance of DM-4103 + DM-4107, assuming that subsequent metabolism of these two chemical species is negligible. In both WT and PCK rat IPLs, CLBL > CLMet > CL03,07 > CLBile suggesting that the formation clearance of metabolites is more substantial than the biliary clearance of TVP, but not as extensive as the basolateral clearance of TVP.
Pharmacokinetic Simulations.
The model described in Fig. 1 also was used to simulate concentrations of TVP and DM-4103 + DM-4107 in the hepatocellular space of PCK and WT IPLs (Fig. 7). The model predicted similar hepatocellular concentrations of TVP at the end of the 30-min perfusion with TVP-containing buffer in the PCK and WT IPLs, whereas TVP concentrations were predicted to be higher in PCK compared with WT IPLs throughout the washout phase and beyond. Based on the simulations, hepatocellular concentrations of DM-4103 + DM-4107 were higher in PCK compared with WT IPLs throughout the 60-min perfusion and beyond.
Discussion
The current ex situ IPL study revealed that the hepatobiliary disposition of TVP and two metabolites was altered in a rodent model of ADPKD. Significantly more DM-4107 (percentage of TVP dose) was recovered in the outflow perfusate from PCK compared with WT IPLs. In contrast, significantly less DM-4103 (percentage of TVP dose) was recovered in the bile of PCK compared with WT IPLs. A pharmacokinetic model was developed to describe the disposition of TVP and metabolites in perfusate and bile and to estimate relevant clearance processes.
The estimated biliary clearance value for TVP (CLBile) was significantly decreased in PCK compared with WT IPLs, suggesting that this biliary excretion process is impaired in PCK rat livers. However, in both groups, basolateral efflux was a more prominent clearance mechanism for TVP (CLBL > CLBile). The formation clearance of DM-4103 + DM-4107 and other metabolites from TVP (CL03,07 and CLMet, respectively) was not significantly different between PCK and WT IPLs; however, in both groups CLMet > CL03,07, indicating that a substantial amount of other metabolites is formed. For example, DM-4103 and DM-4107 are formed through separate sequential metabolic pathways (Furukawa et al., 2014; Mazzarino et al., 2017). However, based on our data set, a more complex pharmacokinetic model could not be justified.
Using the final clearance parameter estimates, simulations revealed that TVP and metabolite concentrations at the end of the perfusion were increased in hepatocytes of PCK rats. These findings are significant because TVP, DM-4107, and especially DM-4103 inhibit hepatic efflux transporters (Slizgi et al., 2016). Inhibition of hepatic efflux transporters may be greater in PCK disease compared with healthy controls due to elevated hepatocellular TVP and metabolite concentrations, resulting in accumulation of hepatotoxic endogenous molecules (e.g., lithocholate, glycine-conjugated deoxycholate, and chenodeoxycholate) (Galle et al., 1990; Fickert et al., 2006; Chatterjee et al., 2014) and leading to liver injury. Therefore, increased hepatocellular concentrations of TVP and DM-4103 + DM-4107 in PCK disease could be related to the higher incidence of hepatotoxic events in patients with ADPKD treated with TVP in the pivotal clinical trial. Bile acid concentrations are increased in PCK rat livers (Munoz-Garrido et al., 2015; Brock et al., 2018), and elevated hepatocyte concentrations of TVP and metabolites may exacerbate this imbalance. Although the IC50 values obtained from previous studies with human hepatic transporters (Slizgi et al., 2016) cannot be compared directly to the simulated hepatocellular TVP and DM-4103 + DM-4107 concentrations in rat IPLs due to species differences, the simulations illustrate the possibility that these chemical species, which are capable of inhibiting hepatic transporters, could accumulate to a greater degree in patients with ADPKD compared with healthy controls. Hepatocellular accumulation of TVP and metabolites may take time to reach inhibitory concentrations, which could be related to ADPKD progression, and this could be one explanation for the delay in TVP toxicity in patients with ADPKD (Watkins et al., 2015).
All rats used in the current study were 20- to 21-week-old males. Although liver function is thought to be preserved in polycystic liver disease (Cnossen and Drenth, 2014), hepatic cyst volume in PCK rats increases with age (Brock et al., 2018), potentially resulting in age-dependent changes in hepatobiliary disposition of compounds. A definitive correlation between hepatic cyst involvement and clinically observed TVP-mediated hepatotoxicity cannot be made. However, of those subjects that were adjudicated as highly likely or probably related to study drug, 76% had some level of hepatic cyst burden (Watkins et al., 2015). Furthermore, in TEMPO 3:4 and 4:4, all three Hy’s law cases and 18 out of 30 subjects with hepatotoxicity attributable to TVP with a probable (or higher) likelihood, were women (Watkins et al., 2015). However, sex differences in TVP-mediated hepatotoxicity could not be concluded.
PCK rat livers weighed significantly more than WT rat livers, consistent with prior studies (Katsuyama et al., 2000; Lager et al., 2001). These differences may be due to the cysts in PCK rat livers (Brock et al., 2018). Interestingly, bile flow in PCK IPLs was slightly lower than in WT IPLs during perfusion. One possible explanation may be that cystic bile ducts in PCK rat livers are dilated (Masyuk et al., 2004). Similarly, renal tubules in PKD are dilated by cyst formation and proliferation (Igarashi and Somlo, 2007). The resulting structural compression in the tubules hinders urine flow (Grantham et al., 2011). Following this logic, it is conceivable that bile flow in the cystic bile ducts of PCK rats is reduced.
Considerable variability was observed in the TVP and metabolite excretion rates in outflow perfusate and bile between individual perfused livers in each group of rats (Supplemental Table 1). Liver function may be inherently different from rat to rat, especially in the case of PCK rats where hepatic cyst burden can differ markedly at a specific age (Brock et al., 2018). Despite variability between individual IPLs, the average data for WT and PCK rat IPLs revealed distinct differences and trends between the two groups.
In conclusion, the hepatobiliary disposition of TVP, DM-4103, and DM-4107 was modified in PCK compared with WT IPLs. The altered disposition resulted in higher simulated concentrations of these chemical species in PCK rat livers, which could impair transport of endogenous compounds such as hepatotoxic bile acids. We hypothesize that this is a contributing mechanism in relation to the increase in liver injury observed in some patients with ADPKD who are administered TVP.
Acknowledgments
We gratefully acknowledge Certara (Princeton, NJ) for providing academic research licenses for the Phoenix software as part of the Center of Excellence program for educational institutions. We thank Dr. Yang Lu for conducting the pilot study, Dr. Brian C. Cooley for surgical expertise, and Dr. Noora Sjöstedt for support with the visual abstract.
Authorship Contributions
Participated in research design: Beaudoin, Roth, Brock, Brouwer.
Conducted experiments: Beaudoin, Mizuno.
Performed data analysis: Beaudoin, Bezençon, Cao, Roth, Brock, Mizuno, Brouwer.
Wrote or contributed to the writing of the manuscript: Beaudoin, Bezençon, Brouwer.
Footnotes
- Received August 10, 2018.
- Accepted November 28, 2018.
This work was supported, in part, by Otsuka Pharmaceutical Development & Commercialization, Inc., and the National Institutes of Health National Institute of General Medical Sciences [Grants R01GM041935 and R35GM122576].
This work was presented, in part, at the 21st North American Meeting of the International Society for the Study of Xenobiotics in 2017.
Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the authors and do not necessarily reflect the views of the Otsuka Pharmaceutical Group or the National Institutes of Health.
↵This article has supplemental material available at dmd.aspetjournals.org.
Abbreviations
- ADPKD
- autosomal dominant polycystic kidney disease
- CLBile
- biliary clearance of tolvaptan
- CLBile,03,07
- biliary clearance of DM-4103 + DM-4107
- CLBL
- basolateral clearance of tolvaptan
- CLBL,03,07
- basolateral clearance of DM-4103 + DM-4107
- CLMet
- formation clearance of other metabolites from tolvaptan
- CLUP
- uptake clearance of tolvaptan
- CL03,07
- formation clearance of DM-4103 + DM-4107 from tolvaptan
- EC
- extracellular
- ESI
- electrospray ionization
- IC
- intracellular
- IPL
- isolated perfused liver
- IS
- internal standard
- KLag,Bile
- transit rate constant for the biliary excretion of tolvaptan
- KLag,Bile,03,07
- transit rate constant for the biliary excretion of DM-4103 + DM-4107
- LDH
- lactate dehydrogenase
- m/z
- mass-to-charge ratio
- PCK
- polycystic kidney
- PKD
- polycystic kidney disease
- TEMPO
- Tolvaptan Efficacy and Safety in Management of Autosomal Dominant Polycystic Kidney Disease and Its Outcomes
- TVP
- tolvaptan
- WT
- wild-type
- Copyright © 2019 by The American Society for Pharmacology and Experimental Therapeutics