The pharmacokinetics of recombinant human growth hormone (rhGH) and its
effects on the induction of insulin-like growth factor I (IGF-I) were
studied in juvenile rhesus monkeys. Disposition profiles of rhGH from
two short-term i.v. infusion studies were described by a
two-compartment model yielding a clearance of 16.1 ml/min and
T1/2 of 2.0 h. Four rhGH treatment groups
were included in this study: group A, ProLease rhGH (24 mg), a
sustained-release microsphere formulation; group B, a single s.c.
injection plus an implanted osmotic pump (24.4 mg); group C, a single
s.c. injection (25.9 mg); group D, daily 0.86-mg s.c. injection for 28 days. Their rhGH input profiles were analyzed by a numerical
deconvolution method. ProLease and osmotic pump provided zero-order
inputs of rhGH and maintained the serum rhGH concentrations around 9 to 13 ng/ml for 16 (group A) and 30 days (group B). For s.c. injections, rhGH underwent first-order absorption. An indirect response model was
applied based on use of a Hill function for stimulation of IGF-I
production. Parameter values obtained included
Smax = 2.2, SC50 = 6.5 ng/ml,
and
(slope coefficient) = 6.8, which were applicable to all
treatments. The area under effect curve showed group B to be most
effective for IGF-I induction, whereas group A produced the highest
peak level in 16 days. Group C had the lowest induction among the four
groups, despite being given the highest dose. Group D had modest IGF-I
induction, but the pulsatile rhGH input is less effective than
continuous input provided by ProLease. Our
pharmacokinetic/pharmacodynamic model demonstrates that ProLease and
osmotic pump delivery were best able to maintain rhGH level above the
s.c.50 value, which provided more effective IGF-I induction
compared with the single or daily subcutaneous injections in solution.
 |
Introduction |
Growth
hormone (GH), a single polypeptide chain of 191 amino acid residues
with a molecular mass of 22 kDa, is a somatotropic hormone
secreted from the anterior pituitary gland. The release of GH is
regulated by growth hormone-releasing hormone (GHRH) secreted from the
hypothalamus (Ascoli and Segaloff, 1996
) in an episodic and pulsatile
fashion (Brook and Hindmarsh, 1992
). GH is used in replacement
treatment for growth hormone-deficient children with short stature.
Because polypeptides generally have low bioavailabilities and short
half-lives, GH replacement is normally given by daily s.c. injections
(Albertsson-Wikland et al., 1986
).
The pharmacodynamic responses to GH form two categories. Immediate
responses include production of insulin-like growth factors I and II
(IGF-I and IGF-2), stimulation of triglyceride hydrolysis in adipose
tissues, and stimulation of hepatic glucose output. Intermediate
responses of GH are anabolism and growth promotion mediated by IGF-I,
including chondrogenesis, skeletal growth, and growth of soft tissues
(Ascoli and Segaloff, 1996
). The liver is the major organ producing
IGF-I, and hepatic IGF-I induction is dependent on GH exposure by
increasing transcription of IGF-I mRNA (Mathews et al., 1986
). The
major form of IGF-I in the circulation is a 150-kDa complex associated
with a 40-kDa insulin-like growth factor binding protein 3 (IGFBP-3)
and an 85-kDa acid-labile subunit (Baxter et al., 1989
). This complex
circulating in the plasma serves to transport IGF-I from liver into
effect sites for growth promotion. Overall, GH plays a major role in
the regulation of IGF-I, IGFBP-3, and acid-labile subunit (Kupfer et
al., 1993
).
The effectiveness from continuous infusion compared with pulsatile GH
injections remains a major issue in therapy. In hypophysectomized (HX)
rats, the increase of IGF-I in rat skeletal tissues is more effective
by pulsatile i.v. GH infusion than by continuous infusion (Isgaard et
al., 1988
). However, rhGH continuous infusion induced serum IGF-I and
GH binding protein (GHBP) in a dose-dependent manner, whereas single
s.c. daily injections were ineffective in HX rats (Clark et al., 1995
).
The same study showed that twice-daily s.c. injections of rhGH were
slightly more effective than continuous infusion for body weight gains,
whereas single-daily injections were the least effective. In
GH-deficient human adults, serum IGF-I induction was greater by
continuous infusion compared with the same daily dose as eight i.v.
injections every 3.5 h (Laursen et al., 1994
). These results
suggest that a sustained-release formulation of rhGH may achieve its
pharmacological effects as the conventional daily injections by
maintaining the serum hormone concentrations above a certain level over time.
In this report, a sustained-release formulation ProLease rhGH was
studied in juvenile rhesus monkeys. ProLease is composed of
biodegradable microspheres containing the desired bioactive molecule
incorporated into a matrix of
poly-(DL-lactide-co-glycolide). In vitro and in vivo
characteristics of ProLease rhGH have been reported elsewhere (Johnson
et al., 1996
). To test the effectiveness of rhGH delivered by ProLease,
four different treatment groups of monkeys were given similar total
rhGH doses (24.0 to 25.9 mg). They included rhGH in ProLease by
s.c. injection, single s.c. injection plus osmotic pump delivery of
rhGH (to mimic the initial release and continuous delivery of rhGH by
ProLease), single s.c. injection of rhGH, and daily s.c. injections of
rhGH (Johnson et al., 1996
; Lee et al., 1997
).
Pharmacokinetic/pharmacodynamic (PK/PD) modeling was enacted to
describe the pharmacokinetic profiles of rhGH for these four treatment
groups and the relationships between rhGH serum concentrations and
IGF-I induction. The disposition of rhGH in monkeys by two short-term
i.v. infusions of the hormone yielded its primary pharmacokinetic
parameters. Numerical deconvolution method (Cutler, 1978
;
Veng-Pedersen, 1988
) was then applied to assess the input profiles of
the four treatments using i.v. data as the disposition function. These
profiles were used to identify the input functions for the
pharmacokinetic analysis. For IGF-I induction, an indirect response
model (Dayneka et al., 1993
) was applied. The SC50 rhGH
potency factor for IGF-I induction was estimated.
 |
