Population pharmacokinetics meta-analysis of recombinant human erythropoietin in healthy subjects

A Full Target-Mediated Drug Disposition (TMDD) Model to Explain the Changes in Recombinant Human Erythropoietin (rhEpo) Pharmacokinetics in Patients with Different Bone Marrow Integrity Following Hematopoietic Transplantation

Our aim was to build a mechanistic full target-mediated drug disposition (TMDD) model for rhEpo to better understand rhEpo disposition, Epo receptor (EpoR) synthesis, and degradation in hematopoietic transplant patients with four distinct bone marrow conditions. All PK data were analyzed simultaneously using the nonlinear mixed effect modeling approach with NONMEM. The final model was a two-compartmental full TMDD model, which adequately characterizes rhEpo PK in patients and provides insight into the dynamics of free EpoR, rhEpo-EpoR, and total EpoR. The model predicted association rate constant (kon), dissociation rate constant (koff), and internalization rate constant (kint) were 0.0276 pM-1h-1, 0.647 h-1, and 0.255h-1, respectively, which were supported by experimental data. Also, the EpoR degradation rate constant (kdeg) was estimated to be 0.461 h-1. EpoR production rate was estimated to be 37.5 pM/h for adults at pre-ablation baseline and 5.91 pM/h, and 4.19 pM/h in the early post-transplant post-engraftment, and late post-transplant full engraftment. Our model provides extensive information on the dynamics of free EpoR, total EpoR and rhEpo-EpoR, and proven to be more robust and can provide more physiologically relevant binding parameters than previous models.

Population Pharmacodynamic Modeling of Epoetin Alfa in End-Stage Renal Disease Patients Receiving Maintenance Treatment Using Bayesian Approach

The ability to control dosage regimens of erythropoiesis-stimulating agents (ESAs) to maintain a desired hemoglobin (HGB) target is still elusive. We utilized a Bayesian approach and informative priors to characterize HGB profiles, using simulated drug concentrations, in patients with end-stage renal disease receiving maintenance doses of epoetin alfa. We also demonstrated an adaptive Bayesian method, applied to individual patients, to improve the accuracy of HGB predictions over time. The results showed that sparse HGB data from daily clinical practice were characterized successfully. The adaptive Bayesian method effectively improved the accuracy of HGB predictions by updating the individual model with new data accounting for within-subject changes over time. The Bayesian approach presented leverages existing knowledge of the model parameters and has a potential utility in clinical practice to individualize dosage regimens of epoetin alfa and ESAs to achieve target HGB. Further studies are warranted to develop an application for practical use.

The Utility of Pharmacometric Models in Clinical Pharmacology Research in Infants

Purpose of commentary

Acquiring knowledge on drug disposition and action in infant is challenging because of the problem of sparse and unbalanced data obtained for each individual infant due to the limited blood volume as well as the issue of extensive inter-subject and intra-subject variability in drug exposure and response due to the fast growth and dynamic maturation changes in infants. This commentary highlights the importance of using population-based pharmacometric models to improve knowledge on drug disposition and action in infants.

Recent findings

Pharmacometric modeling remains to be critical in clinical pharmacology research in infants. Many pediatric covariate models developed for scaling of drug clearance use a combination of allometric weight scaling to account for size change and a sigmoid function of antenatal development and postnatal maturation to characterize the age-related maturation. To expedite the development of safe and effective dosing regimens in infants, a number of strategies have been proposed recently, including the use of pediatric covariate model obtained from one drug for extrapolation to other drugs undergoing similar elimination pathways, as well as the combination of opportunistic clinical studies and population-based pharmacometrics models.

Summary

Population-based pharmacometric modeling plays a pivotal role in clinical pharmacology research in infants. Most of the covariate models reported so far focus on antibiotics undergoing renal elimination. Novel modeling strategies have been proposed recently to facilitate clinical pharmacology research and expedite the dose optimization process in infants.

Dose Correction for a Michaelis-Menten Approximation of a Target-Mediated Drug Disposition Model with a Multiple Intravenous Dosing Regimens

This study aimed to develop a method for implementing dose correction in a Michaelis-Menten (M-M) approximation of a target-mediated drug disposition (TMDD) model with multiple intravenous (IV) bolus administrations. We derived the formula of a correction factor (Fcorr) for each dose in a multiple IV bolus dosing regimens for M-M model. Fcorr depends on the residual free drug amount prior IV bolus dosing event and dose amount. We conducted a stochastic simulation and estimation (SSE) exercise based on therapeutic antibody PK parameters to evaluate the effect of Fcorr on parameter estimation. Previously published clinical PK data of recombinant human erythropoietin (rHuEPO) from four clinical trials in healthy subjects receiving multiple IV bolus doses were analyzed by both M-M model with and without dose correction (MMC and MMNC) as well as the rapid-binding/quasi-steady-state (RB/QSS) TMDD models. Our results showed that MMNC introduced bias to fixed-effect parameter estimates and overestimated random-effect variables. Compared with MMC, MMNC was not able to adequately characterize the nonlinearity in the PK data of antibody and rHuEPO. The MMC-based simulation demonstrated that thricely weekly 10 IU/kg rHuEPO dosing regimen yielded Fcorr = 0.5. This result suggested that the lower-than-expected exposure for rHuEPO at low dose is due to target binding. An M-M approximation of the TMDD model should include a dose correction to avoid model misfitting and potential bias in the estimated PK parameters.

