Clinical dose-response for a broad set of biological products: A model-based meta-analysis

Characterizing clinical dose-response is a critical step in drug development. Uncertainty in the dose-response model when planning a dose-ranging study can often undermine efficiency in both the design and analysis of the trial. Results of a previous meta-analysis on a portfolio of small molecule compounds from a large pharmaceutical company demonstrated a consistent dose-response relationship that was well described by the maximal effect model. Biologics are different from small molecules due to their large molecular sizes and their potential to induce immunogenicity. A model-based meta-analysis was conducted on the clinical efficacy of 71 distinct biologics evaluated in 91 placebo-controlled dose-response studies published between 1995 and 2014. The maximal effect model, arising from receptor occupancy theory, described the clinical dose-response data for the majority of the biologics (81.7%, n = 58). Five biologics (7%) with data showing non-monotonic trend assuming the maximal effect model were identified and discussed. A Bayesian model-based hierarchical approach using different joint specifications of prior densities for the maximal effect model parameters was used to meta-analyze the whole set of biologics excluding these five biologics ( n = 66). Posterior predictive distributions of the maximal effect model parameters were reported and they could be used to aid the design of future dose-ranging studies. Compared to the meta-analysis of small molecules, the combination of fewer doses, narrower dosing ranges, and small sample sizes further limited the information available to estimate clinical dose-response among biologics.

A quantitative systems pharmacology model of blood coagulation network describes in vivo biomarker changes in non-bleeding subjects

Essentials Baseline coagulation activity can be detected in non-bleeding state by in vivo biomarker levels. A detailed mathematical model of coagulation was developed to describe the non-bleeding state. Optimized model described in vivo biomarkers with recombinant activated factor VII treatment. Sensitivity analysis predicted prothrombin fragment 1 + 2 and D-dimer are regulated differently.

SUMMARY

Background Prothrombin fragment 1 + 2 (F1 + 2 ), thrombin-antithrombin III complex (TAT) and D-dimer can be detected in plasma from non-bleeding hemostatically normal subjects or hemophilic patients. They are often used as safety or pharmacodynamic biomarkers for hemostatis-modulating therapies in the clinic, and provide insights into in vivo coagulation activity. Objectives To develop a quantitative systems pharmacology (QSP) model of the blood coagulation network to describe in vivo biomarkers, including F1 + 2 , TAT, and D-dimer, under non-bleeding conditions. Methods The QSP model included intrinsic and extrinsic coagulation pathways, platelet activation state-dependent kinetics, and a two-compartment pharmacokinetics model for recombinant activated factor VII (rFVIIa). Literature data on F1 + 2 and D-dimer at baseline and changes with rFVIIa treatment were used for parameter optimization. Multiparametric sensitivity analysis (MPSA) was used to understand key proteins that regulate F1 + 2 , TAT and D-dimer levels. Results The model was able to describe tissue factor (TF)-dependent baseline levels of F1 + 2 , TAT and D-dimer in a non-bleeding state, and their increases in hemostatically normal subjects and hemophilic patients treated with different doses of rFVIIa. The amount of TF required is predicted to be very low in a non-bleeding state. The model also predicts that these biomarker levels will be similar in hemostatically normal subjects and hemophilic patients. MPSA revealed that F1 + 2 and TAT levels are highly correlated, and that D-dimer is more sensitive to the perturbation of coagulation protein concentrations. Conclusions A QSP model for non-bleeding baseline coagulation activity was established with data from clinically relevant in vivo biomarkers at baseline and changes in response to rFVIIa treatment. This model will provide future mechanistic insights into this system.

Quantitative Systems Pharmacology: A Case for Disease Models

Quantitative systems pharmacology (QSP) has emerged as an innovative approach in model‐informed drug discovery and development, supporting program decisions from exploratory research through late‐stage clinical trials. In this commentary, we discuss the unique value of disease‐scale “platform” QSP models that are amenable to reuse and repurposing to support diverse clinical decisions in ways distinct from other pharmacometrics strategies.

A mechanistic PK/PD model for two anti-IL13 antibodies explains the difference in total IL-13 accumulation observed in clinical studies

