A total of 36 cycles of chemotherapy from 17 patients were evaluated in the present study. Patient characteristics and demographics are summarized in
Table 1. For developing the pharmacokinetic model, we used a total of 131 melphalan plasma concentrations. A two-compartment model with first-order elimination was used and adequately fitted the data (see
Fig. 3A for a representative patient receiving the drug to one eye only). In children receiving tandem therapy, we noted that the interval of time marking the end of infusion of the first eye to the beginning of the infusion into the second eye had an impact on the drug disposition in plasma as shown in
Figure 3B. In delayed tandem therapy with that long interval, a bimodal profile would be observed. This was in contrast to a classical constant infusion one-peak profile found in short intereye times where elimination from the body is limited. As melphalan was administered to the same patient in multiple occasions, interoccasion variability was also incorporated in the base model. The population pharmacokinetic parameters obtained for the final model are presented in
Table 2. A large interindividual variability was observed accounting for 56% and 52% for CL and
V c, respectively. The covariate analysis showed that interindividual variability in CL and
V c could be explained in part by age, weight, and body surface area (
P < 0.05). For the final model, only weight was considered because it explained 81% and 75% of the IIV in CL and
V c, respectively. A relationship between systemic exposure to melphalan (AUC) and dose corrected by weight was observed as depicted in
Figure 4. Interestingly, as shown in
Figure 5A, we observed a significantly higher AUC in children receiving more than 0.48 mg/kg of melphalan (
P < 0.05). In addition, melphalan plasma concentrations after 2 and 4 hours of starting drug infusion significantly differed between children receiving SSOAI to one eye compared with those receiving tandem therapy, as represented in
Figures 5B,
5C (
P < 0.05).