Understanding Inter-Individual Variability in Monoclonal Antibody Disposition

Monoclonal antibodies (mAbs) are currently the largest and most dominant class of therapeutic proteins. Inter-individual variability has been observed for several mAbs; however, an understanding of the underlying mechanisms and factors contributing to inter-subject differences in mAb disposition is still lacking. In this review, we analyze the mechanisms of antibody disposition and the putative mechanistic determinants of inter-individual variability. Results from in vitro, preclinical, and clinical studies were reviewed evaluate the role of the neonatal Fc receptor and Fc gamma receptors (expression and polymorphism), target properties (expression, shedding, turnover, internalization, heterogeneity, polymorphism), and the influence of anti-drug antibodies. Particular attention is given to the influence of co-administered drugs and disease, and to the physiological relevance of covariates identified by population pharmacokinetic modeling, as determinants of variability in mAb pharmacokinetics.

Target-mediated disposition population pharmacokinetics model of erythropoietin in premature neonates following multiple intravenous and subcutaneous dosing regimens

Routine erythropoietin (Epo) therapy for neonatal anemia is presently controversial due to its modest response. We speculate that an important contributor to this modest response is that previous clinical study designs were not driven by rigorous mechanistic and kinetic insights into the complex pharmacokinetics (PK) and pharmacodynamics (PD) of Epo in this population. To address this therapeutic opportunity, we conducted a prospective clinical study to investigate the PK of Epo in very-low-birth-weight (VLBW) premature neonates using a unique Epo dosing algorithm that accounts for complex neonatal erythropoietic physiology. Twenty-seven subjects received up to 10 intravenous or subcutaneous exogenous doses of Epo (600 or 1200 U/kg) during the first 4 weeks of life. Subjects were administered two doses of Epo 1200 U/kg on days 2 and 16, and eight doses of Epo 600 U/kg on days 4, 5, 6, 7, 9, 14, 15, and 28 following birth. We have developed for the first time a mechanistic, target-mediated disposition model that provides novel insights into the mechanisms driving Epo PK in VLBW neonates. Epo association rate, k_on, was estimated to be 0.00610 pM^-1h^-1, and the dissociation rate k_off was 0.112 h^-1. Internalization of the Epo-target complex (k_int) and the total receptor concentration (R_max) were estimated to be 0.118 h^-1 and 133 pM, respectively. Following s.c. administration, the absorption rate (k_a) of Epo was 0.0738h^-1 and bioavailability was 78.0%. Our mechanism-based population pharmacokinetic analysis provided quantitative insight into Epo kinetics in VLBW neonates; the information gained will assist in deriving dosing strategies for neonatal anemia and for neuroprotection efficacy studies.

A Mechanism-Based Population Pharmacokinetics Model of Erythropoietin in Premature Infants and Healthy Adults Following Multiple Intravenous Doses

The objective of the current study was to develop a population pharmacokinetics (PK) model for erythropoietin (Epo) in premature infants and healthy adults to characterize the variation in PK, and to study the differences in Epo PK in these 2 populations. Thirteen very low-birth-weight premature infants (<1500 g at birth), and 10 healthy adults received up to 4 intravenous doses of Epo that ranged from 10 to 500 U/kg. The final model had a target-mediated saturable, nonlinear, elimination pathway that incorporated the mechanism of Epo binding to its receptors along with a parallel linear, central elimination pathway. Epo clearance was found to be significantly higher in preterm infants compared to adults. Epo clearance via the nonlinear pathway was found to be much higher in infants; they had an Epo receptor capacity of 133 pM vs 86.6 pM in adults, which is most likely due to the higher erythroid progenitor cell mass per kilogram of body weight in infants. The parallel linear elimination was found to be more dominant in adults, reaching 91% of the total clearance with a 500-U/kg dose compared to just 6.1% of the total clearance following the same dose in preterm infants. Thus, this mechanism-based population PK model revealed that receptor-based nonlinear elimination is the dominant Epo elimination pathway in premature infants, and parallel linear elimination is dominant in adults.

Semi-mechanistic Pharmacokinetic/Pharmacodynamic model of three pegylated rHuEPO and ior®EPOCIM in New Zealand rabbits

