Department of Pharmaceutical Sciences, Faculty of Pharmacy (E.T.,
K.S.P.) and Department of Pharmacology, Faculty of Medicine (T.L.,
K.S.P.), University of Toronto, Toronto, Ontario, Canada
 |
Introduction |
Hepatic
drug clearance is regulated by hepatic blood flow, vascular and tissue
binding, transport, metabolism, and biliary excretion. Futile cycling,
the metabolic interconversion of two substrates involving different
enzymes, is an additional process that influences drug and metabolite
clearances. Futile cycling has been noted between 4-methylumbelliferone
(4-MU) and 4-methylumbelliferyl sulfate (4-MUS) (Ratna et al.,
1993
) and methylprednisone and methylprednisolone
(Ebling and Jusko, 1986
). The futile cycling of estrone
sulfate (E1S) and estrone (E1), which
represents a pharmacologically important biocycle that conserves and
regulates endogenous estrogens, however, has not been thoroughly investigated.
Since hepatic processing is a distributed-in-space phenomenon with
uptake, metabolism, and efflux occurring repeatedly in hepatocytes
along the direction of flow, it is expected that the futile cycling of
E1S and E1 in liver would be affected by a dual set of transporters and enzymes and their associated zonal
heterogeneities. It is expected that E1 would rapidly
diffuse through the cell membrane due to high lipophilicity, as found
in recent isolated rat hepatocyte studies (Tan and Pang,
2001
). The uptake of E1S into the rat liver
involves transport by members of the organic anion transporting
polypeptide superfamily, Oatp1, Oatp2 and Oatp4 (Jacquemin et
al., 1994
; Noé et al., 1997
;
Cattori et al., 2000
), the multispecific organic anion
transporter 3, Oat3 (Kusuhara et al., 1999
), and the
sodium-dependent taurocholate cotransporting polypeptide, Ntcp
(Hagenbuch et al., 1991
). However, a lack of acinar
heterogeneity was observed for the transport of E1S
(Tan et al., 1999
) and E1 (Tan and
Pang, 2001
) in rat liver.
In rat hepatocytes, E1S and E1 are highly bound
to liver tissue (Tan and Pang, 2001
). Consequently, the
high and nonlinear binding reduces the cellular unbound concentrations
of E1S and E1 for both metabolism and
excretion. Estrone sulfate is mainly deconjugated by arylsulfatase C, a
microsomal enzyme, to E1, which can be further metabolized
to estrone glucuronide (E1G), estradiol (E2),
estriol, their glucuronide and sulfate conjugates, and other minor
metabolites (Roy et al., 1987
). Although arylsulfatase C is found to be evenly dispersed in the liver acinus (El Mouelhi and Kauffman, 1986
), estrogen sulfotransferase, which is
responsible for the sulfation of E1, is predominantly
localized in the perivenous region (Tosh et al., 1996
;
Tan and Pang, 2001
). Both UDP-glucuronosyltransferase-1 (UGT1) and 2 (UGT2) are found to glucuronidate estrone (Tukey and Strassburg, 2000
), and the UGTs are predominantly localized in the perivenous region (Tosh and Burchill, 1996
). The
sulfo- and glucuronide conjugates are found in rat bile, with the
multidrug resistance-associated protein Mrp2 (or cmoat) mediating the
biliary excretion of E1G (Takikawa et al.,
1996
). But the transporter(s) involved in E1S
excretion into bile is unknown. Mrp2 does not appear to be involved
(personal communications with Dr. Dietrich Keppler, University of
Heidelberg, Germany).
Drugs bound to vascular components, namely plasma proteins and red
blood cells (RBCs), are expected to reduce drug clearances (Pang
and Rowland, 1977
; Pang et al., 1995
). Estrone
is bound to both plasma protein and erythrocytes. Consequently, only
3% estrone exists in the unbound form in human blood (Koefoed
and Brahm, 1994
). In addition, E1 is also
metabolized by 17
-hydroxysteroid dehydrogenase, a cytosolic enzyme
in erythrocytes of animals (Challis et al., 1973
;
Tsang, 1976
) and human (Mulder et al.,
1972
) to E2. By contrast, E1S is highly
bound to human plasma protein (1.6% unbound in plasma;
Rosenthal et al., 1972
) but not to erythrocytes. The
determination of the vascular binding and metabolism of E1 and of protein binding of E1S is another important aspect
toward the understanding of factors impacting the hepatic clearances of
E1 and E1S.
In this communication, the metabolic disposition of simultaneously
delivered tracers, [3H]E1S and
[14C]E1, was investigated with the
recirculating, perfused rat liver preparation. Use of dual
radiolabeling of the precursor-product pair allowed for full
characterization of the differential metabolism of
[3H]E1S and
[14C]E1. The strategy was suitable for
investigating the effects of vascular and tissue binding, RBC
metabolism of E1, transport, metabolism, and the various
zonal aspects on the futile cycling between E1S and
E1 in the liver, especially when the RBC distribution and
metabolism of E1 were fully characterized. Finally, a
series-compartment liver model that embodied zonal and subcellular
distribution of metabolic enzymes was developed to interpret the
perfusion results.
 |
