Pharmacokinetic-pharmacodynamic modeling of testosterone and luteinizing hormone suppression by cetrorelix in healthy volunteers

Cetrorelix (CET), a potent luteinizing hormone-releasing hormone (LH-RH) antagonist, was recently approved for the prevention of premature ovulation in patients undergoing a controlled ovarian stimulation (COS), followed by oocyte pickup and assisted reproductive techniques (ART), and is currently under clinical trials for benign prostate hyperplasia, endometriosis, and tumors sensitive to sex hormones. CET suppresses luteinizing hormone (LH), follicle-stimulating hormone (FSH), and testosterone (T) in men. The purpose of this study was to evaluate the pharmacokinetics and absolute bioavailability of 3 mg intravenously and subcutaneously administered CET in healthy male and female volunteers and to develop a pharmacokinetic-pharmacodynamic (PK-PD) model to link the plasma concentrations of CET to the T and LH suppression in males. Following intravenous (IV) (n = 5) and subcutaneous (SC) (n = 6) administration of CET acetate, CET and hormone plasma levels were measured by radioimmunoassay (RIA) and enzyme immunoassay (EIA) methods, respectively. Pharmacokinetics of CET was explained by a three-compartment model for the IV route and by a two-compartment model with first-order absorption for the SC route. Average absolute bioavailability after SC administration was 85%. There were no differences in the pharmacokinetics between male and female subjects (ANOVA, p > 0.05). Single IV and SC doses of CET caused immediate and distinct suppression of LH, FSH, and T levels by 80%, 45% and 95% of their baseline levels, respectively. The duration of hormone suppression was longer for the SC route. An indirect-response PK-PD Emax model was developed to link the measured CET plasma concentrations with the respective T or LH levels. In addition, the circadian rhythm of T levels was accounted by including a cosine function in a second separate PD model. The PD model with cosine function was applied to T baseline levels as well as to the suppressed concentrations after CET dosing. The two models adequately described the PK-PD relationship between plasma levels of CET and T suppression following IV and SC dosing. The EC50 values (mean +/- SD) for the suppression of T were similar (p > 0.05) between the two routes of administration and the two models.     

Pharmacokinetic and pharmacodynamic modeling of recombinant human erythropoietin after single and multiple doses in healthy volunteers

This study describes a pharmacokinetic (PK) model to account for serum recombinant human erythropoietin (rHuEpo) concentrations in healthy volunteers following intravenous (IV) and subcutaneous (SC) dosing; it also characterizes the pharmacodynamics (PD) of SC rHuEpo effects on reticulocytes, red blood cells (RBC), and hemoglobin (Hb) in blood. Data were obtained from 4 clinical studies carried out in healthy volunteers. Epoetin alfa (rHuEpo) was administered as 5 single IV doses ranging from 10 to 500 IU/kg, as 8 single SC doses ranging from 300 to 2400 IU/kg, and as 2 multiple SC dosage regimens (150 IU/kg/3 times a week [tiw] and 600 IU/kg/wk). A dual-absorption rate model (fast zero-order and slow first-order inputs) with nonlinear disposition characterized the PK of SC rHuEpo. A high K(m) value was obtained indicating that clearance was mildly nonlinear. Absorption was slow (t(max) approximately 24 hours), and the bioavailability of SC rHuEpo increased with dose (ranging from 46%-100%). A catenary cell production and loss model with a feedback down regulation component was used to fit the reticulocyte data yielding estimates of the stimulatory capacity (S(max)), sensitivity (SC(50)), and life span parameters. These parameters were used for simulations of RBC and Hb profiles. An SC(50) of 27 to 61 IU/L was estimated indicating that low physiological plasma rHuEpo concentrations were sufficient to produce pharmacological effects. No marked sex-dependent differences in clinical responses to rHuEpo therapy were found despite baseline differences. Realistic pharmacokinetic and physiological models accounted for clinical responses from a wide array of dosing conditions with rHuEpo. The rationale for greater efficacy of SC administration of rHuEpo compared to IV was ascertained.     