Marketed formulations of erythropoietin (EPO) ior®EPOCIM, MIRCERA® and two newly developed pegylated-EPO analogues (PEG-EPO 32 and 40 kDa) formulations were intravenously administered to New Zealand rabbits. A semi-mechanistic Pharmacokinetic/Pharmacodynamic (PK/PD) model describing in a simultaneous and integrated form the time course of reticulocytes, red blood cells and hemoglobin was built to account for the time course of hematopoiesis stimulation after erythropoietin administration. Data analysis was performed based on the population approach with the software NONMEM version 7.3. Erythropoietin disposition of each of the administered formulations was best described with a two compartment model and linear elimination. Different formulations show different clearance and apparent volume of distribution of the central compartment but share estimates of inter-compartmental clearance and apparent peripheral volume of distribution. A semi-mechanistic model including cell proliferation, maturation, and homeostatic regulation provided a good description of the data regardless the type of erythropoietin formulation administered. The system-, and drug-related parameters showed consistency and differed across formulations, respectively. A single IV administration of PEG-EPO 32 and 40 kDa formulations in New Zealand rabbits achieves a median change of 27% and 22% on RET levels, and of 47% and 63% on RBC and HGB levels, respectively compared to MIRCERA®. The administration of new branched PEG-chains formulations improves PK and PD properties of EPO, in terms of increasing elimination half-lives and pharmacological activity on RET, RBC and HGB compared to commercially available formulations (ior®EPOCIM and MIRCERA®).

Pharmacokinetics and Tolerability of a Dual Variable Domain Immunoglobulin ABT-981 Against IL-1α and IL-1β in Healthy Subjects and Patients With Osteoarthritis of the Knee

The interleukin (IL)-1 family of proinflammatory cytokines are thought to play a significant role in the structural progression of osteoarthritis and its associated symptoms. IL-1α and IL-1β are 2 distinct cytokines found in the cartilage, synovial membrane, and synovial fluid of patients with osteoarthritis. The aim of these studies was to evaluate the pharmacokinetics of ABT-981, a dual variable domain immunoglobulin (DVD-Ig) capable of simultaneously binding IL-1α and IL-1β, in healthy subjects and patients with osteoarthritis of the knee. Fifty-six healthy adult subjects were randomized to receive single doses of ABT-981 intravenously (0.3, 1, 3, or 10 mg/kg), subcutaneously (0.3, 1, 3 mg/kg), or matching placebo in a 3:1 ratio. Thirty-six patients with osteoarthritis of the knee were randomized to receive 4 subcutaneous ABT-981 doses of 0.3, 1, or 3 mg/kg administered every 2 weeks, 3 subcutaneous doses of ABT-981 3 mg/kg every 4 weeks, or matching placebo in a 7:2 active:placebo ratio. ABT-981 behaved similarly to conventional monoclonal antibodies following single or multiple doses with mean maximum serum concentrations 2 to 9 days after subcutaneous doses, mean terminal half-lives of 10 to 14 days, and an absolute subcutaneous bioavailability of 46%. Exposure of ABT-981 was approximately linear following single or multiple doses every 2 weeks with monoexponential decline of terminal-phase concentrations. The most common adverse events associated with ABT-981 were diarrhea and headache in healthy subjects and injection site erythema in subjects with osteoarthritis of the knee. Decreased absolute neutrophil counts were observed in response to ABT-981 administration.

Population pharmacokinetic drug-drug interaction pooled analysis of existing data for rifabutin and HIV PIs

OBJECTIVES

Extensive but fragmented data from existing studies were used to describe the drug-drug interaction between rifabutin and HIV PIs and predict doses achieving recommended therapeutic exposure for rifabutin in patients with HIV-associated TB, with concurrently administered PIs.

METHODS

Individual-level data from 13 published studies were pooled and a population analysis approach was used to develop a pharmacokinetic model for rifabutin, its main active metabolite 25-O-desacetyl rifabutin (des-rifabutin) and drug-drug interaction with PIs in healthy volunteers and patients who had HIV and TB (TB/HIV).

RESULTS

Key parameters of rifabutin affected by drug-drug interaction in TB/HIV were clearance to routes other than des-rifabutin (reduced by 76%-100%), formation of the metabolite (increased by 224% in patients), volume of distribution (increased by 606%) and distribution to the peripheral compartment (reduced by 47%). For des-rifabutin, clearance was reduced by 35%-76% and volume of distribution increased by 67%-240% in TB/HIV. These changes resulted in overall increased exposure to rifabutin in TB/HIV patients by 210% because of the effects of PIs and 280% with ritonavir-boosted PIs.

CONCLUSIONS

Given together with non-boosted or ritonavir-boosted PIs, rifabutin at 150 mg once daily results in similar or higher exposure compared with rifabutin at 300 mg once daily without concomitant PIs and may achieve peak concentrations within an acceptable therapeutic range. Although 300 mg of rifabutin every 3 days with boosted PI achieves an average equivalent exposure, intermittent doses of rifamycins are not supported by current guidelines.

Erythropoietin Stimulates Tumor Growth via EphB4

