3.1. Patient Characteristics:
Twenty-nine patients enrolled in the study, and their baseline characteristics are in Table 1. Eight blood samples were collected from each patient for PK assessment. As some patients expired before all the blood samples could be collected, 227 pegasparaginase blood levels were available for pharmacokinetic analysis.
Table 1: Baseline demographics of patients
Baseline characteristics
|
All patients (N=29)
|
Gender (M:F)
|
19:10
|
Age (years)
|
9.7 ± 2.7
|
Weight (kg)
|
30.1 ± 12.4
|
Height (cm)
|
135.3 ±17.6
|
BMI (kg/m2)
|
16.28 ± 3.98
|
BSA (m2)
|
1.06 ± 0.26
|
Type of ALL (n)
B-ALL
T-ALL
|
25
4
|
Fibrinogen (mg/dL)
|
381.43 ± 170.38
|
Antithrombin III (%)
|
118.28 ± 12.77
|
Serum ALT (U/L)
|
47.28 ± 49.42
|
Serum AST (U/L)
|
40.38 ± 33.43
|
Serum ALP (U/L)
|
171.10 ± 70.78
|
Haemoglobin (g/dL)
|
10.02 ± 2.3
|
Platelet count (x109/L)
|
54688.89 ± 91711.44
|
S. Creatinine (mg/dL)
|
0.56 ± 0.10
|
S. Bilirubin (mg/dL)
|
0.54 ± 0.25
|
ALL: Acute Lymphoblastic Leukemia; ALP: serum alkaline phosphatase;
ALT: alanine transaminase, AST: serum aspartate transaminase;
All continuous data are shown as Mean±SD
3.2. Population Pharmacokinetic Model of PEG-Asparaginase:
A one-compartment model with first-order absorption and elimination adequately described pegasparaginase pharmacokinetics. The model was parameterized as absorption rate constant (ka), apparent Clearance (CL), and apparent Volume of distribution (V). The additive error model was superior to other error models. Between-subject variability was quantified using an exponential model.
Figure 1 shows the diagnostic plots for the goodness of fit for the final model. Observed (DV) vs. Individual predicted concentrations (IPRED) showed good agreement, which shows the appropriateness of the structural model (Fig.1-A). The lack of any trend in the conditional weighted residuals vs. the population-predicted concentration plot suggests the adequacy of the residual error model (Fig.1-B). Inter-individual variability was observed in the population-predicted plot (Fig.1-C).
3.3. Relationship between predicted pharmacokinetic parameters and covariates:
The inclusion of age as a covariate on the Volume of distribution (Age on V) in the base model resulted in a significant decrease in the OFV (ΔOFV = -34) (Table 2). Shrinkage was evident in the visual inspection of the observed concentration vs. predicted concentration plot (Figure 1d). Moreover, inter-individual variability (% CV) for the Volume of distribution decreased from 30% in the base model to less than 1%.
Table 2: Stepwise and statistical values used for discrimination
Model No.
|
Model Description
|
OFV (-2 log-likelihood)
|
ΔOFV
|
p
|
Forward-inclusion
|
1
|
Base model
|
2569
|
|
|
2
|
Age (months) on Volume of distribution
|
2535
|
-34
|
<0.05
|
3
|
Age (months) on Clearance
|
2557
|
-12
|
<0.05
|
4
|
BSA on Clearance
|
2558
|
-11
|
<0.05
|
5
|
Age (months) and BSA on Clearance
|
2556
|
-13
|
<0.05
|
6
|
Age (months) on Volume of distribution and Clearance
|
2528
|
-41
|
<0.05
|
Backward-elimination
|
7
|
Removed BSA from model 5
|
2557
|
+1
|
>0.05
|
8
|
Removed Age (Months) from model 6
|
2536
|
+8
|
<0.05
|
ΔOFV: change in Objective Function Value
Applying age as a covariate on Clearance (Age on CL) also significantly decreased OFV (ΔOFV = -12). A similar reduction in OFV was also observed with BSA as a covariate on Clearance (BSA on CL). However, the redundancy of both covariates became apparent when both were included simultaneously, as OFV did not decrease further. As the current cohort of patients was already dosed according to the existing practice of BSA-based dosing, we decided to keep age as the sole covariate for V and CL in the final model. Other tested covariates – baseline lab values (S. bilirubin, S. creatinine, S. albumin, ALP, AST, ALT) and baseline anti-asparaginase antibody status did not significantly decrease in OFV. Equations for the final covariate model are shown in Eqs. 1 and 2:
CL = 0.0117972* Age (months)^ 1.59 ------------ (1)
V = 0.00329742 * Age (months)^ 1.28 ------------ (2)
The parameter estimates from the final model and the associated 95% confidence intervals are provided in Table 3.
Table 3: The parameter estimates from the final model
Parameter
|
Estimate
|
SE
|
CV%
|
2.5% CI
|
97.5% CI
|
tvKa (1/hr)
|
0.027
|
0.003
|
11.494
|
0.021
|
0.033
|
tvV (L)
|
0.003
|
0.001
|
36.817
|
0.001
|
0.006
|
tvCl (L/hr)
|
0.012
|
0.001
|
8.211
|
0.010
|
0.014
|
dVd AGE
|
1.281
|
0.075
|
5.837
|
1.134
|
1.428
|
dCL dAGE
|
1.598
|
0.323
|
20.212
|
0.962
|
2.235
|
SD
|
48.15
|
3.70
|
7.69
|
40.85
|
55.44
|
CI: Confidence interval; CV: coefficient of variability; dCl dAGE = fixed parameter coefficient of age on Clearance; dVd AGE = fixed parameter coefficient of age on Volume of distribution; SD= standard deviation; SE: standard error; tvCL = typical value of Clearance; tvka = typical value of absorption rate constant; tvV = typical value of Volume of distribution
3.4. Model validation:
The model's performance was validated using 1000 replicates generated from the original dataset to evaluate the stability of the final model. The mean values of pharmacokinetic parameters were within the 95% CI of the bootstrap values, indicating that all pharmacokinetic parameters can be estimated with acceptable precision (Table 4).
Table 4: Estimates of population pharmacokinetic parameters obtained by fitting the final model to bootstrap samples
Parameter
|
Mean
|
SE
|
CV%
|
Median
|
2.50% CI
|
97.50% CI
|
tvKa (1/hr)
|
0.034
|
0.004
|
12.767
|
0.036
|
0.025
|
0.042
|
tvV (L)
|
0.525
|
0.797
|
151.830
|
0.003
|
0.001
|
1.741
|
tvCl (L/hr)
|
0.012
|
0.002
|
15.905
|
0.011
|
0.009
|
0.014
|
dCl dAGE
|
0.988
|
0.793
|
80.176
|
1.171
|
0.000
|
2.203
|
dV dAGE
|
0.971
|
0.652
|
67.106
|
1.298
|
0.000
|
1.701
|
SD
|
48.50
|
4.04
|
8.32
|
47.71
|
38.64
|
56.54
|
CI: Confidence interval; CV: coefficient of variability; dCl dAGE = fixed parameter coefficient of age on Clearance; dV dAGE = fixed parameter coefficient of age on Volume of distribution; SD = standard deviation; SE: standard error; tvCL = typical value of Clearance; tvka = typical value of absorption rate constant; tvV = typical value of Volume of distribution
The visual predictive checks show that the observed concentrations correspond with the 90% prediction intervals of the 5th, 50th, and 95th percentiles, calculated from 1000 simulated datasets, as shown in Figure 2.
3.5 Effect of age on Pharmacokinetic parameters:
The impact of a patient's age on pegasparaginase pharmacokinetic parameters was assessed. The volume of distribution increased as the patients' age increased. However, the patient's age had no significant influence on drug clearance. Cmax and AUC were compared in different age groups (Table 5). Cmax in patients over 120 months was 28 % lower than in patients under 120 months (p<0.05). A similar reduction was observed in AUClast (22%) between two age groups, but it was not statistically significant. By stratifying patients into four age groups (75-100 months, 101-125 months, 126-150 months, and 151 – 200 months), we found that Cmax and AUC decrease with age (Figure 3).
Table 5: Observed and Predicted PK parameters
|
Observed$
|
Predicted NMD*
|
Predicted MD#
|
|
Cmax
(IU/L)
|
AUClast
(IU.hr)/L
|
Cmax
(IU/L)
|
AUClast
(IU.hr)/L
|
AUCinf
(IU.hr)/L
|
Cmax
(IU/L)
|
AUCinf
(IU.hr)/L
|
Patient age ≤10 yrs (n=15)
|
503.0
± 23.3
|
83762 ± 4088.6
|
438 ± 19.5
|
77587 ± 6268.9
|
92932 ± 10480.6
|
475.8 ± 20.7
|
100859.9 ± 11170.9
|
Patient age >10 yrs (n=14)
|
358.9
± 41.8
|
64620 ± 8500.7
|
362 ± 25.8
|
64193 ± 6729.5
|
78469 ± 10880.1
|
475.4 ± 33.5
|
103280.3 ± 14625.0
|
All
(n=29)
|
433.4
± 26.8
|
74521 ± 4877.3
|
401 ± 17.3
|
71121 ± 4681.2
|
85950 ± 7537
|
475.6 ± 19.1
|
102028.4 ± 8959.1
|
All continuous data is shown as Mean ± SE;
$: obtained by Noncompartmental analysis (NCA)
*NMD: Non-Modified Dose- 1000 IU/ m2 BSA
#MD: Modified Dose – according to formula (4*age in months+715)/m2 BSA
3.6. Simulation:
The final PopPK model was used to predict the concentration-time profiles of pegasparaginase in the same set of patients in whom the study was conducted (Table 5 and Figure 3). Concentrations predicted from non-modified doses (NMD) were comparable with those observed in the study, indicating that the developed model could accurately predict the levels. Furthermore, predicted Cmax and AUCinf of the modified dose (MD) were uniform across the different age groups.
3.7. Effect of pegasparaginase on lab parameters:
Changes in the safety lab parameters (like blood glucose, antithrombin III, and S. fibrinogen) post-treatment with pegasparaginase showed no differential effects when patients were stratified by age.
After the induction phase, the patient's Minimal Residual Disease (MRD) status was assessed; five out of 29 patients remained MRD positive, with more patients in the age group of more than 120 months (p-ns). The average Cmax in MRD negative and positive patients was 437.25 IU/L and 373.73 IU/L, respectively (p-ns). The average Cmax and AUClast in the patients aged ≤ 120 months and who were MRD negative was 470.9 ± 27.6 IU/L and 102871.3 ± 7681.9 IU.hr/L, respectively.