IMA-638 and IMA-026 are humanized IgG1 monoclonal antibodies (mAbs) that target non-overlapping epitopes of IL-13. Separate first-in-human single ascending dose studies were conducted for each mAb. These studies had similar study designs, but mild to moderate asthmatics were recruited for the IMA-638 study and healthy subjects were recruited for the IMA-026 study. IMA-638 and IMA-026 showed similar pharmacokinetic (PK) profiles, but very different total IL-13 (free and drug bound IL-13) profiles; free IL13 was not measured. IMA-026 treatment induced a dose-dependent accumulation of total IL-13, while IMA-638 treatment led to a much smaller accumulation without any clear dose-response. To understand the differences between the two total IL-13 profiles and to predict the free IL-13 profiles for each mAb treatment, a mechanistic PK/pharmacodynamic model was developed. PK-related parameters were first fit to the mean PK profiles of each mAb separately; thereafter, the target-related parameters were fit to both total IL-13 profiles simultaneously. The IL-13 degradation rate was assumed to be the same for asthma patients and healthy subjects. The model suggests that an approximately 100× faster elimination of IL-13-IMA-638 complex than IL-13-IMA-026 complex could be responsible for the differences observed in total IL-13 profiles for the two mAbs. Furthermore, the model predicts that IMA-638 administration results in greater and more prolonged free IL-13 inhibition than equivalent dosing of IMA-026 despite similar binding KD and PK profile. In conclusion, joint modeling of two similar molecules provided mechanistic insight that the elimination rate of mAb-target complex can regulate the degree of free target inhibition.

Using a Systems Pharmacology Model of the Blood Coagulation Network to Predict the Effects of Various Therapies on Biomarkers

A number of therapeutics have been developed or are under development aiming to modulate the coagulation network to treat various diseases. We used a systems model to better understand the effect of modulating various components on blood coagulation. A computational model of the coagulation network was built to match in-house in vitro thrombin generation and activated Partial Thromboplastin Time (aPTT) data with various concentrations of recombinant factor VIIa (FVIIa) or factor Xa added to normal human plasma or factor VIII-deficient plasma. Sensitivity analysis applied to the model revealed that lag time, peak thrombin concentration, area under the curve (AUC) of the thrombin generation profile, and aPTT show different sensitivity to changes in coagulation factors' concentrations and type of plasma used (normal or factor VIII-deficient). We also used the model to explore how variability in concentrations of the proteins in coagulation network can impact the response to FVIIa treatment.

Predictability of intranasal pharmacokinetics in man using pre-clinical pharmacokinetic data with a dopamine 3 receptor agonist, PF-219061

(R)-3-(4-propylmorpholin-2-yl) phenol (PF-219061) is a potent, selective agonist of the dopamine 3 receptor for the treatment of female sexual dysfunction. In vivo, PF-219061 exhibits liver blood flow clearance in both rat and dog. Oral bioavailability was 0.7% in dog and less than 5% in rat. Intranasal dosing was investigated to improve bioavailability. Pre-clinical assessments in rat and dog demonstrated intranasal bioavailabilities of 16-38% in rat and 54-61% in dog with very rapid absorption. It was predicted that an intranasal dose in man would give approximately 25-50% bioavailability. The clinical data verified the preclinical predictions demonstrating rapid absorption and approximately dose-proportional increases in exposure. The intranasal bioavailability in man was estimated to be 26-38%. These findings indicate the potential utility of intranasal dosing as a route that circumvents the first-pass effects for PF-219061 resulting in high exposures.

Development and validation of a quantitative cell-based bioassay for comparing the pharmacokinetic profiles of two recombinant erythropoietic proteins in serum

An in vitro cell-based bioassay was developed and validated to assess the pharmacokinetic profiles of two novel therapeutic recombinant proteins (EP1 and EP2) with erythropoiesis stimulating properties in Sprague-Dawley rats. While immunoassays are the standard choice for evaluating the pharmacokinetic parameters of drugs, no immunoassay was available for EP2, necessitating the need for a quantitative bioassay capable of measuring both EP1 and EP2 separately so that appropriate comparisons could be made. The bioassay described here utilizes a sub clone of the murine 32D cell line transfected with the gene encoding for the human erythopoietin (HuEPO) receptor. Erythropoietin (EPO), EP1 and EP2 exert their proliferative effect on the cell line by signaling through the HuEPO receptor. The proliferation induced by the erythropoietic proteins was measured by [methyl-(3)H]thymidine incorporation into the cellular DNA. The assay was conducted in 96-well microtiter plates and had relatively high throughput. The Guidelines of the International Conference on Harmonization (ICH) were followed for the validation of the different assay parameters including robustness, linearity, accuracy, precision, limit of quantitation (LOQ) and specificity. The robustness of the bioassay is demonstrated by the lack of an effect of age of the 32D cell culture on the performance of the EP2 bioassay. The bioassay demonstrated good linearity, yielding a coefficient of determination of 0.99 or higher for both EP1 and EP2. The assay showed reproducible dose-response curves for EP1 in the range of 0.039-2.5 ng/mL and for EP2 in the range of 0.125-8 ng/mL. The accuracy estimates ranged between 98% and 108% for EP1 and between 90% and 110% for EP2 in the reproducible range mentioned above. Intermediate precision (within-plate R.S.D.) in the same range was within 26% and 17% for the EP1 and EP2 bioassays, respectively. The validated bioassays for EP1 and EP2 were utilized to quantitatively analyze serum samples from a pharmacokinetic study conducted to compare the profiles of the two compounds in Sprague-Dawley rats.