While recombinant human erythropoietin (rhEpo) has been widely used to treat anemia in cancer patients, concerns about its adverse effects on patient survival have emerged. A lack of correlation between expression of the canonical EpoR and rhEpo's effects on cancer cells prompted us to consider the existence of an alternative Epo receptor. Here, we identified EphB4 as an Epo receptor that triggers downstream signaling via STAT3 and promotes rhEpo-induced tumor growth and progression. In human ovarian and breast cancer samples, expression of EphB4 rather than the canonical EpoR correlated with decreased disease-specific survival in rhEpo-treated patients. These results identify EphB4 as a critical mediator of erythropoietin-induced tumor progression and further provide clinically significant dimension to the biology of erythropoietin.

Quantitative pharmacology of denosumab in patients with bone metastases from solid tumors

Denosumab (XGEVA®) is a recombinant, fully human IgG2 monoclonal antibody directed against the receptor activator of nuclear factor kappa-B ligand (RANKL) that prevents differentiation of osteoclast precursors into mature osteoclasts and acceleration of bone resorption, resulting in the inhibition of osteoclast activation. Denosumab is indicated for the prevention of skeletal-related events (SREs) in adult patients with bone metastases from solid tumors at the dose of 120 mg administered subcutaneously (SC) every 4 weeks. This review is focused on describing its target-mediated disposition and direct inhibitory effect on bone resorption, as well as the modeling and simulation techniques used to integrate the PKPD information collected during clinical development of denosumab. In addition, this review further discusses the clinical relevance of patient covariate effects on denosumab systemic exposure, target engagement and downstream pharmacodynamics biomarkers, and the rationale for dosing regimen selection for Phase 3 studies. Phase 3 clinical studies demonstrated that denosumab was superior to zoledronic acid in inhibiting bone resorption and, consequently, delaying the time to first SRE by a median of 8.2 months in patients with bone metastases from solid tumors. Thus, denosumab may be considered a better alternative treatment than zoledronic acid for the prevention of SRE in patients with bone metastases from solid tumors.

A quantitative approach to analysing cortisol response in the horse

The cortisol response to glucocorticoid intervention has, in spite of several studies in horses, not been fully characterized with regard to the determinants of onset, intensity and duration of response. Therefore, dexamethasone and cortisol response data were collected in a study applying a constant rate infusion regimen of dexamethasone (0.17, 1.7 and 17 μg/kg) to six Standardbreds. Plasma was analysed for dexamethasone and cortisol concentrations using UHPLC-MS/MS. Dexamethasone displayed linear kinetics within the concentration range studied. A turnover model of oscillatory behaviour accurately mimicked cortisol data. The mean baseline concentration range was 34-57 μg/L, the fractional turnover rate 0.47-1.5 1/h, the amplitude parameter 6.8-24 μg/L, the maximum inhibitory capacity 0.77-0.97, the drug potency 6-65 ng/L and the sigmoidicity factor 0.7-30. This analysis provided a better understanding of the time course of the cortisol response in horses. This includes baseline variability within and between horses and determinants of the equilibrium concentration-response relationship. The analysis also challenged a protocol for a dexamethasone suppression test design and indicated future improvement to increase the predictability of the test.

A Bottom-Up Whole-Body Physiologically Based Pharmacokinetic Model to Mechanistically Predict Tissue Distribution and the Rate of Subcutaneous Absorption of Therapeutic Proteins

The ability to predict subcutaneous (SC) absorption rate and tissue distribution of therapeutic proteins (TPs) using a bottom-up approach is highly desirable early in the drug development process prior to clinical data being available. A whole-body physiologically based pharmacokinetic (PBPK) model, requiring only a few drug parameters, to predict plasma and interstitial fluid concentrations of TPs in humans after intravenous and subcutaneous dosing has been developed. Movement of TPs between vascular and interstitial spaces was described by considering both convection and diffusion processes using a 2-pore framework. The model was optimised using a variety of literature sources, such as tissue lymph/plasma concentration ratios in humans and animals, information on the percentage of dose absorbed following SC dosing via lymph in animals and data showing loss of radiolabelled IgG from the SC dosing site in humans. The resultant model was used to predict t_max and plasma concentration profiles for 12 TPs (molecular weight 8–150 kDa) following SC dosing. The predicted plasma concentration profiles were generally comparable to observed data. t_max was predicted within 3-fold of reported values, with one third of the predictions within 0.8–1.25-fold. There was no systematic bias in simulated C_max values, although a general trend for underprediction of t_max was observed. No clear trend between prediction accuracy of t_max and TP isoelectric point or molecular size was apparent. The mechanistic whole-body PBPK model described here can be applied to predict absorption rate of TPs into blood and movement into target tissues following SC dosing.

Assessment of hemoglobin responsiveness to epoetin alfa in patients on hemodialysis using a population pharmacokinetic pharmacodynamic model