Experimental Procedures |
Materials.
[6,7-3H]E1S (ammonium
salt, specific activity, 53 Ci/mmol),
[6,7-3H]E1 (specific activity, 40.6 Ci/mmol),
and [4-14C]E1 (specific activity, 56.6 Ci/mol) were purchased from PerkinElmer Life Sciences (Boston, MA). All
radiochemical purities found by high performance liquid chromatography
(HPLC) or thin-layer chromatography (TLC) were greater than 95%.
E1S, E1, E2, and bovine serum
albumin (BSA) were obtained from Sigma Chemical Co. (St. Louis, MO).
All other reagents were of the highest grade available.
Protein Assay and Hematocrit Count.
In all preparations,
protein was determined by the method of Lowry et al.
(1951)
, with bovine serum albumin as the standard. The
hematocrit was measured by capillary centrifugation in a
microhematocrit centrifuge (IEC MB Centrifuge, Damon, Fisher
Scientific, Mississauga, ON, Canada).
BSA Binding of tracer E1, E2, and
E1S.
BSA binding of E1, E2,
and E1S was investigated using a commercially available
ultrafiltration kit (Centricon 3; Amicon Inc., Beverly, MA). Tracer
[3H]E1 (3.7 ± 0.1 × 106 dpm/ml), [3H]E1S (1.7 ± 0.1 × 106 dpm/ml), or
[3H]E2 (1.1 ± 0.1 × 106 dpm/ml) was added to 4% BSA (v/v) in
Krebs-Henseleit-bicarbonate buffer (pH 7.4). After incubating the
mixture for 10 min at 37°C, an aliquot (2 ml) was removed into a
Centricon tube and centrifuged at 2500g for 20 min. The
radioactivities in the original mixture (0.2 ml) and the resulting
ultrafiltrate (0.2 ml) were quantified by liquid scintillation
spectrometry (model LS6800; Beckman Instruments Canada, Mississauga,
ON). Leakage of BSA into the ultrafiltrate was less than 1% of the
original protein concentration.
Distribution and Metabolism of E1 in
Erythrocytes.
Bovine erythrocytes (a generous gift from Ryding
Regency Meat Packers Ltd., Toronto, ON) were washed three times with
saline and twice with lactated Ringer's solution (Baxter Corporation, Toronto, ON). The distribution and metabolism of tracer
[3H]E1 (2.5 ± 0.4 × 105 dpm/ml or 27 ± 4.2 nM) were studied with
perfusion media of different compositions: 20 and 60% RBC, in the
absence and presence of 4% BSA. Erythrocytes (20 or 60% v/v) and
plasma (Krebs-Henseleit-bicarbonate buffer at pH 7.4 containing 0 or
4% BSA and [3H]E1) were mixed and incubated
under oxygenation (carbogen, 95% oxygen and 5% carbon dioxide; Canox
Gas, Mississauga, ON) at 37°C in the reservoir of the commercially
available TWO-TEN Perfuser. Blood samples were taken at 1, 30, 60, 120, and 180 min, and the hematocrit was measured. The 20 and 60% RBC
yielded hematocrits (HCTs) of 0.15 ± 0.01 and 0.5 ± 0.03, respectively.
Recirculating Rat Liver Perfusion.
Male Sprague-Dawley rats
(290-330 g; Charles River Canada, St. Constant, QC, Canada), which
were fed ad libitum, were used for perfusion at 10 ml/min. The
temperature of the liver was maintained at 37°C with a heating lamp.
Surgery was performed under pentobarbital anesthesia (50 mg/kg,
intraperitoneal injection), and the surgical procedure and the
perfusion apparatus were identical to those described by
deLannoy et al. (1993)
. The perfusion medium consisted of 20% washed bovine erythrocytes, 4% BSA, and 300 mg/dl glucose (50% dextrose injection USP; Travenol Canada, Mississauga, ON) in
Krebs-Henseleit-bicarbonate solution (pH 7.4). The rat liver was
recirculated with blank medium for 20 min during the equilibration period, followed by perfusion with medium containing
[3H]E1S (initial concentration of 2.85 ± 0.23 × 105 dpm/ml or 2.4 ± 0.20 nM) and
[14C]E1 (initial concentration of 1.04 ± 0.12 × 105 dpm/ml or 848 ± 94 nM) from a
second reservoir (200 ml). Reservoir perfusate (1-2 ml) was sampled at
0, 2.5, 10, 30, 60, 90, 110, 130, and 150 min. The total volume removed
from the reservoir was 7% (14 ml) of the initial volume (200 ml), and
no attempt was made to correct for the loss in volume. Bile was
collected at 5- and 10-min intervals so the mid-time of the interval
coincided with the sampling time of the reservoir.
Extraction and TLC Assays of [3H]E1 and
[3H]E2 for the RBC Metabolism Studies.
The blood and its derived plasma obtained by instantaneous
centrifugation (1.5 ml each) were immediately extracted into ethyl acetate (1:2, v/v). One aliquot (1 ml) of the ethyl acetate extract was
directly subjected to liquid scintillation counting, and the total
count of the sample was determined against a calibration curve
constructed of standards containing varying known counts of
[3H]E1 in perfusate and processed in the same
fashion. Since [3H]E1 and
[3H]E2 were completely extracted into ethyl
acetate (>99%), the extraction method furnished a mixture of
[3H]E1 and [3H]E2
in each sample except for time zero, when only
[3H]E1 was present. The ratio of
[3H]E1/[3H]E2 was
further given by TLC described below. A second aliquot (1 ml) of the
ethyl acetate was spotted onto the Silica Gel GF (250 µm) TLC plate
(Analtech, Newark, DE), which had been preloaded with E1
and E2 at the origin to separate
[3H]E1 and [3H]E2.
The plates were developed in a system of toluene:ethanol (9:1, v/v).
Regions for E1 (Rf = 0.76) and
E2 (Rf = 0.57) were visualized
under UV light and scraped into minicounting vials. After the addition
of water (0.5 ml) and liquid scintillation fluor (5 ml, Ready Safe,
Beckman Instruments, Canada) into minicounting vials, the radioactivity
was quantified by liquid scintillation spectrometry (model LS6800;
Beckman Instruments, Palo Alto, CA). Hence, the amounts of
[3H]E1 and [3H]E2
in plasma and blood perfusate were quantified by the combined extraction-TLC method. The amounts of E1 and E2
in RBC were, however, calculated by difference between the quantities
in plasma and blood perfusate of known hematocrit. The radioactivities
were expressed as a percentage of the initial concentration of
[3H]E1 used.
HPLC Assays for Quantitation of E1G, E1S,
and E1 in the Liver Perfusion Studies.
Acetonitrile,
which contained 4 µM danazol (the internal standard), was used to
terminate any metabolic reactions, with 1:4 (v/v) volume ratio. All
perfusate samples (1-2 ml) were immediately transferred to tubes
containing acetonitrile (4-8 ml). Contents of the deproteinized
samples were dried under nitrogen (Canox Gas) and analyzed by HPLC as
described by Tan and Pang (2001)
. Standards of the
calibration curve prepared with samples containing varying known counts
of [3H]E1S and
[14C]E1 were processed in the same fashion.
Bile samples were diluted 1:1 (v/v) with water, and 20-µl aliquots
were directly counted. A portion of the diluted bile (20 µl) was
subjected to HPLC with internal standardization. The radioactivities in
bile and from HPLC radioelution were quantified by liquid scintillation
spectrometry. Eluted radioactivities of less than three times the
background counts were treated as zeroes. All 3H- and
14C-radioactivities quantified in the samples were higher
than 3000 and 1000 dpm, respectively.
Kinetic Modeling of [3H]E1 Metabolism
in Erythrocytes and Fitting.
Various cellular models were tested
for their abilities to predict the disposition of
[3H]E1 and [3H]E2
in erythrocytes. The cellular kinetic model that included plasma
protein binding and red cell binding and metabolism of E1
(Fig. 1) best described the kinetics of
E1 and E2. Mass balanced rate equations (eqs.
1-4) were written to describe the RBC distribution and metabolism of
E1 and E2 (Fig. 1). Oxidation of E2
to E1 was not included since preliminary study revealed
less than 1% metabolism of E2 to E1 over
3 h. The same was observed by Tsang (1976)
.