Receptor-mediated pharmacokinetic/pharmacodynamic model of interferon-beta 1a in humans

Purpose. An integrated receptor-based pharmacokinetic/pharmaco- dynamic (PK/PD) model of interferon- 1a (IFN- 1a) previously developed for monkeys was used to capture the time-course of drug and induced neopterin concentrations after intravenous (IV) and subcutaneous (SC) dosing in humans. Methods. Data were extracted from the literature by digitalization. Single-dose (3 IV doses and 1 SC dose) PK/PD profiles were simultaneously fitted using the basic model and the ADAPT II computer program. Additional submodels incorporating neutralizing antibody formation and negative feedback inhibition were applied to account for drug accumulation and lower than expected neopterin concentrations encountered after multiple-dosing (1 SC dose every 48 hs). Results. The basic model jointly-captured the nonlinear PK behavior of the drug and induced neopterin concentrations after all single doses. Slow and incomplete absorption (F = 0.33) of the SC dose resulted in prolonged drug concentrations reflective of flip-flop kinetics. Despite lower drug concentrations, SC dosing produced a similar neopterin profile as compared with the IV doses; however, with a longer time to peak effect and slightly higher neopterin concentrations at later time points. The PD component of the model represents a modified precursor-dependent indirect response model driven by the amount of internalized drug-receptor complex. The latter stimulated a 6-fold increase in the production of the neopterin precursor (S_max = 5.89). Drug accumulation and lower than expected neopterin concentrations after multiple dosing were also captured after the inclusion of the submodels. Conclusions. The present integrated PK/PD model for IFN- 1a is mechanistic in nature with receptor-mediated disposition and dynamics and was successfully applied to human clinical data.     

Modeling pharmacokinetics/pharmacodynamics of abatacept and disease progression in collagen-induced arthritic rats: a population approach

The PK/PD of abatacept, a selective T cell co-stimulation modulator, was examined in rats with collagen-induced arthritis (CIA) using a nonlinear mixed effect modeling approach. Male Lewis rats underwent collagen induction to produce rheumatoid arthritis. Two single-dose groups received either 10 mg/kg intravenous (IV) or 20 mg/kg subcutaneous (SC) abatacept, and one multiple-dose group received one 20 mg/kg SC abatacept dose and four additional 10 mg/kg SC doses. Effects on disease progression (DIS) were measured by paw swelling. Plasma concentrations of abatacept were assayed by enzyme-linked immunosorbent assay. The PK/PD data were sequentially fitted using NONMEM VI. Goodness-of-fit was assessed by objective functions and visual inspection of diagnostic plots. The PK of abatacept followed a two-compartment model with linear elimination. For SC doses, short-term zero-order absorption was assumed with F = 59.2 %. The disease progression component was an indirect response model with a time-dependent change in paw edema production rate constant (k _in ) that was inhibited by abatacept. Variation in the PK data could be explained by inter-individual variability in clearance and central compartment volume (V ₁ ), while the large variability of the PD data may be the result of paw edema production (k _in ⁰ ) and loss rate constant (k _out ). Abatacept has modest effects on paw swelling in CIA rats. The PK/PD profiles were well described by the proposed model and allowed evaluation of inter-individual variability on drug- and DIS-related parameters.     

Population PK and IgE Pharmacodynamic Analysis of a Fully Human Monoclonal Antibody Against IL4 Receptor

Purpose

For AMG 317, a fully human monoclonal antibody to interleukin receptor IL-4R??, we developed a population pharmacokinetic (PK) model by fitting data from four early phase clinical trials of intravenous and subcutaneous (SC) routes simultaneously, investigated important PK covariates, and explored the relationship between exposure and IgE response.

Methods

Data for 294 subjects and 2183 AMG 317 plasma concentrations from three Phase 1 and 1 Phase 2 studies were analyzed by nonlinear mixed effects modeling using first-order conditional estimation with interaction. The relationship of IgE response with post hoc estimates of exposure generated from the final PK model was explored based on data from asthmatic patients.

