METHODS: Pkdhps were amplified and sequenced from 28 P. knowlesi samples collected in 2008 and 2020 from nine provinces across Thailand. Combining pkdhfr sequencing data from previous work with pkdhps data to analyze polymorphisms of pkdhfr and pkdhps haplotype. Protein modeling and molecular docking were constructed using two inhibitors, sulfadoxine and sulfamethoxazole, and further details were obtained through analyses of protein-ligand interactions by using the Genetic Optimisation for Ligand Docking program. A phylogenetic tree cluster analysis was reconstructed to compare the P. knowlesi Malaysia isolates.

RESULTS: Five nonsynonymous mutations in the pkdhps were detected outside the equivalence of the binding pocket sites to sulfadoxine and sulfamethoxazole, which are at N391S, E421G, I425R, A449S, and N517S. Based on the modeling and molecular docking analyses, the N391S and N517S mutations located close to the enzyme-binding pocket demonstrated a different docking score and protein-ligand interaction in loop 2 of the enzyme. These findings indicated that it was less likely to induce drug resistance. Of the four haplotypes of pkdhfr-pkdhps, the most common one is the R34L pkdhfr mutation and the pkdhps quadruple mutation (GRSS) at E421G, I425R, A449S, and N517S, which were observed in P. knowlesi in southern Thailand (53.57%). Based on the results of neighbor-joining analysis for pkdhfr and pkdhps, the samples isolated from eastern Thailand displayed a close relationship with Cambodia isolates, while southern Thailand isolates showed a long branch separated from the Malaysian isolates.