View larger version (17K):
[in this window]
[in a new window]
|
Fig. 1.
The cellular kinetic model for the description of
E1 and E2 in blood perfusate. Plasma protein
and red cell binding of E1 and E2 and
metabolism of E1 were included in the model. The transport
of E1 and E2 was denoted by the bidirectional
linear transmembrane clearances and the metabolism of E1
was described by the RBC metabolic intrinsic clearance.
|
|
Binding to plasma and RBC is expressed as unbound fractions. The
unbound fraction in blood (fblood) is related to
that in plasma (fp) and the plasma
(CP) and blood (Cblood)
concentrations (Pang and Rowland, 1977
).
|
(1)
|
For species such as E1 and E2 but not
E1S that distribute into RBC, we assumed that protein and
red cell binding and debinding of E1 and E2
occurred almost instantaneously. The binding to plasma albumin results
in plasma unbound fractions of E1
(f
) and
E2 (f
) whereas their distributions in erythrocytes (RBC) yielded unbound fractions of E1
(f
) and
E2
(f
). For these
species, the unbound fraction in blood is related to the HCT and the
unbound fractions in plasma and red blood cell (frbc).
|
(2)
|
For estimation of the bound concentrations of E1 in
plasma ([E1,bound]p), we assume that
K
(the binding dissociation
constant of E1) > [E1,unbound]p, and this exists for the tracer
condition (Tan and Pang, 2001
). The Langmuir binding
isotherm simplifies to
|
(3)
|
where
(nE1[Ptotal]p)
is effective binding concentration, with
nE1 being the number of the binding sites
on BSA and [Ptotal]p, the total
protein concentration in plasma.
The unbound fraction of E1 in plasma
(f
) is derived from
substitution of eq. 3 into eq. 4
|
(4)
|
|
(5)
|
The fixed protein concentration (e.g., 4% BSA) in blood
perfusate in the incubation study will change when the BSA
concentration is expressed in relation to the plasma when the
hematocrit (0, 20, and 60% RBC) is modified. To relate to the in vitro
binding data of 4% BSA in plasma to those for the red cell-albumin
binding study, the following correction was made. The plasma protein
concentration ([Ptotal]p) is
related to the protein concentration in blood perfusate, [Ptotal]perfusate, as given below.
|
(6)
|
Analogously, the corrected unbound fraction of E1 in
perfusate plasma is related to the plasma unbound fraction assessed in
vitro (obtained from substitution of eq. 6 into eq. 5)
|
(7)
|
and the corrected, unbound fraction of E2 in
perfusate plasma was described by Eq. 8.
|
(8)
|
For the 20% RBC and 60% RBC, albumin-free perfusate, the
unbound plasma fractions of E1 and E2 equal
unity. In the presence of 4% BSA, the plasma unbound fractions of
E1 and E2 are expressed in relation to the
corrected, unbound plasma fractions of E1 and E2
(f
and f
) as described by eqs. 7 and 8, respectively.
Clearance terms were normalized to the hematocrit for purposes of
comparison since different hematocrits were used for study. Analogously, the transmembrane clearances of E1 and
E2 described by Koefoed and Brahm (1994)
were also normalized to the hematocrit to provide