Materials and Methods |
Study Design
For the rhGH disposition study, four juvenile rhesus monkeys
(body weight, b.wt. = 2.3 ± 0.05 kg) were given two different doses of rhGH, 0.10 mg/kg and 1.42 mg/kg, by 5-min short-term infusions. These studies were separated by a washout period of 7 days.
Serum rhGH concentrations were determined at 16 times up to 12.08 h
after each dose. Mean values at each time point from these four monkeys
were used in the pharmacokinetic analysis.
Group A (ProLease rhGH).
One hundred sixty micrograms of
microspheres (24.0 mg of rhGH) were given by s.c. injection to each of
the four monkeys (b.wt. = 3.3 ± 0.2 kg) on day 0. Serum rhGH
concentrations were determined at nine times up to 12 h of the
first day, and on 24 days up to day 75. Serum IGF-I concentrations were
determined at six times up to 12 h of the first day, and on 23 days up to day 60.
Group B (Single s.c. Plus Pump).
A single s.c. injection of
3.6 mg of rhGH was given, and an osmotic pump containing 20.8 mg of
rhGH was surgically implanted in each of four monkeys (b.wt. = 3.2 ± 0.2 kg) on day 0. Pumps were removed on day 30. Serum rhGH and IGF-I
concentrations were determined at essentially the same times as in
group A.
Group C (Single rhGH s.c.).
A single s.c. injection of 25.9 mg of rhGH was given to each of four monkeys (b.wt. = 3.3 ± 0.2 kg) on day 0. Serum rhGH concentrations were determined at the same
time as group B, whereas serum IGF-I concentrations were determined at
fewer times.
Group D (Daily rhGH s.c.).
A single s.c. injection of 0.86 mg of rhGH was given to each of four monkeys (b.wt. = 3.0 ± 0.2 kg) daily from day 0 to day 27 (total of 28 days). Serum rhGH
concentrations were determined at nine times up to 12 h of the
first day; predose and 2 h postdose on days 1, 2, 3, 5, 7, 13, and
19; predose and seven times postdose on day 27 and days 28, 34, 38, 40, 47, 49, 52, 55, and 74. Serum IGF-I concentrations were determined at
seven times up to 12 h of the first day, 2 h postdose on 8 days up to day 27, and then on 6 days up to day 55.
Serum rhGH concentrations were determined by immunoradiometric assay
(Radim Group, Rome, Italy). IGF-I was separated from the binding
proteins by acid/ethanol extraction, and the serum concentrations were
determined by radioimmunoassay (Diagnostic System Laboratories, Inc.,
Webster, TX). Details of the methodology are reported elsewhere
(Johnson et al., 1996
; Lee et al., 1997
).
Pharmacokinetic Analysis
Disposition of rhGH.
The biexponential disposition of rhGH
serum concentrations from two short-term infusion studies can be
described by a two-compartment model (shown in Fig. 1) as follows:
|
(1)
|
|
(2)
|
|
(3)
|
where Ap is the amount of rhGH in the
central compartment, Cp is the serum rhGH
concentration, At is the amount of rhGH in the
peripheral compartment, CL is clearance,
Vc is the central volume of distribution,
k12 and k21 are the
distribution rate constants between the central and peripheral
compartments, Ri(0) is the zero-order infusion
rate into the central compartment, and Base1 is
the baseline serum rhGH concentration. Equations 1 to 3 were fitted
simultaneously for the two dose levels to estimate CL,
Vc, k12, and k21
values by the maximum likelihood method using the ADAPT II program
(D'Argenio and Schumitzky, 1997
). The estimated parameters were then
used as constants in the pharmacokinetic analysis for the four
treatment groups.
The elimination rate constant from the central compartment
(kel) is CL/Vc.
Considering the 5-min i.v. infusion is practically a bolus dose in a
12-h study, two slope coefficients for the biexponential equation
(C = A · e
1·t + B · e
2·t) were
calculated from kel, k12, and
k21. The terminal phase half-life (T1/2) was estimated by 0.693/
2.
Volume of distribution at the steady state (Vss)
was calculated by Vc · (1 + k12/k21). Mean residence time was
estimated by Vss/CL. For the distribution
clearance of rhGH between compartments, CLd = k12 · Vc = k21 · (Vss
Vc).
Deconvolution Analysis for Input Rates.
PCDCON, a
deconvolution program for pharmacokinetic applications (Gillespie,
1991
), was applied. The principles for the numerical deconvolution
method can be found elsewhere (Cutler, 1978
; Veng-Pedersen, 1988
).
Briefly, when the plasma concentrations versus time profiles for an
extravascular administration and an i.v. bolus (dose-normalized as the
unit disposition function) of the drug are known, the rate at which
drug reaches the system circulation (input rate or absorption rate) for
the extravascular administration can be solved numerically by explicit
curve-fitting method encoded in the PCDCON program.
The rhGH concentrations from the high-dose short-term infusion were
used to obtain the unit disposition function. The input functions
representing the absorption profiles of rhGH treatments were then
determined. Two profiles are provided. One is "input rate of rhGH
versus time" from which the absorption profiles for each treatment
can be determined; the other one is "cumulative amount of rhGH versus
time" from which the total amount of rhGH absorbed can be assessed.
This information was applied to the subsequent pharmacokinetic analysis.
Pharmacokinetic Analysis of Treatment Groups.
Two-compartment models with different input functions were applied.
Depending on the absorption profiles obtained from the numerical
deconvolution, zero-order inputs to the central compartment and/or an
absorption compartment with a first-order absorption rate constant into
the central compartment were assigned (as shown in Fig. 1). The details
are as follows.
Group A (ProLease rhGH).
The results for the analysis of
rhGH absorption rates (Fig. 3) showed that the release of rhGH from
ProLease has four stages: approximately 4.5 mg/day from day 0 to
24 h, 0.22 mg/day from 24 h to day 17, 0.1 mg/day from day 17 to 60, and no rhGH released afterward. The two-compartment model (Fig.
1) was applied as follows:
|
(4)
|
|
(5)
|
|
(6)
|
where R*(0) represents the zero-order input rates,
Cp,A is the serum rhGH concentration, and
Base2 is the baseline value. Four different stages of
inputs were defined as follows: when 0 < t
TRA1, then R*(0) = RA1; when
TRA1 < t
TRA2, then
R*(0) = RA2; when TRA2 < t
TRA3, then R*(0) = RA3; when t > TRA3, then
R*(0) = 0. Equations 4 to 6 were fitted simultaneously to estimate TRA1, TRA2, TRA3,
RA1, RA2, RA3, and
Base2 values by the maximum likelihood method using the
ADAPT II program (D'Argenio and Schumitzky, 1997
).
Group B (Single s.c. Plus Pump).
The results of the
absorption profile (Fig. 4) show that two different rhGH input rates
were involved. In the first 12 h, first-order absorption from the
single s.c. injection was observed. A zero-order input of about 0.3 mg
of rhGH/day delivered by the osmotic pump was seen up to day 30. The
equations used were as follows:
|
(7)
|
|
(8)
|
|
(9)
|
|
(10)
|
where ka,B is the first-order absorption
rate constant, Cp,B is the serum rhGH
concentration, and Base3 is its baseline value. The rhGH
input rate from the pump was defined as follows: when 0 < t
30 days, then R*(0) = RB1; when
t > 30 days, then R*(0) = 0 (pumps were
removed on day 30). The estimates for ka,B,
Aa,B (0) (initial condition of eq. 7; total amount of
rhGH absorbed from the s.c. injection), RB1, and
Base3 were obtained from computer fitting as for group A. The bioavailability of the subcutaneous rhGH injection can be
calculated as Aa,B(0)/injection dose.
Group C (Single rhGH s.c.).
The absorption profile shown in
Fig. 5 indicated that first-order absorption was involved in the single
rhGH s.c. injection. The equations used were as follows:
|
(11)
|
|
(12)
|
|
(13)
|
|
(14)
|
where ka,C is the first-order absorption
rate constant, Cp,C is the serum rhGH
concentration, and Base4 is the baseline value. The
ka,C, Aa,C(0) (initial condition of
eq. 11; total amount of rhGH absorbed from the subcutaneous injection),
and Base4 were estimated by computer fitting as for group
A. The bioavailability of the subcutaneous rhGH injection can be
calculated as Aa,C(0)/injection dose.
Group D (Daily rhGH s.c.).
The rhGH absorption rate profile
for the first injection (Fig. 6) showed that the absorption is a
first-order process. Assuming the pharmacokinetics of rhGH are linear
and stationary, the serum rhGH concentrations with repeated doses were
described as follows:
|
(15)
|
|
(16)
|
|
(17)
|
|
(18)
|
where ka,D is the first-order rate
constant, FD is the bioavailability of rhGH for
each injection, Cp,D is the serum rhGH concentration, and Base5 is its baseline value. The
multiple-dosing schedule was defined in the data file for the ADAPT II
program, and ka,D, FD, and
Base5 were estimated by computer fitting as for group A.
Area Analysis.
The area under the rhGH serum concentration
curve (AUC) from each animal was calculated using the trapezoidal rule.
The AUCtotal, the mean values of AUC from each treatment
group, have been reported elsewhere (Lee et al., 1997
). The
AUCnet, the area under the concentration curve above the
baseline, was calculated by (AUCtotal
AUCbase). Baseline values were obtained from computer
fittings as described previously.
Pharmacodynamic Analysis: Total IGF-I Induction by rhGH
For the data presented herein, IGF-I concentrations were
considered as "total IGF-I" because an acid/ethanol extraction was applied to separate IGF-I from its binding proteins for the assay. We
assumed that rhGH can increase the formation rate of protein-bound IGF-I in the serum, whereas the elimination rate of the complex remains
unchanged. Indirect response model III (Dayneka et al., 1993
) based on
stimulation of the formation rate for IGF-I was applied as follows
(Fig. 1):
|
(19)
|
where IGF1 is the total IGF-I concentration,
CGH is the rhGH serum concentration, which is a
forcing function defined by the pharmacokinetic analysis,
kin is the formation rate of IGF-I, Smax is the maximum stimulation of
kin, SC50 is the rhGH concentration producing 50% of maximum stimulation of the IGF-I formation rate, kout is the elimination rate of IGF-I, and
is the slope coefficient for the Hill function. Assuming that the IGF-I
baseline level is maintained by growth hormone baseline concentration
GH0, and will be restored by a certain time after the
dosing, the IGF-I formation rate was defined as follows:
|
(20)
|
where IGF10 is the initial condition of
eq. 20 (IGF-I predose level) and CGH,0 is the
baseline rhGH concentration. Data from group A (ProLease rhGH) and
group B (single s.c. plus pump) were fitted simultaneously. Considering
that IGF-I predose levels are different between groups A and B, we
assumed that the elimination rate kout is the
same but that formation rates kin differ between these groups. Pharmacodynamic parameters
(kout, Smax,
SC50,
, IGF10) were estimated by
the maximum likelihood method using the ADAPT II program (D'Argenio
and Schumitzky, 1997
).