Development and evaluation of a population pharmacokinetic-pharmacodynamic model of darbepoetin alfa in patients with nonmyeloid malignancies undergoing multicycle chemotherapy

Anemia is frequently observed in patients undergoing chemotherapy. Administration of darbepoetin alfa, a recombinant erythropoiesis-stimulating agent that has longer residence time than endogenous erythropoietin, to patients with chemotherapy-induced anemia (CIA) increases mean hemoglobin concentration, reduces risk of red blood cell transfusions, and improves patient-reported outcomes. A pharmacokinetic/pharmacodynamic (PkPd) model was developed using data from patients with nonmyeloid malignancies and CIA who were receiving darbepoetin alfa. A 2-compartment Pk model with linear elimination described the Pk data obtained in 140 CIA patients after intravenous and subcutaneous (s.c.) doses of 2.25 microg/kg every week and s.c. doses of 6.75 microg/kg every 3 weeks. The population typical values of key Pk parameters were clearance, 2010 mL/day; steady-state volume of distribution, 3390 mL; and bioavailability, 44.3%. A modified indirect response model, wherein serum concentrations stimulated the production of hemoglobin through an Emax-type equation, described the hemoglobin levels after s.c. doses of 0.5 microg/kg every week to 15 microg/kg every 3 weeks in 573 CIA patients. The estimated incremental maximum stimulation of hemoglobin production was 43.7% and darbepoetin alfa serum concentration at half-maximal stimulation was 3.68 ng/mL. The impact of covariates (body weight and platinum-containing chemotherapy) on the PkPd response was evaluated based on point and interval estimates of parameters, rather than through stepwise hypothesis testing. The final PkPd model adequately predicted hemoglobin response in a test data set, thereby confirming the predictive capability of the model. Based on simulations, it was not possible to categorize the influence of any covariate as clinically important.

The absorption of darbepoetin alfa occurs predominantly via the lymphatics following subcutaneous administration to sheep

PURPOSE

To determine the contribution of the lymphatics to the systemic availability of darbepoetin alfa (DA) using an established sheep model.

MATERIALS AND METHODS

DA was administered either by intravenous (IV) injection (0.2, 0.5 or 2 microg/kg) or by subcutaneous (SC) administration (2 microg/kg) into the interdigital space of the hind leg. A SC control group was used to determine the absolute bioavailability (F (sys)). Cannulation of the peripheral lymphatics in a parallel SC group allowed the continuous collection of lymph draining the injection site and determination of the cumulative amount of DA absorbed via the lymphatics. Serum and lymph concentrations of DA were determined by ELISA. The fraction of the dose absorbed into the lymphatics (F (lymph)) relative to the fraction absorbed directly into the blood (F (blood)) was determined using a compartmental approach.

RESULTS

Dose-linear pharmacokinetics was observed within the dose range investigated. The bioavailability was virtually complete following SC injection into the interdigital space (88.4 +/- 15.7%). A high proportion of the administered dose was recovered in peripheral lymph (90.2 +/- 4.4%) resulting in a substantial reduction in the systemic availability in lymph cannulated animals (3.7%).

CONCLUSION

The high recovery of DA in the peripheral lymph demonstrated near complete absorption of this recombinant protein via the lymphatics in a lymph cannulated sheep model.

Pharmacokinetics of Darbepoetin Alfa after Intravenous or Subcutaneous Administration in Patients with Non-myeloid Malignancies Undergoing Chemotherapy

BACKGROUND AND OBJECTIVE