A population pharmacokinetic pharmacodynamic (PK/PD) model describing the effect of epoetin alfa on hemoglobin (Hb) response in hemodialysis patients was developed. Epoetin alfa pharmacokinetics was described using a linear 2-compartment model. PK parameter estimates were similar to previously reported values. A maturation-structured cytokinetic model consisting of 5 compartments linked in a catenary fashion by first-order cell transfer rates following a zero-order input process described the Hb time course. The PD model described 2 subpopulations, one whose Hb response reflected epoetin alfa dosing and a second whose response was unrelated to epoetin alfa dosing. Parameter estimates from the PK/PD model were physiologically reasonable and consistent with published reports. Numerical and visual predictive checks using data from 2 studies were performed. The PK and PD of epoetin alfa were well described by the model.

A Tutorial on Target-Mediated Drug Disposition (TMDD) Models

Target-mediated drug disposition (TMDD) is the phenomenon in which a drug binds with high affinity to its pharmacological target site (such as a receptor) to such an extent that this affects its pharmacokinetic characteristics.1 The aim of this Tutorial is to provide an introductory guide to the mathematical aspects of TMDD models for pharmaceutical researchers. Examples of Berkeley Madonna2 code for some models discussed in this Tutorial are provided in the Supplementary Materials.

Effect of Type 2 Diabetes Mellitus and Diabetic Nephropathy on IgG Pharmacokinetics and Subcutaneous Bioavailability in the Rat

The objective of this research was to assess the effects of type 2 diabetes mellitus (T2DM) and diabetic nephropathy (DN) on the pharmacokinetics of human IgG (hIgG), an antibody isotype, in Zucker diabetic fatty (ZDF) rats. Furthermore, the specific role of T2DM in the altered disposition of hIgG was evaluated by treating diabetic rats with pioglitazone, while the role of chronic kidney disease (CKD) was assessed using 5/6 nephrectomized Sprague Dawley rats. ZDF male (lean non-diabetic control and obese diabetic) and pioglitazone-treated ZDF rats were studied at ages 12-13 weeks (only DM was present), and at ages 29-30 weeks (progression to DN). All animals were dosed with 1 mg/kg of hIgG intravenously (IV) or subcutaneously (SC). ZDF rats had significantly higher blood glucose concentrations and urinary albumin excretion compared to control rats. Significant increases in total clearance (2.5-fold) and renal clearance (100-fold) of hIgG were observed; however the major increase in total clearance was due to increased non-renal clearance. Greater changes in urinary albumin excretion and total and renal clearances of IgG (3.5-fold and 300-fold, respectively) were observed with progression to DN. SC bioavailability of hIgG in all animal groups was similar (>84%). With pioglitazone-treatment, diabetic animals remained euglycemic and treatment was able to reverse the clearance changes, although incompletely. In the CKD group, no difference in hIgG clearance was observed when compared with controls. In conclusion, the increased clearance of hIgG in ZDF diabetic animals, reversal by pioglitazone treatment and lack of effect of CKD, demonstrate the influence of T2DM on hIgG pharmacokinetics.

Clinical pharmacokinetics and pharmacodynamics of erythropoiesis-stimulating agents

The cloning of the EPO gene in the early 1980s allowed for the development of recombinant erythropoietins and analogues [erythropoiesis-stimulating agents (ESAs)], offering an alternative to transfusion as a method of raising haemoglobin (Hb) levels, which have been used for more than 20 years to treat anaemia in millions of anaemic patients. There are now a number of ESAs available worldwide for the treatment of anaemia, approved for different routes of administration (intravenous and subcutaneous) and dosing intervals (three times weekly, weekly, biweekly and monthly). In this review, we discuss the pharmacokinetic characteristics, including absorption, distribution and elimination processes, across the different ESAs. Incomplete and slow lymphatic absorption, with limited extravascular distribution, and minor contributions of the target-mediated drug disposition to the overall elimination are the common characteristics across the marketed ESA. Additionally, we assess the similarities and differences of ESAs related to pharmacodynamics in the context of the different biomarkers used to monitor the magnitude and duration of the effect, and introduce the concept of the minimum effective concentration of the ESA. The relationship between the minimum effective concentration and the half-life suggests that the time during which drug concentrations are above the minimum effective concentration is the main determinant of ESA efficacy in increasing Hb levels. The tolerance phenomenon and its physiological mechanism and implications for ESA dosing are discussed. Finally, the areas of future clinical pharmacology research are envisioned.

A Population Pharmacokinetic and Pharmacodynamic Analysis of Peginesatide in Patients with Chronic Kidney Disease on Dialysis

Peginesatide (OMONTYS®) is an erythropoiesis-stimulating agent that was indicated in the United States for the treatment of anemia due to chronic kidney disease in adult patients on dialysis prior to its recent marketing withdrawal by the manufacturer. The objective of this analysis was to develop a population pharmacokinetic and pharmacodynamic model to characterize the time-course of peginesatide plasma and hemoglobin concentrations following intravenous and subcutaneous administration. Plasma samples (n = 2,665) from 672 patients with chronic kidney disease (on or not on dialysis) and hemoglobin samples (n = 18,857) from 517 hemodialysis patients (subset of the 672 patients), were used for pharmacokinetic-pharmacodynamic model development in NONMEM VI. The pharmacokinetic profile of peginesatide was best described by a two-compartment model with first-order absorption and saturable elimination. The relationship between peginesatide and hemoglobin plasma concentrations was best characterized by a modified precursor-dependent lifespan indirect response model. The estimate of maximal stimulatory effect of peginesatide on the endogenous production rate of progenitor cells (Emax) was 0.54. The estimate of peginesatide drug concentration required for 50% of maximal response (EC50) estimates was 0.4 µg/mL. Several significant (P<0.005) covariates affected simulated peginesatide exposure by ≤36%. Based upon ≤0.2 g/dL effects on simulated hemoglobin levels, none were considered clinically relevant.