and

, respectively.
The hematocrit normalized intrinsic clearance of E1
(
) is expressed as follows.
|
(9)
|
The hematocrit normalized permeation clearances of
E1 and E2 across the red cell membrane are
|
(10)
|
|
(11)
|
The hematocrit normalized clearances were multiplied back to the
hematocrit to yield the corresponding diffusion and intrinsic clearance
for perfusates of varying composition, namely 20 and 60% RBC (see eq.
12 through eq. 15). The equations that describe the rates of change of
E1 and E2 in the plasma (P) are
|
(12)
|
|
(13)
|
The equations that describe the changes of E1 and
E2 in the RBC space (rbc) are
|
(14)
|
|
(15)
|
The RBC unbound fractions of E1
(f
) and
E2 (f
) and the metabolic intrinsic clearances were estimated by the
least-squares fitting procedure (SCIENTIST version 2; MicroMath
Scientific Software, Salt Lake City, UT) with the weighting schemes of
unity. The goodness of fit was viewed with respect to the coefficient
of variation (standard deviation of parameter estimate/parameter
value), the residual plot and the model selection criterion (MSC).
Kinetic Modeling of E1S and E1
Disposition in the Recirculating Rat Liver Preparation.
A
series-compartment, liver model containing two units representing the
periportal (PP) and perivenous (PV) regions of the liver, is the
minimalized model for purposes of fitting that best predicted the
disposition of E1S and E1 in the recirculating
liver preparation (Fig. 2). In this
model, a reservoir compartment was included for recirculation of
the perfusate. The flow of substrates occurs unidirectionally from the
periportal to the perivenous region, and exchange occurs in the
tranverse and not the longitudinal direction. Linear (nonsaturable)
transport and metabolic intrinsic clearances prevail in view of the
tracer condition studied, and the assumption is justified based on the
observed Km values (in µM) (Tan and
Pang, 2001
) being in excess of the tracer concentrations (in
nM) studied. Species such as E1S, E1, and
E1G that were quantified were modeled. Other metabolites
formed from E1 and E1S (E2 and estriol [E3] and their glucuronide and sulfate conjugates
such as E2S, E2-3S-17G, E3S, and
E3-3S-16G) were collectively represented by M'.