The pharmacokinetics of darbepoetin alfa after intravenous (IV) administration in the oncology setting have not been previously reported. The objective of this study was to evaluate the pharmacokinetics of IV or subcutaneous (SC) darbepoetin alfa in patients with non-myeloid malignancies undergoing multicycle chemotherapy.

METHODS

Fifty-six patients (haemoglobin <or=13.0 g/dL) received weekly darbepoetin alfa 2.25 microg/kg administered either IV (n=27) or SC (n=29) during up to three cycles of chemotherapy. Noncompartmental pharmacokinetic analysis was performed, including analysis of intensive pharmacokinetic profiles collected over 168 hours during week 1 of both the first and third cycles of chemotherapy.

RESULTS

Darbepoetin alfa serum concentrations exhibited a biphasic profile (a rapid distributive phase followed by a slower terminal elimination phase) after IV administration, whereas darbepoetin alfa was slowly absorbed after SC administration. Darbepoetin alfa exhibited limited extravascular distribution after IV administration, with both initial and steady-state mean volumes of distribution (36.1 mL/kg and 55.2 mL/kg, respectively, after a single IV dose) approximating the plasma volume. After a single IV dose, darbepoetin alfa exhibited a mean clearance of 1.05 mL/h/kg, with a mean terminal half-life of 38.8 hours. Similar pharmacokinetic results were observed after single and multiple doses of darbepoetin alfa, for both SC and IV administration.

CONCLUSION

Darbepoetin alfa is cleared slowly after IV administration to patients with cancer receiving chemotherapy, resulting in a terminal half-life of 38.8 hours. No evidence of accumulation and no changes in pharmacokinetic profiles after repeated administration were observed in cancer patients undergoing cyclic chemotherapy, for both IV and SC dosing.

A randomized, multicenter study of subcutaneous and intravenous darbepoetin alfa for the treatment of chemotherapy-induced anemia

BACKGROUND

This randomized, open-label study evaluated the efficacy, safety and pharmacokinetics of darbepoetin alfa administered intravenously (i.v.) or subcutaneously (s.c.) in chemotherapy-induced anemia.

PATIENTS AND METHODS

Patients received darbepoetin alfa i.v. (n=59) or s.c. (n=59) at a dose of 4.5 mug/kg once weekly for 6 weeks (correction phase) followed by 4.5 mug/kg once every 3 weeks for the remainder of the 18-week treatment period (maintenance phase).

RESULTS

During the correction phase, the mean [95% confidence interval (CI)] change in hemoglobin (intention-to-treat) was 1.1 (0.6-1.5) g/dl in the i.v. group and 1.3 (0.9-1.7) g/dl in the s.c. group; using available data, the mean change was 1.4 (1-1.9) g/dl and 1.6 (1.2-2) g/dl, respectively. The percentage (95% CI) of patients maintaining hemoglobin (i.e. average decrease < or =0.5 g/dl) during the maintenance phase was similar between the i.v. (82%; 95% CI 66% to 92%) and s.c. (80%; 95% CI 66% to 90%) groups. Thirty-five per cent (95% CI 20% to 50%) of patients in the i.v. group and 32% of patients in the s.c. group (95% CI 18% to 45%) received red blood cell transfusions during week 5 to the end of the treatment period. Darbepoetin alfa was well tolerated in both groups. No significant difference (P=0.36) in weekly darbepoetin alfa serum concentrations was observed between groups.

CONCLUSIONS

Darbepoetin alfa can be administered i.v. or s.c. at equal doses for the treatment of anemia in this setting.

Passage of erythropoietic agents across the blood-brain barrier: a comparison of human and murine erythropoietin and the analog darbepoetin alfa

Studies have suggested that erythropoietin (EPO) may be used to treat stroke in both animals and humans. It is thought to exert its effects directly on the brain and studies with therapeutic doses have shown that it can cross the blood-brain barrier. Here, we compared in a blinded fashion the ability of three erythropoietic agents (murine erythropoietin, human erythropoietin, and darbepoetin alfa, an analog of human erythropoietin in clinical use) to cross the blood-brain barrier of the mouse. High-performance liquid chromatography (HPLC) results showed that all three erythropoietic agents were enzymatically resistant in brain and blood. The unidirectional blood-to-brain influx rates (Ki) as measured by multiple-time regression analysis showed that all the erythropoietic agents crossed the blood-brain barrier at about the same rate as albumin, suggesting that they cross the blood-brain barrier by way of the extracellular pathways. No saturable component to influx was found, but indirect evidence suggested a brain-to-blood efflux system. The percent of the intravenously injected dose taken up per gram of brain (%Inj/g) ranged from 0.05 to 0.1 %Inj/g among the three erythropoietic agents and peaked about 3 h after IV injection. For other substances, this range of %Inj/g is known to produce direct effects on brain function. We conclude that erythropoietic agents cross the blood-brain barrier by way of the extracellular pathways in amounts that are likely sufficient to explain their neuroprotective effects.

