Abstract
Plasmodium falciparum thymidylate kinase (PfTMK) showed structural and catalytic distinctions from the host enzyme rendering it a hopeful antiprotozoal drug target. Despite the comprehensive enzymologic, structural, inhibitory and chemical synthesis approaches targeting this enzyme, the elucidation of the exact mechanism underlying the recognition of the atypical purine substrates remains to be determined. In this study, molecular dynamics (MD) simulation of a broad range of substrates and inhibitors as well as the inhibitory properties of deoxyguanosine (dG) derivatives were used to assess the PfTMK substructure molecular rearrangements. The estimated changes during the favourable binding of high affinity substrate (TMP) include lower interaction with P-loop, free residue fluctuations of the lid domain and the average RMSD value. The RMSD of TMP complex was higher and more rapidly stabilized than the dGMP complex. The lid domain flexibility is severely affected by dGMP and beta-thymidine derivatives, while being partially fluctuating with other thymidine derivatives. The TMK-purine (dGMP) complex was slowly and gradually stabilized with lower over all structure flexibility and residue fluctuations especially at the lid domain, which closes the active site during its catalytic state. Thymidine derivatives allow structure flexibility of the lid domain being highly fluctuating in alpha- and beta-thymidine derivatives and TMP. dG derivatives remains less efficient than thymidine derivatives in inhibiting TMK. The variations in the structural dynamics of the P-loop and lid domain in response to TMP or dGMP might favour thymidine-based compounds. The provided MD simulation strategy can be used for predicating structural changes in PfTMK during lead optimization.