Results

The best structural model was a two-compartment quasi-steady-state target-mediated drug disposition model with linear and non-linear clearances. For a typical 80-kg, 40-year subject, linear clearance was 35.0?mL/hr, central and peripheral volumes of distribution were 1.78 and 5.03?L, respectively, and SC bioavailability was 24.3%. Body weight was an important covariate on linear clearance and central volume; age influenced absorption rate. A significant treatment effect was observable between the cumulative AUC and IgE response measured.

Conclusion

The population PK model adequately described AMG 317 PK from IV and SC routes over a 60-fold range of doses with two dosing strengths across multiple studies covering healthy volunteers and patients with mild to severe asthma. IgE response across a range of doses and over the sampling time points was found to be related to cumulative AMG 317 exposure.     

The application of mechanism-based PK/PD modeling in pharmacodynamic-based dose selection of muM17, a surrogate monoclonal antibody for efalizumab

muM17 is an anti-mouse CD11a monoclonal antibody (mAb) developed as a surrogate molecule for assessing potential reproductive toxicities of efalizumab, an anti-human CD11a mAb approved for treatment of chronic moderate to severe plaque psoriasis. This article shows the use of a mechanism-based PK/PD model for muM17 to further support the determination of dose equivalency of muM17 in the mouse and efalizumab in humans based on CD11a expression on T-lymphocytes (PD). Patients in clinical studies received 1 mg/kg/week efalizumab subcutaneously for 12 weeks. In the mouse model, a single IV dose of 1 or 10 mg/kg or a single SC dose of 3, 5, or 10 mg/kg muM17 was administered. Drug concentrations and PD were quantitated using ELISA and flow cytometry (FACS) analyses, respectively. The PK/PD model of muM17 in mice was developed and was validated using sparse data from a separate multiple dose PK/PD study. The model was next used to simulate PD profiles with multiple dosing regimens mimicking those of the clinical dose of efalizumab. The model showed that 3 mg/kg/week SC administration of muM17 in mice is the minimum dose that can produce PD effects similar to those produced following 1 mg/kg/week SC of efalizumab in humans.     

Population Pharmacokinetic and Pharmacodynamic Model-Based Comparability Assessment of a Recombinant Human Epoetin Alfa and the Biosimilar HX575

The aim of this study was to develop an integrated pharmacokinetic and pharmacodynamic (PK/PD) model and assess the comparability between epoetin alfa HEXAL/Binocrit (HX575) and a comparator epoetin alfa by a model-based approach. PK/PD data-including serum drug concentrations, reticulocyte counts, red blood cells, and hemoglobin levels-were obtained from 2 clinical studies. In sum, 149 healthy men received multiple intravenous or subcutaneous doses of HX575 (100 IU/kg) and the comparator 3 times a week for 4 weeks. A population model based on pharmacodynamics-mediated drug disposition and cell maturation processes was used to characterize the PK/PD data for the 2 drugs. Simulations showed that due to target amount changes, total clearance may increase up to 2.4-fold as compared with the baseline. Further simulations suggested that once-weekly and thrice-weekly subcutaneous dosing regimens would result in similar efficacy. The findings from the model-based analysis were consistent with previous results using the standard noncompartmental approach demonstrating PK/PD comparability between HX575 and comparator. However, due to complexity of the PK/PD model, control of random effects was not straightforward. Whereas population PK/PD model-based analyses are suited for studying complex biological systems, such models have their limitations (statistical), and their comparability results should be interpreted carefully.     

Pharmacodynamics-mediated drug disposition (PDMDD) and precursor pool lifespan model for single dose of romiplostim in healthy subjects

