Skip to main content
Advertisement

Main menu

  • Home
  • Articles
    • Current Issue
    • Fast Forward
    • Latest Articles
    • Archive
  • Information
    • Instructions to Authors
    • Submit a Manuscript
    • FAQs
    • For Subscribers
    • Terms & Conditions of Use
    • Permissions
  • Editorial Board
  • Alerts
    • Alerts
    • RSS Feeds
  • Virtual Issues
  • Feedback
  • Other Publications
    • Drug Metabolism and Disposition
    • Journal of Pharmacology and Experimental Therapeutics
    • Molecular Pharmacology
    • Pharmacological Reviews
    • Pharmacology Research & Perspectives
    • ASPET

User menu

  • My alerts
  • Log in
  • My Cart

Search

  • Advanced search
Drug Metabolism & Disposition
  • Other Publications
    • Drug Metabolism and Disposition
    • Journal of Pharmacology and Experimental Therapeutics
    • Molecular Pharmacology
    • Pharmacological Reviews
    • Pharmacology Research & Perspectives
    • ASPET
  • My alerts
  • Log in
  • My Cart
Drug Metabolism & Disposition

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Fast Forward
    • Latest Articles
    • Archive
  • Information
    • Instructions to Authors
    • Submit a Manuscript
    • FAQs
    • For Subscribers
    • Terms & Conditions of Use
    • Permissions
  • Editorial Board
  • Alerts
    • Alerts
    • RSS Feeds
  • Virtual Issues
  • Feedback
  • Visit dmd on Facebook
  • Follow dmd on Twitter
  • Follow ASPET on LinkedIn
Rapid CommunicationAccelerated Communication

PBPK Modeling to Unravel Nonlinear Pharmacokinetics of Verapamil to Estimate the Fractional Clearance for Verapamil N-Demethylation in the Recirculating Rat Liver Preparation