View larger version (25K):
[in this window]
[in a new window]
|
Fig. 2.
Schematic representation of the liver by a
series-compartment, liver model that embodied acinar and
subcompartmentalization and subcellular distribution of metabolic
enzymes. The bottom chart conveyed the ascinar distributions of
enzymatic activities of estrone sulfotransferase, estrone
UDP-glucuronosyltransferase, and estrone sulfatase for metabolism of
E1S and E1, expressed as percentages of the
total metabolic intrinsic clearances.
|
|
A new feature of the extended model was the addition of an endoplasmic
reticulum compartment, as proposed by Tirona and Pang (1996)
. This was necessary since the elimination profiles of
E1 differed subsequent to the administration of tracer
[3H]E1S and
[14C]E1. The added compartment segregates the
cytosol from the endoplasmic reticulum where microsomal enzymes are
found. Estrone sulfotransferase is placed in the cytosolic compartment,
whereas estrone sulfatase and UDP-glucuronosyltransferase are placed in
the endoplasmic reticulum compartment. In light of the known, zonal
distributions of estrone sulfotransferase and estrone sulfatase
(Tan and Pang, 2001
) and of UDP-glucuronosyltransferase
(Tosh and Burchill, 1996
) activities in the liver, their
enriched zonal metabolic activities were calculated from the previously
obtained in vitro data (Fig. 2, bottom panel; see description to
follow), as described in previous studies (Abu-Zahra and Pang,
2000
).
The assigned volume of sinusoid (Vs), cytosol
(Vc), endoplasmic reticulum
(Ver), and biliary compartment
(Vbile) were 1.4 ml (Schwab et al.,
1990
), 7.3 ml (Pang et al., 1988
), 0.2 ml
(Tirona et al., 1996
), and 0.07 ml, respectively. The
volume of the biliary compartment was the summation of the biliary
volume (0.044 ml; Reichen and Paumgartner, 1980
) and the
void volume (about 0.026 ml) in the bile-duct cannula (PE50, Becton
Dickinson, Sparks, MD). The apparent biliary excretion clearance of
E1S or E1G was calculated as the biliary
excretion rate divided by the midpoint reservoir concentration of each
respective species.
Fitting of Data to the Series-Compartment Liver Model.
Mass balanced rate equations (see Appendix) were written to
describe events of the series-compartment liver model (Fig. 2). The
amounts of drug and metabolite in both perfusate and bile were
normalized by the dose. Binding to red cell and albumin was assumed to
be rapidly equilibrative such that use of on- and off-rate constants
was not necessary. Under this instance, the unbound concentrations of
E1S and E1 in whole blood perfusate equal those in plasma and in RBC. The unbound fraction in blood may then be calculated from either eq. 1 or eq. 2. The unbound fractions of E1S
(f
) and
E1 (f
) in the liver cytosol were taken from Tan and Pang
(2001)
.
The clearance of E1 in erythrocytes
(CL
) was determined as
dose/area from the in vitro RBC metabolism study. The in vitro
Vmax (nmol/min/mg of S9 or cytosolic protein)
and the intrinsic clearance (ml/min/mg of S9 or cytosolic protein) for
liver metabolism were scaled up with factors,
and
(
= 0.8 mg of S9 protein/106 cells;
= 0.5 mg of
cytosolic protein/106 cells, with 125 × 106 cells/g of liver; Mahler and Cordes,
1966
; Lin et al., 1980
). The
Vmax values for transport (nmol/min/mg of
protein) of E1S and E1 obtained from Tan
and Pang (2001)
were scaled-up with the factor, 200 mg of
protein/g of liver (Mahler and Cordes, 1966
; Lin
et al., 1980
). The uptake clearance was expressed as the sum of
all of the Vmax/Km for
the saturable transport processes and Pdiff, the
linear transport clearance; ml/min/g of liver) and the metabolic
intrinsic clearance was
Vmax/Km (ml/min/g of
liver) in view of the tracer conditions employed; the transport
clearances were perceived as bidirectional. The various binding,
metabolic, and transport parameters obtained experimentally that are
necessary for fitting were assigned as constants (Table 6).
Fitting was performed by a software package SCIENTIST (version 2;
MicroMath Scientific Software). Transport parameters
the sinusoidal
bidirectional transmembrane clearance of E1G
(CL
), the endoplasmic reticulum
influx (CL
) and efflux
(CL
) clearances of
E1, the uptake clearances of E1S
(CL
) and E1G (CL
) for the endoplasmic
reticulum compartment were obtained by the fitting procedure.
Similarly, the biliary intrinsic clearances of E1S
(CL
) and E1G
(CL
), the sulfation
(CL
) and
glucuronidation
(CL
), and formation intrinsic clearance for the pooled metabolites of
E1 (CL
), and the
desulfation (CL
) and
formation intrinsic clearance of other metabolites of E1S (CL
) were optimized by least-square fitting. Data for each experiment for
[3H]E1S and its metabolites,
[3H]E1 and [3H]E1G,
and for [14C]E1 and its metabolites,
[14C]E1S and
[14C]E1G, given as bolus doses into the
reservoir, were fitted simultaneously to the series-compartment, liver
model (Fig. 2). The means and standard deviation of the parameter
estimates of four experiments (n = 4) are summarized
(Table 6). Appropriate weighting schemes of 1/observation (for data of
higher values) and 1/observation2 (for data of lower
values) were used. The goodness of fit was viewed with respect to the
coefficient of variation (standard deviation of parameter
estimate/parameter value), the residual plot and the MSC.
Statistical Analysis.
All data were presented as the
mean ± standard deviation, and the means were compared by use of
ANOVA or the paired t test, with the level of significance
set at 0.05. The MSC and the Akaike Information Criteria
(Akaike, 1974
; Ludden et al., 1994
) were used to select the appropriate model(s).
 |
