Characterization of anti-drug antibody dynamics using a bivariate mixed hidden-markov model by nonlinear-mixed effects approach.

Biological therapies may act as immunogenic triggers leading to the formation of anti-drug antibodies (ADAs). Population pharmacokinetic (PK) models can be used to characterize the relationship between ADA and drug disposition but often rely on the ADA bioassay results, which may not be sufficiently sensitive to inform on this characterization.In this work, a methodology that could help to further elucidate the underlying ADA production and impact on the drug disposition was explored. A mixed hidden-Markov model (MHMM) was developed to characterize the underlying (hidden) formation of ADA against the biologic, using certolizumab pegol (CZP), as a test drug. CZP is a PEGylated Fc free TNF-inhibitor used in the treatment of rheumatoid arthritis and other chronic inflammatory diseases.The bivariate MHMM used information from plasma drug concentrations and ADA measurements, from six clinical studies (n = 845), that were correlated through a bivariate Gaussian function to infer about two hidden states; production and no-production of ADA influencing PK. Estimation of inter-individual variability was not supported in this case. Parameters associated with the observed part of the model were reasonably well estimated while parameters associated with the hidden part were less precise. Individual state sequences obtained using a Viterbi algorithm suggested that the model was able to determine the start of ADA production for each individual, being a more assay-independent methodology than traditional population PK. The model serves as a basis for identification of covariates influencing the ADA formation, and thus has the potential to identify aspects that minimize its impact on PK and/or efficacy.

Population pharmacokinetics and pharmacodynamics of efmarodocokin alfa (IL-22Fc).

Efmarodocokin alfa (IL-22Fc) is a fusion protein of human IL-22 linked to the crystallizable fragment (Fc) of human IgG4. It has been tested in multiple indications including inflammatory bowel disease (IBD). The purposes of the present analyses were to describe the population pharmacokinetics (PK) of efmarodocokin alfa and perform pharmacodynamic (PD) analysis on the longitudinal changes of the PD biomarker REG3A after efmarodocokin alfa treatment as well as identify covariates that affect efmarodocokin alfa PK and REG3A PD. The data used for this analysis included 182 subjects treated with efmarodocokin alfa in two clinical studies. The population PK and PD analyses were conducted sequentially. Efmarodocokin alfa concentration-time data were analyzed using a nonlinear mixed-effects modeling approach, and an indirect response model was adopted to describe the REG3A PD data with efmarodocokin alfa serum concentration linked to the increase in REG3A. The analysis software used were NONMEM and R. A 3-compartment model with linear elimination best described the PK of efmarodocokin alfa. The estimated population-typical value for clearance (CL) was 1.12 L/day, and volume of central compartment was 6.15 L. Efmarodocokin alfa CL increased with higher baseline body weight, C-reactive protein, and CL was 27.6% higher in IBD patients compared to healthy subjects. The indirect response PD model adequately described the longitudinal changes of REG3A after efmarodocokin alfa treatment. A popPK and PD model for efmarodocokin alfa and REG3A was developed and covariates affecting the PK and PD were identified.

Sensitivity of Pegfilgrastim Pharmacokinetic and Pharmacodynamic Parameters to Product Differences in Similarity Studies.

In this work, a previously developed pegfilgrastim (PG) population pharmacokinetic-pharmacodynamic (PKPD) model was used to evaluate potential factors of importance in the assessment of PG PK and PD similarity. Absolute neutrophil count (ANC) was the modelled PD variable. A two-way cross-over study was simulated where a reference PG and a potentially biosimilar test product were administered to healthy volunteers. Differences in delivered dose amounts or potency between the products were simulated. A different baseline absolute neutrophil count (ANC) was also considered. Additionally, the power to conclude PK or PD similarity based on areas under the PG concentration-time curve (AUC) and ANC-time curve (AUEC) were calculated. Delivered dose differences between the products led to a greater than dose proportional differences in AUC but not in AUEC, respectively. A 10% dose difference from a 6 mg dose resulted in 51% and 7% differences in AUC and AUEC, respectively. These differences were more pronounced with low baseline ANC. Potency differences up to 50% were not associated with large differences in either AUCs or AUECs. The power to conclude PK similarity was affected by the simulated dose difference; with a 4% dose difference from 6 mg the power was approximately 29% with 250 subjects. The power to conclude PD similarity was high for all delivered dose differences and sample sizes.

Model-based evaluation of low-dose factor VIII prophylaxis in haemophilia A.

INTRODUCTION: The optimal treatment modality for haemophilia A is lifelong prophylaxis which is expensive and may not be implementable everywhere where factor VIII (FVIII) availability is limited. A less costly alternative to prophylaxis is low-dose prophylaxis (LDP) which was compared to conventional prophylaxis in this model-based simulation study. AIM: To explore whether LDP is motivated where standard prophylaxis is not implementable, including evaluating LDP efficacy compared to high-dose prophylaxis and investigating the potential economic benefit of individualized dosing. METHODS: For a virtual adult haemophilia A population, FVIII activity levels were simulated following alternative treatment regimens, based on a published population PK model. The regimens included very LDP, LDP and conventional prophylaxis twice and thrice weekly. The annual probability of bleeding was predicted based on the weekly time spent below 1 IU/dL, using a previously published relationship. Additionally, PK-based dose individualization was evaluated to determine FVIII savings using Bayesian forecasting. RESULTS: A treatment regimen of 10 IU/kg administered thrice weekly cost 75% less than a standard high-dose regimen and was predicted to have a 5% higher median probability of annual bleeds. PK-based dose individualization may result in further cost-savings, but implementation needs benefit versus feasibility consideration. CONCLUSION: Based on simulations, a promising LDP regimen was identified that decreased treatment costs compared with standard high-dose prophylaxis at a small increase in bleeding risk. The results indicate that LDP is advocated where the standard-of-care is on-demand treatment; however, the results should be considered in the context of any limitations of the applied models.

A Population Pharmacokinetic-Pharmacodynamic Model of Pegfilgrastim.

Neutropenia and febrile neutropenia (FN) are serious side effects of cytotoxic chemotherapy which may be alleviated with the administration of recombinant granulocyte colony-stimulating factor (GCSF) derivatives, such as pegfilgrastim (PG) which increases absolute neutrophil count (ANC). In this work, a population pharmacokinetic-pharmacodynamic (PKPD) model was developed based on data obtained from healthy volunteers receiving multiple administrations of PG. The developed model was a bidirectional PKPD model, where PG stimulated the proliferation, maturation, and margination of neutrophils and where circulating neutrophils in turn increased the elimination of PG. Simulations from the developed model show disproportionate changes in response with changes in dose. A dose increase of 10% from the 6 mg therapeutic dose taken as a reference leads to area under the curve (AUC) increases of ~50 and ~5% for PK and PD, respectively. A full random effects covariate model showed that little of the parameter variability could be explained by sex, age, body size, and race. As a consequence, little of the secondary parameter variability (Cmax and AUC of PG and ANC) could be explained by these covariates.