Qi Joy Yang, Luqin Si, Hui Tang, Helle H. Sveigaard, Edwin C. Y. Chow and K. Sandy Pang
Drug Metabolism and Disposition April 2015, 43 (4) 631-645; DOI: https://doi.org/10.1124/dmd.114.062265
Qi Joy Yang
Leslie Dan Faculty of Pharmacy, University of Toronto, Toronto, Ontario, Canada (Q.J.Y., L.S., H.T., H.H.S., E.C.Y.C., K.S.P.); School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, P.R. China (L.S.); School of Pharmacy, Shihezi University, Shihezi, Xinjiang, P.R. China (H.T.); and Denmark Pharmaceutical University, Copenhagen, Denmark (H.H.S.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Luqin Si
Leslie Dan Faculty of Pharmacy, University of Toronto, Toronto, Ontario, Canada (Q.J.Y., L.S., H.T., H.H.S., E.C.Y.C., K.S.P.); School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, P.R. China (L.S.); School of Pharmacy, Shihezi University, Shihezi, Xinjiang, P.R. China (H.T.); and Denmark Pharmaceutical University, Copenhagen, Denmark (H.H.S.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hui Tang
Leslie Dan Faculty of Pharmacy, University of Toronto, Toronto, Ontario, Canada (Q.J.Y., L.S., H.T., H.H.S., E.C.Y.C., K.S.P.); School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, P.R. China (L.S.); School of Pharmacy, Shihezi University, Shihezi, Xinjiang, P.R. China (H.T.); and Denmark Pharmaceutical University, Copenhagen, Denmark (H.H.S.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Helle H. Sveigaard
Leslie Dan Faculty of Pharmacy, University of Toronto, Toronto, Ontario, Canada (Q.J.Y., L.S., H.T., H.H.S., E.C.Y.C., K.S.P.); School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, P.R. China (L.S.); School of Pharmacy, Shihezi University, Shihezi, Xinjiang, P.R. China (H.T.); and Denmark Pharmaceutical University, Copenhagen, Denmark (H.H.S.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Edwin C. Y. Chow
Leslie Dan Faculty of Pharmacy, University of Toronto, Toronto, Ontario, Canada (Q.J.Y., L.S., H.T., H.H.S., E.C.Y.C., K.S.P.); School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, P.R. China (L.S.); School of Pharmacy, Shihezi University, Shihezi, Xinjiang, P.R. China (H.T.); and Denmark Pharmaceutical University, Copenhagen, Denmark (H.H.S.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
K. Sandy Pang
Leslie Dan Faculty of Pharmacy, University of Toronto, Toronto, Ontario, Canada (Q.J.Y., L.S., H.T., H.H.S., E.C.Y.C., K.S.P.); School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, P.R. China (L.S.); School of Pharmacy, Shihezi University, Shihezi, Xinjiang, P.R. China (H.T.); and Denmark Pharmaceutical University, Copenhagen, Denmark (H.H.S.)
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF + SI
  • PDF
Loading

This article has a correction. Please see:

  • Correction to: “PBPK Modeling to Unravel Nonlinear Pharmacokinetics of Verapamil to Estimate the Fractional Clearance for Verapamil N-Demethylation in the Recirculating Rat Liver Preparation” - July 01, 2015

Abstract

We applied physiologically based pharmacokinetic (PBPK) modeling to study the dose-dependent metabolism and excretion of verapamil and its preformed metabolite, norverapamil, to unravel the kinetics of norverapamil formation via N-demethylation. Various initial verapamil (1, 50, and 100 μM) and preformed norverapamil (1.5 and 5 μM) concentrations, perfused at 12 ml/min, were investigated in the perfused rat liver preparation. Perfusate and bile were collected over 90 minutes, and livers were harvested at the end of perfusion for high-performance liquid chromatography analysis. After correction for the adsorption of 10%–25% dose verapamil and norverapamil onto Tygon tubing and binding to albumin and red blood cell, fitting of verapamil and formed and preformed norverapamil data with ADAPT5 revealed nonlinearity for protein binding, N-demethylation ( Embedded Image nmol/min; Embedded Image μM), formation of other metabolites ( Embedded Image nmol/min; Embedded Image μM), as well as biliary excretion ( Embedded Image nmol/min; Embedded Image μM). The hepatic clearance of verapamil ( Embedded Image ) decreased with the dose (8.16–10.2 ml/min), with values remaining high relative to perfusate blood flow rate among the doses. The hepatic clearance of preformed norverapamil (11 ml/min) remained unchanged for the concentrations studied and approximated perfusate blood flow rate, suggesting a high norverapamil extraction ratio. The fractional formation of norverapamil and biliary excretion of verapamil based on fitted constants were 31.1% and 0.64% of Embedded Image , respectively. Enantiomeric disposition and auto-inhibition of verapamil failed to perturb these estimaties according to PBPK modeling, due to the low values of the Michaelis-Menten constant, K m, and inhibition parameter, k I.

Introduction

Verapamil, a calcium channel blocker (Fleckenstein, 1977) used for the treatment of cardiac arrhythmias and hypertension (McTavish and Sorkin, 1989), is given as a racemic mixture. Much is known about the stereoselective disposition of the more potent S- versus R-verapamil in human and rats (Bhatti and Foster, 1997; Busse et al., 2006). Binding of verapamil to bovine serum albumin (BSA) revealed that the unbound plasma fraction of R-verapamil (∼0.65) exceeds that of S-verapamil (∼0.55) (Mehvar and Reynolds, 1996). The reverse is true for binding to human serum albumin and α 1-acid glycoprotein (Mehvar and Reynolds, 1996; Hanada et al., 1998), whereas the plasma bound fraction of verapamil in human is concentration independent within the therapeutic range and unaltered by the presence of norverapamil (Keefe et al., 1981). The distribution of verapamil and norverapamil into red blood cells (RBCs) is also stereoselective [human (S > R) versus rat (R > S)] (Robinson and Mehvar, 1996).

Much species difference exists for the metabolism of verapamil. First-pass removal of verapamil is attributed to both the intestine and liver in rats (Hanada et al., 2008) and man (Fromm et al., 1998), although intestinal elimination is lacking in dogs (Lee et al., 2001). At least 25 Phase I and 14 Phase II metabolites have been identified in the rat (Walles et al., 2003). First-pass metabolism in rats is saturable and stereoselective: hepatic bioavailability for R-verapamil is higher than that for S-verapamil, but the opposite was observed for intestinal bioavailability, resulting in a higher systemic bioavailability for R-verapamil following oral dosing (Hanada et al., 2008). Extensive first-pass removal and liver clearance (CL) have been reported, with verapamil being mostly metabolized via N- and O-dealkylation (98%) in both rats and humans (Eichelbaum et al., 1979; Woodcock et al., 1981; Choi and Burm, 2008), with norverapamil and D617 as major metabolites formed via N-dealkylation, predominantly by CYP3A4, CYP3A5, CYP2C8, and CYP2E1 (Eichelbaum et al., 1979; Kroemer et al., 1992; Tracy et al., 1999; Sun et al., 2004). Nonlinear metabolism was observed in man following long-term, repeated oral dosing (1–1.6 mg/kg over 1 month), leading to prolongation of half-life and decreased apparent clearance (CL/F, or Dosepo/AUCpo, where AUC denotes area under the curve) (Freedman et al., 1981; Shand et al., 1981). High doses of verapamil (6–19 mg/kg) infused intravenously in man resulted in saturable kinetics and decreased systemic clearance with dose escalation (Toffoli et al., 1997). In addition, verapamil and norverapamil were reported to be mechanism-based inhibitors of cytochrome P450 and P-gp, which resulted in auto-inhibition upon long-term use clinically (Lemma et al., 2006; Wang et al., 2013). Verapamil and norverapamil were present at significantly higher concentrations in human intestine lumen and bile compared with plasma, suggesting excretion mediated by transporters, probably P-gp (von Richter et al., 2001). Norverapamil is transported by P-gp in Caco-2 cells and P-gp-overexpressing L-MDR1 cells, but to a lesser extent than other verapamil metabolites, D-617 and D-620 (Pauli-Magnus et al., 2000).

The need for understanding metabolite formation kinetics is paramount, especially when the metabolite in question is active or toxic (Baillie et al., 2002; Pang, 2009). The purpose of this investigation is to develop a strategy to study dose-dependent verapamil metabolism and biliary excretion and describe metabolite formation kinetics on the importance of a given metabolic pathway (N-demethylation) relative to other elimination pathways (metabolism or biliary excretion) with use of verapamil and its active metabolite, norverapamil. To this end, Mehvar et al. (1994) estimated the extent of formation of norverapamil from verapamil in the perfused rat liver preparation and examined removal of verapamil from plasma. Since the K m value for metabolism ranges from 60 to 140 μM (Hanada et al., 2008), we employed initial concentrations of 1–100 μM verapamil to revisit the problem of metabolite kinetics in the recirculating perfused rat liver preparation. Mechanism-based auto-inhibition for verapamil that exists in man (Wang et al., 2004, 2013) was not observed in rat liver microsomes, hepatocytes, and precision-cut liver slices (Obach, 1999; Shibata et al., 2002; Axelsson et al., 2003; Guo et al., 2007), except with gel entrapment after prolonged exposure with verapamil (Yin et al., 2011). Hence, auto-inhibition is not expected to occur within the short time frame for liver perfusion studies. We applied physiologically based pharmacokinetic (PBPK) modeling, to account for tubing adsorption of verapamil and norverapamil and stereoselective vascular binding and saturable metabolism and excretion.

Materials and Methods

Materials

Verapamil (racemic, R-, and S-; C27H38N2O4) and norverapamil (C26H36N2O4) were purchased from Sigma-Aldrich (Mississauga, Canada). BSA was obtained from Sigma-Aldrich; dextrose (50%) was purchased from Abbott Laboratories (Montreal, Canada). High-performance liquid chromatography (HPLC) grade acetonitrile, methanol, and ethyl acetate were obtained from Sigma-Aldrich. Male Sprague-Dawley rats (324 ± 31.1 g) were supplied by Charles River Laboratories (St. Constant, Canada). [3H]Verapamil, produced via general tritium exchange (specific activity 1 mCi/ml; >97% radiochemical purity by HPLC), was obtained from American Radiolabeled Chemicals Inc. (St. Louis, MO).

Albumin Binding and RBC Partitioning of Verapamil

A mixture of labeled and unlabeled verapamil was used to provide for sensitivity in the binding studies. The total verapamil concentration was given as the sum of the labeled (from the specific activity) and unlabeled verapamil, and the dpm/ml in plasma was the equivalent of the assayed HPLC concentration in plasma. Binding of [3H]verapamil to BSA (2%) was measured in equilibrium dialysis half-cells that were separated by a semipermeable membrane (molecular cutoff 12,000–14,000; Spectrum Laboratories Inc., Rancho Dominguez, CA). Next, 1 ml 2% BSA in Krebs-Henseleit buffer (KHB) (pH 7.4) containing various concentrations of verapamil (0.4–140 μM) was added to the protein-side half-cell, and then 1 ml of KHB was introduced to the other (buffer side) half-cell. The cell was incubated at 37°C using rotating water bath. After 5 hours of equilibration, the time predetermined for equilibration, samples (100 μl) from both sides were removed for protein determination and liquid scintillation counting (using an LS 5801 Counter; Beckman Coulter Canada, Mississauga, Canada). Preservation of protein and water volumes was checked at the end of the experiment. The dpm values in the protein and buffer sides were counted. The plasma unbound fraction of verapamil (f P) was calculated as the ratio of verapamil concentrations in the buffer side to that in the protein side, at the plasma concentration determined at the end of the binding study.

The distribution of verapamil into RBCs was studied by mixing [3H]verapamil and unlabeled (1–440 μM) verapamil in 2% BSA plasma perfusate. Upon mixing this protein solution with an equal volume of blood perfusate containing 40% washed bovine erythrocytes, RBC (or 2× the normal 20% RBC, a kind gift from Ryding-Regency Meat Packers Ltd., Toronto), a blood perfusate of composition identical to that used for perfusion was obtained. Aliquots of the resultant labeled and unlabeled verapamil concentration in blood perfusate (C B) (one-half of the concentration of the initial plasma perfusate), were removed for assay by HPLC, and for incubation at 37°C for 5 minutes. Then, 200 μl blood perfusate was removed at 0.25, 0.5, 0.75, 1, 1.5, 2, 3, 4, and 5 minutes after incubation, and hematocrit (Hct) was measured at 0.25 and 5 minutes. The total verapamil concentration in plasma (C P) was estimated using liquid scintillation counting; the concentration in RBCs (C RBC) was calculated as Embedded Image . Data of verapamil binding to BSA and distribution to RBC were then fitted simultaneously according to the binding model with one class of binding sites (see Appendix A for the mass balance equations) to furnish constants such as the binding association constant (KA ) and modified rate constants for partitioning into ( Embedded Image )and out of ( Embedded Image ) RBCs based on the total drug concentration using ADAPT, version 5 (Biomedical Simulations Resource, University of Southern California, Los Angeles, CA). The fraction unbound in plasma (f P) was predicted according to best-fitted binding parameters.

Recirculating Rat Liver Perfusion and Drug Adsorption

The perfusion apparatus, TWO/TEN Perfuser (MX International, Aurora, CO), was used for recirculation of verapamil and norverapamil at the flow rate of 12 ml/min and 37°C. Perfusate consisted of 1% BSA, 0.6% dextrose, and 20% washed bovine erythrocyte in oxygenated (95% O2 and 5% CO2 at 1 l/min) KHB, pH 7.4. The surgical and liver perfusion procedures were identical to those described previously (Tan and Pang, 2001). Male Sprague-Dawley rats (300–365 g) were anesthetized with a mixture of ketamine (90 mg/kg) and xylazine (10 mg/kg), and recirculating liver perfusion was conducted after cannulation of the bile duct and portal and hepatic veins for the inflow and outflow, respectively; the hepatic artery was ligated. Drug-free blank perfusate was first used to recirculate the liver for 20 minutes for equilibration, followed by perfusion with drug-containing perfusate from a second reservoir (200 ml) consisting of the designated concentration of verapamil (1, 50, or 100 μM, n = 4) or preformed norverapamil (1.5 μM, n = 4, or 5 μM, n = 3) for the next 90 minutes. Reservoir perfusate (1 ml) was removed at 2.5, 7.5, 12.5, 17.5, 22.5, 27.5, 35, 45, 55, and 70 minutes, and 2 ml was removed at 0 and 90 minutes. Bile was sampled at 5 minute intervals between 0 and 30 minutes, or at 10 or 20 minute intervals thereafter. At the end of perfusion (90 minutes), the liver was flushed with 50 ml cold blank KHB to remove any residual blood perfusate, and the liver was weighed, minced, blast frozen with liquid nitrogen, and stored at −80°C until analysis by HPLC.

Since adsorption of verapamil and norverapamil to tubing was found, Tygon tubing of constant length (obtained from Saint-Gobain Performance Plastics, Valley Forge, PA), which exhibited less binding to the compounds, was used for perfusion. For characterization of binding, liver perfusion with drug-free (20 minutes) and then verapamil- or norverapamil-containing (90 minutes) blood perfusate was conducted in the absence of the rat liver to characterize for the loss of drug due to binding to tubing.

HPLC

For Measuring CB and Distribution of Verapamil into RBC.

A HPLC method was used for the separation of verapamil and its internal standard (50 μl of 100 μg/ml diltiazem) at the wavelength of 230 nm (Garcia et al., 1997). The HPLC system consisted of a Shimadzu 6A UV spectrophotometric detector, LC-6AD liquid chromatograph, SIL-9A autoinjector, CR-4A chromatopac (Mandel Scientific Company, Guelph, ON), and a Waters 15–20 μm μBondapak C18 reverse column (3.9 × 300 mm, Waters Limited, Mississauga, ON, Canada). The mobile phase [40% 0.04 M ammonium acetate: triethylamine (2:0.04 v/v) + 60% methanol] was delivered at a flow rate of 0.9 ml/min. The typical retention times for verapamil and diltiazem were 18 and 23 minutes, respectively. The extraction procedure of verapamil for blood perfusate samples was as follows: to 1 ml sample in a 12 ml polypropylene tube, 50 μl of 100 μg/ml diltiazem, and 75 μl of 1 M NaOH were added and mixed thoroughly, and then 4 ml of a hexane and 2-propanol mixture (2:1 v/v) was added. The content was mixed for 2 minutes, followed by centrifugation at 3000 rpm. After removal and drying of the extract under N2, the residue was reconstituted in 200 μl acetontrile and water (1:1 v/v) for HPLC injection.

Verapamil and Norverapamil Determination in Tubing Adsorption and Liver Perfusion Studies.

In these studies, the Shimadzu 10A HPLC system was used to separate verapamil and norverapamil at the detection wavelength of 278 nm. The system consisted of a Shimadzu SCL-10A system controller, a SPD-10A UV-visible detector, a SIL-10XL automatic injector, a FCV-10AL solvent delivery unit and a LD-10AT liquid chromatograph. A 10 μm C18 reverse column (4.6 × 250 mm, Altech Associates, Deerfield, IL) was used to separate verapamil, norverapamil, and the internal standard (diazepam). The mobile phase consisted of 55% of 0.02 M ammonium acetate buffer (0.4% triethylamine, pH adjusted to 6 with acetic acid) and 45% of acetonitrile, and was maintained for 30 minutes for each injection at a flow rate of 1 ml/min. The calibration curves for verapamil and norverapamil were linear over the ranges of 0.25–60, and 0.125–30 nmol, respectively, in 2 ml perfusate, with correlation coefficients of ∼0.999 (n = 5). The limits of quantitation were 0.25 nmol for verapamil and 0.125 nmol for norverapamil. The typical retention times for verapamil, norverapamil, and diazepam were 15, 12, and 23 minutes, respectively.

Blood perfusate samples (0.25–1.8 ml) were made up to a final volume of 2 ml with blank perfusate, mixed with 50 μl internal standard (50 μg/ml diazepam) and 50 μl of 1 M NaOH, and extracted against 6 ml ethyl acetate. After repeated, 3 minute vigorous mixing, the mixture was centrifuged at 3000 rpm for 10 minutes and the supernatant was transferred and dried under nitrogen gas. The dried residue that was reconstituted in the 100 μl mobile phase was centrifuged, and 50 μl was injected. Liver tissue was obtained from homogenization in 1.5× volume of ice-cold saline for 30 seconds 3×. Then, 0.2 ml was removed, and 1 ml ethyl acetate was added and mixed for 3 minutes, followed by centrifugation at 3000 rpm for 10 minutes. The organic layer was removed and dried under nitrogen gas and the residue was reconstituted with the 200 μl mobile phase, and 150 μl was injected into the HPLC. Calibration curves were prepared based on a set of standards containing varying known concentrations of verapamil and norverapamil in perfusate blood or tissue, prepared under identical conditions as the samples.

Bile samples were made up to 60 μl with H2O, 10 μl acetonitrile and 10 μl internal standard (50 μg/ml diazepam) were added, and then 1 ml ethyl acetate was extracted. The nitrogen-dried supernatant was reconstituted in the 200 μl mobile phase and 90 μl was injected into the HPLC. For the norverapamil perfusion studies, bile collected during the entire 90 minutes of perfusion was pooled, mixed with 50 μl acetonitrile, 50 μl of 50 μg/ml diazepam, and 50 μl NaOH (1 M), and then extracted with 3 ml ethyl acetate. After the supernatant was nitrogen dried, the residue was reconstituted in the 200 μl mobile phase, and 50 μl was injected into the HPLC. The concentrations of verapamil and norverapamil were determined with the use of calibration curves, prepared under identical conditions.

PBPK Modeling

We applied a PBPK model to fit the perfusate, liver, and bile data of verapamil and preformed and formed norverapamil simultaneously (Fig. 1). The model was constructed under a number of assumptions and using the mass balance equations given in the Appendices A–C. The tubing adsorption rate constants, estimated for reservoir for verapamil ( Embedded Image or Embedded Image ) and norverapamil ( Embedded Image or Embedded Image ) via fitting of the perfusion data in absence of liver (sham experiments), were viewed to be identical to other perfusion experiments with liver. These tubing adsorption rate constants, along with other physiologic parameters such as liver volume (V L) and hepatic blood perfusate flow (Q L), were assigned for modeling of the liver perfusion data. The recirculating liver preparation was constructed as two compartments, liver blood (LB) and liver tissue (L), which were connected to the reservoir (R) by the hepatic perfusate blood flow (Q L). A flow-limited distribution of verapamil and preformed and formed norverapamil was assumed because these drugs are lipophilic compounds ( Embedded Image and 3.3, respectively), and therefore can diffuse freely across the membrane. The influx clearance ( Embedded Image or Embedded Image ) and efflux clearance ( Embedded Image or Embedded Image ) were assumed to be the same, and were assigned as 5× the value of perfusate blood flow (Q L) as initial estimates for the influx/efflux clearances in the PBPK model. Since only unbound verapamil and norverapamil were able to be transported or eliminated, the unbound fractions of drug in blood perfusate, f B (measured experimentally and incorporated in the model, see Appendices A–C), and liver (f L; fitted in the model) were considered. Verapamil was metabolized to norverapamil or other metabolites, with the metabolic intrinsic clearances, Embedded Image and Embedded Image [V max/K m and V max/(C L,u + K m) for linear and nonlinear conditions, respectively, where C L,u denotes the unbound liver concentration]. The assumption that stereoselective metabolism of verapamil under our high concentrations being unimportant was made. Norverapamil was further metabolized with the intrinsic metabolic clearance or Embedded Image (or Embedded Image ). Verapamil and norverapamil were secreted into bile with secretory intrinsic clearances, Embedded Image (or Embedded Image ) or Embedded Image (or Embedded Image ), respectively. Transit compartments (A tr) containing time delays (represented as the τ function) were incorporated to describe the delayed biliary secretion of verapamil and norverapamil (see Appendix C for the mass balance equations).

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

PBPK model of verapamil (racemic) and norverapamil in the recirculating perfused rat liver. Perfusate blood flow rate (Q L) interconnects the reservoir with liver blood compartment. The compartments are denoted as the reservoir (R), Tygon tubing (t), liver blood (LB), liver tissue (L), and bile. Enzymes (Cyps) and transporters (P-gp) are considered in the model. Both verapamil and norverapamil are adsorbed onto the Tygon tubing (top), and concentrations of unbound verapamil and norverapamil in reservoir blood (denoted by Embedded Image VERLB and Embedded Image NORLB, respectively) reach liver blood and rapidly diffuse into and out of the liver (via CLin and CLef, respectively, the transport clearances that equal each other). Unbound verapamil in liver (of concentration Embedded Image VERL) is metabolized to norverapamil via N-demethylation, with the intrinsic clearance, Embedded Image . VER is also eliminated in liver to form other metabolites or is excreted into bile, with intrinsic clearances of Embedded Image and Embedded Image , respectively. Similarly, unbound norverapamil (of concentration Embedded Image NORL), the formed metabolite, is further metabolized ( Embedded Image ) and secreted ( Embedded Image ). Since all these metabolic and secretory pathways are saturable, the intrinsic clearance parameters may be expressed in terms of V max and K m as well as the unbound liver concentration (see Appendix C for the mass balance equations). Two transit compartments (of amounts Embedded Image and Embedded Image ; Embedded Image and Embedded Image ) were defined in the model to account for the time delay (τ) in biliary secretion.

Our strategy was to first consider linearity in the removal processes for verapamil metabolism and excretion, and then consideration of each pathway as nonlinear. Michaelis-Menten equations, with the appropriate maximum velocity of the enzyme or transporter, V max, and Michaelis-Menten constant, K m, terms were used in lieu of CLint for metabolism and biliary secretion: the metabolic intrinsic clearance of the liver, CLint,met, or the excretory intrinsic clearance of the liver, CLint,sec, is equal to Embedded Image , namely, the intrinsic clearance terms under nonlinear conditions were given as V max divided by the sum of the unbound substrate concentration in liver (C L,u) and the K m term. In parallel, data for preformed norverapamil metabolism and excretion were fitted, with and without nonlinearity. The final fit consisted of simultaneously fitting of both verapamil (1, 50, and 100 μM) and formed and preformed norverapamil (1.5 and 5 μM) data.

The model fitting criterion was evaluated based on the Akaike information criteria (AIC) score, or Embedded Image , where k is the number of parameters in the model and likelihood is the maximum value of the likelihood function for the model. The AIC score evaluates the goodness-of-fit with consideration of the penalty for adding additional parameters to fit the data set; the lower the AIC score, the better is the goodness-of-fit. Also, the correlation between predicted and observed concentrations and weighted residuals, calculated as (observations − predictions)/observations versus time plots were used to examine systematic trends and extreme values. The sum of squared weighted residuals (SSWR) was calculated for each of the candidate models to define which model was most appropriate to describe the pharmacokinetics of verapamil and formed and preformed norverapamil.

Statistical Analysis

The estimated parameters were shown as the mean ± S.D. (CV%). The CV% was calculated as S.D./mean × 100%. The two-tailed student t test with unequal variance was used to compare the means between two groups, and a P value of less than 0.05 was considered as significant. One-way analysis of variance and the post-hoc Tukey honest significance difference test were used to evaluate the means in three dosing groups (1, 50, and 100 μM) using SPSS, version 22 (IBM Corp., New York). The F test was used to assess model improvement (Boxenbaum et al., 1974). The error proportional method was used to calculate uncertainties: Embedded Image was used when adding or subtracting two or more values with the S.D.; Embedded Image was used when multiplying or dividing two or more values with the S.D.

Results

Protein Binding and RBC Partitioning.

Binding of verapamil (0.4–140 μM) to BSA (2%) in plasma was found to be concentration dependent, with the unbound fraction in plasma (f P) varying from 0.4 to 0.6 within the concentration range studied (Fig. 2A). The ratio of verapamil concentration in perfusate blood to plasma (C B/C P) was found to be concentration dependent, with the value increasing from 1.11 ± 0.07 to 1.29 ± 0.05 (P < 0.001) within the concentration range of 1–441 μM examined. An average unbound fraction of verapamil in blood perfusate (0.45 ± 0.09) may be estimated, based on the average value of f P (f B = C p f p/C B) (Pang and Rowland, 1977). The ratio of verapamil concentrations in RBCs and plasma (C RBC/C P) was around 2–3 within the concentration range studied when the Hct was 0.134 ± 0.014, showing that verapamil can distribute into RBCs. Upon fitting the data to the model (Appendix A), the nonlinear protein binding and distribution into RBC were adequately described by the model with one class of binding sites (Fig. 2B), yielding fitted values of modified rate constants Embedded Image (13.2 ± 4.4 min−1) and Embedded Image (19.5 ± 8.7 min−1), the number of protein binding sites (n = 0.301 ± 0.160), and binding association constant (KA = 0.0091 ± 0.0057 μM−1). Since the turnover time ( Embedded Image ) of verapamil in red cells was 3.1 seconds—a value much lower than the mean transit time of blood in the liver, or 14–15 seconds (Pang et al., 1995)—verapamil in plasma and RBC would rapidly achieve equilibrium and will not affect verapamil removal. Nonlinear protein binding due to one single class of binding sites and concentration-dependent C B/C P were incorporated in the PBPK modeling (Appendices B and C). Protein binding according to two classes of binding sites was also considered for fitting, although the fit was suboptimal (data not shown). Binding would not materially affect the clearance of verapamil since the drug is highly cleared, and the trend was verified, namely, small changes in f B (with values fluctuating by +0.1 or −0.1), would not materially alter subsequent metabolism or excretion (simulations not shown).

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Nonlinear plasma protein binding of verapamil (racemic) to BSA (2%) (A), and distribution of verapamil into RBCs (B). (A) The unbound fraction of verapamil in plasma (f P) was plotted against log(CP) to reveal concentration-dependent binding within the range of 0.4–140 μM. The sold line represents f P values predicted according to one class of binding site (fitted constants, n = 0.301 ± 0.160 and K A = 0.0099 ± 0.0057 μM−1) (Appendix A). (B) Observed (circle) and fitted (line) are CRBC and CP over time (n = 6) in the RBC distribution study (see Materials and Methods for details). After the model was fitted to the data of verapamil in RBC and plasma, the transfer rate constants, Embedded Image and Embedded Image , were found to be 19.5 ± 8.7 and 13.2 ± 4.4 min−1, respectively (Appendix A). The numbers next to the fitted curves are the total initial CP (1–441 μM in 2% BSA) used for admixture with 40% RBC containing blood perfusate (in 2% BSA).

Tubing Adsorption.

The adsorption profiles of verapamil and norverapamil to Tygon tubing are shown in Fig. 3. The rate constants for binding (k R→t) to tubing and debinding (k t→R), obtained from fitting the adsorption data without rat liver, are summarized in Table 1. Although the absolute amounts of drug adsorbed to Tygon tubing increased with dose, the %dose of verapamil and preformed norverapamil adsorbed decreased with increasing concentrations perfused. The lowest adsorbed amount (50 nmol or 25% dose) occurred for the lowest concentration (for the initial concentration of 1 μM), whereas about 10% dose (1 or 2 μmol) was adsorbed for higher doses (for initial concentrations of 50 and 100 μM). Similarly, the amount of norverapamil adsorbed onto the Tygon tubing increased from 0.06 to 0.1 μmol after perfusion with initial concentrations of 1.5 and 5 μM norverapamil. The rate constants for binding to Tygon tubing ( Embedded Image and Embedded Image ) decreased as the dose escalated, suggesting that tubing binding sites might be saturated at higher drug concentrations. There was no change in the return rate constants ( Embedded Image and Embedded Image ). The tubing rate constants (Table 1) were assigned as constants into the PBPK model (see Appendix C for the mass balance equations).

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

The amounts (percentage of dose) of verapamil (A) and preformed norverapamil (B) adsorbed onto Tygon tubing of the perfusion apparatus. The solid symbols denote the measured amounts (percentage of dose), and the solid lines denote the fitted amounts over time (see Appendix C for the mass balance equations). Different colored symbols and lines denote the different initial concentrations of verapamil (racemic) and norverapamil (racemic) used for perfusion.

View this table:
  • View inline
  • View popup
TABLE 1

Binding (k R→t) and debinding (k t→R) constants for verapamil and norverapmil to the tubing of the perfusion apparatus, obtained by fitting the perfusate data in the absence of a rat liver in sham-liver perfusion studies (Fig. 1)

Nonlinear Kinetics of Verapamil and Norverapamil.

The physiologic and experimentally derived estimates for the verapamil (Supplemental Table 1; Table 2) and preformed norverapamil (Supplemental Table 2; Table 3) perfusion studies are summarized. Preformed norverapamil was investigated to provide information on the sequential handling of norverapamil formed from verapamil. The temporal profiles of verapamil and its preformed and formed metabolite norverapamil in perfusate, liver, and biliary secretion are shown in Fig. 4. Volume recovery was ∼85%, showing small experimental error. In reservoir perfusate, the verapamil concentrations decreased rapidly following an apparent log-linear decay, and the terminal elimination rate constant (calculated from the slope of elimination phase, β) was significantly decreased from 0.0536 ± 0.0008 to 0.0386 ± 0.0075 min−1 for the high-dose groups (50 and 100 μM) compared with the low-dose group (1 μM; P = 0.0007), suggesting that the apparent half-life was prolonged with dose, exhibiting concentration-dependent or nonlinear kinetics. Similarly, the hepatic clearance ( Embedded Image ) for verapamil decreased from 10.2 to 8.75 ml/min (P = 0.04), although the values remained high and close to the perfusate blood flow rate; Embedded Image or the extraction ratio of Embedded Image was high (0.8). Concentrations of formed norverapamil in perfusate were very low (< 4 μM), and no apparent trend was observed over time. The dose-normalized Embedded Image and Embedded Image increased with increasing input concentrations, suggesting accumulation of drug and metabolite in perfusate with dose (Supplemental Table 1).

View this table:
  • View inline
  • View popup
TABLE 2

Experimentally derived parameters with different verapamil initial concentrations delivered to the recirculating, RBC-perfused rat liver preparation at 12 ml/min

View this table:
  • View inline
  • View popup
TABLE 3

Experimentally derived parameters with dosing of preformed norverapamil at initial concentrations of approximately 1.5 and 5 μM

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Temporal profiles of fitted and observed verapamil and formed norverapamil in the perfusate (red), liver (green), and the cumulative secreted amounts in bile (blue) at initial concentrations of 1, 50, and 100 μM of racemic verapamil (n = 4) in recirculating rat liver preparations. The lines denote predictions according to the final fit with the PBPK model (Model H) describing various nonlinear pathways ( Embedded Image , Embedded Image , Embedded Image , and Embedded Image ) (see Appendix C for the mass balance equations). The solid and open symbols are the observed concentrations or amounts of verapamil and formed norverapamil, respectively, in liver and bile.

Preformed norverapamil decayed rapidly and monoexponentially in reservoir perfusate (Fig. 5), yielding similar terminal elimination rate constants (0.0484 ± 0.0022 versus 0.0612 ± 0.008 min−1; P = 0.15) for the 1.5 and 5 μM input concentration groups. Biliary excretion of norverapamil was not remarkable. The estimated extraction ratio of preformed norverapamil ( Embedded Image ) was high (0.9). The fraction of hepatic clearance of verapamil forming norverapamil ( Embedded Image ), based on the assumption of linearity for the lowest verapamil concentration (1 μM) and the equation comparing the AUC of formed versus preformed norverapamil, normalized to the availability of norverapamil, Embedded Image , was 23.2% (Pang and Kwan, 1983). Biliary excretion clearance of verapamil was 1% of Embedded Image for the lowest input concentration of 1 μM (Table 2); the biliary secretion of preformed norverapamil was only 0.3%–0.4% of Embedded Image .

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Temporal profiles of observed (open symbols) and fitted preformed norverapamil in the reservoir perfusate (red), liver (green), and cumulative excreted amounts in bile (blue) at initial norverapamil (racemic) concentrations of 1.5 (A) (n = 4) and 5 (B) (n = 3) μM in recirculating rat liver preparations. The solid lines denote the fitted values according to the final PBPK model (Model H) with nonlinear Embedded Image (see Appendix C for the mass balance equations).

PBPK Modeling.

The data for verapamil and formed and preformed norverapamil were fitted simultaneously using ADAPT5 (Biomedical Simulations Resource). The strategy was to first consider all metabolism and excretion as first-order processes; nonlinearity in metabolism or excretion was then included stepwise into the fitting routine. Physiologic parameters such as volume of tissue (V) and perfusate blood flow rates (Q L) were assigned (see Table 5). Temporal profiles of fitted verapamil and formed and preformed norverapamil in perfusate, liver, and bile (Figs. 4 and 5) for the best model included nonlinearity in N-demethylation, formation of other verapamil metabolites, and biliary excretion. Evidence of improvement of fit is summarized in Table 4. According to the F test, incorporation of nonlinearity in Embedded Image , Embedded Image , Embedded Image , and Embedded Image to the base model significantly improved the goodness-of-fit, evidenced by the lowest SSWR and AIC score (Table 4). With the nonlinear pathways, the predicted concentrations matched the measured concentrations well and fell on the line of identity; the predicted versus observed and weighted residuals revealed a random distribution around 0 (Supplemental Figs. 1 and 2).

View this table:
  • View inline
  • View popup
TABLE 4

F tests comparing various PBPK models (A–H) with different combinations of linear and nonlinear elimination pathways

In the final model (Model H) that best fit the data, the CV% was between 5% and 48% of all parameter estimates (Table 5). The unbound fraction in perfusate blood estimated for norverapamil was 0.44, a value similar to that for verapamil, whereas liver tissue binding to verapamil and norverapamil were high. Reasonable transfer clearances ( Embedded Image and Embedded Image ) as well as the V max and K m parameters with low CV% were obtained. The K m values for the secretion of verapamil and norverapamil were around 5 μM, which suggests that the P-gp pump could be readily saturated. The Embedded Image or fraction of hepatic clearance of verapamil that forms norverapamil via N-demethylation according to the fitted constants, Embedded Image , was 0.311 ± 0.189, and 68.3% of Embedded Image was conducive to formation of other verapamil metabolites; biliary secretion of verapamil was only 0.64% Embedded Image . Based on the ratio of the fitted V max and K m, the recalculated Embedded Image (31.1%) was around 10% higher than that based on the lowest dose under the linear condition (or 23.3%), and was significantly higher than that of 12% obtained by Mehvar et al. (1994).

View this table:
  • View inline
  • View popup
TABLE 5

Fitted estimates obtained with fitting the final model (Model H) to liver perfusion data (Table 4), with incorporation of nonlinear protein binding and nonlinear metabolism and biliary secretion of verapamil and norverapamil (see Appendix C for the mass balance equations)

PBPK Modeling on Stereoselective Disposition and Autoinhibition of Verapamil Metabolism.

One limitation of the present study, however, was that stereoselectivity in verapamil and norverapamil disposition has not been fully addressed. In terms of stereoselective binding, it has been reported that R-VER has a higher (∼1.6×) unbound fraction in plasma (f P) over S-VER at both low and high concentrations (Mehvar et al., 1994; Mehvar and Reynolds, 1995; Robinson and Mehvar, 1996), but the C B/C P ratio for R-VER is ∼1.5 times that of S-VER (Hanada et al., 2008). Since the unbound fraction in blood f B = f P × C P/C B, we conclude that the unbound fractions in blood for S-VER and R-VER are similar, as are the blood profiles of (R)- and (S)-verapamil. Therefore, it is reasonable to assume that the f B value for the racemic mixture was the same as those for S-VER and R-VER. We further found that small changes in f B would not materially affect the clearances of S-VER, R-VER, or racemic VER (simulations not shown). This is reasonable since the compounds are all highly cleared.

We had assumed that stereoselectivity in metabolism was not an important factor under our dosing input concentrations, which were high in relation to the K m values (initial concentrations of 50 and 100 μM), and stereoselectivity effects in verapamil metabolism could be negligible (Mehvar and Reynolds, 1995; Hanada et al., 2008). As shown in the simulations based on the literature data on the K m values for the R- and S-VER pathways for the other metabolite formation and N-demethylation (see Supplemental Table 3), the S-VER/R-VER K m ratios were 1.9 and 1.34, respectively. The K m estimates obtained for racemic or (R,S)-VER for the other metabolites or N-demethylation pathways (Table 5) were 10.4 and 14.1 μM, respectively. These averaged K m values for (R,S)-verapamil corresponded to (1.9x + 1x)K m,R-VER/2 and (1.3x + 1x)K m,R-VER/2, respectively, for the other metabolites and N-demethylation pathways. From our calculations, the K m values for R- and S-VER for N-demethylation are 9.04 and 11.8 μM, whereas those for the other metabolites pathways are 9.72 and 18.5 μM, respectively (Supplemental Table 3). Simulations with these K m constants obtained for R-VER and S-VER showed that the difference in stereoisomeric kinetics was minor (Fig. 6), and existed only for the early-in-time liver profiles at the lowest concentration. Also, the summed simulated R- and S-VER values taken to represent racemic (R,S) profiles were similar to those observed (see Fig. 7).

Fig. 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 6.

Simulated concentrations of verapamil enantiomers (S- and R-verapamil) and their corresponding formed metabolites, S- and R-norverapamil, in reservoir perfusate, and amounts in liver and bile at initial concentrations of 0.5 μM (A), 25 μM (B), and 50 μM (C) (one-half of the concentrations of racemic verapamil used experimentally). Stereoselective elimination, with the K m values shown in Supplemental Table 3, was considered for the simulations (with the equations shown in Appendix C). The black arrow in (A) shows that concentrations of R-VER in liver were slightly lower than those of S-VER at comparable perfusion time, due to the lower K m value for faster metabolism of R-VER at the lowest dose.

Fig. 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 7.

Simulated, concentrations or amounts of racemic verapamil (S-VER + R-VER) and formed norverapamil (S-NOR + R-NOR) in reservoir perfusate and liver (A, C, and E) and excreted amounts in bile (B, D, and F) (data from Fig. 6, with consideration of stereoselective elimination; black solid or dotted lines) were presented against experimental (colored symbols) and fitted results on racemic verapamil and formed norverapamil from PBPK modeling (colored lines, also shown in Fig. 4), without consideration of stereoselective elimination. The observed, fitted (lack of consideration of stereoselective elimination) and simulated data (with consideration of stereoselective elimination) of summed enantiomeric verapamil and norverapamil concentrations/amounts were virtually identical, showing that stereoselective elimination was unimportant under the conditions examined.

In addition, inhibitory effects of verapamil on its own P-gp-mediated biliary secretion and metabolite formation, events observed in humans, should be negligible during the short 90 minute perfusion, inasmuch as lack of precedence in rat liver for auto-inhibition (Obach, 1999; Shibata et al., 2002; Axelsson et al., 2003; Guo et al., 2007; Yin et al., 2011). Moreover, incorporation of time-dependent inhibition components (e −kI t, where k I is the in vivo inhibition rate constant, which equals the in vitro inactivation rate constant, k obs) (Wang et al., 2004) to the various verapamil eliminatory processes (or CLint) failed to significantly improve the goodness-of-fit of the models (see Supplemental Table 4). The fitted k I values for metabolic formation ( Embedded Image and Embedded Image ) and biliary secretions ( Embedded Image ) were low, ranging from 0.002 to 0.14 min−1 (data not shown). Fitted results, with and without auto-inhibition, support the notion that time-dependent inactivation of various verapamil elimination pathways would not perturb the pharmacokinetics of verapamil in the perfused rat liver.

Discussion

Verapamil has been considered to be a good candidate to study transporter-enzyme interplay, due to the fact that it is a substrate of both CYP450 and P-gp. Under linear kinetic conditions, inhibition or induction of metabolism versus excretion has been studied, since induction of the competing pathway reduces the intracellular substrate concentration for the given pathway (Pang et al., 2009). Data in the perfused rat liver preparation not only revealed nonlinear metabolism and secretion of verapamil (1, 50, and 100 μM) but also a poor extent of biliary excretion (Table 5), suggesting that verapamil is a poor candidate for the examination of transporter-enzyme interplay between CYP450 and P-gp in the rat liver. Additionally, the fact that fitted K m values for P-gp mediated excretion (4.75 and 5.36 μM) are lower than those for the metabolic pathways, including N-demethylation or formation of other metabolites (10.4 and 14.1 μM), suggests that transporter-enzyme interplay kinetics would be extremely complex under these nonlinear conditions. Induction of liver metabolism would not greatly perturb the extent of secretion since biliary excretion of verapamil was low, despite that verapamil would serve as a better candidate for the study of transporter-enzyme interplay in the intestine since there is more verapamil secreted (Johnson et al., 2001). However, complexity in data interpretation still applies when nonlinearity metabolism and P-gp-mediated secretion exists, and when the secreted substrate can undergo reabsorption (Tam et al., 2003). Under nonlinear conditions, induction of intestinal P-gp that increases secretion could result in an increased extent of metabolism upon reabsorption of the secreted drug to furnish a more efficient metabolic system, as suggested by Pang and colleagues (Tam et al., 2003; Pang et al., 2009; Fan et al., 2010).

By contrast, verapamil is ideal for the consideration of metabolite (norverapamil) kinetics since there are competing metabolic pathways. Thus, we have used this precursor-product pair to develop strategies for the estimation of formation clearances of metabolites under linear and nonlinear conditions, for a drug that is bound to both RBC and plasma proteins. After consideration of the substantial loss due to tubing adsorption, we found that the extraction ratios of verapamil and preformed norverapamil were both high in the rat liver preparation, and the measured verapamil (∼7%–8% dose) and formed norverapamil (∼3%–5% dose) recovered in the liver tissue at the end of perfusion were similar to those in previous perfusion studies involving a high dose of verapamil administered (Mehvar et al., 1994). However, our reported average unbound fraction of verapamil in the blood perfusate ( Embedded Image ) was higher (0.45 for racemic-VER versus 0.12 for S-VER and 0.22 for R-VER) (Mehvar and Reynolds, 1996), and the total intrinsic hepatic clearance of verapamil ( Embedded Image ) was lower (30 ml/min versus about 1475 ml/min for S- and R-VER) compared with those reported by Mehvar et al. (1994); the product, Embedded Image was also smaller (16 ml/min versus 234 to 491 ml/min) (Pang and Yang, 2015). The likely explanation is due to the consideration of verapamil binding to RBCs and tubing as well as nonlinear metabolism in our study. That verapamil elimination is dose dependent, with Embedded Image increasing and Embedded Image decreasing at higher dosing levels, confirmed nonlinearity in verapamil disposition, an observation that is in agreement with that inferred for the clinical pharmacokinetics of verapamil (Freedman et al., 1981; Shand et al., 1981; Anderson et al., 1982; Tartaglione et al., 1983; Toffoli et al., 1997; Maeda et al., 2011).

Our strategy to study the kinetics of verapamil in the determination of whether formation of norverapamil is the major clearance pathway rests on capacity-limiting PBPK models, since these models are excellent for including sequential metabolism and drug- or metabolite-specific transporters and enzymes, distinguishing between kinetics of formed and preformed metabolites (Pang and Durk, 2010). In this study, we considered the nonlinear binding to plasma proteins and distribution into red cells and accounted for binding to tubing to fit the liver perfusion data of verapamil (1, 50, and 100 μM) and preformed (1.5 and 5 μM) and formed norverapamil simultaneously. For this reason, metabolite modeling of formed norverapamil was achieved with greater assurance. Additionally, the amount of verapamil and norverapamil binding to Tygon tubing was not negligible (∼10%–25% of total dose), suggesting the necessity to account for drug loss due to tubing adsorption to avoid overestimation of hepatic clearance (calculated from Embedded Image ). We further emphasized that plasma protein binding estimated in a test tube must be extended to include the total plasma space (plasma + Disse space) within the liver vasculature to consider the equilibrium of the drug between plasma and RBCs (see Appendix B). Different PBPK models with increasing degrees of nonlinearity were then built to account for linear versus nonlinear metabolism and biliary excretion by defining intrinsic clearance CLint in terms of V max and K m, and the best-fitted model was determined according to the AIC score, SSWR (F test), and diagnostic plots (prediction versus observation and weighted residuals versus time plots).

Our final PBPK model that showed the most significant improvement, with minimum AIC and SSWR, incorporated multiple nonlinear factors describing metabolic and biliary secretion pathways of verapamil and norverapamil. PBPK modeling showed that stereoselective binding, metabolism, and auto-inhibition were unimportant during our dosing conditions (see Figs. 6 and 7; Supplemental Tables 3 and 4). Although PBPK modeling of verapamil and norverapamil that included nonlinear metabolic clearance and biliary excretion has been previously published (Neuhoff et al., 2013; Wang et al., 2013), to our knowledge, we are the first one to describe sequential metabolism and nonlinear kinetics of verapamil by fitting verapamil and its preformed and formed metabolite levels in blood perfusate and bile simultaneously. The fitted results have shown that the rate of verapamil and formed and preformed norverapamil diffuse into the hepatocytes ( Embedded Image and Embedded Image ) at rates that are much greater than the perfusate blood flow (Q L), implying that these compounds did equilibrate rapidly between perfusate and liver. This finding is reasonable because verapamil and norverapamil are both lipophilic compounds. Since the efflux process outcompeted the biliary secretion ( Embedded Image ; Embedded Image ), verapamil and norverapamil were expected to accumulate in the reservoir, with limited amounts appearing in bile, events that are consistent with our observations. Our fitted unbound fraction of verapamil in liver tissue ( Embedded Image ) matched literature values in rats (0.02) (Yamano et al., 2000), but was lower than that in healthy subjects (0.09) (Giacomini et al., 1984).

The K m values determined in the final model (H) (Table 4) are quite low relative to perfusate concentration, ranging from 5 to 14 μM, suggesting at the high concentrations (50 and 100 μM), most of the enzymes and P-gp are saturated, and nonlinear kinetics will likely occur. The fitted K m values for metabolizing verapamil to norverapamil and other metabolites ( Embedded Image and Embedded Image , respectively) are in agreement with the literature values (Tracy et al., 1999; Hanada et al., 2008) after adjusting for Embedded Image or the unbound fraction in liver. The reported K m values for P-gp-mediated biliary secretion of verapamil ( Embedded Image ) are similar to that found in a previous study (K m = 4.1 μM), where the biliary excretion of verapamil mediated by P-gp (∼8 μM) was measured using human MDR1 reconstituted in liposomes (Kimura et al., 2007).

Proper estimation of norverapamil formation was adequately addressed in the study. With the assumption of linear pharmacokinetics, the fraction ( Embedded Image ) of hepatic clearance ( Embedded Image ) that forms norverapamil can be calculated using the formula Embedded Image (Table 2), where Embedded Image represents the fraction of Embedded Image that furnishes the formed norverapamil to the circulation, and Embedded Image denotes the hepatic availability of norverapamil, since Embedded Image needs to account for the sequential metabolism of locally formed norverapamil, which never reaches the systemic circulation (Pang and Gillette, 1979; Pang and Kwan, 1983). The Embedded Image value (0.19–0.23) from Mehvar et al. (1994) was similar to the value based on the lowest dose (0.23) estimated by AUC comparisons assuming linear kinetics. This value, however, remained an underestimation due to prevailing nonlinearity in metabolism compared with the value obtained upon comparison of the formation of intrinsic clearance normalized to the total intrinsic clearance (0.31).

In conclusion, our comprehensive analysis with PBPK modeling revealed that the metabolic pathways converting verapamil to norverapamil and to other metabolites as well as the biliary excretion of verapamil and norverapamil are concentration dependent in verapamil and norverapamil liver perfusion data using PBPK modeling. Stereoselective binding and metabolism (Figs. 6 and 7) and auto-inhibition (Supplemental Table 4) were unimportant. With the careful consideration of nonlinearity stepwise in PBPK modeling, we are able to unravel which elimination pathways are nonlinear and demonstrate a strategy to assess the fractional pathway of hepatic clearance responsible for either excretion of drug into bile or in formation of metabolites, after consideration of sequential metabolism of the metabolite and nonlinear protein binding of drug.

Acknowledgments

The authors thank Dr. Jubo Liu for conducting BSA and RBC binding experiment. The authors also thank Lundbeck, Denmark, for providing [3H]verapamil for the binding studies.

Appendix A

In Vitro Distribution of Verapamil into RBC and Plasma Protein Binding

Scheme 1 shows the exchange of verapamil between RBC and plasma and for binding of verapamil to BSA. Here, f RBC is the unbound fraction of drug in RBC and f RBC C RBC is the unbound drug concentration in RBC; correspondingly, f P is the unbound fraction of drug in plasma and f P C P is the unbound drug concentration in plasma. The rates of change of verapamil in RBC and plasma in in vitro red cell distribution studies, where verapamil in plasma (2% BSA) was incubated with blank blood perfusate (40% RBC and 2% BSA) to result in a perfusate composition of 20% RBC and 2% BSA, appear in the equations shown below with rate constants, Embedded Image ( Embedded Image ) and Embedded Image ( Embedded Image ), based on total concentrations, since f RBC is unknown Embedded Image (A1) Embedded Image (A2)

Scheme 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Scheme 1.

The binding of verapamil to one class of BSA binding site is given by

Embedded Image (A3)

where KA is the binding association constant, which equals k 1/k 2 or the ratio of the on- and off-rate constants for protein binding; n is the number of binding sites; [Pt] is the total BSA concentration; f P is the unbound fraction in plasma; and C p,u is the unbound concentration of verapamil in plasma. Equations (A1)–(A3) were used simultaneously to fit data obtained in the RBC distribution and protein binding studies, and the fitted results are shown in Fig. 2.

Appendix B

Incorporation of Verapamil Distribution into RBC and Plasma Protein Binding for Modeling Liver Perfusion Data

In perfused liver experiments, drug in sinusoidal RBC equilibrates with drug in the total sucrose space or V plasma,tot (sinusoidal plasma and sucrose Disse space). The mass balance equation at equilibrium is summarized as follows:

Embedded Image (B1)

Upon substitution, given that V plasma,tot is (1 + γ Suc)V P, where γ Suc is the interstitial or sucrose Disse space/plasma space,V p, and knowing that V p = V B(1-Hct), we obtain (Pang et al., 1988; Liu et al., 2005)

Embedded Image (B2)

From the mass balance, the sinusoidal blood (C B), plasma (C P), and red blood cell (C RBC) concentrations are expressed as follows:

Embedded Image (B3)

Upon substitution of eq. (B2) into eq. (B3), one obtains

Embedded Image (B4)

Nonlinearity in protein binding will result in a changing C B/C P ratio, and this could be accommodated by eq. (B4), where γ Suc equals 0.52 (Pang et al., 1988).

The fitted Embedded Image and Embedded Image were used to calculate the ratio of drug concentrations in RBCs and plasma (C RBC/C p) [eq. (B2)] and C B from (C B/C p) [eq. (B4)], the ratio of drug concentrations in perfusate blood to plasma. However, for the sake of simplicity, binding of verapamil to liver tissue and norverapamil to liver and plasma proteins was assumed to be linear.

Appendix C

PBPK Modeling of Verapamil and Formed and Preformed Norverapamil, with Consideration of Tubing Adsorption in Perfusion Experiments

In the following equations, superscripts VER, NOR, and NOR p denote the parent drug (verapamil) and its formed and preformed metabolite (norverapamil), respectively. Here, A and C denote the amount and concentration, respectively; Q L, V, and f represent the perfusate blood flow rate, volume, and the unbound fraction, respectively; subscripts t→R, R→t, in, and ef denote transfer from the tubing to the reservoir, from the reservoir to the tubing, influx from liver blood (LB) into liver, and efflux from liver back to liver blood, respectively. Other pharmacokinetic parameters such as metabolic clearance and biliary secretion have been described previously in the Materials and Methods section.

Tubing adsorption was first considered in the model by incorporating dose-specific adsorption rate constants for VER and NOR p . Then, the influx clearance ( Embedded Image or Embedded Image ) and efflux clearance ( Embedded Image or Embedded Image ) were assumed to be the same because verapamil and norverapamil are lipophilic drugs ( Embedded Image and 3.3, respectively) and can diffuse freely across the membrane. Transit compartments were used to account for the time delay (τ) of verapamil and norverapamil biliary secretion, and two transit compartments (Transit 1 and Transit 2, denoted with amounts, A tr1 and A tr2, respectively) best described the bile data (data not shown). Additionally, Embedded Image , Embedded Image , Embedded Image , and Embedded Image , and Embedded Image and Embedded Image were assigned as the same constants for preformed and formed norverapamil. Data of three input concentrations of verapamil, two input concentrations of preformed norverapamil, and the corresponding formed norverapamil from verapamil in perfusate, liver, and bile were fitted simultaneously using the PBPK model (Fig. 1).

For rates of change of verapamil in perfusate (B) adsorbed onto Tygon tubing (t) during verapamil perfusion without liver

Embedded Image (C1) Embedded Image (C2)

For rates of change of norverapamil (both preformed and formed) in perfusate (B) and adsorbed onto Tygon tubing (t) during norverapamil perfusion without liver

Embedded Image (C3) Embedded Image (C4)

Fitted values of Embedded Image , Embedded Image , Embedded Image , and Embedded Image are summarized in Table 1, and these values were fixed for later modeling with verapamil and preformed and formed norverapamil data in perfused liver preparations.

For rates of change of verapamil in reservoir perfusate (B), tubing, liver blood (LB, whole liver blood; LRBC, liver RBCs; LP, liver plasma), and liver (L) and bile compartments

Embedded Image (C5) Embedded Image (C6) Embedded Image (C7) Embedded Image (C8)

For modeling of the perfusion data, the unbound concentration of verapamil in blood perfusate (C B,u) was assumed to equal C P,u based on one class of binding sites

Embedded Image (C9)

We then converted the C P to C B based on the C B/C P ratio (R) in eq. (B4) according to the fitted Embedded Image and Embedded Image and γ Suc values, and express the unbound concentration in blood perfusate as

Embedded Image (C10)

The fitted constants obtained from the in vitro binding and distribution studies for n = 0.301, the binding association constant KA (0.0091 μM−1) and [Pt], and the total BSA concentration in plasma (1512 μM) were used (Appendix A).

For biliary excretion

Embedded Image (C11)

For the time delay of verapamil biliary secretion

Embedded Image (C12) Embedded Image (C13)

For the rates of formed norverapamil in the reservoir, tubing, liver blood perfusate, liver, and blood

Embedded Image (C14) Embedded Image (C15) Embedded Image (C16) Embedded Image (C17) Embedded Image (C18)

For the time delay of norverapamil biliary secretion

Embedded Image (C19) Embedded Image (C20)

For saturable N-demethylation, Embedded Image is

Embedded Image (C21)

For the saturable metabolic pathway, verapamil to other metabolites, Embedded Image is

Embedded Image (C22)

For the saturable verapamil and norverapamil biliary excretion, Embedded Image and Embedded Image are

Embedded Image (C23) Embedded Image (C24)

For the rates of change of preformed norverapamil in the perfusate (B), tubing, liver blood (LB), liver (L), and bile compartments

Embedded Image (C25) Embedded Image (C26) Embedded Image (C27) Embedded Image (C28) Embedded Image (C29)

For the time delay of preformed norverapamil biliary secretion

Embedded Image (C30) Embedded Image (C31)

For the nonlinear biliary excretion of norverapamil, Embedded Image is

Embedded Image (C32)

where Embedded Image may also be expressed as Embedded Image ; however, in the preliminary study, incorporation of this nonlinear term failed to significantly improve the goodness-of-fit when norverapamil data alone was fitted (data not shown), only first-order conditions were assumed, and nonlinearity was not considered when all data were fitted simultaneously.

Authorship Contributions

Participated in research design: Si, Sveigaard, Pang.

Conducted experiments: Si, Yang, Tang, Chow, Sveigaard, Pang.

Performed data analysis: Yang, Si, Tang, Pang.

Wrote or contributed to the writing of the manuscript: Yang, Si, Pang.

Footnotes

    • Received November 20, 2014.
    • Accepted February 3, 2015.
  • ↵ 1 Current affiliation: School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, P.R. China.

  • ↵ 2 Current affiliation: Pharmacy School, Shihezi University, Shihezi City, Xinjiang, P.R. China.

  • Q.J.Y. and L.S. contributed equally to this work.

  • L.S. and H.T. were supported by the China Scholarship Council (CSC), and K.S.P. is supported by Canadian Institute of Heath Research (CIHR).

  • The work was presented in part, at the 2014 AAPS Annual Meeting; San Diego, CA.

  • dx.doi.org/10.1124/dmd.114.062265.

  • ↵ Embedded ImageThis article has supplemental material available at dmd.aspetjournals.org.

Abbreviations

AIC
Akaike information criterion
AUC
area under the curve
BSA
bovine serum albumin
CL
clearance
Hct
hematocrit
HPLC
high-performance liquid chromatography
KHB
Krebs-Henseleit buffer
PBPK
physiologically based pharmacokinetic
RBC
red blood cell
SSWR
sum of squared weighted residuals
  • Copyright © 2015 by The American Society for Pharmacology and Experimental Therapeutics

References

  1. ↵
    1. Anderson P ,
    2. Bondesson U , and
    3. Sylvén C
    (1982) Clinical pharmacokinetics of verapamil in patients with atrial fibrillation. Eur J Clin Pharmacol 23 : 49–57.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Axelsson H ,
    2. Granhall C ,
    3. Floby E ,
    4. Jaksch Y ,
    5. Svedling M , and
    6. Sohlenius-Sternbeck AK
    (2003) Rates of metabolism of chlorzoxazone, dextromethorphan, 7-ethoxycoumarin, imipramine, quinidine, testosterone and verapamil by fresh and cryopreserved rat liver slices, and some comparisons with microsomes. Toxicol In Vitro 17 : 481–488.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Baillie TA ,
    2. Cayen MN ,
    3. Fouda H ,
    4. Gerson RJ ,
    5. Green JD ,
    6. Grossman SJ ,
    7. Klunk LJ ,
    8. LeBlanc B ,
    9. Perkins DG , and
    10. Shipley LA
    (2002) Drug metabolites in safety testing. Toxicol Appl Pharmacol 182 : 188–196.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Bhatti MM and
    2. Foster RT
    (1997) Pharmacokinetics of the enantiomers of verapamil after intravenous and oral administration of racemic verapamil in a rat model. Biopharm Drug Dispos 18 : 387–396.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Boxenbaum HG ,
    2. Riegelman S , and
    3. Elashoff RM
    (1974) Statistical estimations in pharmacokinetics. J Pharmacokinet Biopharm 2 : 123–148.
    OpenUrlCrossRefPubMed
  6. ↵
    1. Busse D ,
    2. Templin S ,
    3. Mikus G ,
    4. Schwab M ,
    5. Hofmann U ,
    6. Eichelbaum M , and
    7. Kivistö KT
    (2006) Cardiovascular effects of (R)- and (S)-verapamil and racemic verapamil in humans: a placebo-controlled study. Eur J Clin Pharmacol 62 : 613–619.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Choi JS and
    2. Burm JP
    (2008) Effect of pioglitazone on the pharmacokinetics of verapamil and its major metabolite, norverapamil, in rats. Arch Pharm Res 31 : 1200–1204.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Eichelbaum M ,
    2. Ende M ,
    3. Remberg G ,
    4. Schomerus M , and
    5. Dengler HJ
    (1979) The metabolism of DL-[14C]verapamil in man. Drug Metab Dispos 7 : 145–148.
    OpenUrlAbstract
  9. ↵
    1. Fan J ,
    2. Maeng HJ , and
    3. Pang KS
    (2010) Interplay of transporters and enzymes in the Caco-2 cell monolayer: I. effect of altered apical secretion. Biopharm Drug Dispos 31 : 215–227.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Fleckenstein A
    (1977) Specific pharmacology of calcium in myocardium, cardiac pacemakers, and vascular smooth muscle. Annu Rev Pharmacol Toxicol 17 : 149–166.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Freedman SB ,
    2. Richmond DR ,
    3. Ashley JJ , and
    4. Kelly DT
    (1981) Verapamil kinetics in normal subjects and patients with coronary artery spasm. Clin Pharmacol Ther 30 : 644–652.
    OpenUrlCrossRefPubMed
  12. ↵
    1. Fromm MF ,
    2. Dilger K ,
    3. Busse D ,
    4. Kroemer HK ,
    5. Eichelbaum M , and
    6. Klotz U
    (1998) Gut wall metabolism of verapamil in older people: effects of rifampicin-mediated enzyme induction. Br J Clin Pharmacol 45 : 247–255.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Garcia MA ,
    2. Aramayona JJ ,
    3. Bregante MA ,
    4. Fraile LJ , and
    5. Solans C
    (1997) Simultaneous determination of verapamil and norverapamil in biological samples by high-performance liquid chromatography using ultraviolet detection. J Chromatogr B Biomed Sci Appl 693 : 377–382.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Giacomini KM ,
    2. Massoud N ,
    3. Wong FM , and
    4. Giacomini JC
    (1984) Decreased binding of verapamil to plasma proteins in patients with liver disease. J Cardiovasc Pharmacol 6 : 924–928.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Guo Y ,
    2. Wang H , and
    3. Zhang C
    (2007) Establishment of rat precision-cut fibrotic liver slice technique and its application in verapamil metabolism. Clin Exp Pharmacol Physiol 34 : 406–413.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Hanada K ,
    2. Akimoto S ,
    3. Mitsui K ,
    4. Mihara K , and
    5. Ogata H
    (1998) Enantioselective tissue distribution of the basic drugs disopyramide, flecainide and verapamil in rats: role of plasma protein and tissue phosphatidylserine binding. Pharm Res 15 : 1250–1256.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Hanada K ,
    2. Ikemi Y ,
    3. Kukita K ,
    4. Mihara K , and
    5. Ogata H
    (2008) Stereoselective first-pass metabolism of verapamil in the small intestine and liver in rats. Drug Metab Dispos 36 : 2037–2042.
    OpenUrlAbstract/FREE Full Text
  18. ↵
    1. Johnson BM ,
    2. Charman WN , and
    3. Porter CJ
    (2001) The impact of P-glycoprotein efflux on enterocyte residence time and enterocyte-based metabolism of verapamil. J Pharm Pharmacol 53 : 1611–1619.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Keefe DL ,
    2. Yee YG , and
    3. Kates RE
    (1981) Verapamil protein binding in patients and in normal subjects. Clin Pharmacol Ther 29 : 21–26.
    OpenUrlCrossRefPubMed
  20. ↵
    1. Kimura Y ,
    2. Kioka N ,
    3. Kato H ,
    4. Matsuo M , and
    5. Ueda K
    (2007) Modulation of drug-stimulated ATPase activity of human MDR1/P-glycoprotein by cholesterol. Biochem J 401 : 597–605.
    OpenUrlAbstract/FREE Full Text
  21. ↵
    1. Kroemer HK ,
    2. Echizen H ,
    3. Heidemann H , and
    4. Eichelbaum M
    (1992) Predictability of the in vivo metabolism of verapamil from in vitro data: contribution of individual metabolic pathways and stereoselective aspects. J Pharmacol Exp Ther 260 : 1052–1057.
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Lee YH ,
    2. Perry BA ,
    3. Lee HS ,
    4. Kunta JR ,
    5. Sutyak JP , and
    6. Sinko PJ
    (2001) Differentiation of gut and hepatic first-pass effect of drugs: 1. Studies of verapamil in ported dogs. Pharm Res 18 : 1721–1728.
    OpenUrlCrossRefPubMed
  23. ↵
    1. Lemma GL ,
    2. Wang Z ,
    3. Hamman MA ,
    4. Zaheer NA ,
    5. Gorski JC , and
    6. Hall SD
    (2006) The effect of short- and long-term administration of verapamil on the disposition of cytochrome P450 3A and P-glycoprotein substrates. Clin Pharmacol Ther 79 : 218–230.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Liu L ,
    2. Mak E ,
    3. Tirona RG ,
    4. Tan E ,
    5. Novikoff PM ,
    6. Wang P ,
    7. Wolkoff AW , and
    8. Pang KS
    (2005) Vascular binding, blood flow, transporter, and enzyme interactions on the processing of digoxin in rat liver. J Pharmacol Exp Ther 315 : 433–448.
    OpenUrlAbstract/FREE Full Text
  25. ↵
    1. Maeda K ,
    2. Takano J ,
    3. Ikeda Y ,
    4. Fujita T ,
    5. Oyama Y ,
    6. Nozawa K ,
    7. Kumagai Y , and
    8. Sugiyama Y
    (2011) Nonlinear pharmacokinetics of oral quinidine and verapamil in healthy subjects: a clinical microdosing study. Clin Pharmacol Ther 90 : 263–270.
    OpenUrlCrossRefPubMed
  26. ↵
    1. McTavish D and
    2. Sorkin EM
    (1989) Verapamil. An updated review of its pharmacodynamic and pharmacokinetic properties, and therapeutic use in hypertension. Drugs 38 : 19–76.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Mehvar R and
    2. Reynolds J
    (1995) Input rate-dependent stereoselective pharmacokinetics. Experimental evidence in verapamil-infused isolated rat livers. Drug Metab Dispos 23 : 637–641.
    OpenUrlAbstract
  28. ↵
    1. Mehvar R and
    2. Reynolds J
    (1996) Reversal of stereoselectivity in the hepatic availability of verapamil in isolated perfused rat livers. Role of protein binding. Drug Metab Dispos 24 : 1088–1094.
    OpenUrlAbstract
  29. ↵
    1. Mehvar R ,
    2. Reynolds JM ,
    3. Robinson MA , and
    4. Longstreth JA
    (1994) Enantioselective kinetics of verapamil and norverapamil in isolated perfused rat livers. Pharm Res 11 : 1815–1819.
    OpenUrlCrossRefPubMed
  30. ↵
    1. Neuhoff S ,
    2. Yeo KR ,
    3. Barter Z ,
    4. Jamei M ,
    5. Turner DB , and
    6. Rostami-Hodjegan A
    (2013) Application of permeability-limited physiologically-based pharmacokinetic models: part II - prediction of P-glycoprotein mediated drug-drug interactions with digoxin. J Pharm Sci 102 : 3161–3173.
    OpenUrlCrossRefPubMed
  31. ↵
    1. Obach RS
    (1999) Prediction of human clearance of twenty-nine drugs from hepatic microsomal intrinsic clearance data: An examination of in vitro half-life approach and nonspecific binding to microsomes. Drug Metab Dispos 27:1350–1359.
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Pang KS
    (2009) Safety testing of metabolites: Expectations and outcomes. Chem Biol Interact 179 : 45–59.
    OpenUrlCrossRefPubMed
  33. ↵
    1. Pang KS ,
    2. Barker F ,
    3. Simard A ,
    4. Schwab AJ , and
    5. Goresky CA
    (1995) Sulfation of acetaminophen by the perfused rat liver: the effect of red blood cell carriage. Hepatology 22 : 267–282.
    OpenUrlCrossRefPubMed
  34. ↵
    1. Pang KS and
    2. Durk MR
    (2010) Physiologically-based pharmacokinetic modeling for absorption, transport, metabolism and excretion. J Pharmacokinet Pharmacodyn 37 : 591–615.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Pang KS and
    2. Gillette JR
    (1979) Sequential first-pass elimination of a metabolite derived from a precursor. J Pharmacokinet Biopharm 7 : 275–290.
    OpenUrlCrossRefPubMed
  36. ↵
    1. Pang KS and
    2. Kwan KC
    (1983) A commentary: methods and assumptions in the kinetic estimation of metabolite formation. Drug Metab Dispos 11 : 79–84.
    OpenUrlPubMed
  37. ↵
    1. Pang KS ,
    2. Lee WF ,
    3. Cherry WF ,
    4. Yuen V ,
    5. Accaputo J ,
    6. Fayz S ,
    7. Schwab AJ , and
    8. Goresky CA
    (1988) Effects of perfusate flow rate on measured blood volume, disse space, intracellular water space, and drug extraction in the perfused rat liver preparation: characterization by the multiple indicator dilution technique. J Pharmacokinet Biopharm 16 : 595–632.
    OpenUrlCrossRefPubMed
  38. ↵
    1. Pang KS ,
    2. Maeng HJ , and
    3. Fan J
    (2009) Interplay of transporters and enzymes in drug and metabolite processing. Mol Pharm 6 : 1734–1755.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Pang KS and
    2. Rowland M
    (1977) Hepatic clearance of drugs. I. Theoretical considerations of a “well-stirred” model and a “parallel tube” model. Influence of hepatic blood flow, plasma and blood cell binding, and the hepatocellular enzymatic activity on hepatic drug clearance. J Pharmacokinet Biopharm 5 : 625–653.
    OpenUrlCrossRefPubMed
  40. ↵
    1. Pang KS and
    2. Yang QJ
    (2015) Letter in Response to Rebuttal from Resa Mehvar: Fractional clearance for Verapamil N-demethylation in the isolated rat liver preparation; published ahead of print May15 2015; DOI:10.1124/dmd.113.052225
  41. ↵
    1. Pauli-Magnus C ,
    2. von Richter O ,
    3. Burk O ,
    4. Ziegler A ,
    5. Mettang T ,
    6. Eichelbaum M , and
    7. Fromm MF
    (2000) Characterization of the major metabolites of verapamil as substrates and inhibitors of P-glycoprotein. J Pharmacol Exp Ther 293 : 376–382.
    OpenUrlCrossRefPubMed
  42. ↵
    1. Robinson MA and
    2. Mehvar R
    (1996) Enantioselective distribution of verapamil and norverapamil into human and rat erythrocytes: the role of plasma protein binding. Biopharm Drug Dispos 17 : 577–587.
    OpenUrlCrossRefPubMed
  43. ↵
    1. Shand DG ,
    2. Hammill SC ,
    3. Aanonsen L , and
    4. Pritchett EL
    (1981) Reduced verapamil clearance during long-term oral administration. Clin Pharmacol Ther 30 : 701–706.
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Shibata Y ,
    2. Takahashi H ,
    3. Chiba M , and
    4. Ishii Y
    (2002) Prediction of hepatic clearance and availability by cryopreserved human hepatocytes: an application of serum incubation method. Drug Metab Dispos 30 : 892–896.
    OpenUrlPubMed
  45. ↵
    1. Sun L ,
    2. Zhang SQ , and
    3. Zhong DF
    (2004) In vitro identification of metabolites of verapamil in rat liver microsomes. Acta Pharmacol Sin 25 : 121–128.
    OpenUrlAbstract/FREE Full Text
  46. ↵
    1. Tam D ,
    2. Sun H , and
    3. Pang KS
    (2003) Influence of P-glycoprotein, transfer clearances, and drug binding on intestinal metabolism in Caco-2 cell monolayers or membrane preparations: a theoretical analysis. Drug Metab Dispos 31 : 1214–1226.
    OpenUrlAbstract/FREE Full Text
  47. ↵
    1. Tan E and
    2. Pang KS
    (2001) Sulfation is rate limiting in the futile cycling between estrone and estrone sulfate in enriched periportal and perivenous rat hepatocytes. Drug Metab Dispos 29 : 335–346.
    OpenUrlPubMed
  48. ↵
    1. Tartaglione TA ,
    2. Pieper JA ,
    3. Lopez LL , and
    4. Mehta J
    (1983) Pharmacokinetics of verapamil and norverapamil during long-term oral therapy. Res Commun Chem Pathol Pharmacol 40 : 15–27.
    OpenUrlCrossRefPubMed
  49. ↵
    1. Toffoli G ,
    2. Robieux I ,
    3. Fantin D ,
    4. Gigante M ,
    5. Frustaci S ,
    6. Nicolosi GL ,
    7. De Cicco M , and
    8. Boiocchi M
    (1997) Non-linear pharmacokinetics of high-dose intravenous verapamil. Br J Clin Pharmacol 44 : 255–260.
    OpenUrlCrossRefPubMed
  50. ↵
    1. Tracy TS ,
    2. Korzekwa KR ,
    3. Gonzalez FJ , and
    4. Wainer IW
    (1999) Cytochrome P450 isoforms involved in metabolism of the enantiomers of verapamil and norverapamil. Br J Clin Pharmacol 47 : 545–552.
    OpenUrlCrossRefPubMed
  51. ↵
    1. von Richter O ,
    2. Greiner B ,
    3. Fromm MF ,
    4. Fraser R ,
    5. Omari T ,
    6. Barclay ML ,
    7. Dent J ,
    8. Somogyi AA , and
    9. Eichelbaum M
    (2001) Determination of in vivo absorption, metabolism, and transport of drugs by the human intestinal wall and liver with a novel perfusion technique. Clin Pharmacol Ther 70 : 217–227.
    OpenUrlCrossRefPubMed
  52. ↵
    1. Walles M ,
    2. Thum T ,
    3. Levsen K , and
    4. Borlak J
    (2003) Metabolism of verapamil: 24 new phase I and phase II metabolites identified in cell cultures of rat hepatocytes by liquid chromatography-tandem mass spectrometry. J Chromatogr B Analyt Technol Biomed Life Sci 798 : 265–274.
    OpenUrlCrossRefPubMed
  53. ↵
    1. Wang J ,
    2. Xia S ,
    3. Xue W ,
    4. Wang D ,
    5. Sai Y ,
    6. Liu L , and
    7. Liu X
    (2013) A semi-physiologically-based pharmacokinetic model characterizing mechanism-based auto-inhibition to predict stereoselective pharmacokinetics of verapamil and its metabolite norverapamil in human. Eur J Pharm Sci 50 : 290–302.
    OpenUrlAbstract/FREE Full Text
  54. ↵
    1. Wang YH ,
    2. Jones DR , and
    3. Hall SD
    (2004) Prediction of cytochrome P450 3A inhibition by verapamil enantiomers and their metabolites. Drug Metab Dispos 32 : 259–266.
    OpenUrlCrossRefPubMed
  55. ↵
    1. Woodcock BG ,
    2. Schulz W ,
    3. Kober G , and
    4. Rietbrock N
    (1981) Direct determination of hepatic extraction of verapamil in cardiac patients. Clin Pharmacol Ther 30 : 52–56.
    OpenUrlAbstract/FREE Full Text
  56. ↵
    1. Yamano K ,
    2. Yamamoto K ,
    3. Kotaki H ,
    4. Takedomi S ,
    5. Matsuo H ,
    6. Sawada Y , and
    7. Iga T
    (2000) Quantitative prediction of metabolic inhibition of midazolam by erythromycin, diltiazem, and verapamil in rats: implication of concentrative uptake of inhibitors into liver. J Pharmacol Exp Ther 292 : 1118–1126.
    OpenUrlCrossRefPubMed
  57. ↵
    1. Yin J ,
    2. Meng Q , and
    3. Dong X
    (2011) Auto-inhibition of verapamil metabolism in rat hepatocytes of gel entrapment culture. Biomed Pharmacother 65 : 328–333.
    OpenUrlCrossRefPubMed
PreviousNext
Back to top

In this issue

Drug Metabolism and Disposition: 43 (4)
Drug Metabolism and Disposition
Vol. 43, Issue 4
1 Apr 2015
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Editorial Board (PDF)
  • Front Matter (PDF)
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for sharing this Drug Metabolism & Disposition article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
PBPK Modeling to Unravel Nonlinear Pharmacokinetics of Verapamil to Estimate the Fractional Clearance for Verapamil N-Demethylation in the Recirculating Rat Liver Preparation
(Your Name) has forwarded a page to you from Drug Metabolism & Disposition
(Your Name) thought you would be interested in this article in Drug Metabolism & Disposition.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Rapid CommunicationAccelerated Communication

Nonlinear Metabolite Kinetics of Verapamil

Qi Joy Yang, Luqin Si, Hui Tang, Helle H. Sveigaard, Edwin C. Y. Chow and K. Sandy Pang
Drug Metabolism and Disposition April 1, 2015, 43 (4) 631-645; DOI: https://doi.org/10.1124/dmd.114.062265

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Rapid CommunicationAccelerated Communication

Nonlinear Metabolite Kinetics of Verapamil

Qi Joy Yang, Luqin Si, Hui Tang, Helle H. Sveigaard, Edwin C. Y. Chow and K. Sandy Pang
Drug Metabolism and Disposition April 1, 2015, 43 (4) 631-645; DOI: https://doi.org/10.1124/dmd.114.062265
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Acknowledgments
    • Appendix A
    • Appendix B
    • Appendix C
    • Authorship Contributions
    • Footnotes
    • Abbreviations
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF + SI
  • PDF

Related Articles

Cited By...

More in this TOC Section

  • Use of CYP Transgenic Mouse Models
  • Production of Recombinant hsa-mir-27b
Show more Accelerated Communication

Similar Articles

  • Home
  • Alerts
Facebook   Twitter   LinkedIn   RSS

Navigate

  • Current Issue
  • Fast Forward by date
  • Fast Forward by section
  • Latest Articles
  • Archive
  • Search for Articles
  • Feedback
  • ASPET

More Information

  • About DMD
  • Editorial Board
  • Instructions to Authors
  • Submit a Manuscript
  • Customized Alerts
  • RSS Feeds
  • Subscriptions
  • Permissions
  • Terms & Conditions of Use

ASPET's Other Journals

  • Journal of Pharmacology and Experimental Therapeutics
  • Molecular Pharmacology
  • Pharmacological Reviews
  • Pharmacology Research & Perspectives
ISSN 1521-009X (Online)

Copyright © 2021 by the American Society for Pharmacology and Experimental Therapeutics