Results |
Plasma Binding of E1S, E1, and
E2.
The unbound fraction of E1S in 4% BSA
plasma was 0.03 ± 0.01 (n = 3), whereas those for
E1 and E2 were 0.05 ± 0.01 and 0.04 ± 0.01, respectively (see Table 1,
n = 3). The unbound fractions of E1 and
E2 in perfusate of different compositions as determined by
eq. 1 are summarized in Table 1.
Incubation of a Tracer Dose of [3H]E1 in
Blood Perfusate.
The time courses for
[3H]E1 and [3H]E2
in erythrocytes are shown in Fig. 3.
Different areas under the concentration-time curves (AUC) for
[3H]E1 were noted in the presence and absence
of BSA, and similar observations were found for
[3H]E2 (Table
2). The RBC clearance of E1
(CL
or 0.035 ± 0.02 ml/min) in 60% RBC blood-perfusate was higher than that of the
20% RBC blood-perfusate (0.0092 ± 0.006 ml/min), and these rates
for E2 formation decreased in the presence of 4% BSA
(Table 2, Fig. 3). Upon normalization to the hematocrit, normalized
values of the RBC clearance of E1
(
) in
20 and 60% RBC perfusates became similar (0.061 ± 0.04 and 0.069 ± 0.04 ml/min/HCT, respectively). These were dramatically reduced to 0.0031 ± 0.001 and 0.0024 ± 0.001 ml/min/HCT,
respectively, in the presence of 4% BSA. The RBC to plasma
partitioning ratios of E1 and E2 in the 20 and
60% RBC albumin-free perfusate were higher than those in the presence
of 4% BSA. In the presence of 4% albumin in perfusate, values of
approximately unity were obtained for the RBC to plasma partitioning
ratio for E1 and E2 (Fig.
4). Moreover, the ratios reached their
equilibrium values almost immediately.

View larger version (25K):
[in this window]
[in a new window]
|
Fig. 3.
Time-dependent profiles for tracer
[3H]E1 incubated with perfusates of various
composition at 37°C: Amounts of E1 in blood perfusate
(A), E1 in plasma (B), E1 in erythrocytes (C),
E2 in blood perfusate (D), E2 in plasma (E),
and E2 in erythrocytes (F). The experiment was conducted
with perfusate of four different compositions: 20% RBC ( ), 60% RBC
( ) at 0% BSA, and 20% RBC ( ) and 60% RBC ( ) in 4% BSA. All
were means ± S.D. of five experiments, and the lines were the
fitted lines based on eqs. 12 to 15 with the optimum weighting scheme
of unity.
|
|
View this table:
[in this window]
[in a new window]
|
TABLE 2
CL and AUC of E1 and E2 in plasma and RBC after
incubation with a tracer concentration of [3H]E1 in
blood perfusates
|
|

View larger version (15K):
[in this window]
[in a new window]
|
Fig. 4.
RBC to plasma partitioning of E1 (A) and
E2 (B) in perfusate of different compositions: 20% RBC
( ) and 60% RBC ( ) in the absence of BSA, and 20% RBC ( ) and
60% RBC ( ) in 4% BSA. All data were means ± S.D. of five
experiments.
|
|
Fitted Results for the Kinetic Model of E1 and
E2 in Erythrocytes.
Upon examination of the composite
data for plasma and RBC, the results showed that E1 and
E2 rapidly reached equilibrium in less than a minute (Fig.
4). The same observation was found by Koefoed and Brahm
(1994)
. The fitted RBC unbound fractions of E1
(f
) and E2, obtained upon the simultaneously fitting of the
composite data, were 0.073 ± 0.032 and 0.10 ± 0.07, respectively, showing that both E1 and E2 were
highly bound to erythrocytes. The fitted RBC metabolic intrinsic
clearance of E1
(
) was 0.11 ± 0.07 ml/min/HCT. Good fits were obtained although high coefficients of variation were found associated with the fitted parameters. The best fit to the model that considered red blood cell
and plasma binding E1 and E2 and metabolism of
E1 is presented in Fig. 3, and the optimized parameters of
five experiments are summarized in Table
3.
View this table:
[in this window]
[in a new window]
|
TABLE 3
Assigned and fitted parameters for the cellular kinetic model that
described the distribution and metabolism of E1 and E2
in blood
|
|
Metabolism of the Tracer Dose of [3H]E1S
in the Perfused Rat Liver Preparation.
From the plasma unbound
fraction (Table 3), values of the unbound fraction of E1S
in the blood perfusate
(f
) were
0.027 ± 0.004 according to eq. 1. The apparent hepatic clearance of [3H]E1S was 5.8 ± 0.9 ml/min in the
recirculating rat liver perfusion (Table
4). Within 150 min of recirculation, a
rapid monoexponential decline of [3H]E1S
(t1/2 = 27 ± 1 min) to around 1% of its
initial concentration was observed (Fig.
5A). The accumulation of
[3H]E1G was higher than that of
[3H]E1 in perfusate within the first hour,
followed by a gradual descent (Table 4). The decay half-life of
[3H]E1 paralleled that of its precursor,
[3H]E1S.
View this table:
[in this window]
[in a new window]
|
TABLE 4
Hepatic CLs and AUCs of [3H]E1S,
[14C]E1, and their metabolites in the recirculating
rat liver preparation
|
|