View larger version (17K):
[in this window]
[in a new window]
|
Fig. 1.
Pharmacokinetic model for rhGH and indirect response
model for IGF-I induction by rhGH. Symbols are defined in the text.
|
|
Simulations for the total IGF-I induction by rhGH for group C (single
rhGH s.c.) and group D (daily rhGH s.c.) were also performed using the
proposed pharmacodynamic model (eqs. 19 and 20) and the estimated
parameters. Values of IGF10 were assigned from
the predose levels. Total area under the effect curve
(AUECtotal) for total IGF-I inductions were calculated for
each group using the trapezoidal rule. The AUECnet, area
under the effect curve above the IGF-I baseline value, was calculated
by AUECtotal
AUECbase. The IGF-I baseline
levels were from the least-squares estimated values (groups A and B) or
assigned predose levels (groups C and D).
 |
Results |
Pharmacokinetics
Disposition of rhGH.
The serum rhGH concentration profiles
from two short-term rhGH infusion studies are shown in Fig.
2. For the high dose, the rhGH
concentrations decline biexponentially with a fairly steep distribution
phase followed by a shallow elimination phase. For the low dose, there
is a steep disposition slope parallel to that of the high dose.
However, the second phase is not very obvious. Both profiles have a
baseline rhGH of about 3 ng/ml. The estimated parameters from the
computer fitting (eqs. 1-3) are listed in Table 1. Our CL, Vc, and
Vss values for rhGH in monkeys are very close to
those of Mordenti et al. (1991)
. The distribution clearance (CLd = 0.35 ml/min) is much lower than the
elimination clearance (CL = 16.1 ml/min), indicating
that the steep distribution phase is due to rapid elimination rather
than distribution from plasma.

