Background: To develop a population pharmacokinetic model for intravenous omeprazole in critically ill children.
Methods: One hundred eighty-six omeprazole concentration-time data from 40 critically ill children were analyzed using the nonlinear mixed-effects approach with the nonlinear mixed-effects modeling software, version 7.2 software. Patients were randomized into 2 groups and received intravenous omeprazole at a dose of 0.5 or 1 mg/kg twice daily. Blood samples were drawn at 0.5, 2, 6, 12, 24, and 48 hours after the first infusion.
Results: The pharmacokinetic profile was best described by a 2-compartment model with a first-order elimination process. Between-patient variability could only be associated with plasma clearance (CL). The typical values for plasma CL were 24.9 L·h·70 kg (10.08%), with a distributional clearance of 53.9 L·h·70 kg (11.00%) and central and peripheral compartment distribution volumes of 4.23 L/70 kg (19.62%) and 674 L/70 kg (0.89%), respectively. Allometric size models seemed to predict changes adequately in all the pharmacokinetic parameters. High values of between-patient variability of CL [75.50% (2.60%)] and residual variability [130.0% (5.26%)] were still found in the final model. Model-based simulations suggested that the most suitable dose was 1 mg/kg because this yielded similar exposure (defined by the area under the concentration-time curve) to that obtained in adults after a 20-mg dose of omeprazole intravenously.
Conclusions: An allometric size model allows changes to be predicted in all the pharmacokinetic parameters, making dose adjustment by body weight important to achieve the most effective omeprazole exposure. This is the first step toward a population pharmacokinetic study, including more data to develop a predictable model to be used during therapeutic drug monitoring.