Population pharmacokinetic meta-analysis of denosumab in healthy subjects and postmenopausal women with osteopenia or osteoporosis

BACKGROUND AND OBJECTIVE

Inhibition of the receptor activator of nuclear factor κ-B ligand (RANKL) is a therapeutic target for treatment of bone disorders associated with increased bone resorption, such as osteoporosis. The objective of this analysis was to characterize the population pharmacokinetics of denosumab (AMG 162; Prolia®), a fully human IgG2 monoclonal antibody that binds to RANKL, in healthy subjects and postmenopausal women with osteopenia or osteoporosis.

METHODS

A total of 22944 serum free denosumab concentrations from 495 healthy subjects and 1069 postmenopausal women with osteopenia or osteoporosis were pooled. Denosumab was administered as either a single intravenous dose (n = 36), a single subcutaneous dose (n = 469) or multiple subcutaneous doses (n = 1059), ranging from 0.01 to 3 mg/kg (or 6-210 mg as fixed mass dosages), every 3 or 6 months for up to 48 months. An open, two-compartment pharmacokinetic model with a quasi-steady-state approximation of the target-mediated drug disposition model was used to describe denosumab pharmacokinetics, using NONMEM Version 7.1.0 software. Subcutaneous absorption was characterized by the first-order absorption rate constant (k(a)), with constant absolute bioavailability over the range of doses that were evaluated. Clearance and volume of distribution parameters were scaled by body weight, using a power model. Model evaluation was performed through visual predictive checks.

RESULTS

The subcutaneous bioavailability of denosumab was 64%, and the k(a) was 0.00883 h-1. The central volume of distribution and linear clearance were 2.49 L/66 kg and 3.06 mL/h/66 kg, respectively. The baseline RANKL level, quasi-steady-state constant and RANKL degradation rate were 614 ng/mL, 138 ng/mL and 0.00148 h-1, respectively. Between-subject variability in model parameters was moderate. A fixed dose of 60 mg provided RANKL inhibition similar to that achieved by equivalent body weight-based dosing. The effects of age and race on the area under the serum concentration-time curve of denosumab were less than 15% over the range of covariate values that were evaluated.

CONCLUSIONS

The non-linearity in denosumab pharmacokinetics is probably due to RANKL binding, and denosumab dose adjustment based on the patient demographics is not warranted.

Pharmacokinetic similarity of biologics: analysis using nonlinear mixed-effects modeling

Our objective was to show, using two examples, that a pharmacokinetic (PK) similarity analysis can be performed using nonlinear mixed-effects models (NLMEM). We used two studies that compared different biosimilars: a three-way crossover trial with somatropin and a parallel-group trial with epoetin-α. For both data sets, the results of NLMEM-based analysis were compared with those of noncompartmental analysis (NCA). For the latter analysis, we performed an NLMEM-based equivalence Wald test on secondary parameters of the model: the area under the curve and the maximal concentration. Somatropin PK was described by a one-compartment model and epoetin-α PK by a two-compartment model with linear and Michaelis-Menten elimination. For both studies, similarity of PK was demonstrated by means of both NCA and NLMEM, and both methods led to similar results. Therefore, for establishing similarity, PK data can be analyzed by either of the methods. NCA is an easier approach because it does not require data modeling; however, NLMEM leads to a better understanding of the underlying biological system.

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.

Exposure-response modeling of darbepoetin alfa in anemic patients with chronic kidney disease not receiving dialysis

A population pharmacokinetic and pharmacodynamic model (PK/PD) of darbepoetin alfa following intravenous (IV) or subcutaneous (SC) administration in participants with chronic kidney disease (CKD) was developed. Darbepoetin alfa concentrations from 96 CKD participants, who received IV or SC darbepoetin alfa, and Hgb concentration from 332 CKD participants not on dialysis, who received SC doses of darbepoetin alfa, were used to develop the PK/PD model. An open 2-compartment model with sequential zero- and first-order absorption was used to characterize darbepoetin alfa pharmacokinetics. Darbepoetin alfa was assumed to trigger concentration-dependent stimulation of production of progenitor cells of red blood cells (RBCs) in bone marrow, which become red blood cells and died after life span expiration. Model evaluation was performed through nonparametric bootstrap and posterior predictive checks. Absolute bioavailability, total mean absorption time, clearance, and volume of distribution were estimated to be 44%, 52 h, 3.4 L/d/70 kg, and 5.9 L/70 kg, respectively. The estimates of drug potency, efficacy, and RBC life span were 0.41 ng/mL, 64%, and 77 days, respectively. Pharmacokinetic or pharmacodynamic parameters of darbepoetin alfa were not affected by age and sex. The qualified model supports the use of darbepoetin alfa administered biweekly (SC) in CKD patients for anemia correction and monthly (SC) for hemoglobin maintenance. In addition, the model is deemed appropriate to conduct simulations to support dose selection for additional clinical studies.

