Bioequivalence risk assessment of oral formulations containing racemic ibuprofen through a chiral physiologically based pharmacokinetic model of ibuprofen enantiomers

The characterization of the time course of ibuprofen enantiomers can be useful in the selection of the most sensitive analyte in bioequivalence studies. Physiologically based pharmacokinetic (PBPK) modelling and simulation represents the most efficient methodology to virtually assess bioequivalence outcomes. In this work, we aim to develop and verify a PBPK model for ibuprofen enantiomers administered as a racemic mixture with different immediate release dosage forms to anticipate bioequivalence outcomes based on different particle size distributions. A PBPK model incorporating stereoselectivity and non-linearity in plasma protein binding and metabolism as well as R-to-S unidirectional inversion has been developed in Simcyp®. A dataset composed of 11 Phase I clinical trials with 54 scenarios (27 per enantiomer) and 14,452 observations (7129 for R-ibuprofen and 7323 for S-ibuprofen) was used. Prediction errors for AUC0-t and Cmax for both enantiomers fell within the 0.8-1.25 range in 50/54 (93 %) and 42/54 (78 %) of scenarios, respectively. Outstanding model performance, with 10/10 (100 %) of Cmax and 9/10 (90 %) of AUC0-t within the 0.9-1.1 range, was demonstrated for oral suspensions, which strongly supported its use for bioequivalence risk assessment. The deterministic bioequivalence risk assessment has revealed R-ibuprofen as the most sensitive analyte to detect differences in particle size distribution for oral suspensions containing 400 mg of racemic ibuprofen, suggesting that achiral bioanalytical methods would increase type II error and declare non-bioequivalence for formulations that are bioequivalent for the eutomer.

Quantitative assessment of the exposure-efficacy relationship of glucocerebrosidase using Markovian elements in Gaucher patients treated with enzyme replacement therapy

AIMS

The aims of this study are (i) to develop a population pharmacokinetic model of enzyme activity in Gaucher-type 1 (GD1) patients after intravenous administration of enzyme replacement therapy (ERT) and, (ii) to establish an exposure-efficacy relationship for bone marrow infiltration to propose dose adjustments according to patient covariate values.

MATERIALS AND METHODS

A prospective follow-up, semi-experimental multi-centre study was conducted in four hospitals to evaluate the pharmacokinetics, efficacy and safety of ERT in GD1 patients. 25 individuals with 266 glucocerebrosidase (GCase) observations in plasma and leukocytes and 14 individuals with 68 Spanish Magnetic Resonance Imaging (S-MRI) observations were enrolled.

RESULTS

A two concatenated compartments with zero-order endogenous production and first-order distribution (CL₁ =3.85 x10^-1 L/d) and elimination (CL₂ = 1.25 L/d) allowed to describe GCase observations in plasma and leucocytes, respectively. An exponential time-dependency (k_T =6.14 x10^-1 d^-1 ) effect on CL₁ was incorporated. The final exposure-efficacy model was a longitudinal logistic regression model with a first-order Markov element. An Emax function (EC₅₀ =15.73 U/L and Emax=2.33) linked steady-state concentrations of GCase in leucocytes to the probability of transition across the different S-MRI stages.

CONCLUSION

A population pharmacokinetic model successfully characterized the leukocyte activity-time profiles of GCase following intravenous administration of ERT in GD1 patients together with an exposure-efficacy relationship in bone marrow using markovian elements. The information obtained from this study could be of high clinical relevance in individualization of ERT in GD1 patients, as this could lead to anticipate decision-making regarding clinical response in bone and optimal dosing strategy. NONSTANDARD ABBREVIATIONS: -2LL: -2xlog(likelihood); ERT: enzyme replacement therapy; GCase: glucocerebrosidase activity; GD1: Gaucher disease type 1; GOF: goodness-of-fit plots; IIV: inter-individual variability; NLME: non-linear mixed effects modelling; OFV: objective function value; pc-VPC: prediction-corrected visual predictive check; PK: pharmacokinetic; RSE: relative standard error; RUV: residual unexplained variability, S-MRI: Spanish Magnetic Resonance Imaging, TDM: therapeutic drug monitoring.

