A population pharmacokinetics analysis will be performed using a non-linear model of mixed effects with the software Nonlinear Mixed-effects Model (NONMEM), R and Piranha. Parasite burden will be determined by 7SLRNA Quantitative Polymerase Chain Reaction (qPCR) of Leishmania from swab samples of lesions and extralesional tissues before and at the end of treatment. The relationship between pharmacokinetics and parasite persistence/burden will be determined by correlation analysis and pharmacodynamic modeling.