Target-mediated drug disposition model: approximations, identifiability of model parameters and applications to the population pharmacokinetic-pharmacodynamic modeling of biologics

Models for drugs exhibiting target-mediated drug disposition (TMDD) describe biological processes in which drug-target binding significantly influences both pharmacodynamics (PD) and pharmacokinetics (PK). TMDD models are often over-parameterized and their parameters are difficult to estimate based on available data. Approximations of the general model have been suggested, but even these simpler forms can be over-parameterized when, for example, target and drug-target complex concentrations are not available. This work i) reviews TMDD equations, their approximations and methods to study identifiability of model parameters; ii) reviews the publications that used TMDD equations to describe PK and PD of biologics; and iii) discusses issues of identifiability of the TMDD model parameters related to study design and data analysis. Examples demonstrate that use of the TMDD equations for the population PK and PD modeling is most successful when the target and drug-target complex concentrations are available in addition to the drug concentration data. TMDD parameter estimates can be trusted only when they are identifiable, that is, can be estimated from the available data with sufficient precision. Parameter identifiability analysis should be an integral part of the TMDD system investigation. It also should be used prospectively for optimal study design.

Pharmacodynamic analysis of recombinant human erythropoietin effect on reticulocyte production rate and age distribution in healthy subjects

OBJECTIVE

To evaluate the effect of recombinant human erythropoietin (rHuEPO) on the reticulocyte production rate and age distribution in healthy subjects.

METHODS

Extensive pharmacokinetic and pharmacodynamic data collected from 88 subjects who received a single subcutaneous dose of rHuEPO (dose range 20-160 kIU) were analysed. Four nonlinear mixed-effects models were evaluated to describe the time course of the percentage of reticulocytes and their age distribution in relation to rHuEPO pharmacokinetics. Model A accounted for stimulation of the production of progenitor cells in bone marrow, and model B implemented shortening of differentiation and maturation times of early progenitors in bone marrow. Model C was the combination of models A and B, and model D was the combination of model A with an increase in the maturation times of the circulating reticulocytes. Model evaluation was performed using goodness-of-fit plots, a nonparametric bootstrap and a posterior predictive check.

RESULTS

Model D was selected as the best model, and evidenced accurate and precise estimation of model parameters and prediction of the time course of the percentage of reticulocytes. At baseline, the estimated circulating reticulocyte maturation time was 2.6 days, whereas the lifespan of the precursors in the bone marrow was about 5 days. The rHuEPO potency for the stimulatory effect (7.61 IU/L) was higher than that for the increase in reticulocyte maturation times (56.3 IU/L). There was a significant 1- to 2-day lag time in the reticulocyte response. The effect of rHuEPO on the reticulocyte age distribution consisted of a transient increase in the reticulocyte maturation time from baseline up to 6-7 days, occurring 1 day after administration. The dose-dependent amplitude of the changes in the age distribution lasted for 12-14 days. The model-predicted peak increase in the reticulocyte release rate ranged from 140% to 160% of the baseline value and was maximal on days 7-8 following rHuEPO administration.

CONCLUSIONS

A semiphysiological model quantifying the effect of rHuEPO on the reticulocyte production rate and age distribution was developed. The validated model predicts that rHuEPO increases the reticulocyte production rate and modifies the reticulocyte age distribution in a dose-dependent manner.

Basic pharmacodynamic models for agents that alter the lifespan distribution of natural cells

A new class of basic indirect pharmacodynamic models for agents that alter the loss of natural cells based on a lifespan concept are presented. The lifespan indirect response (LIDR) models assume that cells (R) are produced at a constant rate (k(in)), survive during a certain duration T(R), and finally are lost. The rate of cell loss is equal to the production rate but is delayed by T(R). A therapeutic agent can increase or decrease the baseline cell lifespan to a new cell lifespan, T(D), by temporally changing the proportion of cells belonging to the two modes of the lifespan distribution. Therefore, the change of lifespan at time t is described according to the Hill function, H(C(t)), with capacity (E(max)) and sensitivity (EC(50)), and the pharmacokinetic function C(t). A one-compartment cell model was examined through simulations to describe the role of pharmacokinetics, pharmacodynamics and cell properties for the cases where the drug increases (T(D) > T(R)) or decreases (T(D) < T(R)) the cell lifespan. The area under the effect curve (AUCE) and explicit solutions of LIDR models for large doses were derived. The applicability of the model was further illustrated using the effects of recombinant human erythropoietin (rHuEPO) on reticulocytes. The cases of both stimulation of the proliferation of bone marrow progenitor cells and the increase of reticulocyte lifespans were used to describe mean data from healthy subjects who received single subcutaneous doses of rHuEPO ranging from 20 to 160 kIU. rHuEPO is about 4.5-fold less potent in increasing reticulocyte survival than in stimulating the precursor production. A maximum increase of 4.1 days in the mean reticulocyte lifespan was estimated and the effect duration on the lifespan distribution was dose dependent. LIDR models share similar properties with basic indirect response models describing drug stimulation or inhibition of the response loss rate with the exception of the presence of a lag time and a dose independent peak time. The current concept can be applied to describe the pharmacodynamic effects of agents affecting survival of hematopoietic cell populations yielding realistic physiological parameters.

