(-)-Ampelopsin F from Dryobalanops oblongifolia against Plasmodium falciparum Dihydrofolate Reductase Thymidylate Synthase – Quadruple Mutant: A Guide to In Silico Drug Design

Axl Laurens Lukas Windah1,2

Muhammad Ikhlas Abdjan3

Nanik Siti Aminah2,4,Email

Imam Siswanto2,5

Alfinda Novi Kristanti2,4

Theint Su Wai1,2

Geraldo Figo Hariyanto2

Yoshiaki Takaya6

1Doctoral Program of Mathematics and Natural Sciences, Faculty of Science and Technology, Universitas Airlangga, Surabaya, 60115, Indonesia
2Departement of Chemistry, Faculty of Science and Technology, Universitas Airlangga, Surabaya, 60115, Indonesia 
3Department of Chemistry, Faculty of Mathematics and Natural Science Universitas Negeri Surabaya, Surabaya, 60213, Indonesia
4Biotechnology of Tropical Medicinal Plants Research Group, Universitas Airlangga, Surabaya, 60115, Indonesia 
5Bioinformatic Laboratory, UCoE Research Center for Bio-Molecule Engineering Universitas Airlangga, Surabaya, 60115, Indonesia 
6Faculty of Pharmacy, Meijo University, 150 Yagotoyama, Tempaku, Nagoya, 468-8503, Japan

Abstract

Malaria is a deadly, life-threatening disease that poses significant global health challenges in terms of mortality and morbidity. The development of antimalarial drugs is essential to address this issue. This study aims to break down the inhibition molecular mechanism of (-)-ampelopsin F (AMP), one of the dimer stilbenoids isolated from Dryobalanops oblongifolia that have been previously demonstrated in our research to exhibit potent inhibition of the Plasmodium falciparum strain 3D7 (Pf3D7). A combination of computational chemistry methods, including pharmacokinetic analysis, molecular docking, and molecular dynamics simulations, was employed. This study focuses on the quadruple mutant (N51I+C59R+S108N+I164L) Plasmodium falciparum dihydrofolate reductase-thymidylate synthase (PfDHFR-TS) protein to evaluate the stability and molecular behavior of AMP as an inhibitor. Based on this research, we present a novel finding that the -OH functional groups in the AMP structure may enhance solvent accessibility within the PfDHFR active site. Additionally, the prediction of binding free energy (∆Gbind) using the quantum mechanics/molecular mechanics-generalized born surface area (QM/MM-GBSA) method indicates that AMP-PfDHFR (-16.78 kcal.mol-1) exhibits better inhibitory potential compared to PC6 (pyrimethamine)-PfDHFR (-4.47 kcal.mol-1) and 6,6-dimethyl-1-[3-(2,4,5-trichlorophenoxy)propoxy]-1,6-dihydro-1,3,5-triazine-2,4-diamine(WRA)-PfDHFR (-12.37 kcal.mol-1). Moreover, the interaction energy analysis of each system shows that the ∆Gbind and ∆Gexp values are strongly correlated, with AMP-PfDHFR demonstrating the most stable interaction.