View larger version (19K):
[in this window]
[in a new window]
|
Fig. 2.
Serum concentrations of rhGH after short-term i.v.
infusions of 1.42 mg/kg ( ) and 0.10 mg/kg ( ) of rhGH. Data points
are mean values (±S.D.) from four monkeys. The solid lines are the
fittings for eqs. 1 to 3.
|
|
Four rhGH Treatments
ProLease rhGH Dose.
The rhGH concentration versus time and the
rhGH absorption profiles are shown in Fig.
3. ProLease delivered rhGH by three different rates and stopped around 60 days postdosing. The cumulative amount of rhGH absorbed at day 60 is about 12 mg, indicating
bioavailability of approximately 50%. That is close to results from
AUC analysis where the overall bioavailability is 0.47 and the net
bioavailability is 0.51 (Table 3). The parameters estimated by computer
fitting (eqs. 4-6) are listed in Table
2. During days 0 to 1.3, rhGH was released at a high rate of 4.077 mg/day. During days 1.3 to 16, a
moderate rhGH release rate of 0.184 mg/day kept the serum rhGH concentration at a steady state around 9.2 ng/ml. Finally, during days
16 to 62, a low rhGH release rate of 0.063 mg/day maintained the serum
rhGH concentration at a steady state about 4.0 ng/ml. After day 62, the
rhGH concentration dropped to its baseline value (estimated as 1.3 ng/ml). These results were approximated by the model based on the
assumption that rhGH was absorbed by three different constant rates
over three different periods. By calculating the sum of (input
rate · input duration) for these three zero-order inputs, the
total amount of rhGH absorbed is about 10.9 mg (bioavailability = 10.9 mg/24 mg = 0.45). It is consistent with the results from the
AUC and deconvolution analyses.

