CIC-0901 Inserm Necker-Cochin (Assistance Publique-Hopitaux de Paris), Paris EA-3620, France

Université Paris Descartes Sorbonne Paris Cité et Institut IMAGINE, Paris, France

Pharmacologie, Hôpital Cochin, Assistance Publique-Hopitaux de Paris, Paris, France

Unité d’Immuno-Hématologie et Rhumatologie pédiatriques, Hôpital Necker-Enfants Malades, Assistance Publique-Hopitaux de Paris, Paris, France

Abstract

Background

Anakinra pharmacokinetics and pharmacodynamics were investigated in children and adolescents treated for systemic-onset juvenile idiopathic arthritis (SJIA) and autoinflammatory syndromes.

Methods

Anakinra was given subcutaneously at doses between 2 and 10 mg/kg (maximum 100 mg) per day. Anakinra concentrations were recorded in patients, as well as C-reactive protein (CRP) levels, on different occasions. The data were fitted to a pharmacokinetic-pharmacodynamic model via a population approach using Monolix.

Results

A total of 87 children and adolescents, 8 months to 21 years old, were available for pharmacokinetic evaluation. A one compartment model with linear absorption and elimination described the pharmacokinetics. Taking into account bodyweight to explain variations in apparent clearance (CL/F) and distribution volume (V/F) significantly reduced the associated between-subject and between-occasion variabilities. The final estimates were 6.24 L/h/70 kg and 65.2 L/70 kg for CL/F and V/F respectively. A mixture pharmacodynamic model described the CRP level change during anakinra treatment for the SJIA patients with 2 subpopulations, patients with high baseline and large CRP decrease and patients with low baseline and small CRP decrease followed by a re-increase in CRP levels. There was no significant effect of the combined anti-inflammatory treatment. The proportion of patients for which the development of a resistance to treatment was significant was 62% and the corresponding time was approximately 60 days.

Conclusions

Based on effects in SJIA, a prospective dosage adjustment was proposed based on a 0.4 mg/L Css target in order to obtain a CRP decrease to 10 mg/L or below.

Background

Anakinra, a recombinant nonglycosylated homolog of human IL-1 receptor antagonist, competitively inhibits the binding of IL-1α and IL-1β to the IL-1 receptor and thus inhibits the effects of this pro-inflammatory cytokine. Il-1 blockade using selective IL-1β blockade (canakinumab) or IL-1 α and β blockade (anakinra, rilonacept) has proven efficacy in cryopyrin-associated periodic syndromes (CAPS)

The aim of the present study was to investigate anakinra pharmacokinetics in children and young adult patients with SJIA and autoinflammatory conditions. The effects of anakinra on CRP was also modelled in relation to the anakinra pharmacokinetics in children with SJIA.

Methods

Patients

Following IRB (CPP Paris V) approval and patients and/or parents (for children) written informed consent, pharmacokinetic data were obtained from a phase IIB trial testing anakinra in SJIA patients (ANAJIS)(5) and from patients subsequently treated for diverse autoinflammatory conditions. In ANAJIS trial all the patients received anakinra once-a-day at the dose of 2 mg/kg subcutaneously, maximum 100 mg. In patients with autoinflammatory syndromes treated afterwards, anakinra was given at doses ranging between 2 and 10 mg/kg daily (the highest doses were in low-weight CAPS patients who had failed to respond to lower doses). In the patients who took part to ANAJIS trial, CRP was recorded at each visit and retained for the pharmacodynamic modeling.

Anakinra plasma determination

Anakinra plasma determinations were performed on blood taken at one or repeated occasions depending on the study group. Whole blood samples were collected into tubes containing heparin. Plasma was separated immediatly after sampling and frozen at −20°C. Concentrations of anakinra in plasma samples were determined using the antibody (Ab) ELISA purchased from R&D Systems (Minneapolis, Minnesota, USA). Briefly, samples and quality controls were diluted with buffered animal serum and added to a microtiter plate which have been pre-coated with a monoclonal antibody specific for IL-1ra. An enzyme-linked polyclonal antibody specific for IL-1ra (horseadish peroxydase) was added to the wells. Following a wash, a substrate solution was added for color development. Reaction was stopped with sulfuric acid. Optical density was determined using a microplate reader set to 450 nm. Anakinra concentration were calculated for each sample by log-log curve fitting of the plate standards dilutions. The lower limit of quantification (LOQ) for anakinra concentrations in plasma samples was 40 ng/mL.

Pharmacokinetic modelling

Pharmacokinetic data was ascribed to an open one-or two-compartment models with linear absorption. Zero-order absorption and absorption with transit compartments models, as well as the possibility of non-linear elimination, were also considered.

Pharmacodynamic modelling

The CRP levels as a function of time and drug treatment were ascribed to a indirect response model. In this model, anakinra is thought to inhibit the response production rate, k_{TR}*R0 (transit time rate constant multiplied by response at baseline). The model equation was then