Semi-Mechanistic Pharmacokinetic Model to Guide the Dose Selection of Nimotuzumab in Patients with Autosomal Dominant Polycystic Kidney Disease

Autosomal dominant polycystic kidney disease (ADPKD) is a genetic disease characterized by an overexpression of epidermal growth factor receptor (EGFR). Nimotuzumab is a recombinant humanized monoclonal antibody against human EGFR. The aim of this study was to develop a population pharmacokinetic model for nimotuzumab and to identify demographic and clinical predictive factors of the pharmacokinetic variability. The population pharmacokinetics (PopPK) of nimotuzumab was characterized using a nonlinear mixed-effect modeling approach with NONMEM^®. A total of 422 log-transformed concentration-versus-time datapoints from 20 patients enrolled in a single-center phase I clinical trial were used. Quasi steady state approximation of the full TMDD (target-mediated drug disposition) model with constant target concentration best described the concentration-time profiles. A turnover mediator was included which stimulates the non-specific clearance of mAb in the central compartment in order to explain the reduced levels at higher doses. Covariates had no influence on the PK (pharmacokinetics) parameters. The model was able to detect that the maximum effective dose in ADPKD subjects is 100 mg. The developed PopPK model may be used to guide the dose selection for nimotuzumab during routine clinical practice in patients with polycystic kidney disease. The model will further support the ongoing investigations of the PK/PD relationships of nimotuzumab to improve its therapeutic use in other disease areas.

Evaluation of ABC gene polymorphisms on the pharmacokinetics and pharmacodynamics of capecitabine in colorectal patients: Implications for dosing recommendations

AIMS

The aims are to develop a population pharmacokinetic model of capecitabine (CAP) and its main metabolites after the oral administration of CAP in colorectal cancer patients with different polymorphisms of the ATP-binding cassette (ABC) gene and a population pharmacokinetic/pharmacodynamic model capable of accounting for the neutropenic effects, and to optimize the dosing strategy based on the polymorphisms of the ABC gene and/or the administration regimen as a single agent or in combination.

METHODS

Forty-eight patients diagnosed with colorectal cancer were included, with 432 plasma levels of CAP, 5'-desoxi-5-fluorouridine (5'-DFUR) and 5-fluorouracil (5-FU), and 370 neutrophil observations. Capecitabine doses ranged from 1250 to 2500 mg/m² /24 h. Plasma measurements of CAP, 5'-DFUR and 5-FU were obtained at 1, 2 and 3 hours post administration. Neutrophil levels were measured between day 15 and day 24 post administration.

RESULTS

The pharmacokinetic model incorporates oxaliplatin as a covariate on absorption lag time, rs6720173 (ABCG5 gene) on clearance of 5'-DFUR (182% increase for mutated rs6720173) and rs2271862 (ABCA2 gene) on clearance of 5-FU (184% increase for mutated rs2271862). System- (Circ₀ = 3.54 × 10⁹ cells/mL, MTT = 204 hours and γ = 6.0 × 10^-2 ) and drug-related (slope [SLP] = 3.1 × 10^-2 mL/mg). Co-administration of oxaliplatin resulted in a 2.84-fold increase in SLP. The predicted exposure thresholds to G3/4 neutropenia in combination and monotherapy were 26 and 70 mg·h/L, respectively.

CONCLUSIONS

The population pharmacokinetic/pharmacodynamic model characterized the time course of capecitabine and its metabolites in plasma. Dose recommendations of capecitabine in patients with mutated and wild allele for single nucleotide polymorphisms rs2271862 of ≤3000 and ≤2400 mg/m² /24 h in monotherapy and ≤1750 and ≤600 mg/m² /24 h in combination with oxaliplatin, respectively, have been proposed.