View larger version (22K):
[in this window]
[in a new window]
|
Fig. 3.
Pharmacokinetic profiles for group A (ProLease rhGH).
a, rhGH absorption rate versus time. b, cumulative amount of rhGH
absorbed versus time (obtained by the numerical deconvolution). c,
serum rhGH concentrations (inset, early times) after s.c. injection of
160 mg of ProLease (24.0 mg of rhGH). Data points are mean values
(±S.D.) from four monkeys. Solid line is the fitting for eqs. 4 to
6.
|
|
Single s.c. Dose Plus Pump.
The rhGH concentration
versus time and the rhGH absorption profiles are shown in Fig.
4. The rhGH s.c. injection on day 0 had a
monoexponential decline, indicating that first-order absorption was
involved. Subsequently, the implanted pump delivered rhGH at a constant
rate until the pump was surgically removed (day 30). The parameters
estimated by computer fitting (eqs. 7-10) are listed in Table 2. The
first-order absorption rate constant for the s.c. injection is 7.75 days
1. The estimated release rate from the pump is 0.295 mg/day, which maintained the serum rhGH concentration at a steady state
of about 13.7 ng/ml for 30 days. After that, the rhGH concentration
declined to its baseline value (estimated as 0.97 ng/ml). The
bioavailability of the single s.c. dose can be estimated by
Aa,B/Injection dose = 2.433/3.6 = 0.68. The overall bioavailability of single s.c. plus pump rhGH
administration estimated by the AUC method is 0.50, and the net
bioavailability is 0.49 (Table 3).

View larger version (22K):
[in this window]
[in a new window]
|
Fig. 4.
Pharmacokinetic profiles for group B (single s.c.
plus pump). a, rhGH absorption rate versus time. b, cumulative amount
of rhGH absorbed versus time profile (obtained by the numerical
deconvolution). c, serum rhGH concentrations (inset, early times) for
s.c. injection of rhGH (3.6 mg) and implanted osmotic pump (20.8 mg of
rhGH) for 30 days. Data points are mean values (±S.D.) from four
monkeys. Solid line is the fitting for eqs. 7 to 10.
|
|
Single rhGH s.c. Dose.
The rhGH concentration versus time and
the rhGH absorption profiles are shown in Fig.
5. The absorption rate profile showed a
monoexponential decline on the first day, and the cumulative amount of
rhGH absorbed versus time curve exhibited a typical first-order
absorption process. The estimated parameters by computer fitting (eqs.
11-14) are listed in Table 2. After the s.c. dose, the peak
concentration (Cmax) of about 2700 ng/ml
appeared shortly after the dosing (Tmax is about
1 h). The concentration returned to its baseline level (1.8 ng/ml)
in 2 days. Bioavailability of rhGH can be estimated by
Aa,C/Injection dose = 15.07/25.9 = 0.58. The overall and net bioavailabilities by the AUC method are 0.63 and 0.71 (Table 3).

View larger version (21K):
[in this window]
[in a new window]
|
Fig. 5.
Pharmacokinetic profile for group C (single rhGH
s.c.). a, rhGH absorption rate versus time. b, cumulative amount of
rhGH absorbed versus time profile (obtained by the numerical
deconvolution). c, serum rhGH concentrations (inset, early times) after
subcutaneous injection of rhGH (25.9 mg). Data points are mean values
(±S.D.) from four monkeys. Solid line is the fitting for eqs.
11-14.
|
|
Daily rhGH s.c. Dosing.
The serum rhGH concentration versus
time and the rhGH absorption profiles are shown in Fig.
6. For the first dose, the results showed
that first-order absorption was involved. Estimations of parameters by
computer fitting (eqs. 15-18) are listed in Table 2. Each s.c. dose
had a peak concentration of about 200 ng/ml at 1 h postdose. The
accumulation effect of multiple dosing is not seen. After the last dose
on day 27, rhGH concentration declined to its baseline value estimated
as 2.52 ng/ml. The estimated bioavailability for each dose from
computer fitting is 0.86, which is very close to the results from the
AUC analysis done for the first dose, where overall and net
bioavailabilities are 0.85 and 0.88 (Table 3).

