![]() |
|
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Department of Pharmacology, Faculty of Medicine, the Chinese University of Hong Kong, Shatin, N.T., Hong Kong, SAR (R.Y., N.L.K., S.-L.L., G.L.); and Sinoveda Canada Inc., Edmonton, Alberta, Canada (Y.K.T.)
(Received July 20, 2007; Accepted November 21, 2007)
| Abstract |
|---|
|
|
|---|
In the past decades, the phytochemistry of Chuanxiong has been studied extensively. Nineteen compounds from L. chuanxiong have been reported to be pharmacologically active (Naito et al., 1995
; Ko et al., 1997
, 1998
, 2002
; Matsumoto et al., 1998
; Chong and Feng, 1999
; Shih et al., 2002
). In particular, phthalides are present in large quantities in Chuanxiong essential oil; these compounds have been shown to have in vitro and in vivo pharmacological activities (Naito et al., 1995
; Ko et al., 1998
, 2002
; Matsumoto et al., 1998
; Chong and Feng, 1999
; Cao et al., 2006
; Chan et al., 2006
, 2007
; Du et al., 2006
; Lu et al., 2006
). Ligustilide (Z-3-butylidene-4,5-dihydroisodenzofuranone) (Fig. 1), a phthalide derivative, is the most abundant constituent in the herb (Li et al., 2002
, 2006
; Yan et al., 2005
). Vasodilatation, antiplatelet aggregation, antithrombotic, serotonergic activity, and antiproliferative properties of ligustilide have been well documented (Naito et al., 1995
; Cao et al., 2006
; Deng et al., 2006
; Du et al., 2006
; Lu et al., 2006
; Chan et al., 2007
). Hence, ligustilide is considered to be a key active ingredient in Chuanxiong.
|
The present studies investigated pharmacokinetics of ligustilide in its pure form and in a Chuanxiong extract and also evaluated the effects of administration routes on pharmacokinetics of ligustilide. Furthermore, in vitro studies were conducted to study ligustilide metabolism in greater detail.
| Materials and Methods |
|---|
|
|
|---|
|
In Vivo Protocols. Male Sprague-Dawley (SD) rats (200–230 g) were supplied by the Laboratory Animal Service Center, The Chinese University of Hong Kong. Animals were housed under standard conditions of temperature, humidity, and light. Food and water were provided ad libitum. A well established in vivo pharmacokinetic model in the rat (Su et al., 2006
; Yan et al., 2007
) was adopted in the present study. The day before administration, a jugular vein cannula was implanted for blood sampling under light anesthesia with diethyl ether. Rats were divided randomly into five groups (five animals in each group). For animals receiving i.v. dosing, a femoral vein cannula was also implanted. After cannulation, rats were allowed to recover and fasted overnight with free access to water. For i.v. administration, the "purified" ligustilide sample (ligustilide dose, 15.6 mg/kg) or Chuanxiong extract (100 mg/kg, ligustilide dose, 14.9 mg/kg) formulated in Pharmatek formulation-6 (Pharmatek, San Diego, CA) was injected into the left femoral vein through the catheter. A single dose of ligustilide formulated in normal saline containing 3% Tween 80 was given p.o. (100, 360, or 500 mg/kg) or i.p. (26 or 52 mg/kg) to individual groups of rats. Blood (
0.3 ml/sample) was collected into heparinized Eppendorf tubes at appropriate time intervals over a 48-h period as illustrated in Fig. 2. After each collection, 0.3 ml of saline containing 40 IU of heparin/ml was injected to compensate for blood loss. Blood samples were centrifuged at 3000g for 10 min, and plasma samples were harvested. Each plasma sample (160 µl) was mixed with 5 µl of internal standard (
-naphthoflavone, 40 µg/ml) and 165 µl of acetonitrile, vortexed, and centrifuged at 16,000g for 5 min. The supernatant was filtered with a 0.45-µm syringe filter, and 100 µl of the filtrate was subjected to HPLC analysis.
|
Stability of Ligustilide in Simulated Gastric and Intestinal Fluids. One volume of ligustilide (1 mg/ml) in 3% Tween 80 saline solution was spiked into three volumes of simulated gastric fluid containing no pepsin (0.05 M sodium chloride adjusted to pH 1.5 using HCl) (The United States Pharmacopeial Convention, 2006
) or simulated intestinal fluid containing no pancreatin (0.05 M sodium dihydrogen phosphate buffer adjusted to pH 6.8 using NaOH) (The United States Pharmacopeial Convention, 2006
). The mixtures were incubated for 1 h at 37°C. Aliquots (100 µl) of incubates were collected at 0, 5, 15, 30, and 60 min. Each aliquot was immediately extracted with ethyl acetate (200 µl), followed by centrifugation at 500g for 5 min. After removing the organic layer, the sample was extracted again with the same procedure. The organic layers were combined, and an aliquot (150 µl) was mixed with 150 µl of methanol containing 1.75 µgof
-naphthoflavone. The resultant solution was filtered, and an aliquot (100 µl) of the filtrate was subjected to HPLC analysis. Each test was conducted in triplicate. Degradation of ligustilide in the simulated gastric and intestinal fluids was determined by comparing peak area ratios (ligustilide/
-naphthoflavone) of the samples with that at time 0.
Identification of Ligustilide and Its Metabolites. Ligustilide and its metabolites senkyunolide I (L1), senkyunolide H (L2), and butylidenephthalide (L13) in biological samples were unambiguously identified by comparing their retention time, MS, and UV spectra with those of authentic compounds. Identities of other metabolites were tentatively elucidated using an on-line HPLC/UV/MS method as described blow.
The chromatographic condition was the same as described under Quantitative Analysis. The online HPLC/UV/MS analysis was performed on a PerkinElmer series 200 liquid chromatograph connected to a PerkinElmer Sciex API-2000 triple-quadrupole mass spectrometer (PerkinElmerSciex Instruments, Boston, MA) equipped with electrospray ionization interface. Positive electrospray ionization/MS was conducted under the following conditions: nebulizing gas, 40 psi; auxiliary gas, 60 psi; curtain gas, 30 psi; turbo-ion spray temperature, 400°C; declustering potential, 61 V; focusing potential, 380 V; and ionization potential, 5500 V. Full scan mass spectrum was obtained over a range of m/z 150 through 700. Mass spectrum was recorded with a PE Sciex MassChrom data system (version 1.1.1).
Quantitative Analysis. HPLC/UV method. The HPLC system consisted of an Agilent series 1100 HPLC/diode array detector (Agilent Technologies, Palo Alto, CA) and a Waters (Milford, MA) Symmetry C18 column (5 µm, 150 x 4.6 mm) coupled with a guard column (Waters Spherisorb S5 ODS2, 10 x 4.6 mm). The samples were eluted using a mobile phase containing A (0.25% aqueous acetic acid, v/v) and B (methanol). A previously reported gradient system was adopted for the analysis of both in vivo and in vitro samples (Yan et al., 2005
). Ligustilide was monitored at a wavelength of 284 nm.
Calibration curves of ligustilide. A methanol stock solution of ligustilide was diluted to the desired concentrations with methanol. Aliquot (5 µl) of each diluted solution was spiked into pooled plasma (three untreated male SD rats) to give concentrations ranging from 0.2 to 25 mg/l (0.2, 0.8, 2, 4, 8, 16, and 25 mg/l) for the in vivo study, or into Tris buffer to give a concentration ranging from 1.9 to 95 mg/l for the in vitro study. The resultant samples were mixed thoroughly, then treated and analyzed in the same manner as described under In Vivo Protocols. Samples of each concentration were analyzed in triplicate. Calibration curves were derived by plotting the peak area ratios of ligustilide to the internal standard as a function of concentration of ligustilide.
Method validation. Three concentrations of ligustilide at high, medium, and low levels of the corresponding calibration curves were chosen to determine the intraday and interday variability. Samples at each concentration were analyzed in triplicates to examine intraday variability. The interday variability was determined on three separate days over a period of 1 week. Peak area ratio for each concentration was measured, and concentration of ligustilide was calculated from the corresponding calibration curve. Relative standard deviation and percentage difference between amounts spiked and determined (recovery) were taken as measures of precision and accuracy, respectively.
Aliquots of methanol stock solution were also spiked into blank plasma or incubation buffer to determine limits of quantification (LOQ) of ligustilide. The resultant samples were treated and analyzed using the procedures described in this section. LOQ was determined when the concentration of a sample has a signal-to-noise ratio at or higher than 5.
Analysis of Pharmacokinetic Parameters. Pharmacokinetic parameters were calculated from the corresponding plasma concentration-time curves using noncompartmental analysis (WinNonlin 4.0, Pharsight, Mountain View, CA). The area under the plasma concentration-time curve from time 0 to time infinity (AUC0
) was calculated using the linear/log trapezoidal rule-extrapolation method, in which the linear trapezoidal rule was used for the calculation of the area during the ascending phase with the logarithmic trapezoidal rule for the declining phase. Cmax and Tmax were determined from the concentration-time profile. Other pharmacokinetic parameters, including apparent plasma clearance (CL/F), apparent volume of distribution based on the terminal phase (Vd/F), mean residence time (MRT), and absolute bioavailability (F) were calculated using the following equations:
![]() |
z is the terminal rate constant; AUMC0
is the area under the first moment curve from time 0 to infinity; Dosei.v. is the i.v. dose of ligustilide; and Dosep.o./i.p. is the p.o. or i.p. dose of ligustilide. AUCi.v. is the area under the plasma concentration-time curve of ligustilide after i.v. administration; AUCp.o./i.p. is the corresponding area under the plasma concentration-time curve of ligustilide after p.o. or i.p. administration. Statistical Analysis. All the data were expressed as mean ± S.D. Unpaired Student's t test was used to compare results between two groups of animals receiving ligustilide i.p., whereas one-way analysis of variance followed by Bonferroni post hoc test was performed to compare results among different groups. Level of p < 0.05 was considered significant.
| Results |
|---|
|
|
|---|
Pharmacokinetics of Ligustilide. Figure 2 shows plasma concentration-time curves of ligustilide after three routes of administration. The calculated pharmacokinetic parameters are summarized in Table 2. After i.v. administration, ligustilide exhibited extensive distribution in the body (Vd, 3.76 ± 1.23 l/kg) and rapid elimination (t1/2, 0.31 ± 0.12 h) from plasma. The plasma concentration decreased below LOQ after 1.5 h (Fig. 2A). When a low dose of ligustilide (26 mg/kg) was given i.p. (Fig. 2B), it was rapidly absorbed (Tmax, 0.05 ± 0.02 h) and quickly eliminated (t1/2, 0.36 ± 0.05 h). The i.p. bioavailability was estimated to be 52%, indicating an extensive hepatic first-pass metabolism. When a higher i.p. dose (52 mg/kg) of ligustilide was given, a markedly higher Cmax and a significantly lower CL/F were obtained compared with the lower i.p. dose; the absolute bioavailability was approximately 98%, suggesting nonlinear and dose-dependent pharmacokinetics of ligustilide.
|
In the case of p.o. administration of ligustilide, plasma concentration-time profile and pharmacokinetic parameters could only be obtained at the highest dose (500 mg/kg). At lower doses, data were not meaningful because plasma levels of ligustilide were either below detectable limit (at 100 mg/kg) or they were not high enough for accurate pharmacokinetic analysis. At the dose of 500 mg/kg (Fig. 2C; Table 2), ligustilide was rapidly absorbed (Tmax, 0.36 ± 0.19 h) reaching a Cmax of 0.66 ± 0.23 µg/ml. Plasma concentration of ligustilide declined in a multiphase manner. Levels decreased quickly to approximately 120 ng/ml within 4 h, and this level was maintained for another 4 to 8 h before decreasing below LOQ. The longer-lasting terminal phase (4–8 h) significantly increased the retention of ligustilide (MRT, 5.14 ± 1.56 h) when compared with that of i.v. and i.p. data (Table 2). Values of CL/F were significantly higher and dose-normalized AUC lower than those obtained after i.v. and i.p. dosing. Oral bioavailability was estimated to be 2.6% at the 500-mg/kg dose.
The elimination of ligustilide when administered in the form of Chuanxiong extract i.v. was significantly faster than that for its pure form (CL, 20.35 ± 3.05 versus 9.14 ± 1.27 l/h/kg; p < 0.01) (Fig. 2A; Table 2), suggesting interaction between ligustilide and components present in the extract.
Identification of Metabolites of Ligustilide. Fifteen peaks (L1–L15) were found as potential metabolites in vivo (Fig. 1, A–C). Among them, eight peaks (L1–L2, L4–L8, and L13) were also found in the in vitro incubation studies with commercially available pure ligustilide (Fig. 1E). Formation of these eight metabolites required the presence of NADPH-regenerating system, indicating the involvement of NADPH-dependent monooxygenases. The spectra of seven of eight metabolites, except L8, found in both in vivo and in vitro studies were obtained using the online HPLC/MS analysis (Fig. 3).
|
|
Structures of other three metabolites (L4–L6) were tentatively assigned. Mass spectrum of L4 (Fig. 3) exhibited the characteristic ions of [M+Na]+ at m/z 229, [M+H]+ at m/z 207, and [M+H-H2O]+ at m/z 189 as the base peak corresponding to the loss of an H2O molecule from the primary/secondary hydroxyl group. The mass data showed a molecular mass of 206 Da, 16 mass units higher than that of ligustilide, indicating a monooxygenated metabolite of ligustilide. Moreover, the UV spectrum of L4 was similar to that of ligustilide, suggesting the presence of an unchanged chromophore of the parent compound. Therefore, oxidation might occur at the end of 3-substituted side chain of ligustilide, and L4 was tentatively identified as 11-hydroxyligustilide. This metabolite was expected to be generated from direct oxidation of ligustilide (Fig. 4). L5 and L6 showed similar mass spectra with a protonated molecular ion ([M+H]+) at m/z 514, which corresponded to the molecular mass of glutathione (GSH) conjugate of monohydroxylated ligustilide. The fragmentation ions at m/z 385 and 207 represented characteristic ions as a result of the loss of glutamate ([M+H-Glutamate]+) and GSH ([M+H-GSH]+) moiety, respectively. Furthermore, the diagnostic ion resulting from the loss of both GSH and H2O ([M+H-GSH-H2O]+) at m/z 189 indicated the presence of a hydroxyl group in both metabolites. Therefore, these two metabolites were assigned as isomers of glutathione conjugate of monohydroxylated ligustilide (Fig. 4). However, the substituted positions for both GSH and hydroxyl group were unable to be unequivocally assigned. Based on the identified seven metabolites, the metabolic pathways of ligustilide in both in vivo and in vitro studies were proposed as shown in Fig. 4.
In Vitro Metabolic Profiles. Ligustilide was metabolized neither in the rat intestinal S9 fractions nor microsomes. When pure ligustilide (300 µM) was incubated with rat hepatic S9 preparations, approximately 60% ligustilide was metabolized in the presence of NADPH-regenerating system (Fig. 1D), but no significant metabolism was observed in the absence of NADPH-regenerating system (Fig. 1E). A similar metabolic pattern was found when ligustilide was incubated with hepatic microsomes. After ligustilide was incubated with hepatic microsomes in the presence of NADPH-regenerating system for 1 h, more than 98% of the substrate was metabolized. The most abundant metabolite generated was L1 (
21% of ligustilide), whereas L13 was identified as a minor metabolite (0.88% of ligustilide). The metabolic rate was rapid, and significant saturation was not attained at concentrations as high as 800 µM (Fig. 5). However, because of solubility constraint, accurate Michaelis-Menten parameters, such as Vmax and Km, for ligustilide could not be determined. The results indicated that ligustilide mainly underwent NADPH-dependent metabolism in the rat liver, which might be responsible for its extensive hepatic first-pass effect observed in vivo.
|
Stability of Ligustilide in Simulated Gastric and Intestinal Fluids. Because no significant intestinal metabolism of ligustilide was observed, the stability of ligustilide in gastrointestinal tract was examined as a potential cause of low oral bioavailability. No significant degradation of ligustilide occurred in simulated gastric fluid and about 15 ± 8% degraded in simulated intestinal fluid within 60 min (detailed data not shown).
| Discussion |
|---|
|
|
|---|
Chuanxiong is one of the most commonly prescribed CM herbs for the treatment of cardiovascular diseases in China. However, pharmacokinetics of the main bioactive ingredients in this herb is largely unknown. There is only one published study on the oral pharmacokinetics of ligustilide (Shi et al., 2006
). However, it is difficult to compare the reported results directly with our current findings because the published protocol and sample preparations were not described accurately in detail. Moreover, oral bioavailability was not determined (Shi et al., 2006
). In the present study, pharmacokinetics of ligustilide was investigated. The isolated ligustilide sample, which contained 5% butylidenephthalide with the composition ratio of two phthalides similar to that in the extract (approximately 18:1) (Table 1), was used because 1) concentration of commercially available ligustilide was too low for in vivo dosing; and 2) pure ligustilide is light and thermally labile but stable when kept in solvents, media, or mixed with other compounds (Cui et al., 2006
; Li et al., 2007
).
The results showed that ligustilide was rapidly absorbed after i.p. and p.o. dosing and exhibited extensive distribution and rapid elimination after i.v. administration. The terminal phase after p.o. administration was significantly longer (p < 0.001) than that of other two routes (Fig. 2). Given the oral bioavailability is less than 3%, the concentration achieved at the portal system during absorption was unlikely to be high enough to saturate hepatic enzymes; therefore, nonlinear metabolism that leads to slower elimination is unlikely. However, the mechanism behind this observation is not clear. In vitro metabolic studies (Fig. 5) suggested that hepatic elimination could be quantitatively the most important route of elimination and partly contribute to the low oral bioavailability of ligustilide. On the other hand, the role of intestine leading to a low oral bioavailability is uncertain. Our previous study using the Caco-2 cell model revealed high intestinal absorption of ligustilide with an absorptive permeability (Papp, 4.19 x 10–5 cm/s) comparable with that of highly absorbable drugs (Ko et al., 2003
), and the present in vitro studies showed that the intestinal metabolism and degradation during the time frame of absorption are unlikely contributors to the extensive loss of ligustilide in the gastrointestinal tract. Based on the data at hand, a plausible explanation is the rapid metabolism triggered by enzymes present in the intestinal lumen, which include but are not limited to that from intestinal flora. Moreover, lipophilic drugs are known to bind to tissues nonspecifically; this process is extensive and saturable (Hussain et al., 1994
). It is not clear whether ligustilide would extensively bind to intestinal and liver tissues. More studies are warranted.
Dose-dependent pharmacokinetics was observed after i.p. administration of ligustilide. The mechanisms responsible for this observation are not known. Saturation of metabolizing enzymes is unlikely because nonsaturable metabolism was shown at concentrations as high as 800 µM in vitro. Furthermore, the elimination profiles after low and high i.p. dose were not significantly different from each other (Fig. 2B). The total body clearance value was much higher than that of hepatic blood flow, suggesting extrahepatic elimination is involved. However, it is not certain as to how extrahepatic elimination is involved. The disproportionate increase in AUC after a high dose of ligustilide may be because of saturation of nonspecific binding in the liver. More study is required to elucidate this hypothesis.
In this study, constituents in Chuanxiong extract were found to significantly increase clearance of ligustilide (CL, 20.35 ± 3.05 versus 9.14 ± 1.27 l/h/kg; p < 0.01). The elimination of ligustilide is so high that it is likely to be flow-limited. Because Chuanxiong and its products are used for the relief of pain induced by blood stagnation (State Pharmacopoeia Commission of P.R. China, 2005
), it is possible that the increase in ligustilide clearance may be a direct result of hepatic blood flow change. However, to date there is no published evidence supporting this contention. Interestingly, no significant pharmacokinetic interactions between senkyunolide A and components in Chuanxiong were found previously (Yan et al., 2007
). Mechanisms responsible for pharmacokinetic interaction of ligustilide are unknown and worthy of investigation. The results revealed a variable and unpredictable pharmacokinetics of ligustilide, depending on components concurrently administered. Nevertheless, it is common to take more than one ingredient in herbal therapy. Therefore, therapeutic outcomes may vary when different Chuanxiong-based herbal preparations are used.
Eight ligustilide metabolites were identified both in vivo and in vitro. Butylidenephthalide, senkyunolide I, senkyunolide H, and 3-hydroxybutylphthalide were unequivocally characterized; 11-hydroxyligustilide and two isomers of hydroxyligustilide glutathione conjugate were tentatively identified, whereas structure of L8 was not elucidated. Several other peaks, which were only observed in chromatograms of plasmas after dosing but not in the incubated samples and blank plasma (Fig. 1), showed no ligustilide-related ions in their mass spectra and thus were not identified. The proposed metabolic pathways of ligustilide were reported for the first time (Fig. 4). All the metabolites were formed via oxidation pathways mediated by NADPH-dependent monooxygenases, and two were subsequently converted to glutathione conjugates. Butylidenephthalide, senkyunolide I, and senkyunolide H also occur naturally in Chuanxiong (Li et al., 2003
; Yan et al., 2005
) and have been reported to produce various pharmacological activities (Kobayashi et al., 1992
; Liang et al., 2005
; Chan et al., 2006
). Being naturally present in herb and also as metabolites in vivo, it obviously makes overall effects of all the bioactive ingredients and their individual contributions in the body more complicated and difficult to be predicted. Therefore, it should be emphasized that quality control of the bioactive ingredients in Chuanxiong and its products must consider not only herbal contents but also quantities of their active metabolites in vivo.
Although to date there are no approved remedies containing ligustilide as a single active ingredient, large varieties of Chuanxiong-based proprietary products and herb itself, which all contain ligustilide as one of the most abundant constituents, are widely used in China and other communities worldwide. Although most of these products are taken p.o. in CM practice, the present findings for ligustilide together with a low oral bioavailability of senkyunolide A reported previously (Yan et al., 2007
) provided important scientific data to challenge the validity of p.o. administration. The present study, in combination with our previous report on senkyunolide A, provided sound scientific evidence to support that another administrative route, such as sublingual delivery, to avoid extensive first-pass effects in the gut and the liver might be more preferable. Furthermore, pharmacokinetic interactions between ligustilide and other components in the herb showed that pharmacokinetic data obtained from a pure compound might be inadequate to predict its in vivo profile in herbal products.
| Acknowledgments |
|---|
| Footnotes |
|---|
Article, publication date, and citation information can be found at http://dmd.aspetjournals.org.
ABBREVIATIONS: CM, Chinese medicine; HPLC, high-performance liquid chromatography; MS, mass spectrometry; SD, Sprague-Dawley; LOQ, limit of quantification; AUC0
, area under the plasma concentration-time curve from time 0 to time infinity; MRT, mean residence time; GSH, glutathione.
Address correspondence to: Ge Lin, Department of Pharmacology, The Chinese University of Hong Kong, Shatin, Hong Kong, SAR. E-mail: linge{at}cuhk.edu.hk
| References |
|---|
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||