View larger version (23K):
[in this window]
[in a new window]
|
Fig. 5.
Time-dependent profiles of
[3H]E1S and [14C]E1
in the recirculating rat liver preparation. All data
{[3H]E1S ( ),
[3H]E1 ( ),
[3H]E1G ( ),
[14C]E1S ( ),
[14C]E1 ( ), and
[14C]E1G ( )} were mean ± S.D. of
four experiments. The lines were fitted lines based on average
parameters shown in Table 6, with mass balanced rate equations
described under Appendix and appropriate weighting
schemes.
|
|
During recirculation, the excreted amounts of
[3H]E1S and [3H]E1G
in bile increased with time and reached asymptotic levels at 150 min
(Fig. 5B), and the total amounts of [3H]E1S
and [3H]E1G found in bile were 2.5 ± 0.4 and 6.5 ± 0.6 percent dose, respectively (Table
5). However, very little
[3H]E1 was detected in the bile (below the
detection sensitivity). When the biliary excretion clearances for
[3H]E1S and [3H]E1G
were plotted against time, a time-dependent declining profile was
observed for [3H]E1G (Fig. 6B); the bile flow
declined slightly with perfusion time, as expected of the rat liver
upon depletion of bile salts (Fig. 6A).
The excretion clearance of preformed [3H]E1S
reached an asymptotic level by 150 min after reaching distribution equilibrium in the system. At the end of the experiment, the
radioactivities in reservoir, bile, and liver accounted for 3.5 ± 0.4, 54 ± 3, and 43 ± 6 percent dose, respectively.
View this table:
[in this window]
[in a new window]
|
TABLE 5
Biliary excretion of [3H]E1S,
[14C]E1, and their metabolites during simultaneous
delivery to the recirculating perfused rat liver preparation
|
|

View larger version (14K):
[in this window]
[in a new window]
|
Fig. 6.
A, bile flow rates of the rat liver during the
recirculating liver preparation. B, the apparent biliary excretion
clearances of [3H]E1S ( ),
[3H]E1G ( ),
[14C]E1S ( ), and
[14C]E1G ( ). The apparent biliary
excretion clearances of E1S and E1G were
calculated as the biliary excretion rate divided by the midpoint
reservoir concentration of each respective species. The data were
means ± S.D. of four experiments.
|
|
Metabolism of the Tracer Dose of [14C]E1
in the Perfused Rat Liver Preparation.
Estrone was highly cleared
upon the recirculation of [14C]E1 with an
apparent hepatic clearance of 9.4 ± 2.2 ml/min (Table 4), despite
that the unbound fraction of E1 in the blood
(f
) was very
low. The unbound fraction of E1 in blood (0.036 ± 0.006) based on the plasma unbound fraction and blood/plasma
concentration ratio (eq. 1) agreed well with the value (0.053 ± 0.02) estimated according to eq. 2 that utilized the fitted value for
f
(0.073). The
closeness in the values suggest the soundness in the estimation of
f
. The
concentrations of [14C]E1 declined
monoexponentially with a slightly shortened half-life of about 20 ± 1.6 min (Fig. 5C). The elimination profile of the metabolite,
[14C]E1S, was more prolonged and the
half-life was eventually similar to that of
[3H]E1S and [3H]E1
but failed to decay in unison with [14C]E1.
The accumulation of [14C]E1S and
[14C]E1G in the reservoir was comparable and
increased within the first hour, followed by a gradual decline
thereafter. During the time course of the experiment, the
dose-corrected AUC of [3H]E1G, when
extrapolated to time infinity, was not different from that of
[14C]E1G (P > 0.05). Large
variations were, however, observed (Table 4).
The amounts of [14C]E1S and
[14C]E1G excreted in bile increased with time
and reached asymptotic levels at 150 min (Fig. 5D), yielding 1.7 ± 0.01 and 8.2 ± 0.2 percent dose, respectively (Table 5). These
values were not significantly different from those of the
[3H]E1S dose (ANOVA, P > 0.05). Again, little [14C]E1 was
detected in the bile (below the detection limit). When the biliary
excretion clearances of [14C]E1S and
[14C]E1G were plotted against time,
time-dependent declining excretion clearances were observed for both
[14C]E1S and
[14C]E1G (Fig. 6B). At the end of the
experiment, the radioactivities in reservoir, bile, and liver accounted
for 3.3 ± 0.7, 54 ± 6, and 43 ± 8 percent dose.
Fitted Results for the Kinetic Model of E1 and
E1S in the Perfused Liver Preparation.
Upon
simultaneous fitting of perfusate and bile data consisting of
[3H]E1, [3H]E1S,
[3H]E1G, [14C]E1,
[14C]E1S, and
[14C]E1G in each study for the same liver
preparation, good fits were obtained although high coefficients of
variation were found associated with the fitted parameters. Parameter
unidentifiability existed among the fitted parameters due to the
remoteness of the endoplasmic reticulum with respect to the sampling
compartment and the high correlation among parameters. For example, the
endoplasmic reticulum efflux clearance
(CL
), the glucuronidation
intrinsic clearance
(CL
), and the
"pooled" metabolic intrinsic clearance
(CL
) for E1 were
all highly correlated, as were the bidirectional transmembrane clearance (CL
), the desulfation
intrinsic clearance
(CL
), and the pooled
metabolic intrinsic clearance
(CL
) for E1S.
These highly correlated parameters exhibited coefficients of variation
larger than one, and their reliability was much reduced. The other
parameters exhibited lower coefficients of variation since the
parameters were not correlated.
The optimized fit that considered both zonal and subcellular
localization of metabolic enzymes is presented in Fig. 5, and the
assigned parameters and the mean ± S.D. of the optimized
parameters of four experiments are summarized in Table 6. The fitted
sinusoidal bidirectional transmembrane clearance for E1G
(CL
) was 339 ± 22 ml/min.
The endoplasmic reticulum influx
(CL
) and efflux
(CL
) clearances of
E1 were 86 ± 40 and 17 ± 2 ml/min,
respectively, suggesting a 5-fold partitioning of E1 into
the endoplasmic reticulum compartment. The bidirectional transmembrane
clearances of E1S
(CL
) and E1G
(CL
) for the endoplasmic
reticulum were 742 ± 146 and 0.018 ± 0.001 ml/min,
respectively. The biliary intrinsic clearances of E1S
(CL
) and E1G
(CL
) were 8.0 ± 0.1 and
1.8 ± 0.2 ml/min, respectively. The sulfation
(CL
), the
glucuronidation
(CL
) intrinsic clearances of E1, and the desulfation intrinsic clearance
of E1S (CL
) were 318 ± 38, 105 ± 30, and 332 ± 44 ml/min, respectively. Last, the
pooled metabolic intrinsic clearances for E1
(CL
) and E1S
(CL
) in the formation of
all other metabolites were 255 ± 60 and 214 ± 57 ml/min,
respectively. One should further be aware that these sinusoidal
transmembrane, endoplasmic reticulum transmembrane, metabolic
intrinsic, and biliary intrinsic clearances are highly interrelated,
and the set of values is not unique because other combinations could
possibly be consistent with the data.
Inclusion of the endoplasmic reticulum compartment in modeling
was justified since a high partitioning of E1 into the
endoplasmic reticulum space was observed by Zakim and Vessey
(1977)
. Absence of the endoplasmic reticulum compartment
(achieved with high clearances between the endoplasmic reticulum and
cytosol) furnished an inferior fit, predicting a much higher formation
of [3H]E1 (see simulation in Fig.
7).