View larger version (22K):
[in this window]
[in a new window]
|
Fig. 6.
Pharmacokinetic profile for group D (daily rhGH
s.c.). a, rhGH absorption rate versus time. b, cumulative amount of
rhGH absorbed versus time profile (obtained by the numerical
deconvolution for the first dose). c, serum rhGH concentrations (inset,
early times) after daily s.c. injection of rhGH (0.86 mg/day). Data
points are mean values (±S.D.) from four monkeys. Solid line is the
fitting for eqs. 15 to 18.
|
|
Pharmacodynamics
The IGF-I concentrations versus time profiles for group A
(ProLease rhGH) and group B (single s.c. plus pump) are shown in Fig.
7. The parameters estimated for eqs. 19
and 20 by computer fitting are listed in Table
4. For group A, IGF-I has a baseline value estimated as 310.5 ng/ml. After ProLease rhGH was given, IGF-I
gradually increased to 800 ng/ml at about 20 days after dosing and
returned to its predose baseline slowly by about 40 days. For group B,
the baseline for IGF-I is estimated as 245.4 ng/ml. The IGF-I induction
has a maximum value of 750 ng/ml at about 30 days after dosing and
returned to the baseline by about 40 days. The estimated
SC50 value is 6.8 ng/ml. When the serum rhGH concentration
is higher than this level, growth hormone can increase IGF-I formation
rate by as much as 3.2-fold (estimated Smax is
2.2). The elimination rate constant for total IGF-I
(kout) is 0.18 day
1, which means
the half-life of total IGF-I is about 4 days. This value should be much
longer than the half-life of free IGF-I (Guler et al., 1989
; Lieberman
et al., 1992
).

View larger version (19K):
[in this window]
[in a new window]
|
Fig. 7.
Total IGF-I induction by rhGH for group A (ProLease
rhGH) and group B (single s.c. plus pump). Data points are mean values
(±S.D.) from four monkeys in each group. Solid lines are the fittings
for eqs. 19 and 20.
|
|
Simulations for groups C and D were performed using eqs. 19 and 20 and
the predicted parameter values. The superposition of our data and the
simulated results shown in Fig. 8
indicate that the proposed PK/PD model can describe and predict IGF-I
induction under different rhGH-dosing regimens. For group C, the IGF-I
has a low baseline (190.8 ng/ml). The peak IGF-I level appeared at 1 to
2 days after dosing and returned to baseline in 5 days. For group D,
the accumulation of IGF-I reached its steady state (peak concentration
350 ng/ml) about 20 days postdosing. Daily
injection was stopped at day 27; the IGF-I then declined to its predose level.

View larger version (18K):
[in this window]
[in a new window]
|
Fig. 8.
Total IGF-I induction by rhGH for group C (single
rhGH s.c.) and group D (daily rhGH s.c.). Data points are the mean
values (±S.D.) from four monkeys in each group. Solid lines are
simulation values for eqs. 19 and 20 using parameter values listed in
Table 4.
|
|
The AUEC for total IGF-I induction are listed in Table
5. Group B produced the largest net AUEC
value. Group C has the lowest net AUEC value despite being given the
highest rhGH dose. Comparing ProLease versus pulsatile injections,
results for group A show that monkeys given ProLease rhGH can induce
IGF-I more effectively than daily solution injection (group D).
To demonstrate the comparative relationships between rhGH serum
concentrations and IGF-I induction, simulations for the four rhGH
treatment groups based on the proposed PK/PD model are presented in
Figs. 9 and
10. Groups A and B had similar rhGH
input schemes in which rhGH was delivered into the system by continuous
zero-order inputs. As a result, serum rhGH concentrations were
maintained above the SC50 value much longer than the single
(group C) or daily (group D) SC injections. ProLease was
able to produce a surge of rhGH concentration for about 1.3 days and
maintain the level at another steady-state concentration for 16 days.
As shown in Fig. 10, IGF-I induction is very effective in that period
of time as the rhGH concentrations were higher than the estimated SC50 value. For group B, a peak rhGH concentration is
produced by the single s.c. injection followed by a steady-state
concentration maintained by the osmotic pump. For 30 days, the rhGH
level was kept above the SC50 value producing appreciable
IGF-I induction. For group C, although the large single dose of rhGH
produced the highest peak concentration among the four groups, the rhGH
level fell below the SC50 quickly by 2 days after dosing
(Fig. 10). Lowest IGF-I induction was seen in this group (Fig. 10). For
group D, daily injections provided pulsatile growth hormone inputs, but rhGH serum concentrations were below SC50 value for about
one-third of the day during the daily injection period (Fig. 9). As a
result, modest IGF-I induction was observed (Fig. 10). Although each
Cmax produced by the daily injection in group D
is lower than the Cmax produced by a single
injection in group C, group D had a greater IGF-I induction by the end
of the treatment due to its accumulative effects.

