Abstract
Hepatobiliary excretion mediated by transporters, organic anion-transporting polypeptide (OATP) 1B1 and multidrug resistance-associated protein (MRP) 2, is the major elimination pathway of an HMG-CoA reductase inhibitor, pravastatin. The present study examined the effects of changes in the transporter activities on the systemic and liver exposure of pravastatin using a physiologically based pharmacokinetic model. Scaling factors, determined by comparing in vivo and in vitro parameters of pravastatin in rats for the hepatic uptake and canalicular efflux, were obtained. The simulated plasma and liver concentrations and biliary excretion profiles were very close to the observed data in rats under linear and nonlinear conditions. In vitro parameters, determined in human cryopreserved hepatocytes and canalicular membrane vesicles, were extrapolated to in vivo parameters using the scaling factors obtained in rats. The simulated plasma concentrations of pravastatin were close to the reported values in humans. Sensitivity analyses showed that changes in the hepatic uptake ability altered the plasma concentration of pravastatin markedly but had a minimal effect on the liver concentration, whereas changes in the ability of canalicular efflux altered the liver concentration of pravastatin markedly but had a small effect on the plasma concentration. In conclusion, the model allows the prediction of the disposition of pravastatin in humans. The present study suggests that changes in the OATP1B1 activities may have a small and a large impact on the therapeutic efficacy and side effect (myopathy) of pravastatin, respectively, whereas those in the MRP2 activities may have opposite impacts (i.e., large and small impacts on the therapeutic efficacy and side effect).
Predicting the disposition of drugs in humans, particularly in the early stages of drug development, has been a critical issue in selecting the proper candidate drugs because the exposure of drugs to target organs is the major factor determining their pharmacological and/or toxicological activity. Human liver microsomes allow the reliable prediction of the metabolic clearance of drugs in humans (Rane et al., 1977; Iwatsubo et al., 1997; Obach, 1999; Naritomi et al., 2001). Biliary excretion, another hepatic elimination pathway, is the major systemic elimination pathway, particularly for amphipathic anionic drugs such as HMG-CoA reductase inhibitors (statins) and angiotensin II receptor antagonists. Because multiple transporters on the sinusoidal and canalicular membranes are involved, it is necessary to separately determine three kinetic parameters: 1) uptake, 2) sinusoidal efflux, and 3) canalicular efflux, to predict biliary clearance with regard to the plasma concentration (Giacomini and Sugiyama, 2005; Shitara et al., 2006a). The uptake clearance determined in freshly isolated rat hepatocytes correlates well with that determined with the multiple indicator dilution method (Miyauchi et al., 1993), and rat hepatocytes are reported to be a useful tool for predicting the hepatic clearance of drugs with significant hepatic uptake (Soars et al., 2007). Although cryopreserved human hepatocytes and canalicular membrane vesicles (CMVs) are commercially available, their usefulness in predicting in vivo hepatic uptake and canalicular efflux clearance remains to be examined. No method to quantify in vitro sinusoidal efflux has yet been established.
A physiologically based pharmacokinetic (PBPK) model, in which compartments representing tissues are connected with the blood flow, has been used to predict the time profiles of plasma and tissue concentrations (Kawai et al., 1998; Jones et al., 2006). The PBPK model is quite useful for simulating the effects of drug-drug interactions and genetic variations in drug-metabolizing enzymes and transporters on the exposure of drugs to the blood and organs and, ultimately, their effects on the pharmacological actions of drugs (Jones et al., 2006; Shitara and Sugiyama, 2006a). The purpose of this study was to establish a PBPK model to describe the disposition of pravastatin for which transporters are deeply involved in its hepatobiliary transport. Pravastatin, one of the statins used for the treatment of hyperlipidemia, was selected as the model compound in this study. The liver is a target organ for the pharmacological actions of statins, whereas myotoxic adverse effects, sometimes severe, including myopathy or rhabdomyolysis, are associated with the use of statins. Therefore, it is very important to simulate the exposure of statins to the liver and skeletal muscle to predict their pharmacological and toxicological effects. Hepatobiliary transport is the main elimination pathway of pravastatin from the systemic circulation and is mediated by uptake and efflux transporters in the liver (Shitara and Sugiyama, 2006b). The hepatic uptake of pravastatin is mainly mediated by organic anion-transporting polypeptide (OATP) 1B1, and its biliary excretion is predominantly mediated by multidrug resistance-associated protein (MRP) 2 (Yamazaki et al., 1993, 1997; Nakai et al., 2001). Pravastatin undergoes urinary excretion by tubular secretion and by glomerular filtration in humans (Singhvi et al., 1990). Organic anion transporter 3 has been suggested to be responsible for the basolateral uptake of pravastatin in rats and humans (Hasegawa et al., 2002; Nakagomi-Hagihara et al., 2007), whereas the transporter involved in its luminal efflux is yet to be identified.
In this study, in vivo experiments were carried out using male rats to obtain concentration-time profiles of pravastatin in the plasma, liver, kidney, muscle, brain, and lung. The kinetic parameters for the hepatic uptake and canalicular efflux of pravastatin were determined from in vitro transport studies using freshly isolated rat hepatocytes and CMVs, respectively. In vitro-in vivo scaling factors (SFs) were obtained for the hepatic uptake and subsequent canalicular efflux of pravastatin in rats. A PBPK model was constructed to simulate the systemic and liver exposure of pravastatin in rats. Using the PBPK model, the SFs determined in rats and kinetic parameters determined using human materials, the plasma concentration-time curve of pravastatin in humans was also simulated. Finally, the effects of changes in these transporter activities, caused by genetic polymorphisms and drug-drug interactions, on the concentration profiles of pravastatin in plasma and the liver were examined using the PBPK model.
Materials and Methods
Materials
[3H]Pravastatin (45.5 Ci/mmol), unlabeled pravastatin, and a pravastatin analog, R-122798, were provided by Daiichi Sankyo Co., Ltd. (Tokyo, Japan). Cryopreserved human hepatocytes and human liver S9 fractions were purchased from In Vitro Technologies (Baltimore, MD). Human liver S9 fractions were also purchased from XenoTech, LLC (Lenexa, KS) and Tissue Transformation Technology (Edison, NJ). All other chemicals and reagents were of analytical grade and were readily available from commercial sources.
Animals
Male Sprague-Dawley rats (6–7 weeks old) were purchased from Nippon SLC (Hamamatsu, Japan). All animals were maintained under standard conditions with a reversed light/dark cycle and were treated humanely. Food and water were available ad libitum. The studies were carried out in accordance with the guidelines of the Institutional Animal Care Committee, Graduate School of Pharmaceutical Sciences, The University of Tokyo, Tokyo, Japan.
Animal Experiments
Male Sprague-Dawley rats, weighing approximately 250 to 320 g, were used throughout the experiments. Under ether anesthesia, the femoral artery was cannulated with a polyethylene catheter (SP-31) for the collection of blood samples. The bile duct was cannulated with a polyethylene catheter (PE-10) for bile collection, and the bladder was cannulated with a silicon catheter to collect urine. The femoral vein or the duodenum was cannulated with a polyethylene catheter (SP-31) for the administration of pravastatin. Each rat was placed in a Bollman cage and allowed to recover from the anesthesia before the experiments were continued. The rats were given pravastatin intravenously at 0.2, 1, 10, 50, or 200 mg/kg or intraduodenally at 20 mg/kg. Blood samples were collected at the designated times and centrifuged at 1500g for 10 min at 4°C to obtain plasma. Bile and urine samples were collected in preweighed test tubes at the designated intervals throughout the experiment. After the last blood sample had been taken, each rat was killed, and the liver, kidney, brain, lungs, and skeletal muscle were excised immediately for the tissue distribution study. The tissues were weighed and flash frozen in liquid nitrogen. All the samples were stored at −20°C until quantification.
Transport Study Using Human Cryopreserved Hepatocytes
This experiment was performed as described previously (Shitara et al., 2003). In brief, immediately before the study, the hepatocytes were thawed at 37°C. After they had been washed twice with ice-cold Krebs-Henseleit buffer, the cells were resuspended in Krebs-Henseleit buffer to a cell density of 1.0 × 106 viable cells/ml for the uptake study. After preincubation of the cells (1.2 × 105 cells/reaction) at 37°C for 3 min, drug uptake was initiated by the addition of labeled and unlabeled substrates to the cell suspension. The reaction was terminated after 0.5 or 2 min by separating the cells from the substrate solution. For this purpose, an aliquot of 100 μl of incubation mixture was placed in a centrifuge tube (450 μl) containing 50 μl of 2 N NaOH under a layer of 100 μl of oil (density = 1.015, a mixture of silicone oil and mineral oil; Sigma-Aldrich, St. Louis, MO). The sample tube was centrifuged for 10 s in a tabletop centrifuge (10,000g; Beckman Microfuge E; Beckman Coulter, Fullerton, CA). After overnight incubation in alkali to dissolve the hepatocytes, the centrifuge tube was cut, and each phase was transferred to a scintillation vial. The phase containing the dissolved cells was neutralized with 50 μl of 2 N HCl, mixed with scintillation cocktail, and its radioactivity was measured in a liquid scintillation counter (LS6000SE; Beckman Coulter). The time course for the uptake of [3H]pravastatin into hepatocytes was expressed as the uptake volume (microliters per 106 viable cells) of the radioactivity taken up into the cells (disintegrations per minute per 106 cells) divided by the concentration of radioactivity in the incubation buffer (disintegrations per minute per microliter). The initial uptake velocity of [3H]pravastatin was calculated from the slopes of the uptake volume versus time plots obtained at 0.5 and 2 min and expressed as the uptake clearance (microliters per minute per 106 cells).
Metabolism Study Using the Liver S9 Fraction
It has been reported that pravastatin is metabolized by sulfotransferase in male rats (Kitazawa et al., 1993). Therefore, we used the liver S9 fraction as the enzyme source. The liver S9 fraction was prepared from four rats using standard procedures and stored at −80°C until use. The protein concentration was determined by the Lowry method, using bovine serum albumin as the standard. Pravastatin was incubated with a reaction mixture consisting of rat liver S9 fraction (final concentration, 8 mg/ml), NADPH-generating system (0.8 mM NADP+, 8 mM glucose 6-phosphate, 1 U/ml glucose-6-phosphate dehydrogenase, and 3 mM MgCl2), and 3′-phosphoadenosine 5′-phosphosulfate (final concentrations, 0.5 and 5 mM for low and high pravastatin concentrations, respectively) in the presence of 100 mM phosphate buffer, pH 7.4. After preincubation at 37°C for 10 min, pravastatin (final concentration, 0.1–500 μM) was added to initiate the enzyme reaction. At the designated time, the reactions were terminated by mixing them with equal volumes of methanol containing R-122798, an analytical internal standard, followed by centrifugation at 15,000g for 10 min at 4°C. The supernatant was subjected to liquid chromatography/mass spectrometry (LC/MS) analysis. In studies with the human liver S9 fraction, the concentrations of pravastatin and 3′-phosphoadenosine 5′-phosphosulfate were 5 μM and 1 mM, respectively; other incubation conditions were the same as in the rat studies.
LC/MS Analysis
Liver, kidney, brain, lung, and skeletal muscle were added to 3 to 5 volumes of physiological saline (w/v) and homogenized. Tissue homogenates and plasma, bile, and urine samples were deproteinated with 2 volumes of methanol containing the internal standard (1 μg/ml R-122798) and centrifuged at 15,000g for 10 min at 4°C. High-concentration samples were diluted appropriately with blank matrix before deproteination. The supernatant was subjected to LC/MS analysis. The appropriate standard curves were prepared in the equivalent blank matrix and used for each analysis.
The LC/MS consisted of an Alliance HT 2695 separation module with an autosampler (Waters, Milford, MA) and a Micromass ZQ mass spectrometer with an electron ion spray interface (Waters). The optimum operating conditions used were as follows: electrospray probe (capillary) voltage, 3.2 kV; sample cone voltage, 20 V; and source temperature, 100°C. The spectrometer was operated at a drying desolvation gas flow rate of 350 l/h. The mass spectrometer was operated in the selected ion monitoring mode using the respective MH− ions, m/z 423.3 for pravastatin and m/z 409.3 for the internal standard. The mobile phase used for high-performance liquid chromatography was acetonitrile/ammonium acetate buffer (10 mM), pH 4 = 7:3 (v/v), and the flow rate was 0.3 ml/min. Chromatographic separation was achieved on a C18 column (Inertsil ODS-3 column, 50 × 2.1 mm; particle size, 3 μm) (GL Sciences, Tokyo, Japan).
Data Analysis of Metabolic Clearance in Liver S9
The metabolic velocity was calculated from the slope of the natural log (concentration)-time plot. Because the Eadie-Hofstee plot showed curvature, the kinetic parameters were obtained using eq. 1:
where v is the initial velocity (picomoles per minute per milligram of protein), S is the substrate concentration (micromolar), Vmax1 and Vmax2 are the maximum velocities (picomoles per minute per milligram of protein), and Km1 and Km2 are the Michaelis constants (micromolar). Fitting was performed with the nonlinear least-squares method using the MULTI program (Yamaoka et al., 1981). The input data were weighted as the reciprocals of the observed values, and the Damping Gauss-Newton algorithm was used for fitting.
Model Development
The PBPK model was constructed to describe the pharmacokinetics of pravastatin in rats and humans (Fig. 1). The key features of this model are as follows. 1) Active uptake (PSinf) and passive diffusion clearances (PSdif) on the sinusoidal membrane, and biliary clearance (PSbile) on the canalicular membrane in the liver are incorporated. 2) The liver compartment consists of five units of extracellular and subcellular compartments, connected by blood flow in tandem, to fit the hepatic disposition to the “dispersion” model. Because the hepatic elimination of pravastatin in rats is blood flow limited, the dispersion model is the appropriate model for the hepatic elimination of such high-clearance drugs (Roberts and Rowland, 1986; Iwatsubo et al., 1997; Naritomi et al., 2001). The number of liver compartments was determined by comparing the hepatic availability (Fh,n) and Fh predicted using the dispersion model. Fh,n is the product of the availability in the liver compartments (eq. 2).
where n represents the number of compartments. The integer number n, which gave the Fh,n value closest to that in the dispersion model, was selected. 3) The brain and muscle, target tissues for the adverse effects of statins, were included. 4) Although urinary excretion is a minor elimination pathway in male rats, the kidney was included because the kidney/blood concentration ratio for pravastatin is high in male rats, probably because of the efficient uptake and/or reabsorption of pravastatin. The renal clearance of pravastatin in male rats was lower than the glomerular filtration rate corrected by the blood unbound fraction. In contrast, renal clearance must be taken into consideration in humans. Because this study focused on hepatobiliary transport, renal elimination occurs from the systemic compartment. 5) The rapid equilibrium distribution of pravastatin between the blood and tissues other than the liver was assumed. 6) The initial distribution volume, estimated by fitting the plasma concentration time profiles of pravastatin in rats after the intravenous administration of 0.2 and 1 mg/kg to a two-compartment model, was used as the volume of the rapid equilibrium compartment including the blood compartment, and it was assumed that there is no interspecies difference in the initial distribution volume. The differential equations are shown in Appendix I, and all simulations were performed with SAAM II (SAAM Institute, Seattle, WA).
Schematic diagram of the PBPK model predicting the concentration-time profiles of pravastatin. The liver compartment was divided into five compartments to mimic the dispersion model. Indicated are the blood flow (Q), the active hepatic uptake clearance (PSinf), the passive diffusion clearance (PSdif), the biliary clearance (PSbile), and the metabolic clearance (CLmet,int), human (H), and rat (R). The enterohepatic circulation was incorporated in the case of humans.
Estimation of Kinetic Parameters Used in the Simulation
In Vitro Parameters (Rats and Humans).
The in vitro active uptake (PSinf,vitro) and passive diffusion clearances (PSdif,vitro) of pravastatin on the sinusoidal membrane were determined from the uptake studies using isolated hepatocytes. The parameters for rats were taken from previous reports (Yamazaki et al., 1993; Ishigami et al., 1995), and those for humans were determined in the present study. PSinf,vitro and PSdif,vitro were regarded as the saturable and nonsaturable components, respectively, in the uptake clearance into hepatocytes. A physiological scaling factor of 1.2 × 108 cells/g liver was used for scaling up to the organ level (Iwatsubo et al., 1997). The in vitro biliary clearance (PSbile,vitro) of pravastatin was calculated from the ATP-dependent uptake clearance into the CMVs using eq. 3 (Niinuma et al., 1999):
where Vinitial represents the velocity of the initial ATP-dependent uptake by CMVs corrected by medium concentration (6.08 μl/min/mg protein for rats and 1.90 μl/min/mg protein for humans), R represents the recovery of liver homogenate protein (174 mg homogenate protein/g liver for rats and 133 mg homogenate protein/g liver for humans), E represents the enrichment of the CMV fraction (70.4 for rats and 61.8 for humans), and IO represents the population of inside-out CMVs (0.347 for rats and 0.555 for humans).
In Vivo Parameters (Rats).
The in vivo intrinsic biliary clearance (PSbile,vivo) at the canalicular membrane was calculated by dividing the biliary excretion rate by the hepatic unbound concentration at steady state (Yamazaki et al., 1996b, 1997). Systemic elimination other than biliary excretion was regarded as the hepatic metabolism because renal elimination in male rats is negligible. Thus, in vivo intrinsic metabolic clearance (CLmet,int,vivo) was obtained with eq. 4:
The in vivo passive diffusion clearance on the sinusoidal membrane was assumed to be the same as PSdif,vitro. The in vivo active uptake clearance (PSinf,vivo) was estimated using eq. 5:
where CLint,all represents the overall hepatic intrinsic clearance estimated from the hepatic availability using the dispersion model, with a dispersion number of 0.17, which was obtained by dividing the bioavailability by the fraction absorbed (Komai et al., 1992), assuming negligible metabolism in the small intestine. The average of the tissue/blood concentration ratios at 30, 60, and 90 min after the intravenous administration at 10 mg/kg pravastatin were used as the tissue/blood partition coefficient (Kp), assuming a pseudo-steady state (Table 2). Actually, the tissue/blood concentration ratios at 30, 60, and 90 min were similar (muscle, 0.28, 0.21, and 0.18; brain, 0.045, 0.029 and 0.034; kidney, 13, 14, and 15; lung, 0.76, 0.67, and 0.77 at 30, 60, and 90 min, respectively). The absorption rate constants were estimated by noncompartment analysis using the plasma concentration data.
Results
In Vivo Pharmacokinetics of Pravastatin in Rats.
Figure 2 shows time profiles of the plasma concentration of pravastatin after its intravenous (0.2 mg/kg) and intraduodenal (20 mg/kg) administration and the cumulative amount of pravastatin excreted into the bile. The total blood clearance (CLtot,B) was similar to the hepatic blood flow rate. The bioavailability of pravastatin after intraduodenal administration was calculated to be 0.0087 by comparing the AUC for pravastatin after intravenous and intraduodenal administration. Forty-six percent of the dose was recovered in the bile as the parent compound after intravenous administration, whereas the amount excreted into the urine was less than 4% of the dose. Even after intraduodenal administration, 33% was recovered in the bile.
Simulated and observed plasma concentrations and biliary excretion rates for pravastatin in rats after intravenous (●, 0.2 mg/kg) or intraduodenal (○, 20 mg/kg) administration. The symbols and solid lines represent experimentally observed and simulated values, respectively. Each point represents the mean ± S.E. (n = 3).
The nonlinearity of the disposition of pravastatin was examined. The plasma concentrations and cumulative amounts excreted into the bile after its intravenous administration were determined at doses ranging from 0.2 to 200 mg/kg (Fig. 3). CLtot,B was independent of the dose up to 50 mg/kg but decreased to 27 ml/min/kg at 200 mg/kg pravastatin. The cumulative biliary excretion increased slightly from 46 to 60% at doses above 0.2 mg/kg and was significantly delayed at 200 mg/kg.
Simulated and observed plasma concentrations and biliary excretion rates for pravastatin in rats after the intravenous administration of various doses. The symbols and lines represent experimentally observed and simulated values, respectively. Each point represents the mean ± S.E. (n = 3). ———- ●, 0.2 mg/kg; - - - - ○, 1 mg/kg; · · · · · · · · · · , 10 mg/kg; —- · —- · —- ▵, 50 mg/kg; —- · · —- · ♦, 200 mg/kg.
Hepatic Metabolism of Pravastatin in Rats.
The metabolism of pravastatin in the liver was examined using S9 fractions prepared from rat liver. It exhibited biphasic kinetics with high-affinity (Km1, 0.846 ± 0.403 μM; Vmax1, 4.47 ± 1.92 pmol/mg/min) and low-affinity (Km2, 80.3 ± 12.6 μM; Vmax2, 240 ± 16.2 pmol/mg/min) components (mean ± S.D.). The sum of the in vitro metabolic clearance for the high- and low-affinity components, corrected with the physiological scaling factor of 96.1 mg protein/g liver, was used as the in vitro metabolic clearance (CLmet,int,vitro), which was 0.793 ml/min/g liver (Table 1).
Kinetic parameters for hepatic intrinsic clearance
Active hepatic uptake and passive diffusion clearances on the sinusoidal membrane, biliary clearance on the canalicular membrane, and metabolic clearance were estimated by both in vitro and in vivo experiments. The details of these estimations are described in the text. Values within parentheses indicate the Km value (micromolar) for each clearance.
Simulation of Concentration-Time Profiles of Pravastatin in Rats.
All parameters used in the simulation are summarized in Tables 1 and 2. The initial distribution volume was estimated to be 0.393 l/kg from the plasma concentration profile after the intravenous administration of 0.2 mg/kg pravastatin, which was used as the volume of the rapid equilibrium compartment in the model. Figures 2 and 3 show the simulated plasma concentrations and biliary excretion time profiles for pravastatin, together with the observed data after its intravenous and intraduodenal administration. To reproduce the in vivo pharmacokinetic profiles using in vitro parameters, SFs were necessary. For the in vitro-in vivo extrapolation of the transporter-mediated clearances, the ratio of the in vivo/in vitro intrinsic clearances of each process in rats was given as the SF (Table 1). Furthermore, using the Km values, the simulated plasma concentration and biliary excretion time profiles gave similar values to the observed data, even under nonlinear conditions. This model reproduced the time profiles for pravastatin in the liver and other peripheral tissues after its intravenous administration (Figs. 4 and 5). In particular, the nonlinearity of the liver concentration-time profiles could also be simulated (Fig. 4). Although the model reasonably describes the experimental data, the simulated lines showed some deviation from the observed data at the terminal phase in Figs. 2 (left) and 5. This may be caused by the lack of a compartment corresponding to the organ, which is associated with the terminal phase of systemic pravastatin. Moreover, the simulation results of biliary excretion of pravastatin administered at 1 and 10 mg/kg showed some deviation from the observed data (Fig. 3, right). Because hepatic metabolism is saturated at these doses, this may be attributed to the deviation of the Km value for metabolism.
Physiological and kinetic parameters for modeling in rats and humans
Simulated and observed liver concentration profiles for pravastatin in rats after intravenous administration. Dashed and solid lines, simulated plasma and liver concentrations, respectively. Open and closed symbols, experimentally observed plasma and liver concentrations, respectively (circles, 10 mg/kg; triangles, 200 mg/kg). Each point represents the mean ± S.E. (n = 3).
Simulated and observed tissue concentration profiles for pravastatin in rats after intravenous administration at 10 mg/kg. Symbols and solid lines, experimentally observed and simulated values, respectively. Each point represents the mean ± S.E. (n = 3).
Prediction of Pharmacokinetics in Humans.
The uptake clearance determined using eight lots of human cryopreserved hepatocytes was 4.5 ± 2.9 μl/min/106 cells at 1 μM pravastatin and 0.77 ± 0.63 μl/min/106 cells at 100 μM pravastatin (mean ± S.D.). Using the physiological scaling factor of 1.2 × 108 cells/g liver, PSinf,vitro,human and PSdif,vitro,human were calculated to be 0.448 and 0.0924 ml/min/g liver, respectively (Table 1). Unlike the rat liver S9 fraction, no metabolism of pravastatin was observed up to 180 min in the human liver S9 fractions purchased from three different vendors. Thus, the hepatic metabolism of pravastatin might be negligible in the human liver. The in vivo kinetic parameters for pravastatin in humans were predicted by multiplying the corresponding in vitro parameters obtained using human materials by the SF obtained from rat studies. For PSbile, the saturable (ATP-dependent) biliary clearance in humans was predicted as described above (eq. 3), and the nonsaturable component of the biliary clearance in humans was assumed to be the same as that in rats. Thus, the predicted PSbile, vivo, human was 0.388 ml/min/g liver. Assuming that the distribution of pravastatin to the tissues, except the liver, occurs by passive diffusion, the tissue/blood partition coefficient (Kp) was calculated by the following equation: Kp = fB/fT, where fB and fT represent the blood unbound fraction (=plasma unbound fraction/blood-to-plasma concentration ratio) and the unbound fraction in the tissues, respectively. It was assumed that there is no species difference in fT between rats and humans based on the previous report by Sawada et al. (1985). The estimated or reported physiological, anatomical, and kinetic parameters for humans used in the simulation are shown in Tables 1 and 2. Using these parameters, the plasma concentration-time profiles for pravastatin in humans after intravenous or oral administration were predicted. A lag time of 17 min was taken into consideration in the simulation of oral administration. The predicted concentration-time profiles were similar to the observed data (Fig. 6).
Predicted and observed plasma concentration profiles for pravastatin in humans. Closed and open symbols, reported plasma concentrations after intravenous (9.9 mg) and oral (19.2 mg) administration, respectively (Singhvi et al., 1990). Solid lines, simulated values using the parameters shown in Tables 1 and 2.
Effect of Transporter Activity on Systemic and Target Exposure.
Sensitivity analyses were performed to understand the effects of the changes in transporter activities on the time profiles for the plasma and liver (a target organ) concentrations of pravastatin in humans. The plasma and liver concentrations after the oral administration (40 mg) of pravastatin were simulated using the PBPK model constructed in this study, with varying hepatic transport activities over a range of 0.33 to 3.0 times the initial value. The simulated concentration-time profiles and the changes in the AUC are shown in Fig. 7 and Table 3, respectively. Changes in the active hepatic uptake ability affected the plasma concentration profiles dramatically but did not greatly affect the liver concentration profiles. On the contrary, changes in the ability of canalicular efflux altered the liver concentration of pravastatin markedly but had a small effect on the plasma concentration. Changes in the passive diffusion clearance hardly affect the plasma and the liver concentration profiles.
Effects of changes in transporter activity on the time profiles of plasma and liver (target organ) concentrations of pravastatin in humans. Plasma and liver concentrations after oral administration (40 mg) were simulated using the PBPK model with varying hepatic transport activities over a 1/3- to 3-fold range of the initial values shown in Table 1. (———-, initial; – – –, ×1/3; · · · · · · · , ×3)
Changes in the AUC (percentage of the control) for plasma and liver concentrations of pravastatin after its oral administration when the transporter function changes
Discussion
It is now well recognized that drug transporters play important roles in the processes of absorption, distribution, and excretion (Giacomini and Sugiyama, 2005; Shitara et al., 2006a). The purpose of this study was to construct a PBPK model to evaluate the concentration-time profiles for drugs in the plasma and peripheral organs in humans using physiological parameters, SFs, and drug-related parameters (unbound fraction and metabolic and membrane transport clearances extrapolated from in vitro experiments). The principle of the prediction was as follows. First, SFs were obtained by comparing in vitro and in vivo parameters in rats. Then, the in vitro human parameters were extrapolated in vivo using the SFs obtained in rats (Naritomi et al., 2001). Pravastatin was selected as the model compound because many studies have investigated the mechanisms involved in the drug disposition in rodents, and clinical data after intravenous and oral administration are available.
Consistent with a previous report (Yamazaki et al., 1996a), the hepatic elimination of pravastatin is blood flow limited. Considering that the maximum amount of intact pravastatin excreted into the bile was 50%, it is likely that pravastatin undergoes hepatic metabolism in rats because pravastatin is excreted negligibly in the urine. Incubating pravastatin with the rat liver S9 fractions caused a reduction in intact pravastatin with time and consisted of two different mechanisms with high- and low-affinity sites. The kinetic parameters related to hepatic clearance (PSinf, PSdif, PSbile, and CLmet,int) were estimated from various in vivo experiments and were incorporated into the PBPK model. As a result, plasma concentration and biliary excretion-time profiles for pravastatin were successfully reproduced (Figs. 2⇑⇑–5). Moreover, nonlinear pharmacokinetics were also reproduced using the Km values for hepatic uptake, biliary excretion, and metabolic clearances (Fig. 3). The liver concentrations of pravastatin were similar to the observed data, even under nonlinear conditions (200 mg/kg) (Fig. 4). These results suggest that the PBPK model constructed in this study is appropriate for describing the pharmacokinetics of pravastatin in rats.
The kinetic parameters PSinf, PSdif, and PSbile were also determined in vitro using rat hepatocytes and CMVs to obtain the relevant SFs (Table 2). The corresponding parameters were also determined using human cryopreserved hepatocytes and CMVs. These parameters were extrapolated in vivo using the SFs determined in rats. Because there is no evidence that active transport mechanisms are involved in the sinusoidal efflux of pravastatin, the clearance corresponding to the nonsaturable component (PSdif) of the uptake was used as the clearance for sinusoidal efflux. Unlike the rat liver S9 fractions, pravastatin was not metabolized in the human liver S9 fractions. Therefore, the metabolic clearance was set to zero in humans. Using the human parameters, simulated plasma concentration-time profiles of pravastatin after the intravenous and oral administration were fairly close to the observed data for humans (Fig. 6), showing that the predicted value was not far from the true value. It should be noted that the sinusoidal efflux clearance (passive diffusion clearance) was lower than the intrinsic biliary clearance with regard to the liver concentration, indicating that the hepatobiliary transport of pravastatin is likely uptake limited and that the hepatic intrinsic clearance can be approximated to PSinf (Shitara et al., 2006a). Therefore, even though the predictability of the absolute values for biliary and sinusoidal efflux clearance is low, the simulated results will be close to the observed data as far as the uptake clearance is correctly predicted. To validate the predictability of those clearances, the liver concentrations must be determined in humans, which should be possible with imaging technologies such as positron emission tomography, single-photon emission computed tomography, and magnetic resonance imaging. Ghibellini et al. (2007) recently developed a methodology for the real-time measurement of the biliary excretion profiles of drugs in humans using a gamma scintigraphy technique. Further efforts are required to use such in vivo imaging technologies to increase the predictability of these pharmacokinetic parameters.
To date, clinical studies have demonstrated that the genetic variations of OATP1B1 and drug-drug interactions involving OATP1B1 are associated with interindividual differences in the systemic exposure of pravastatin and other substrate drugs (Nishizato et al., 2003; Maeda et al., 2006; Niemi et al., 2006; Shitara and Sugiyama, 2006b). Because the pharmacological target of pravastatin is inside the cell, the liver exposure is a critical factor for its pharmacological activities. Based on the pharmacokinetic concepts, the AUC in the liver concentration is governed only by the sequestration clearance from the liver as far as the renal elimination is negligible and is independent of the change in uptake clearance (eq. A4) (see Appendix II). When renal clearance makes a significant contribution, the changes in hepatic uptake activity can affect both the liver and the plasma AUC (Fig. 9; Appendix II). Actually, the renal elimination of pravastatin makes a significant contribution to the total body clearance (47% of the total body clearance) (Singhvi et al., 1990). Therefore, it is possible that the liver concentrations of pravastatin are affected to some extent also by the changes in the hepatic uptake activity. To support this concept, a simulation was performed with different uptake clearances (Fig. 7). Changes in the hepatic uptake clearance had a great impact on the plasma concentrations of pravastatin but less impact on the liver concentrations. In accordance, the effects of the genetic polymorphisms of OATP1B1 on the cholesterol-lowering effects of pravastatin will be small or absent at least at steady state (in other words, after relatively long-term treatment). The alteration of pharmacological effect of pravastatin with its chronic administration has not been observed in subjects with OATP1B1 polymorphisms although alteration of inhibitory effect of HMG CoA reductase activities in short-term treatments was reported (Takane et al., 2006; Kivistö and Niemi, 2007; Zhang et al., 2007). In contrast, changes in the intrinsic canalicular efflux activity should dramatically affect the liver concentration of pravastatin, whereas the plasma concentration is not affected as much by changes in the intrinsic biliary clearance (Fig. 7). Because the biliary excretion of pravastatin is mainly mediated by MRP2, the factors affecting MRP2 function, such as the use of MRP2 inhibitors or the genetic mutations causing Dubin-Johnson syndrome, will affect the pharmacological action of pravastatin. Furthermore, changes in the sinusoidal efflux clearance had only a slight impact on both the plasma and the liver concentrations. This is because, even under these conditions, the uptake process is still the rate-limiting process. Although the predictability of the sinusoidal efflux clearance remains unknown, changes within this range will not affect the simulated results.
One of the serious adverse effects of statins is myopathy (rhabdomyolysis). Because its target organ is the skeletal muscle, the systemic exposure should be the determinant factor of this adverse effect. The sensitivity analyses showed that the changes in the hepatic uptake clearance had a great impact on the systemic exposure of pravastatin, whereas those in the canalicular efflux had a minimal impact (Fig. 7). The results suggest that patients with an impaired OATP1B1 might be more susceptible to pravastatin-induced myopathy than those with normal one. Morimoto et al. (2004) reported that the frequency of the OATP1B1*15 haplotype was significantly higher in patients who experienced myopathy after receiving pravastatin or atorvastatin (which is also an OATP1B1 substrate) than in patients without myopathy, and a genome-wide study elucidated that the variants in OATP1B1 are strongly associated with an increased risk of simvastatin-induced myopathy (Link et al., 2008).
In the present study, a PBPK model, including transporter-mediated membrane transport processes, was constructed, which allows the prediction of the pharmacokinetics of pravastatin in humans. It also extends our understanding of the effects of changes in the transport processes on the pharmacological and adverse effects of drugs by simulating the exposure of the systemic circulation and tissues to them. The present study suggests that changes in the OATP1B1 activities may have a small and a large impact on the therapeutic efficacy and side effect (myopathy) of pravastatin, respectively, whereas those in the MRP2 activities may have opposite impacts (i.e., a large and a small impact on the therapeutic efficacy and side effect).
Acknowledgments
We thank Ayako Takada and Yuji Sekiya for excellent technical assistance.
Appendix I: Differential Mass Balance Equations for the PBPK Model
Nomenclature
General.
Q, blood flow rate; V, tissue weight; C, pravastatin concentration; Kp, tissue/blood partition coefficient; fB, blood unbound fraction; fT, tissue unbound fraction; Vm, maximum transport velocity; Km, half-saturation concentration for transport; CLR, renal clearance; PSinf, intrinsic hepatic uptake clearance; PSdif, passive diffusion clearance on the sinusoidal membrane; PSbile, intrinsic biliary clearance; CLmet, intrinsic metabolic clearance; Fa, fraction absorbed; ka, absorption rate constant.
Subscripts.
B, blood; LU, lung; BR, brain; MU, muscle; R, kidney; GI, gastrointestinal tract; H, liver; HE, liver extracellular space; inf, influx; met, metabolism.
Model Equations
Hepatic uptake, biliary excretion, and metabolic clearances in humans were linear parameters.
Blood pool:
Lung:
Brain, muscle, kidney:
Liver 1 to 5:
Liver extracellular compartment 1:
Liver extracellular compartments 2 to 5:
Bile or gastrointestinal tract:
Appendix II: Effect of Renal Clearance on the Impact of the Change in the Uptake Clearance on the AUC of the Plasma and Liver
Qh and CLr represent the hepatic blood flow and renal clearance, respectively. PSinf, PSeff, and CLint are hepatic uptake, sinusoidal efflux, and intrinsic sequestration clearances, respectively. V and C represent volume and concentration, respectively. Subscripts b, i, and h represent blood, inlet, and liver, respectively. Mass balance differential equations for each compartment in the simple model shown in Fig. 8 are as follows:
Integrating these equations gives:
where AUCb and AUCh represent the area under the concentration-time curve for the blood and liver, respectively. Substituting CLR = 0 yields:
Simple model to analyze the effects of changes in hepatic uptake activity and intrinsic clearance on blood and liver concentrations. Indicated are the hepatic blood flow (Qh), renal clearance (CLr), volume (V), concentration (C), unbound fraction (f), hepatic uptake clearance (PSinf), sinusoidal efflux clearance (PSeff), intrinsic clearance (CLint), blood (b), inlet (i), and liver (h).
Equations A3 and A4 indicate that AUCh depends only on CLint when the renal clearance is negligible. In contrast, AUCb is inversely proportional to PSinf. If the renal clearance is maximal, that is, the renal blood flow (Qr), eq. A2 can be converted to:
where
When the hepatic uptake is the rate-limiting process, so CLint ≫ PSeff, eq. A5 can be converted to:
In accordance, the R value can be higher than CLint by up to Q × CLint/(fb × PSinf). Figure 9 shows the effects of renal clearance on the impact of the changes in hepatic uptake and intrinsic sequestration clearance on the AUC of pravastatin in the plasma and liver. A simulation was performed using eqs. A1 and A2 and the parameters shown in Tables 1 and 2, when the renal clearance was 0, 4 (one quarter of the renal blood flow), and 16 (the renal blood flow) ml/min/kg b.wt.
Effects of renal clearance on the impact of changes in hepatic uptake and intrinsic sequestration clearances on the AUC of pravastatin in the blood and liver. Relative AUC (100% as the initial value) were estimated by varying the uptake and intrinsic sequestration clearances over a 0.1- to 10-fold range of the initial value when the renal clearance was 0 (——), 4 (– – –), and 16 ( · · · · · · ) ml/min/kg.
Footnotes
Article, publication date, and citation information can be found at http://jpet.aspetjournals.org.
doi;10.1124/jpet.108.146647.
-
ABBREVIATIONS:
- CMV
- canalicular membrane vesicle
- PBPK
- physiologically based pharmacokinetic
- OATP
- organic anion-transporting polypeptide
- MRP
- multidrug resistance-associated protein
- SF
- scaling factor
- R-122798
- (3R,5R)-3,5-dihydroxy-7-[(1S,2S,6S,8S,8aR)-6-hydroxy-8-(isobutyryloxy)-2-methyl-1,2,6,7,8,8a-hexahydronaphthalen-1-yl]heptanoic acid
- LC/MS
- liquid chromatography/mass spectrometry
- PS
- permeability surface product
- inf
- influx
- dif
- diffusion
- CL
- clearance
- met
- metabolism
- tot
- total
- B
- blood
- AUC
- area(s) under the concentration-time curve.
- Received September 25, 2008.
- Accepted November 7, 2008.
- Copyright © 2009 by The American Society for Pharmacology and Experimental Therapeutics