The objective of this study was to characterize the pharmacokinetics and pharmacodynamics (PK-PD) of romiplostim after single-dose administration in healthy subjects. The mean serum romiplostim concentrations (PK data) and mean platelet counts (PD data) collected from 32 subjects receiving a single intravenous (0.3, 1 and 10 μg/kg) or subcutaneous (0.1, 0.3, 1, and 2 μg/kg) dose were fitted simultaneously to a mechanistic PK-PD model based on pharmacodynamics-mediated drug disposition (PDMDD) and a precursor pool lifespan concept. The two-compartment PK model incorporated receptor-mediated endocytosis and linear mechanisms as parallel elimination pathways. The maximal concentration of receptors (assumed to be proportional to the platelet count), the equilibrium dissociation constant, and the first-order internalization rate constant for endocytosis of the drug-receptor complex were 0.022 fg/platelet, 0.131 ng/mL, and 0.173 h?1, respectively. Romiplostim concentration stimulates the production of platelet precursors via the Hill function, where the SC?? was 0.052 ng/mL and S (max) was 11.2. The estimated precursor cell and platelet lifespans were 5.9 and 10.5 days, respectively. Model-based simulations revealed that the romiplostim exposure and the platelet response are both dependent on the dose administered and the baseline platelet counts. Also, weekly dosing produced a sustained PD response while dosing intervals ≥2 weeks resulted in fluctuating platelet counts. Thus, the mechanistic PK-PD model was suitable for describing the romiplostim PK-PD interplay (PDMDD), the dose-dependent platelet stimulation, and the lifespans of thrombopoietic cell populations.     

Population Pharmacokinetics (PK) and Pharmacodynamics (PD) Analysis of LY3015014, a Monoclonal Antibody to Protein Convertase Subtilisin/Kexin Type 9 (PCSK9) in Healthy Subjects and Hypercholesterolemia Patients

Purpose

LY3015014 is a humanized immunoglobulin G4 (IgG4) monoclonal antibody that binds to the catalytic domain of PCSK9 and reduce low-density lipoprotein cholesterol (LDL-C) in patients with hypercholesterolemia that is poorly controlled by maximally tolerated statin therapy. The objective of this pharmacokinetic/pharmacodynamics (PK/PD) analysis was to characterize the PK and PD properties of LY3015014 and assess the effect of covariates on the LY3015014 PK-PD profiles.

Methods

Single and multiple dose data from three phase1 studies in healthy subjects (n?=?133), as well as a phase 2 study in hypercholesterolemia patients (n?=?527) were combined into a single dataset for analysis. In this dataset, healthy subjects received single intravenous (IV) doses of 0.03 to 10 mg/kg, or multiple subcutaneous (SC) doses of 1.0 to 3.0 mg/kg, administered every 2 to 4 weeks, while patients received 20 to 300 mg every 4 weeks or 100 to 300 every 8 weeks. PK/PD analysis was performed using NONMEM (ICON, software version 7.0 level 3). PK and PD modeling were conducted sequentially, with PK parameters fixed during the development of the PK/PD model. PD parameters and estimated intersubject and intrasubject variability were obtained based on pharmacological drug exposure-response relationships. Age, baseline total PCSK9, body weight, diabetes diagnosis, hypercholesterolemia disease status, dose, ezetimibe administration, gender, ethnic origin, metabolic syndrome, and satin administration were evaluated as potential covariates in the PK model. Baseline total PCSK9, baseline LDL-C, diabetes diagnosis, disease status, ezetimibe administration, gender, ethnic origin, metabolic syndrome, and Statin administration were evaluated as potential covariates in the PD model.

Results

LY3015014 PK profile was consistent across all the studies and between healthy subjects and hypercholesterolemia patients. The PK time course data were well described by a two compartment PK model with first order absorption, and covariates identified for PK parameters included weight on both clearance (CL) and central volume (V2), dose on CL, race on bioavailability (F), and age on V2. The PD (LDL-C) was described using an indirect response model with LY3015014 acting to stimulate the elimination of LDL-C. Covariates identified to have a statistically significant impact on PD were coadministration of statins, baseline LDL-C, metabolic syndrome status and gender.