View larger version (23K):
[in this window]
[in a new window]
|
Fig. 9.
Simulations of rhGH serum concentrations for the
four rhGH treatments. The dashed lines are the rhGH SC50
value (ng/ml) for IGF-I induction estimated by the indirect response
model (Fig. 1).
|
|

View larger version (19K):
[in this window]
[in a new window]
|
Fig. 10.
Simulations of IGF-I induction by four different
rhGH treatments using the indirect response model (Fig. 1).
|
|
 |
Discussion |
The pharmacokinetics of rhGH was assessed in monkeys following two
5-min short-term infusions. A two-compartment model describes the fast
elimination (CL = 16.1 ml/min) and slow distribution (CLd = 0.35 ml/min) of this hormone. The
pharmacokinetic parameters are similar to values reported by Mordenti
et al. (1991)
who found CL = 14.7 ml/min,
Vc = 199 ml, and Vss = 314 ml. The numerical deconvolution method (Cutler, 1978
;
Veng-Pedersen, 1988
) then allows us to study the input rates for
different rhGH delivery systems. Once the input profiles were
identified, input parameters (such as absorption rate constant, input
rate, and input duration) can be estimated by computer fittings using
fixed pharmacokinetic parameters determined in the disposition study.
ProLease provides three different rhGH zero-order deliver rates in
monkeys in three periods: an initial surge (4.077 days) for 1.3 days as
a loading dose, followed by a moderate rate (0.184 mg/day) for 16 days,
and then a low rate (0.063 mg/day) for up to 62 days. This third phase
of low rhGH input rate may be due to a slow release of the hormone from
the microspheres or to the cross-reactivity with endogenous monkey GH
(Lee et al., 1997
). Because these estimates are model-dependent,
different release rates over different periods for ProLease rhGH would
be obtained if more (or fewer) phases for the rhGH release were assumed
for the analysis. Our results show that this sustained-release
formulation can maintain rhGH concentration above the SC50
value (6.8 ng/ml) for IGF-I induction for 16 days postdosing.
Therefore, the AUEC of IGF-I induction for ProLease rhGH is much higher
than from single and daily s.c. injections. Similar to ProLease, the
osmotic pump delivers rhGH as a zero-order input and induces IGF-I
effectively until removal of the pump (30 days).
Two bioavailability values for each rhGH treatment group were
calculated by the AUC method (Table 3). The overall bioavailability (F), calculated by the AUCtotal values,
represents the bioavailability of the growth hormone without adjusting
its endogenous baseline level. The net bioavailability
(Fnet) was calculated by the AUCnet values, which were produced exclusively by the hormone treatments. Theoretically, Fnet may better estimate the true
rhGH bioavailability; however, the variability and uncertainties in the
estimations for the endogenous baseline levels may create some modest
errors for Fnet.
The indirect response model applied in the PK/PD analysis provides an
excellent characterization of the diverse treatments with a common set
of parameters. Using the pharmacokinetic profiles of rhGH as forcing
functions in the dynamic analysis, IGF-I level is described as a
reversible induction process. The use of the slope coefficient
and
the definition of the initial conditions of the differential equations
allowed the model to better describe the abrupt nature of the response
curve and reflect the baseline hormone levels for the response. This
approach also allows us to estimate the potency factor
SC50. This factor provides important information to
optimize rhGH delivery systems as it is crucial to maintain rhGH serum
concentration above the SC50 value to induce IGF-I
effectively. We showed that both ProLease rhGH and osmotic pump
delivery were able to attain this goal for a lengthy period. Their
IGF-I inductions were much more effective than single and daily s.c.
injections (as shown in Figs. 9 and 10).
Several studies have tested a similar concept by prolonging serum GH
circulating time. When rhGH was coadministered with rhGH-binding protein (rhGHBP) by s.c. injection in HX rats, the GH serum
concentration profile had a wider peak with a similar maximum
concentration compared with GH alone. As a result, rhGH plus rhGHBP
daily injections for 8 days produced better weight gain than GH alone
(Clark et al., 1996a
). Another study showed that polyethylene
glycol-conjugated GH formulations had longer half-lives compared with
unconjugated GH in HX rats and produced similar weight gains when
polyethylene glycol-GH s.c. injections were given twice in 12 days
compared with daily GH injections (Clark et al., 1996b
). Furthermore,
Fielder et al. (1996)
showed that the coadministration of GH and IGF-I in HX rats was more effective for long-term weight gains than GH or
IGF-I alone. Therefore, the continuous exposure of GH may promote
desirable growth effect resulting from the induction of IGF-I and the
presence of the hormone itself. Although our monkey studies showed
promising results for IGF-I induction (Johnson et al., 1996
; Lee et
al., 1997
), the effectiveness of long-term ProLease GH treatment in
GH-deficient children is now under investigation. The proposed PK/PD
model, which characterizes the biological events for GH disposition and
the subsequent IGF-I induction profile, may help us better understand
the role of this hormone in the future.
The numerical deconvolution method (Cutler, 1978
; Veng-Pedersen, 1988
)
is very useful to identify drug input profiles and estimate the
cumulative amount of the drug in the system. There are some assumptions
and limitations of this method. For example, rhGH kinetics are assumed
to be linear and stationary for the range of rhGH doses given in the
study (Table 3). Low interindividual variation of rhGH kinetics is
assumed to apply the same rhGH disposition function among different
groups of monkeys. To interpret the absorption rate and cumulative
amount of GH absorbed, an approximation process was needed to utilize
the information in the compartmental analysis. For instance, ProLease
rhGH has an irregular input profile in the first 2 days (Fig. 3). A
zero-order input process was assumed, and the model predicted steady
state to be reached in about 12 h. This is consistent with the
release profile obtained from the in vitro study (Johnson et al.,
1996
). However, the model overestimates the rhGH concentration between
days 1 and 1.3. An alternative approach is to use a spline-type
polynomial equation to describe the input rate. However, such a method
cannot predict the rhGH kinetic profile within differential equations.
Our current approach is advantageous in the drug development process
because it allows extrapolation.
Short-term i.v. infusion studies clearly showed that rhGH
concentrations declined biexponentially; however, this was not evident when s.c. rhGH injections were given. This is common for extravascular administration when ka is close to
2 (Gibaldi and Perrier, 1982
). Since the
ka values (Table 2) for rhGH s.c. injections are
lower than the disposition rate constants (Table 1), flip-flop kinetics may occur where terminal slopes are not as steep as the slope in the
short-term infusion studies. Protein binding of growth hormone is
another important issue. The primary source of serum GHBP is an
extracellular domain of growth hormone receptor cleaved from cell
membranes (Leung et al., 1987
). Under physiological conditions in
humans, circulating growth hormone is about 39 to 59% bound to GHBP
(Baumann et al., 1988
). In a rhGH replacement study in GH-deficient
children (Tauber et al., 1993
), continuous infusions (0.1 I.U./kg
daily) for 6 months increased the GHBP level from predose 8.6 ± 3.1% to 22.5 ± 9.8%. The normal range for children is reported
as 24.8 ± 1.7% (Tar et al., 1990
). Meanwhile, total IGF-I was
increased from 134 to 329 ng/ml in the same study. As demonstrated by
Clark et al. (1996a)
, GHBP alters the pharmacokinetics of rhGH, and its
dynamic responses need further investigation. Monkey GH and native
human GH differ by four amino acids (Li et al., 1986
); therefore, it is
possible that the tested monkeys will develop antibodies after rhGH is
given. Anti-hGH antibodies were detected in only one of four monkeys
given ProLease rhGH, and the results suggested that rhGH released from
ProLease is no more immunogenic than unformulated protein (Johnson et
al., 1996
; Lee et al., 1997
).
In conclusion, our proposed PK/PD model can describe rhGH
pharmacokinetics and total IGF-I induction in monkeys under diverse dosing conditions. ProLease is an effective drug delivery system for
continuous growth hormone replacement therapy. Continuous rhGH input
delivered by ProLease and osmotic pump is more effective than pulsatile
injections in induction of IGF-I, especially when the delivery systems
can maintain growth hormone concentration above the SC50.
We thank Dr. Jogarao Gobburu (SUNY at Buffalo) for his
suggestions for the modeling and data analysis.
Accepted for publication January 28, 1999.
Received for publication November 10, 1998.
rhGH, recombinant human growth hormone;
IGF-I, insulin-like growth factor I;
PK/PD, pharmacokinetics/pharmacodynamics;
AUC, area under the curve;
AUEC, area under the effect curve;
CL, clearance;
Vc, volume of distribution in the
central compartment, Vss, volume of distribution
at the steady state;
ka, first-order absorption
rate constant, k12 and
k21, distribution rate constants;
1 and
2, slope coefficients;
kel, elimination rate constant, F,
bioavailability;
s.c., subcutaneous;
MRT, mean residence time;
kin, IGF-I formation rate;
kout, IGF-I elimination rate;
Smax, maximum stimulation of IGF-I formation
rate;
SC50, rhGH concentration producing 50% of the
maximum effect;
, slope coefficient for the Hill function;
HX, hypophysectomized;
GHBP, growth hormone-binding protein.