In this work we have collected a set of 30 trypanosomicidal naphthoquinones and developed pharmacophoric and 3D-QSAR models as tools for the design of new potential anti-Chagasic compounds. Firstly, qualitative information was obtained from SAR and pharmacophoric models identifying some fragments around the 2-aryloxynaphthoquinone scaffold important for the antiparasitic activity. Then, 3D-QSAR CoMFA and CoMSIA models were developed. The models showed adequate statistical parameters where the steric, electrostatic, and hydrophobic features explain the trypanosomicidal effect. Therefore, to validate our models, we carried out the design, synthesis, and biological evaluation on T. cruzi epimastigotes of five new compounds (33a-e). According to CoMFA model, three out of five compounds showed pIC50 values within one logarithmic unit of deviation. The two compounds that did not fit the predictions were those with high lipophilicity, which agreed with the SAR and pharmacophore models. Docking and molecular dynamic studies were performed on T. cruzi trypanothione reductase, in a proposed binding site for this type of naphthoquinone. Interestingly, 33a-e showed the same interaction pattern as a naphthoquinone inhibitor (2). Finally, predicted drug-likeness properties indicated that 33a-e have optimal oral bioavailability. Thus, this study provides new in silico models for obtaining novel trypanosomicidal compounds.
Keywords: 3D-QSAR; Trypanosoma cruzi; aryloxy-quinones; pharmacophore analysis; structure-activity relationship; trypanosomicidal activity.