Lymphatic Absorption Is the Primary Contributor to the Systemic Availability of Epoetin Alfa following Subcutaneous Administration to Sheep

The contribution of the lymphatics to the absorption and systemic availability of recombinant human epoetin alfa (rHuEPO) following s.c. injection was examined using a cannulated sheep model. Parallel studies were conducted in sheep where a single bolus dose was administered either by i.v. (10, 100, or 1000 IU/kg) or s.c. (400 IU/kg) injection. The first s.c. group served as a control for the calculation of absolute bioavailability. In the second group, the efferent popliteal lymphatic duct was cannulated and peripheral lymph draining the injection site was continuously collected. In the third group, the thoracic duct was cannulated to allow collection of central lymph just prior to entry into the systemic circulation. Blood was periodically sampled from all animals, and concentrations in serum and lymph were determined by enzyme-linked immunosorbent assay. The cumulative amount of rHuEPO recovered in peripheral and central lymph was 83.9 +/- 6.6% and 75.3 +/- 3.9% of the administered dose, respectively, indicating almost complete absorption from the s.c. injection site and minimal clearance during transit through the lymphatic system. After i.v. administration, the systemic clearance of rHuEPO decreased with increasing dose, reflecting capacity-limited elimination kinetics. A pharmacokinetic model was developed to simultaneously fit experimental data for all treatment groups and estimate bioavailability. The direct measurement of >75% of the dose in peripheral and central lymph independently verifies the calculated bioavailability of 87% and demonstrates the major role of the lymphatic route in the overall s.c. bioavailability of rHuEPO after s.c. administration with this animal model.

Kinetics of haematopoietic recovery after dose-intensive chemo/radiotherapy in mice: optimized erythroid support with darbepoetin alpha

Despite its frequency and impact on clinical outcomes, anaemia in cancer patients remains poorly understood and suboptimally treated. The definition of optimum treatment schedules with erythropoietic agents requires a suitable model of chemotherapy-induced progressive anaemia. This study investigated novel strategies such as once-per-chemotherapy-cycle dosing, synchronization between erythroid supportive care and chemotherapy, and definition of the optimum timing of erythroid support. A murine model of carboplatin chemotherapy/radiotherapy (CRT)-induced anaemia was used, which caused progressive anaemia across multiple cycles. Weekly administration of recombinant human erythropoietin (rHuEPO) was effective, but the longer-acting darbepoetin alpha resulted in superior responses. In all animals, anaemia became progressive and more refractory across cycles because of accumulated bone marrow damage. Exploiting a specific enzyme-linked immunosorbent assay, which could distinguish between darbepoetin alpha and endogenous erythropoietin, the effect of CRT upon the pharmacokinetics of darbepoetin alpha showed that clearance of darbepoetin alpha, and presumably erythropoietin, was at least partially dependent on a chemotherapy-sensitive pathway. Scheduling data suggested that administration of erythropoietic agents prior to chemotherapy was more effective than administration after chemotherapy. There was no evidence that erythropoietic agents exacerbated anaemia, even when administered immediately prior to CRT in an attempt to "prime" erythroid cells for the effects of CRT.

Pharmacokinetics and pharmacodynamics of darbepoetin alfa and epoetin in patients undergoing dialysis

OBJECTIVE

The aim of this multicenter, randomized, open-label study was to compare the pharmacokinetic and pharmacodynamic profiles of darbepoetin alfa, a new erythropoiesis-stimulating protein, and recombinant human erythropoietin (epoetin) after repeated intravenous dosing in patients with chronic kidney disease receiving hemodialysis.

METHODS

Forty-seven patients were randomized to receive darbepoetin alfa administered once weekly (n = 17) or 3 times weekly (n = 15) or epoetin administered 3 times weekly (n = 15) for up to 52 weeks. Pharmacokinetic profiles were measured during weeks 1 and 12 and at hemoglobin steady state (defined as a hemoglobin concentration within the target range for 4 consecutive weeks after week 12 with no change in study drug dose) or between weeks 36 and 40, whichever occurred first.

