Computational exploration of compounds in Xylocarpus granatum as a potential inhibitor of Plasmodium berghei using docking, molecular dynamics, and DFT studies.

ABSTRACT

Malaria remains a global health concern, with the emergence of resistance to the antimalarial drug atovaquone through cytochrome b (cyt b) being well-documented. This study was prompted by the presence of this mutation in cyt b to enable new drug candidates capable of overcoming drug resistance. Our objective was to identify potential drug candidates from compounds of Xylocarpus granatum by computationally assessing their interactions with Plasmodium berghei cyt b. Using computational methods, we modeled cyt b (GenBank AF146076.1), identified the binding cavity, and analyzed the Ramachandran plot against cyt b. Additionally, we conducted drug-likeness and absorption, distribution, metabolism, excretion, and toxicity (ADMET) studies, along with density functional theory (DFT) analysis of the compounds. Molecular docking and molecular dynamics simulation (MDS) were used to evaluate the binding energy and stability of the cyt b-ligand complex. Notably, our investigation highlighted kaempferol as a promising compound due to its high binding energy of 7.67 kcal/mol among all X. granatum compounds, coupled with favorable pharmacological properties (ADMET) and antiprotozoal properties at Pa 0.345 > Pi 0.009 (PASS value). DFT analysis showed that kaempferol has an energy gap of 4.514 eV. MDS indicated that all tested ligands caused changes in bonding and affected the structural conformation of cyt b, as observed before MDS (0 ns) and after MDS (100 ns). The most notable differences were observed in the types of hydrogen bonds between 0 and 100 ns. Nevertheles, MDS results from a 100 ns simulation revealed consistent behavior for kaempferol across various parameters including root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), solvent-accessible surface area (SASA), molecular mechanics-Poisson Boltzmann surface area (MM-PBSA), and hydrogen bonds. The cyt b-kaempferol complex demonstrated favorable energy stability, as supported by the internal energy distribution values observed in principal component analysis (PCA), which closely resembled those of the atovaquone control. Additionally, trajectory stability analysis indicated structural stability, with a cumulative eigenvalue of 24.7 %. Dynamic cross-correlation matrix (DCCM) analysis revealed a positive correlation among catalytic cytochrome residues within the amino acid residues range 119-268. The results of our research indicate that the structure of kaempferol holds promise as a potential candidate against Plasmodium.