Conclusions

The population PK/PD model adequately describes the PK and PD profiles of LY3015014. Identification of clinically significant covariates will support the design and dose selection for the pivotal registration studies, ensuring that patients are dosed appropriately.

     

Population pharmacokinetic/pharmacodynamic (PK/PD) modelling of the hypothalamic-pituitary-gonadal axis following treatment with GnRH analogues

AIMS: To develop a population pharmacokinetic/pharmacodynamic (PK/PD) model of the hypothalamic-pituitary-gonadal (HPG) axis describing the changes in luteinizing hormone (LH) and testosterone concentrations following treatment with the gonadotropin-releasing hormone (GnRH) agonist triptorelin and the GnRH receptor blocker degarelix. METHODS: Fifty-eight healthy subjects received single subcutaneous or intramuscular injections of 3.75 mg of triptorelin and 170 prostate cancer patients received multiple subcutaneous doses of degarelix of between 120 and 320 mg. All subjects were pooled for the population PK/PD data analysis. A systematic population PK/PD model-building framework using stochastic differential equations was applied to the data to identify nonlinear dynamic dependencies and to deconvolve the functional feedback interactions of the HPG axis. RESULTS: In our final PK/PD model of the HPG axis, the half-life of LH was estimated to be 1.3 h and that of testosterone 7.69 h, which corresponds well with literature values. The estimated potency of LH with respect to testosterone secretion was 5.18 IU l(-1), with a maximal stimulation of 77.5 times basal testosterone production. The estimated maximal triptorelin stimulation of the basal LH pool release was 1330 times above basal concentrations, with a potency of 0.047 ng ml(-1). The LH pool release was decreased by a maximum of 94.2% by degarelix with an estimated potency of 1.49 ng ml(-1). CONCLUSIONS: Our model of the HPG axis was able to account for the different dynamic responses observed after administration of both GnRH agonists and GnRH receptor blockers, suggesting that the model adequately characterizes the underlying physiology of the endocrine system.     

Pharmacokinetic/pharmacodynamic modeling of luteinizing hormone (LH) suppression and LH surge delay by cetrorelix after single and multiple doses in healthy premenopausal women

Cetrorelix (CET) is a potent luteinizing hormone-releasing hormone (LH-RH) antagonist and is used to prevent premature ovulation in IVF (in vitro fertilization) procedures. The objective of the present study was to develop a pharmacokinetic/pharmacodynamic (PK/PD) model for the LH suppression and LH surge delay after single doses (SD) and multiple doses (MD) of CET in healthy premenopausal women without ovarian stimulation. CET was given by subcutaneous route (SD, 0.25, 0.5, or 1 mg) on cycle day 3 and as similar multiple once-a-day doses from cycle day 3 to day 16 in two consecutive menstrual cycles. The concentration-time data of CET and LH were used for PK/PD modeling. A two-compartment model described the PK of CET with median terminal half-life estimates of 9.2 and 54.5 hours after SD and MD, respectively. An indirect-response Emax model was used to describe the LH suppression and the LH surge delay. LH suppression was linked to plasma concentrations of CET, while the delay in the LH surge was linked to the PK of CET through a hypothetical effect compartment. Since the SD regimen on day 3 did not cause significant delay, these values were used as controls in the analysis of surge delay in MD data. The IC50 (for suppression) estimate was 0.73 ng/ml for SD, and EC50 (surge delay) was 1.42 ng/ml for MD. The PK/PD model adequately described the LH suppression and the surge delay.     

Modification of the Fc Region of a Human Anti‐oncostatin M Monoclonal Antibody for Higher Affinity to FcRn Receptor and Extension of Half‐life in Cynomolgus Monkeys