RESULTS

At each of the 3 time points evaluated, the terminal half-life of darbepoetin alfa was 2 to 3 times longer and the clearance approximately 4 times slower than those of epoetin. At week 12, the terminal half-life was 23.4 hours with darbepoetin alfa once weekly, 18.3 hours with darbepoetin alfa 3 times weekly, and 8.0 hours with epoetin 3 times weekly. The pharmacokinetics of darbepoetin alfa was not dependent on dose or time. Mean hemoglobin values at steady state were all approximately 11 g/dL, within the target range of 9.0 to 13.0 g/dL. Safety analyses revealed no differences between darbepoetin alfa and epoetin.

CONCLUSIONS

The pharmacokinetic and pharmacodynamic profiles and safety data for darbepoetin alfa demonstrate that it can be administered less frequently than epoetin in patients with chronic kidney disease receiving hemodialysis, thus simplifying anemia management.

Clinical trial simulation of a 200-microg fixed dose of darbepoetin alfa in chemotherapy-induced anemia

Our objective was to assess, using clinical trial simulation, the feasibility of a fixed 200-microg dose of darbepoetin alfa (Aranesp) administered every 2 weeks in chemotherapy-induced anemia. A pharmacokinetic/pharmacodynamic model was developed using clinical data from 547 cancer patients who received darbepoetin alfa at various doses and schedules. Monte Carlo simulations were performed for weight-based (3 microg/kg every 2 weeks) and fixed-dose (200 microg every 2 weeks) regimens and were compared with observed clinical data. Mean hemoglobin changes from baseline to end of treatment were +1.61 g/dL, +1.83 g/dL, and +1.79 g/dL for observed data, the weight-based simulation, and the fixed-dose simulation, respectively. The rates of required transfusions (hemoglobin < or = 8 g/dL) were also similar between groups. For patients between 45 and 95 kg (over 90% of the population), the impact of a fixed dose on mean hemoglobin change was negligible. There was a slight weight effect at body weight extremes (< 45 kg and > 95 kg). Clinical outcomes from simulations of weight-based andfixed dosing of darbepoetin alfa were similar to those of observed weight-based data. Given the weight distribution of a typical cancer population, the majority would be expected to benefit equally from weight-based and fixed-dose darbepoetin alfa in the amelioration of chemotherapy-induced anemia.

Pharmacokinetics of novel erythropoiesis stimulating protein (NESP) in cancer patients: preliminary report

Anaemia is a common occurrence in patients with cancer, and currently can be treated in several ways. Novel erythropoiesis stimulating protein (NESP, darbepoetin alfa) was created using site-directed mutagenesis to have 8 more sialic acid side chains than recombinant human erythropoietin (rHuEPO). The additional sialic acid content has resulted in an approximately 3-fold greater half-life relative to rHuEPO in patients with chronic renal failure. This study evaluates the pharmacokinetic profile of NESP in patients receiving multiple cycles of chemotherapy. Anaemic patients (haemoglobin < or = 11.0 g dl(-1)) who had non-myeloid malignancies received NESP weekly (2.25 mcg kg(-1) wk(-1)) under the supervision of a physician, starting on day 1 of chemotherapy for 3 chemotherapy cycles given at 3-week intervals. Blood samples were collected during chemotherapy cycles 1 and 3 for pharmacokinetic analysis. All patients were followed for 4 weeks after treatment. NESP was well tolerated by all patients. After a single dose during chemotherapy cycle 1, pharmacokinetic parameters (mean (SD), n) for the first 15 patients were: T(max)86.1 (22.8) h (n = 14); C(max)9.0 (5.1) ng ml(-1)(n = 14); t(1/2,z)32.6 (11.8) h (n = 7); CL/F 3.7 (1.0) ml h(-1) kg(-1)(n = 7). The subjects for whom all parameters could be calculated may represent a sub-group of the entire population. Similar results were obtained in cycle 3. In addition, haemoglobin response data suggests that, in this patient population, dosing less frequently than the 3 times weekly doses used for rHuEPO may be possible while improving anaemia.

Percutaneous absorption and disposition of [14C]chlordecone in young and adult female rats