View larger version (17K):
[in this window]
[in a new window]
|
Fig. 7.
Simulated profiles of
[3H]E1S, [14C]E1,
and their metabolites in the recirculating rat liver preparation based
on rapid equilibration of all species between the cytosolic and
endoplasmic reticulum compartments. A, time-dependent profiles for
[3H]E1S ( ),
[3H]E1 (- - -), and
[3H]E1G (......). B, time-dependent
profiles for [14C]E1S ( ),
[14C]E1 (- - -), and
[14C]E1G (......).
|
|
Further Simulation for Understanding the Futile Cycling Kinetics of
E1 and E1S in the Perfused Liver
Preparation.
Simulations were further performed based on the
fitted and assigned parameters shown in Table
6. If rapid equilibration of the species
existed between the cytosolic and endoplasmic reticulum compartments
(high transport clearance of 1000 ml/min for E1), similar
elimination profile for E1S and E1 would result
pursuant to [3H]E1S and
[14C]E1 dosing as expected of futile cycling
(Fig. 7); values for the transport clearances of E1S and
E1G, when increased to 1000 ml/min, failed to further
affect the shapes of the curves. With high exchange of E1
between the cytosolic and endoplasmic reticulum compartments
(inter-compartmental clearance of 1000 ml/min), the observed,
discrepant half-lives of E1 resulting from tracer
[14C]E1 dose and not the
[3H]E1S dose now disappeared.
Effects of the Reversible Pathway.
For understanding the
effects of futile cycling on the clearances of estrone and estrone
sulfate, the sulfation intrinsic clearance of E1
(CL
) was set to
zero to eliminate futile cycling of E1, the formed
metabolite. The result was the accumulation of
[3H]E1 upon elimination of the resulfation
pathway after the [3H]E1S dose (Fig.
8A). The profile of
[3H]E1G remained virtually unchanged because
sulfation is a minor pathway. When the desulfation intrinsic clearance
of E1S
(CL
) was set as zero
to prevent futile cycling of E1S as the formed metabolite
of estrone, a greater accumulation of the
[14C]E1S resulted, and formation of
[14C]E1G was reduced after the
[14C]E1 dose (Fig. 8B).

View larger version (15K):
[in this window]
[in a new window]
|
Fig. 8.
Simulated profiles of E1S,
E1, and E1G following the administration of
[3H]E1S in the absence of sulfation, with the
sulfation intrinsic clearance of E1
(CL ) equal to zero
(A) and following the administration of
[14C]E1 with the desulfation intrinsic
clearance of E1S
(CL ) equal to zero
(B), in the recirculating rat liver preparation. The symbols
[3H]E1S ( ),
[3H]E1 ( ), and
[3H]E1G ( ) in (A) are the mean of four
experiments, and the lines [3H]E1S ( ),
[3H]E1 (- - -), and
[3H]E1G (...) represent the simulated
profiles. The symbols [14C]E1S ( ),
[14C]E1 ( ), and
[14C]E1G ( ) in (B) are the mean of four
experiments and the lines [14C]E1S ( ),
[14C]E1 (- - -), and
[14C]E1G (...) represent the simulated
profiles.
|
|
 |
Discussion |
E1 and E2 are potent estrogens
that are bound tightly to albumin and to red blood cells. The
erythrocyte distribution of E1 and E2 and
metabolism of E1 in the presence and absence of 4% BSA
were cha