where R, C_{50} and C stand for the pharmacodynamic response, 50% inhibitor concentration and drug concentration. When the drug treatment starts, the system is at equilibrium (stable disease) what is defined by the baseline parameter, R = R_{0}.

Some individual time-courses showed that there was a loss of drug effect during the treatment time. To take this into account, an empirical resistance function was as a function of time was defined

where k_{RESI} is a time rate constant of resistance appearance. Then the response model becomes

Therefore, using a mixture model, the patients were ascribed to either equation (1) or (2). The CRP levels were analysed in the 22 SJIA patients from the ANAJIS trial.

Population pharmacokinetic and pharmacodynamic analysis

Pharmacokinetic and pharmacodynamic data were analysed using the nonlinear mixed effect modelling software program Monolix version 3.2 (_{10}–transformed to take into account the wide range of observed data during the analysis. The data were analysed sequentially; the pharmacokinetic estimates were fixed for the pharmacodynamic analysis. Parameters were estimated by computing the maximum likelihood estimator of the parameters without any approximation of the model (no linearization) using the stochastic approximation expectation maximization (SAEM) algorithm combined to a MCMC (Markov Chain Monte Carlo) procedure. The number of MCMC chains was allowed to vary in order to obtain a nice and reliable convergence of the SAEM algorithm. Additive and proportional error models were used to describe the residual variability (ϵ) for the pharmacokinetic and pharmacodynamic data respectively, and the between-subject or between occasion variabilities (η or φ) were generally ascribed to an exponential model, except specific indication. The likelihood ratio test (LRT) including the log-likelihood, the Akaike information criterion (AIC) and the bayesian information criterion (BIC) was used to test different hypotheses regarding the final model, covariate effect on structural parameter(s), residual variability model (proportional versus proportional plus additive error model), structure of the variance-covariance matrix for the BSV parameters. The normalised prediction distribution errors (NPDE) metrics

Results

Population pharmacokinetic modeling

A total of 87 children (32 girls, 52 boys) with 148 anakinra concentrations were available for pharmacokinetic evaluation, four concentrations were observed below the limit of quantification (BLQ) and coded as left censored data for the analysis. The distribution of sampling times can be observed in Figure

Anakinra concentration-time courses for once daily administrations

**Anakinra concentration-time courses for once daily administrations.** Solid and dashed lines, median and 5^{th}/95^{th} percentiles of the final model predictions.

The 22 SJIA patients were 2.26 – 16.8 years old (median 7.6) weighing 10 – 83 kg (median 21). The other patients were 0.73 – 21 years old (median 8) weighing 4.3 – 60 kg (median 21). These included 20 patients with CAPS (14 with the chronic, inflammatory, neurologic, cutaneous and articular (CINCA) syndrome / neonatal onset multisystem inflammatory disease (NOMID) and 6 with Muckle-Wells syndrome), mevalonate kinase deficiency (n = 3), TNF-receptor associated periodic syndromes (n = 2), familial mediterranean fever (n = 1) and patients with genetically undetermined autoinflammatory conditions.

A one-compartment with linear absorption and elimination model adequately described the data. Other candidate models included one-compartment with non-linear elimination and two-compartment models. The parameters of the model were Ka, V/F and CL/F, respectively the absorption rate, the apparent volume of distribution and clearance, F being the unknown bioavailability. The statistical model included a between-subject variability for CL/F, η_{CL}, and a between-occasion variability for V/F, γ_{V/F}. The residual variability was described as an additive component.

The only covariates that influenced the pharmacokinetic parameters were age and bodyweight (BW). CL/F and V/F were then related to BW using an allometric function with estimated power exponents. This decreased the CL/F and V/F variabilities from 0.41 to 0.28 and from 1.34 to 0.475 and the AIC and BIC criteria by more than 15 units. This BW effect finally removed the effect of age. No other covariate effect could be identified, gender or combined use of corticoids or anti-inflammatory drugs, AINS. The final models for CL/F and V/F were then

Figure ^{th}/95^{th} percentiles of the model predictions. Table

**Parameter**

**Estimate**

**Relative standard error (%)**

Key:

CL/F (L h^{-1} 70 kg^{-1})

6.24

8

β_{CL}, TV(CL)∙(BW/70)^{βCL}

0.47

14

V/F (L)

65.2

12

β_{V}, TV(V)∙(BW/70)^{βV}

0.76

16

Ka (h^{-1})

0.38

19

η_{CL/F}

0.28

15

γ_{V/F}

0.47

17

ϵ, mg/L

0.072

10

Pharmacodynamic modeling

An indirect response model assuming that the production rate of the inflammatory process is inhibited by anakinra concentrations described well the CRP time-courses in the 22 ANAJIS trial patients (195 observations were available). The C_{50} parameter was related to the mean steady-state anakinra concentration (Css), allowing the determination of the corresponding dosage to obtain for example 90% of the maximal effect (C_{90} ~ 10 x C_{50}), i.e., dose rate (mg/h) = CL/F*C_{90}. A mixture model with 2 subpopulations, patients with high baseline and large CRP decrease and patients with low baseline and small CRP decrease followed by a slow re-increase in baseline levels (delayed “resistance” to treatment), was finally retained. The corresponding observed baselines (median [range]), 120 [44, 230] and 43 [2.5, 152] mg/L, were significantly different (p = 0.02, Kruskal-Wallis test). However it was impossible to estimate a clear cut-off value between the 2 groups, probably because of the small sample size. There was no significant effect of the combined anti-inflammatory treatment. The proportion of patients for which the development of a “resistance” to treatment was significant was 62% and the corresponding time was approximately 60 days. Table ^{th}/95^{th} percentiles of the model predictions.

**Parameter**

**Estimate**

**Relative standard error (%)**

Key: _{
TR
}, transit time rate constant; _{
50
}, anakinra concentration that induces a 50% decrease of CRP level; _{
RESI
}, time rate constant of resistance appearance;

Baseline (mg/L)

RESP

141

27

RESI

37.9

29

k_{TR} (day^{-1})

0.042

27

C_{50} (mg/L)

0.03

37

k_{RESI} (day^{-1})

0.0048

0.0018

Proportion of RESP

0.37

31

η_{BASELINE.RESI}

0.79

24

η_{KTR}

081

25

ϵ, mg/L (*)

0.39

6

Effect of anakinra on the C-reactive concentration-time courses during the ANAJIS trial

**Effect of anakinra on the C-reactive concentration-time courses during the ANAJIS trial.** The top and bottom plots stand for the 2 response groups : top, responders; bottom, patients with a delayed “resistance” to treatment. Solid and dashed lines, median and 5^{th}/95^{th} percentiles of the final model predictions.

Evaluation and validation

Figures

Dosage simulations

Figure

Model-predicted effect of anakinra on the C-reactive concentration-time courses assuming 0.02 to 0.40 mg/L mean steady-state anakinra concentrations in the 2 subgroups of patients

**Model-predicted effect of anakinra on the C-reactive concentration-time courses assuming 0.02 to 0.40 mg/L mean steady-state anakinra concentrations in the 2 subgroups of patients.** Left, high base level with large CRP decrease. right, moderate base CRP with initial decrease followed by a re-increase in CRP concentrations.

Daily dose of anakinra (mg/kg, thick curve) as a function of bodyweight in order to reach the mean anakinra steady-state concentration of 0.4 mg/L

**Daily dose of anakinra (mg/kg, thick curve) as a function of bodyweight in order to reach the mean anakinra steady-state concentration of 0.4 mg/L.** The 0.4 mg/L target corresponds to a maximal effect on the CRP biological marker of inflammation. Bold text in rectangles defines possible mean dosage recommendation for 3 bodyweight ranges. Text on top side stands for the maximal dosage in the corresponding bodyweight range.

Discussion

This is the first pharmacokinetic study of anakinra in children. In pediatric patients with SJIA and diverse autoinflammatory conditions, the pharmacokinetics of anakinra following sub-cutaneous administration was satisfactorily described by a one compartment model with 1^{st} order absorption. Bodyweight was identified as the sole significant covariate that could finally explain a significant part of the variability and this BW effect was modeled via an allometric scaling of CL/F and V/F. The allometric power values are typically 1 and 0.75 for BW effect on V and CL respectively

Assuming that anakinra inhibits the production of the inflammatory marker CRP via an “Emax” model allowed a satisfactory description of the CRP time-course during anakinra treatment in SJIA patients. As seen in Figure _{50} estimate was low (0.03 mg/L) and nearly all the observed concentrations were above this limit. According to the predicted time-course effects of different mean Css values (Figure

Conclusions

Bodyweight significantly influenced the pharmacokinetics of anakinra. An allometric BW scaling of apparent clearance and distribution volume was sufficient to describe the variation of anakinra pharmacokinetics in these 87 pediatric patients weighing from 4 to 80 kg and 9 months to 21 years old. The effect on the C-reactive protein was adequately described by a turn-over model and allowed to derive some concentration-dose guidelines for SJIA disease.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SU, JMT and PQ wrotre the paper; SU and FF performed the modelling and analysis of the data; PQ, JMT, BBM, CB designed the research; CB was responsible for the drug assay; BBM, RM, SCL, BF, CW, BN, PQ conducted the research; SU, PQ and JMT had primary responsibility for the final content. All authors read and approved the final content.

Acknowledgements

We acknowledge Agnès Mogenet, MD (CIC Hopital Necker), Solimda Sotou-Bere and Yamina Boulahdaj, clinical research assistants, and the nurses for their time and dedication to this study. We also acknowledge the young patients and their families who participated to these trials.

Pre-publication history

The pre-publication history for this paper can be accessed here: