Computer-aided discovery of novel SmDHODH inhibitors for schistosomiasis therapy: Ligand-based drug design, molecular docking, molecular dynamic simulations, drug-likeness, and ADMET studies

Schistosomiasis, also known as bilharzia or snail fever, is a tropical parasitic disease resulting from flatworms of the Schistosoma genus. This often overlooked disease has significant impacts in affected regions, causing enduring morbidity, hindering child development, reducing productivity, and creating economic burdens. Praziquantel (PZQ) is currently the only treatment option for schistosomiasis. Given the potential rise of drug resistance and the limited treatment choices available, there is a need to develop more effective inhibitors for this neglected tropical disease (NTD). In view of this, quantitative structure-activity relationship studies (QSAR), molecular docking, molecular dynamics simulations, drug-likeness, and ADMET predictions were applied to 31 inhibitors of Schistosoma mansoni Dihydroorotate dehydrogenase (SmDHODH). The designed QSAR model demonstrated robust statistical parameters including an R2 of 0.911, R2adj of 0.890, Q2cv of 0.686, R2pred of 0.807, and cR2p of 0.825, confirming its robustness. Compound 26, identified as the most active derivative, emerged as a lead candidate for new potential inhibitors through ligand-based drug design. Subsequently, 12 novel compounds (26A-26L) were designed with enhanced inhibition activity and binding affinity. Molecular docking studies revealed strong and stable interactions, including hydrogen bonding and hydrophobic interactions, between the designed compounds and the target receptor. Molecular dynamics simulations over 100 nanoseconds and MM-PBSA free binding energy (ΔGbind) calculations validated the stability of the two best-designed molecules (26A and 26L). Furthermore, drug-likeness and ADMET prediction analyses affirmed the potential of these designed compounds, suggesting their promise as innovative agents for treating schistosomiasis.