Population Pharmacokinetic and Pharmacodynamic Modeling of Artemisinin Resistance in Southeast Asia

Orally administered artemisinin-based combination therapy is the first-line treatment against uncomplicated P. falciparum malaria worldwide. However, the increasing prevalence of artemisinin resistance is threatening efforts to treat and eliminate malaria in Southeast Asia. This study aimed to characterize the exposure-response relationship of artesunate in patients with artemisinin sensitive and resistant malaria infections. Patients were recruited in Pailin, Cambodia (n = 39), and Wang Pha, Thailand (n = 40), and received either 2 mg/kg/day of artesunate mono-therapy for 7 consecutive days or 4 mg/kg/day of artesunate monotherapy for 3 consecutive days followed by mefloquine 15 and 10 mg/kg for 2 consecutive days. Plasma concentrations of artesunate and its active metabolite, dihydroartemisinin, and microscopy-based parasite densities were measured and evaluated using nonlinear mixed-effects modeling. All treatments were well tolerated with minor and transient adverse reactions. Patients in Cambodia had substantially slower parasite clearance compared to patients in Thailand. The pharmacokinetic properties of artesunate and dihydroartemisinin were well described by transit-compartment absorption followed by one-compartment disposition models. Parasite density was a significant covariate, and higher parasite densities were associated with increased absorption. Dihydroartemisinin-dependent parasite killing was described by a delayed sigmoidal Emax model, and a mixture function was implemented to differentiate between sensitive and resistant infections. This predicted that 84% and 16% of infections in Cambodia and Thailand, respectively, were artemisinin resistant. The final model was used to develop a simple diagnostic nomogram to identify patients with artemisinin-resistant infections. The nomogram showed > 80% specificity and sensitivity, and outperformed the current practice of day 3 positivity testing.