Population pharmacokinetic modeling of epoetin delta in pediatric patients with chronic kidney disease

This analysis quantifies the population pharmacokinetics of subcutaneous and intravenous epoetin delta, an epoetin produced in a human cell line, in pediatric patients with chronic kidney disease and estimates the effects of covariate factors on epoetin delta and epoetin alfa pharmacokinetic parameters. Erythropoietin serum concentration data, taken from a phase III study conducted in 60 patients aged 1 to 17 years, were best described by a 1-compartment model with first-order absorption and elimination. The typical point estimates were clearance (0.268 L/h), central volume of distribution (1.03 L), absorption rate constant (0.0554 h(-1)), and bioavailability (0.708) for a 35-kg male < or = 10 years who was predialysis and on subcutaneous epoetin delta treatment. Erythropoietin pharmacokinetic parameters were similar in pediatric patients as compared with adults when scaled by weight. The subcutaneous administration of epoetin alfa exhibited lower systemic bioavailability than subcutaneous administration of epoetin delta.

Pharmacokinetic and pharmacodynamic profiles of extended dosing of epoetin alfa in anemic patients who have chronic kidney disease and are not on dialysis

BACKGROUND AND OBJECTIVES

Emerging evidence suggests that epoetin alfa can be administered at extended intervals of up to 4 wk. This open-label, randomized study was performed to characterize the pharmacokinetic and pharmacodynamic profiles of four dosing regimens of epoetin alfa administered subcutaneously in anemic patients who had chronic kidney disease and were not on dialysis.

DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS

Thirty-eight patients, enrolled from nine centers in the United States, were > or =18 yr of age and had hemoglobin <11.0 g/dl and GFR 12 to 60 ml/min per 1.73 m(2). Patients received one of four epoetin alfa dosing regimens: 50 IU/kg three times per week, 10,000 IU once weekly, or 20,000 IU every 2 wk for 36 d or 40,000 IU every 4 wk for 64 d. Each regimen provided a similar dosage of epoetin alfa over 4 wk. Dosage adjustments were not permitted.

RESULTS

Drug exposure to epoetin alfa over 4 wk, based on area under the curve, was somewhat higher with the extended interval regimens compared with the three-times-weekly regimen. Mean change in hemoglobin during the study period was similar for all regimens. No patients were transfused. Three patients experienced five serious adverse events, none of which was considered treatment related.

CONCLUSIONS

Extended dosing interval regimens of epoetin alfa yielded modest pharmacokinetic differences but a similar pharmacodynamic response, suggesting that less frequent, higher dosages of epoetin alfa may be as effective as the current three-times-weekly regimen in anemic patients who have chronic kidney disease and are not on dialysis.

Simulating pharmacokinetic and pharmacodynamic fuzzy-parameterized models: a comparison of numerical methods

Statistical techniques have been traditionally used to deal with parametric variation in pharmacokinetic and pharmacodynamic models, but these require substantial data for estimates of probability distributions. In the presence of limited, inaccurate or imprecise information, simulation with fuzzy numbers represents an alternative tool to handle parametric uncertainty. Existing methods for implementing fuzzy arithmetic may, however, have significant shortcomings in overestimating (e.g., conventional fuzzy arithmetic) and underestimating (e.g., vertex method) the output uncertainty. The purpose of the present study is to apply and compare the applicability of conventional fuzzy arithmetic, vertex method and two recently proposed numerical schemes, namely transformation and optimization methods, for uncertainty modeling in pharmacokinetic and pharmacodynamic fuzzy-parameterized systems. A series of test problems were examined, including empirical pharmacokinetic and pharmacodynamic models, a function non-monotonic in its parameters, and a whole body physiologically based pharmacokinetic model. Our results verified that conventional fuzzy arithmetic can lead to overestimation of response uncertainty and should be avoided. For the monotonic pharmacokinetic and pharmacodynamic models, the vertex method accurately predicted fuzzy-valued output while incurring the least computational cost. It turned out that the choice of a suitable method for fuzzy simulation of the non-monotonic function depended on the required accuracy of the results: the vertex method was capable of eliciting an initial approximate solution with few function evaluations; for more accurate results, the transformation method was the most superior approach in terms of accuracy per unit CPU time.