The purpose of this study was to evaluate the pharmacokinetics (PK) of anti‐oncostatin M (OSM) IgG1 monoclonal antibodies, CNTO 1119 and its Fc variant (CNTO 8212), which incorporates the LS(Xtend) mutation to extend terminal half‐life (T_1/2), after a single intravenous (IV) or subcutaneous (SC) administration in cynomolgus monkeys, and to predict human PK. In study 1, single doses of CNTO 1119 and CNTO 8212 were administered IV or SC at 3 mg/kg to cynomolgus monkeys (n = 3 per group). In study 2, single doses of CNTO 8212 were administered IV at 1, 5 or 20 mg/kg, or SC at 5 mg/kg to cynomolgus monkeys (n = 5 per group). Serial blood samples were collected for assessment of serum concentrations of CNTO 1119 and/or CNTO 8212. A two‐compartment population PK model with first‐order elimination was utilized to simultaneously describe the serum concentrations of CNTO 1119 and CNTO 8212 over time after IV and SC administration in cynomolgus monkeys. The typical population PK parameter estimates for CNTO 1119 in cynomolgus monkeys were clearance (CL) = 2.81 mL/day/kg, volume of distribution of central compartment (V₁) = 31.3 mL/kg, volume of distribution of peripheral compartment (V₂) = 23.3 mL/kg, absolute bioavailability (F) = 0.84 and T_1/2 = 13.4 days. In comparison, the typical population PK parameter estimates for CNTO 8212 in cynomolgus monkeys were CL = 1.41 mL/day/kg, V₁ = 39.8 mL/kg, V₂ = 32.6 mL/kg, F = 0.75 and T_1/2 = 35.7 days. The mean CL of CNTO 8212 was ~50% lower compared with that for CNTO 1119 in cynomolgus monkeys. The overall volume of distribution (V₁+V₂) for CNTO 8212 was about 32% larger compared with that for CNTO 1119, but generally similar to the vascular volume in cynomolgus monkeys. The T_1/2 of CNTO 8212 was significantly (p < 0.05) longer by about 2.7‐fold than that for CNTO 1119 in cynomolgus monkeys. Thus, the modification of the Fc portion of an anti‐OSM IgG1 mAb for higher FcRn binding affinity resulted in lower systemic clearance and a longer terminal half‐life in cynomolgus monkeys. CNTO 8212 demonstrated linear PK after a single IV dose (1–20 mg/kg) in cynomolgus monkeys. The predicted human PK parameters suggest that CNTO 8212 is likely to exhibit slow clearance and long terminal half‐life in human beings and may likely allow less frequent dosing in the clinical setting.     

Population Pharmacokinetic Modeling of LY2189102 after Multiple Intravenous and Subcutaneous Administrations

Interleukin-1 beta (IL-1β) is an inflammatory mediator which may contribute to the pathophysiology of rheumatoid arthritis (RA) and type 2 diabetes mellitus (T2DM). Population pharmacokinetics (PK) of LY2189102, a high affinity anti-IL-1β humanized monoclonal immunoglobulin G4 evaluated for efficacy in RA and T2DM, were characterized using data from 79 T2DM subjects (Study H9C-MC-BBDK) who received 13 weekly subcutaneous (SC) doses of LY2189102 (0.6, 18, and 180 mg) and 96 RA subjects (Study H9C-MC-BBDE) who received five weekly intravenous (IV) doses (0.02–2.5 mg/kg). Frequency of anti-drug antibody (ADA) development appears dose-dependent and is different between studies (36.7% in Study H9C-MC-BBDK vs. 2.1% in Study H9C-MC-BBDE), likely due to several factors, including differences in patient population and background medications, administration routes, and assays. A two-compartment model with dose-dependent bioavailability best characterizes LY2189102 PK following IV and SC administration. Typical elimination and distribution clearances, central and peripheral volumes of distribution are 0.222 L/day, 0.518 L/day, 3.08 L, and 1.94 L, resulting in a terminal half-life of 16.8 days. Elimination clearance increased linearly, yet modestly, with baseline creatinine clearance and appears 37.6% higher in subjects who developed ADA. Bioavailability (0.432–0.721) and absorption half-life (94.3–157 h) after SC administration are smaller with larger doses. Overall, LY2189102 PK is consistent with other therapeutic humanized monoclonal antibodies and is likely to support convenient SC dosing.KEY WORDS: interleukin-1 beta (IL-1β), LY2189102, population pharmacokinetics     