The objective of this study was to evaluate the effect of age and dosage on percutaneous absorption and disposition of [14C]chlordecone (Kepone) and to describe results using a physiological based pharmacokinetic (PBPK) model. Female Fischer 344 rats 33 and 82 days old were used as the young and adult animal models, respectively, and were studied over a 10-fold dose range. [14C]Chlordecone (0.286 micromol/cm2) was applied to dorsal skin (2. 3% BSA) and radioactivity was quantified in selected tissues and excreta up to 120 h. Absorption and disposition were also determined at three dose levels up to 2.68 micromol/cm2; fraction absorbed decreased as dose increased. In vitro percutaneous absorption was measured by static and flow-through methods; these yielded similar penetration rates, which were lower than those obtained in vivo. In vivo percutaneous absorption over 120 h was 14.4+/-0.99 and 14.2+/-1. 5% dose in young and adults, respectively. Organ and tissue content increased over time (carcass>liver>kidney), indicating prolonged absorption. Statistical differences between young and old were found for liver, skin, and urine, but not for absorption. Excretion occurred primarily in feces, but also in urine. A biophysically based percutaneous model was fitted to both young and adult in vivo absorption data. This was embedded in a whole body PBPK model which, upon optimization with SAAM II, estimated apparent tissue partition coefficients, urinary and fecal excretion rates, and parameters characterizing hepatic nonlinear uptake of bound chlordecone. The model reasonably predicted tissue chlordecone content at higher doses, when decreased absorption was accounted for.

A pharmacokinetic/pharmacodynamic comparison of SAAM II and PC/WinNonlin modeling software

This paper presents a detailed comparison of the kinetic analysis software packages SAAM II and PCNonlin/WinNonlin, based on benchmark modeling problems reported in "Pharmacokinetic andPharmacodynamic Data Analysis: Concepts and Applications" (Gabrielsson and Weiner, 1994) and seven additional models. For each model, both software packages were presented with identical implementations. Models were initially executed in PCNonlin or WinNonlin and automated comparisons with SAAM II made using Microsoft Test. Models investigated included one- and multicompartment models with nonlinearities, multiple inputs and samples, multiple simultaneous experiments, and linear equations. Maximum number of compartments, data sets, and parameters were 9, 5, and 10, respectively. We compared 88 different models, many of them in different configurations, e.g., different weighting schemes or different parameter limits. The total number of attempted comparisons between SAAM II and PCNonlin was 161, of which 142 executed without problems. Parameter estimates, their precision (standard errors), and model predictions were compared; a difference of 1% or less was considered "agreement". Observed differences, mainly in parameter standard errors, can be accounted for in terms of different optimization algorithms, convergence criteria, and individual capabilities. In general, there was good agreement (<1% difference) between SAAM II and PCNonlin in terms of parameter estimates and model predictions. However, due to differences in the optimization procedure, parameter standard errors showed considerable differences. Additionally, there were differences when multiple data sets were fitted, indicating the importance of different fitting procedures for interpreting multiple kinetic data sets. The full results of the comparison and the model files in SAAM II and PCNonlin/WinNonlin formats are available from the authors.

Pharmacokinetic-pharmacodynamic modeling of stimulatory and sedative effects of alprazolam: timing performance deficits

Alprazolam decreased the reinforcement rate and increased the shorter-response rate of contingency-controlled timing behavior under a differential reinforcement of low-rate schedule (DRL 45-s) in rats. An integrated pharmacokinetic-pharmacodynamic (PK-PD) model was developed to describe and characterize the effects of i.v. and s. c. administration of alprazolam. The onset, peak and disappearance of alprazolam effects were evaluated during a 3-hr session. After s. c. alprazolam administration, two peak increases in shorter-response rate occurred at moderate alprazolam serum levels, first in the ascending and then in the descending limb of the concentration-time profile. We used a stimulation-sedation PD model incorporating two opposing effect-link sigmoidal Emax functions to model the two peaks after s.c. alprazolam administration. The model suggested that alprazolam possesses both stimulatory and sedative effects in a continuous but sequential fashion, which corresponded to low- and high-concentration effects as indicated by the EC50 values of 0.09 and 0.18 microg/ml, respectively. Owing to the rapid onset of i.v. administration, the first peak (a transition phase before the onset of the sedative effect) was absent, with the presence of the second peak again coinciding with the offset of the sedative effect. The reinforcement rate (IC50 = 0.02 microg/ml) characterized by the indirect response model to account for the initial hysteresis is an index for evaluating the deficit in timing performance. Although the effects of alprazolam can be described in behavioral terms, simultaneous PK-PD optimization numerically defines the performance and hypothesizes the coexistence of stimulation and sedation components for alprazolam. The stimulation-sedation model may help in delineating the possible mechanisms for adverse rebound side effects and of tolerance in humans.