Pharmacokinetics and pharmacodynamics of propofol in patients undergoing abdominal aortic surgery

Available propofol pharmacokinetic protocols for target-controlled infusion (TCI) were obtained from healthy individuals. However, the disposition as well as the response to a given drug may be altered in clinical conditions. The aim of the study was to examine population pharmacokinetics (PK) and pharmacodynamics (PD) of propofol during total intravenous anesthesia (propofol/fentanyl) monitored by bispectral index (BIS) in patients scheduled for abdominal aortic surgery. Population nonlinear mixed-effect modeling was done with Nonmem. Data were obtained from ten male patients. The TCI system (Diprifusor) was used to administer propofol. The BIS index served to monitor the depth of anesthesia. The propofol dosing was adjusted to keep BIS level between 40 and 60. A two-compartment model was used to describe propofol PK. The typical values of the central and peripheral volume of distribution, and the metabolic and inter-compartmental clearance were V(C) = 24.7 l, V(T) = 112 l, Cl = 2.64 l/min and Q = 0.989 l/min. Delay of the anesthetic effect, with respect to plasma concentrations, was described by the effect compartment with the rate constant for the distribution to the effector compartment equal to 0.240 min(-1). The BIS index was linked to the effect site concentrations through a sigmoidal E(max) model with EC(50) = 2.19 mg/l. The body weight, age, blood pressure and gender were not identified as statistically significant covariates for all PK/PD parameters. The population PK/PD model was successfully developed to describe the time course and variability of propofol concentration and BIS index in patients undergoing surgery.     

Pharmacokinetic/Pharmacodynamic Modellingof GnRH Antagonist Degarelix: A Comparisonof the Non-linear Mixed-Effects Programs NONMEM and NLME

In this paper, the two non-linear mixed-effects programs NONMEM and NLME were compared for their use in population pharmacokinetic/pharmacodynamic (PK/PD) modelling. We have described the first-order conditional estimation (FOCE) method as implemented in NONMEM and the alternating algorithm in NLME proposed by Lindstrom and Bates. The two programs were tested using clinical PK/PD data of a new gonadotropin-releasing hormone (GnRH) antagonist degarelix currently being developed for prostate cancer treatment. The pharmacokinetics of intravenous administered degarelix was analysed using a three compartment model while the pharmacodynamics was analysed using a turnover model with a pool compartment. The results indicated that the two algorithms produce consistent parameter estimates. The bias and precision of the two algorithms were further investigated using a parametric bootstrap procedure which showed that NONMEM produced more accurate results than NLME together with the nlmeODE package for this specific study.     

Target-mediated pharmacokinetic and pharmacodynamic model of exendin-4 in rats, monkeys, and humans

A mechanism-based pharmacokinetic-pharmacodynamic (PK/PD) model was developed for exendin-4 to account for receptor-mediated endocytosis via glucagon-like peptide 1 receptor (GLP-1R) as the primary mechanism for its nonlinear disposition. Time profiles of exendin-4 concentrations after intravenous, subcutaneous, and continuous intravenous infusion doses in rats, intravenous and subcutaneous doses in monkeys, and intravenous infusion and subcutaneous doses in humans were examined. Mean data for glucose and insulin after glucose challenges during exendin-4 treatment in healthy rats were analyzed. The PK model components included receptor binding, subsequent internalization and degradation, nonspecific tissue distribution, and linear first-order elimination from plasma. The absorption rate constant (k(a)) decreased with increasing doses in all three species. The clearance from the central compartment (CL(c)) (rats, 3.62 ml/min; monkeys, 2.39 ml · min(-1) · kg(-1); humans, 1.48 ml · min(-1) · kg(-1)) was similar to reported renal clearances. Selected PK parameters (CL(c), V(c), and k(off)) correlated allometrically with body weight. The equilibrium dissociation constant (K(D)) was within the reported range in rats (0.74 nM), whereas the value in monkeys (0.12 pM) was much lower than that in humans (1.38 nM). The effects of exendin-4 on the glucose-insulin system were described by a feedback model with a biphasic effect equation driven by free exendin-4 concentrations. Our generalized nonlinear PK/PD model for exendin-4 taking into account of drug binding to GLP-1R well described PK profiles after various routes of administration over a large range of doses in three species along with PD responses in healthy rats. The present model closely reflects underlying mechanisms of disposition and dynamics of exendin-4.     