Clearance of postprandial and lipolytically modified human HDL in rabbits and rats

Triglyceride (TG) enrichment of high density lipoproteins (HDL) in hypertriglyceridemic states renders the particles vulnerable to lipolysis, which reduces their size. In the present study we modified the size and composition of HDL in vivo in hypertriglyceridemic humans by administering a bolus of intravenous heparin, and tested the subsequent clearance of the isolated HDL particles in rabbits and rats. HDL was isolated by ultracentrifugation from 21 moderately hypertriglyceridemic humans, 5 h after ingestion of a high fat meal and then 15 min after an intravenous heparin bolus (60 U/kg). Postprandial large TG-rich preheparin HDL and small, TG-poor postheparin HDL were labeled with either 125I or 131I. The clearance of apoA-I associated with each HDL tracer was determined by injecting the tracers 1) simultaneously (n = 13) and 2) sequentially (n = 8) into male New Zealand White rabbits, an hepatic lipase-deficient animal, and 3) by injecting the tracers simultaneously into male Sprague-Dawley rats (n = 8), an animal that has hepatic lipase. Die-away curves of each radiolabeled tracer were analyzed using a two-pool model that assumes the existence of an intravascular pool in dynamic equilibrium with an extravascular pool. In the rabbit studies, the fractional catabolic rate (FCR) of small, postheparin TG-poor HDL was greater than the FCR of the larger TG-rich HDL (11% greater in the simultaneous study, P < 0.001, and 45% greater in the sequential study, P < 0.001). Opposite results were observed in rats as large TG-rich preheparin particles showed a greater FCR (1.8-fold) than smaller TG-poor postheparin HDL (P < 0.05). These data suggest that although size and composition of HDL can influence its catabolism, the effect is not always in the same direction and depends on other factors present in vivo.

A physiologically based pharmacokinetic model incorporating dispersion principles to describe solute distribution in the perfused rat hindlimb preparation

A physiologically based pharmacokinetic model incorporating dispersion principles has been developed to describe outflow data from the isolated perfused rat hindlimb preparation, for the three reference markers 14C-sucrose, 14C-urea, and 3H-water and three 14C-labeled 5-n-alkyl-5-ethyl barbiturates; the methyl, butyl, and nonyl homologues. Also 51Cr-RBC and 125I-albumin were studied. The model consists of four parallel components representing each of the tissues comprising the hindlimb: skeletal muscle, skin, bone, and adipose. Attempts to simplify the model by using the principle of tissue lumping were made by examining the tissue equilibration rate constant k tau for each of respective tissues for each compound. It was found that simplification was only possible in the case of 3H-water data. The model took into account a possible shunting component in the skin tissue and incomplete mass but not volumetric recovery from the system. The dispersion model characterizes the relative spreading of solute on transit through a tissue bed by a dimension-less parameter DN. The estimated dispersion numbers (DN) obtained were in the region of 2.7-4.72, 8.39-15.54, 0.61-2.74, and 6.02-14.0 for skeletal muscle, skin, bone, and adipose, respectively, and were independent of the compound studied. These values are much larger than the range reported in the literature for hepatic outflow data, DN = 0.2-0.5, and suggest a greater heterogeneity of vascular flow in the different component tissues of the rat hindlimb.

Discrepancies in pharmacokinetic parameter estimation between bolus and infusion studies in the perfused rat hindlimb

Isolated, perfused rat hindlimb consists of skeletal muscle, skin, bone, and adipose. Hence, it is a heterogeneous preparation composed of slowly equilibrating tissues of different characteristics and fractional flow rates. This paper shows how caution should be exercised in interpreting the results following bolus administration and subsequent statistical moment analysis of intravascular markers (51Cr-erythrocytes and 125I-albumin) and lipophilic barbiturates. For the intravascular markers, the events in the hindlimb are overshadowed by events in the connecting tubing and cannulas, due to their comparable volumes. For the barbiturates, these estimates appear to apply to short-term effects as the volume estimates obtained following infusion to steady state are greater than after bolus administration. For the extravascular markers, 14C-sucrose, 14C-urea, and 3H-water, no such time dependency was shown. However, it is only from the outflow profiles following bolus administration that events in the tissue beds can be elucidated.