A pharmacokinetic/pharmacodynamic approach to predict the cumulative cortisol suppression of inhaled corticosteroids

The suppression of endogenous cortisol release is one of the major systemic side effects of inhaled corticosteroids in the treatment of asthma. The circadian rhythm of the endogenous cortisol release and the resulting plasma concentrations as well as the release suppression during corticosteroid therapy could previously be described with an integrated PK/PD model. Based on this model, a PK/PD approach was developed to quantify and predict the cumulative cortisol suppression (CCS) as a surrogate marker for the systemic activity of inhaled corticosteroid therapy. The presented method was applied to predict CCS after single doses and during short-term multiple dosing of the inhaled corticosteroids flunisolide (FLU), fluticasone propionate (FP), and triamcinolone acetonide (TCA), and after oral methylprednisolone as systemic reference therapy. Drug-specific PK and PD parameters were obtained from previous single-dose studies and extrapolated to the multiple-dose situation. For single dosing, a similar CCS within the range of 16-21% was predicted for FP 250 micrograms, FLU 500 micrograms, and TCA 1000 micrograms. For multiple dosing, a respective CCS of 28-33% was calculated for FLU 500 micrograms bid, FP 250 micrograms, bid, and TCA 1000 micrograms bid. Higher cortisol suppression compared to these single and multiple dosing regimens of the inhaled corticosteroids was predicted after oral doses of only 1 mg and 2 mg methylprednisolone, respectively. The predictive power of the approach was evaluated by comparing the PK/PD-based simulations with data reported previously in clinical studies. The predicted CCS values were in good correlation with the clinically observed results. Hence, the presented PK/PD approach allows valid predictions of CCS for single and short-term multiple dosing of inhaled corticosteroids and facilitates comparisons between different dosing regimens and steroids.     

Pharmacokinetics/Pharmacodynamics of Nondepleting Anti-CD4 Monoclonal Antibody (TRX1) in Healthy Human Volunteers

Purpose TRX1 is a nondepleting anti-CD4 monoclonal IgG1 antibody being developed to induce tolerance by blocking CD4-mediated functions. The purpose of this study is to describe the pharmacokinetics (PK) and pharmacodynamics (PD) of TRX1 and to develop a receptor-mediated PK/PD model that characterizes the relationships between serum TRX1 concentration and total and free CD4 expression in healthy male volunteers. Methods Nine subjects from three dosing cohorts in double-blinded, placebo-controlled phase I clinical study was included in the analysis. Serum TRX1 levels were determined using enzyme-linked immunosorbent assay. Blood total and free CD4 receptor levels were determined by using flow cytometric analyses. The receptor-mediated PK/PD model was developed to describe the dynamic interaction of TRX1 binding with CD4 receptors. Results and Conclusions TRX1 displayed nonlinear pharmacokinetic behavior and the CD4 receptors on T cells were saturated and down-modulated following treatment with TRX1. Results from in vitro studies using purified human T cells suggested that CD4-mediated internalization may constitute one pathway by which CD4 is down-modulated and TRX1 is cleared in vivo. The developed receptor-mediated PK/PD model adequately described the data. This PK/PD model was used to simulate PK/PD time profiles after different dosing regimens to help guide the dose selection in future clinical studies.