Next Article in Journal
Glycoproteomic Analysis of Urinary Extracellular Vesicles for Biomarkers of Hepatocellular Carcinoma
Previous Article in Journal
Impact of Ag/ZnO Reinforcements on the Anticancer and Biological Performances of CA@Ag/ZnO Nanocomposite Materials
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 28
Issue 3
10.3390/molecules28031292
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unveiling the Efficacy of Sesquiterpenes from Marine Sponge Dactylospongia elegans in Inhibiting Dihydrofolate Reductase Using Docking and Molecular Dynamic Studies

by
Abdelsattar M. Omar
1,2,3,*,
Khadijah A. Mohammad
1,
Ikhlas A. Sindi
4,
Gamal A. Mohamed
5 and
Sabrin R. M. Ibrahim
6,7,*
1
Department of Pharmaceutical Chemistry, Faculty of Pharmacy, King Abdulaziz University, Jeddah 21589, Saudi Arabia
2
Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Al-Azhar University, Cairo 11884, Egypt
3
Center for Artificial Intelligence in Precision Medicines, King Abdulaziz University, Jeddah 21589, Saudi Arabia
4
Department of Biology, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
5
Department of Natural Products and Alternative Medicine, Faculty of Pharmacy, King Abdulaziz University, Jeddah 21589, Saudi Arabia
6
Preparatory Year Program, Department of Chemistry, Batterjee Medical College, Jeddah 21442, Saudi Arabia
7
Department of Pharmacognosy, Faculty of Pharmacy, Assiut University, Assiut 71526, Egypt
*
Authors to whom correspondence should be addressed.
Molecules 2023, 28(3), 1292; https://doi.org/10.3390/molecules28031292
Submission received: 30 December 2022 / Revised: 25 January 2023 / Accepted: 27 January 2023 / Published: 29 January 2023
Browse Figures
Versions Notes

Abstract

:
Dihydrofolate reductase (DHFR) is a crucial enzyme that maintains the levels of 5,6,7,8-tetrahydrofolate (THF) required for the biological synthesis of the building blocks of DNA, RNA, and proteins. Over-activation of DHFR results in the progression of multiple pathological conditions such as cancer, bacterial infection, and inflammation. Therefore, DHFR inhibition plays a major role in treating these illnesses. Sesquiterpenes of various types are prime metabolites derived from the marine sponge Dactylospongia elegans and have demonstrated antitumor, anti-inflammation, and antibacterial capacities. Here, we investigated the in silico potential inhibitory effects of 87 D. elegans metabolites on DHFR and predicted their ADMET properties. Compounds were prepared computationally for molecular docking into the selected crystal structure of DHFR (PDB: 1KMV). The docking scores of metabolites 34, 28, and 44 were the highest among this series (gscore values of −12.431, −11.502, and −10.62 kcal/mol, respectively), even above the co-crystallized inhibitor SRI-9662 score (−10.432 kcal/mol). The binding affinity and protein stability of these top three scored compounds were further estimated using molecular dynamic simulation. Compounds 34, 28, and 44 revealed high binding affinity to the enzyme and could be possible leads for DHFR inhibitors; however, further in vitro and in vivo investigations are required to validate their potential.

1. Introduction

DHFR (Dihydrofolate reductase) is a substantial enzyme that is accountable for the conversion of DHF (7,8-dihydrofolate) to THF (5,6,7,8-tetrahydrofolate), as well as folate to DHF in the presence of NADPH (nicotinamide adenine dinucleotide phosphate) [1]. This enzyme exists in various organisms, such as humans, plants, animals, and bacteria [2]. It is needed for maintaining THF, which is essential for carbon atom donation in the synthesis of pyrimidines, purines, and amino acids (methionine, glycine, N-formyl-methionyl tRNA, and serine). Thus, it is crucial for proper cellular proliferation and growth. Its prohibition leads to pyrimidine and purine starvation, which consequently stops DNA and RNA synthesis, resulting in cell death [3]. DHFR is a remarkable target for folate antagonists that are known as valuable therapeutic agents for inflammatory, neoplastic, parasitic, and infectious illnesses [4,5,6,7]. For example, it is targeted by anticancer agents such as methotrexate, which is renowned for treating leukemia, rheumatoid arthritis, osteosarcoma, lymphoma, and breast and lung cancers [5,8], as well as by antimalarial (e.g., proguanil and pyrimethamine) [7] and antibacterial (e.g., trimethoprim) drugs [4,9].
Owing to the growing drug resistance to some of the available DHFR inhibitors, research for the discovery of new and selective DHFR inhibitors has been intensely increased, using molecular modeling, synthetic techniques, in vivo and in vitro biological investigation, mechanistic studies, and structure-activity relationships [10,11,12].
Computer-aided drug design (CADD), including molecular design, molecular modeling, computational chemistry, and rational drug design, is a contemporary computational tool that is utilized in drug discovery methods for identifying and developing potential leads in silico, leading to a decrease in the number of metabolites to be experimentally assessed [13]. These methods have made crucial contributions to the development of drugs that are in clinical trials or have clinical uses. In this tool, diverse programs and software are utilized to generate and filter a set of compounds based on specific criteria, predict their physicochemical characters, predict proper targets, and evaluate their binding affinity to the predicted targets. SBDD (structure-based drug design) is a category of CADD that uses the 3D structure of the target to perform molecular dynamic simulation (MD) and docking investigations. Docking assessed the strength of compound-target binding; however, MD estimated the ligand-protein complex behavior and stability in aqueous circumstances to simulate the physiological condition [14].
The marine environment is a promising source of valuable natural metabolites, many of which are launching in the market or being evaluated in clinical trials as new drugs, particularly for cancer treatment and antimicrobials [15,16,17]. Recently, molecular docking has become a substantial tool in marine drug research worldwide for screening various marine metabolites to predict their possible bioactivities and mechanisms of action because of its relatively uncomplicated procedures [18,19].
Among marine sponges, Dactylospongia elegans (D. elegans) has been demonstrated to be a rich source of diverse metabolites with substantial bioactivities, including sesquiterpenes such as hydroquinone, quinone, and tetronic acid derivatives. The sesquiterpene hydroquinone/quinone family of terpenoids possessing a drimane or a rearranged drimane framework is a distinguished type of metabolite that is a result of mixed biosynthesis. They feature a bicyclic sesquiterpene moiety linked to a quinol or quinone moiety [20]. In our continued research interest to discover untapped bioactivities and underpinning mechanisms of the marine reported metabolites using SBDD, 87 sesquiterpene derivatives reported from D. elegans were screened for their DHFR inhibitory potential using molecular docking. Additionally, the highly ranked metabolites were further assessed by MD simulation.

2. Results and Discussion

2.1. Molecular Docking Evaluation

Diverse natural metabolites of various structural classes, including terpenoids, have been assessed for their DHFR inhibition (DHFRI) capacity using in silico or in vitro approaches [21,22,23,24]. For example, cynaropicrin is a sesquiterpene lactone reported from artichoke that has notable DHFRI potential (IC50 7.1 µM) [24]. Herrera-Acevedo et al. reported that the kaurane-type diterpene, 3α-cinnamoyloxy-ent-kaur-16-en-19-oic acid, possessed a high binding affinity to LbDHFR-TS (Leishmania panamensis dihydrofolate-reductase-thymidylate-synthase), which is an important target for L. panamensis [25]. Further, in a docking study by Kumar et al., labdane diterpenoid-derived aulocarpin purified from Afromomum aulocarpus seed displayed Plasmodium wild-type DHFR antagonistic activity (Figure 1) [21].
Furthermore, some studies have reported the DHFRI potential of terpenoids containing the quinone moiety; for example, stachybotrydial is a drimane sesquiterpene containing a dioxgenated isobenzofuranone moiety, as reported by Kwon et al. from the soil derived Stachybotrys sp. FN298. It demonstrated S. aureus`s DHFR (IC50, 41 µM) inhibition potential and prohibited MRSA growth (MIC, 32 µg/mL), suggesting its possible efficacy as an antibacterial agent versus MRSA [26]. Coscinoquinol purified from the Australian Coscinoderma sp. sponge had cytotoxic potential versus P-388, A-549, HT-29, and CV-1 (IC50 values of 0.25, 0.5, 0.25, and 0.5 µg/mL, respectively) cell lines and inhibited DHFR, TOPO II isomerase, and glutathione reductase (IC50, 2.5, 0.5, and 15.0 µg/mL, respectively) [27]. Panicein F1, a sesquiterpene hydroquinone from Reniera mucosa sponge, was found to exhibit DHFRI (IC50 3 µg/mL) (Figure 1) [28].
The inhibitory potential of previously reported quinone-containing terpenoids towards the DHFR suggested that other unexamined and structurally related metabolites, such as sesquiterpenes derived from D. elegans, could also have inhibitory effects on the same protein. Therefore, hDHFR was selected as the study target. About 87 metabolites (Figure 2, Figure 3, Figure 4 and Figure 5) were virtually investigated in the docking studies performed entirely in the Schrödinger program (Schrödinger Release 2022-3, LLC, New York, NY, USA, 2021).
It is noteworthy that sesquiterpenes are the prime constituents separated from D. elegans that revealed cytotoxic, antitumor, anti-inflammatory, and antibacterial capacities (Table 1) [29].
For validation purposes of the docking method, the co-crystallized inhibitor SRI-9662 was redocked inside the active site of the prepared hDHFR, and the original and redocked inhibitors were superimposed. By comparing the two results, the redocked inhibitor produced a nearly identical pose to that of the original crystal structure. The calculated root-mean-square deviation (RMSD) value of the superimposition was 0.2105 Å within the acceptable range (Figure 6).
The docking results of the study compounds were listed in Table 2, where they were ranked according to their gscores from highest to lowest in free energy of binding; the more negative scores imply better binding. The top-ranked compounds were 34, 28, and 44, with gscores of −12.431, −11.502, and −10.62 kcal/mol, respectively. The gscores of these compounds exceeded the value of the reference inhibitor SRI-9662, which had a binding energy of −10.432 kcal/mol (Table 2). The detailed docking parameters and scores are presented in Table S1.
The 2D and 3D structures of the co-crystallized inhibitor SRI-9662 revealed that the 4-amino group in the 5-deazapteridine ring interacted with Val115 and Ile7 backbones through hydrogen bonds (H-bond). Moreover, N-1 and N-8 in the ring are also bound through H-bonds to Glu30 with and without a water bridge, respectively. The pyrimidine ring involved in π-π stacking interaction with Phe34, as well as the rest of the molecule, are hydrophobically bound with the nearby hydrophobic residues in the hDHFR active site (Figure 7).
For compound 34, there was an array of H-bonding interactions between the carbonyl, hydroxyl, and carboxylic groups in the molecule and amino acid residues such as Asn64, Pro66, Gln35, Lys68, and Arg70. Water molecules were involved as well in such interactions (Figure 8). Additionally, Arg70 formed an ionic bond with the carboxylate ion in 34, and the decahydronaphthalene ring was involved in hydrophobic interactions.
The structure of compound 28 was very similar to that of compound 34; hence, it involved similar binding interactions with the active site residues as well as water molecules (Figure 9).
The carboxylic acid was replaced by a phenyl ring that formed a π-π interaction with Phe31 in compound 44 (Figure 10).

2.2. In silico ADMET Properties

The Maestro’s QikProp module in Schrödinger was applied to predict the drug-likeness, ADME properties, and toxicity (ADMET) of the metabolites under investigation [30]. Table 3 displayed the ADMET properties and other descriptors. In general, most of the predicted properties for the compounds were within the recommended ranges. However, some compounds fell beyond the recommended ranges of certain descriptors. The logP (QPlogPo/w) and binding to human serum albumin (QPlogKhsa) were high for compounds 66, 67, 68, 69, and 70. Additionally, the toxicity was evaluated in terms of the number of reactive functional groups (#rtvFG) and HERG K+ channel inhibition (QPlogHERG). It was observed that none of the compounds exceeded the acceptable range (0–2) of reactive groups. However, 11 compounds (35, 36, 37, 43, 44, 56, 66, 67, 68, 69, and 70) were predicted to inhibit the HERG K+ channel. The results suggested that the high lipophilicity of these compounds is a factor that contributes to HERG inhibition and plasma protein binding. Moreover, no CNS activity was predicted for any compound.

2.3. MD Simulation

We performed MD simulation for the top three scoring compounds from the docking study (34, 28, and 44), as well as to the co-crystallized inhibitor SRI-9662, using Desmond software in Schrödinger [31,32]. The RMSD of proteins (blue) and ligands (red) are presented at the left and right Y-axes of the plot, respectively.
The results showed that the hDHFR protein and the co-crystallized inhibitor SRI-9662 were stable during the 100 ns of the simulation run time (Figure 11A). The fluctuations were insignificant and lay within the acceptable range of 1–3 Å (the differences were within 1 and 1.8 Å for the protein and ligand, respectively). This confers a high-binding and stable protein-ligand complex throughout the run.
In the case of compounds 34, the RMSD plot analysis revealed a stable protein-ligand complex, and the fluctuations of both the protein and ligand RMSD charts were within the acceptable range as well (Figure 12A). A small jump in the RMSD of the protein between 55 and 75 ns was observed and then resumed to its normal level until the end of the run. This change was within the acceptable range.
The RMSD of the protein during its interaction with compound 28 was also within range; however, the compound was less stable inside the binding pocket as the RMSD fluctuated slightly beyond the 1–3 Å range (Figure 13A). A possible explanation would be the absence of the 2-hydroxyethyl group, which was the only difference between the 28 and 34 structures (Figure 8A and Figure 9). The chirality of the OH group and adjacent carbon atom may stabilize the carboxylic acid side chain by restricting its free rotation, leading to a decrease in the number of structural conformations available for binding to the protein active site. The RMSD plot analysis of compound 44 was comparable to 34 in terms of complex stability and RMSD fluctuation range (Figure 14A).
The secondary structure elements (SSE), alpha helices in orange and beta strands in light blue, of hDHFR (PDB: 1KMV) complexed with each ligand, were monitored during the 100 ns simulation time. The results showed that the total %SSE was maintained for all ligands and that most of the individual SSEs (by residue index) were stable throughout the run. Minor differences were observed in the turn and loop regions between the reference inhibitor SRI-9662 and the three investigated metabolites (Figure 11B, Figure 12B, Figure 13B, and Figure 14B), as a few parts of these regions were less stable when each of the metabolites interacted with the protein.
Additionally, the MD analysis displayed the specific contact points between the inhibitor and the amino acid residues in the enzyme active site. Figure 15A showed that SRI-9662 formed strong H-bond interactions with Ile7 and Val115, which lasted about 90% of the simulation time (as demonstrated by the interaction fraction of 0.9).
However, the strongest binding interaction was observed with Glu30, which involved multiple H-bonds with the 2-amino, N-1, and N-8 of the 5-deazapteridine ring either directly or through a water bridge (Figure 15B). The combined effect of these interactions was maintained for about 240% of the run time (Figure 15A). An H-bond through a water bridge was also formed with Trp24. Additional significant hydrophobic interactions with values near or above 1 were observed with Phe31 and Phe34, as well as, to a lesser extent, with Pro61. Detailed ligand atom interactions that lasted more than 30% of the simulation time with selected amino acid residues are presented in Figure 16B. The dimethoxyphenyl ring formed a π-π stacking interaction with Phe31, and the amino groups and nitrogen atoms of the 5-deazapteridine ring H-bonded with several residues, as detailed above. Moreover, the total specific interactions were presented in the top panel of Figure 16C, while the bottom plot showed the ligand-protein interactions by residues in each trajectory frame. The darker and more continuous the orange color is, the stronger the binding, which was observed with the key residues of Ile7, Trp24, Glu30, Phe31, Phe34, and Val115.
It was observed that compound 34 was bound to a different set of amino acid residues in the active site of hDHFR compared to the native inhibitor. The differences could be related to the polar nature of 34. The major types of interactions between 34 and the active residues were the direct and water-aided H-bonds with Gln35, Asn64, Lys68, and Arg70, where the latter was the strongest (> 200%) and involved in additional ionic interaction with the carboxylate anion in the ligand (Figure 16A–C). Minor hydrophobic contacts were formed with Phe31, Phe34, Ile60, and Leu67 with <50% (Figure 16A).
Similar binding modes were observed with compound 28 and the key residues mentioned above (Figure 17A–C), owing to the structural similarity between 28 and 34.
Moreover, it was found that compound 44 bound to amino acid residues of Phe31, Phe34, and Glu30 in the active site (Figure 18A) similar to the reference inhibitor (Figure 16A). It might be due to the hydrophobic nature of 44, in which the phenyl ring in the side chain replaced the carboxylic acid in compounds 34 and 28, which formed strong hydrogen and ionic interactions with the protein. Compound 44 formed hydrophobic interactions with Phe31 and Phe34 and H-bonds mostly in the presence of water molecules with Gly20, Leu27, Glu30, and Ser59 (Figure 18A,B). However, the availability of these interactions did not exceed 60% of the simulation time. The inconsistent binding affinity was also demonstrated by the low total number of specific contacts and incomplete orange lines with each residue in Figure 18C, top and bottom panels, respectively.

3. Materials and Methods

3.1. Ligand and Protein Preparation

The docking study was performed with the Schrodinger program (Schrödinger Release 2022-3: Schrödinger, LLC, New York, NY, USA, 2021). The crystal structure of hDHFR complexed with the SRI-9662 inhibitor was downloaded from the protein data bank (PDB; ID: 1KMV) [33]. The protein was prepared using the “Protein Preparation Wizard” tool in Maestro software, where the missing hydrogens were added to the residues, the metal ionization state was corrected, and the water molecules > 5 Å from protein residues were deleted. The protein was then refined by predicting the pKa of the ionizable residues using PROPKA and water molecules >3 Å (not involved in the water bridge), which were removed [34]. Finally, the protein was minimized by applying the OPLS4 force field. All ligands (sesquiterpene metabolites in addition to the co-crystallized inhibitor SRI-9662) were prepared before docking using the “LigPrep” tool, where their 2D structures were converted to 3D and energy-minimized using the OPLS3 force field [35]. The hydrogens were added, and all possible ionization states and tautomeric forms were created at a pH of 7.0 ± 0.2 by Epik; a desalt option was also chosen. The H-bonds were optimized by predicting the pKa of ionizable groups using PROPKA [34].

3.2. Grid Generation and Molecular Docking

Before docking, a grid box was generated around the active site of hDHFR (PDB: 1KMV) containing the co-crystallized inhibitor SRI-9662, aided by Glide’s “Receptor-Grid-Generation” tool in the Schrödinger suite [36]. The box was built around the co-crystallized ligand by selecting the “centroid of workspace ligand” function. The length of the box in each of the X, Y, and Z dimensions was set by default at 10 Å. All compounds under investigation were docked inside the grid box once with standard precision (SP) and three times with extra precision (XP) protocols (for validation purposes), and all other parameters were set to default [37]. The non-polar atoms were set for the VdW radii scaling factor at 1.0, and the partial charge cutoff was 0.25. Finally, the “Ligand Docking” tool was implemented for docking [38]. To further validate the docking method, the co-crystallized inhibitor was re-docked inside the grid box and evaluated. The docking results were assessed in terms of the gscore (ranks different compounds), emodel (ranks different conformers), and XP gscore. Glide uses emodel scoring to select the best poses of the docked compounds; then, it ranks the best poses based on the given gscores. The XP gscore ranks the poses generated by the XP Glide mode. The XP Glide takes into consideration the major driving forces and structural motifs that contributed to protein-ligand binding affinity, as described in the following equations [37]:
XP Glide Score = Ecoul + EvdW + Ebind + Epenalty
Ebind = Ehyd_enclosure + Ehb_nn_motif + Ehb_cc_motif + EPI + Ehb_pair + Ephobic_pair
Epenalty = Edesolv + Eligand_strain
where:
E is energy (calculated for each of the following descriptors); Ecoul is Coulomb energy, EvdW is van der Waal, Ebind is the energy that favors binding, Epenalty is the penalty that disfavors binding, Ehyd_enclosure is hydrophobic enclosure, Ehb_nn_motif is special neutral-neutral hydrgen-bond motifs, Ehb_cc_motif is special charged-charged hydrogen-bond motifs, EPI is pi-cation interactions, Ehb_pair is hydrogen bond pair, Ephobic_pair is lipophilic pair, and Edesolv is desolvation energy [37].
The Molecular Mechanics Generalized Born Surface Area (MMGBSA) method was used to calculate the binding free energy (ΔG) of the protein-ligand complexes in a solvent based on the following equation:
ΔG binding = G complex − (G protein + G ligand)
where G complex is the free energy of the protein-ligand complex, G protein and G ligand are the free energies of unbound protein and ligand in the solvent, respectively [39,40].

3.3. ADME Properties

The ADME properties (absorption, distribution, metabolism, and excretion) and toxicity of compounds under investigation were predicted using the QikProp module of the Schrodinger suite [30]. This module is useful in predicting the physicochemical properties and other descriptors to facilitate the drug discovery and development process by identifying and eliminating non-drug-like compounds from entering the clinical stage and failing. The predicted descriptors were molecular weight (mol_MW), drug-likeness (#Stars), total solvent accessible surface area (SASA), number of hydrogen bond donors and acceptors (donorHB and acceptHB), predicted octanol-water partitioning (QPlogPo/w), estimated binding to human serum albumin (QPlogKhsa), number of the possible metabolites (# metab), predicted blood-brain partitioning (QPlogBB), percentage of human oral absorption, predicted IC50 for inhibiting HERG-K+ channels (QPogHERG), predicted apparent Caco-2 cell permeability in nm/s for gut-blood barrier (QPPCaco), central nervous system activity (CNS), number of reactive functional groups present (#rtvFG), and percent of human oral absorption. The predicted values were then compared to the recommended range derived from values determined for 95% of known drugs.

3.4. MD Simulation

MD simulations were performed using Desmond in the Schrödinger package [31,32]. The software created a simulated environment that resembles the dynamic nature of the molecular system under physiological conditions to evaluate the virtual stability of interactions between the protein and ligand [41]. The RMSD plots assess the stability of protein-ligand complexes by calculating the deviation of the protein and ligand atoms inside the binding pocket at the end of the simulation time of 100 ns and comparing the results to their initial positions before the simulation at 0 ns [42]. The hDHFR protein was first complexed with the desired ligand in the docking experiment. The protein-ligand complex was then tuned through the “System-Builder” tool to generate the solvated system for simulation. The solvent model was set to TIP3P, the selected box shape was orthorhombic, and the box dimensions were 10 Å. Sodium ions (Na+) were added to neutralize the system. The simulation parameters were set up in the “Molecular Dynamic” tool, where the protein-ligand complex was evaluated at pH 7.0 ± 0.2 over a simulation time of 100 ns. The ensemble class was set as NPT to maintain a constant temperature and pressure of 300 K and 1.01325 bar, respectively, throughout the run. The generated results were analyzed at the end of the MD simulation.

4. Conclusion

Our findings suggested that 34, 28, and 44 possessed potent interactions with the active site of hDHFR compared to SRI-9662 (reference standard), as demonstrated by the docking studies. The MD simulation revealed detailed information about protein-ligand complex stability and specific binding contacts between each ligand and the hDHFR binding site. Additionally, the ADMET prediction demonstrated that all physicochemical parameters and ADMET properties are within the satisfactory range described for human treatment for the majority of sesquiterpene metabolites. These metabolites could be possible leads for DHFRI candidates; however, more in vitro, in vivo, and mechanistic investigation, as well as developing semi- and synthetic derivatives to enhance their DHFRI effectiveness, should be the focus of future research.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28031292/s1, Table S1: In-silico docking results of sesquiterpene metabolites with hDHFR (PDB: 1KMV) compared to the reference inhibitor SRI-9662. In addition to Simulation Interactions Diagram Report for the following compounds SRI-9662 (native ref.), 34, 28, and 44.

Author Contributions

Conceptualization, A.M.O., S.R.M.I., and G.A.M.; methodology, A.M.O., S.R.M.I., K.A.M. and G.A.M.; software, A.M.O., I.A.S., and K.A.M.; validation, A.M.O., S.R.M.I., and G.A.M.; investigation, A.M.O. and K.A.M.; resources, A.M.O., S.R.M.I., K.A.M., I.A.S., and G.A.M.; data curation, A.M.O., S.R.M.I., G.A.M. and K.A.M.; writing—original draft preparation, A.M.O., S.R.M.I., G.A.M. and K.A.M.; writing—review and editing, A.M.O., S.R.M.I., G.A.M., I.A.S., and K.A.M.; project administration, A.M.O.; funding acquisition, A.M.O. All authors have read and agreed to the published version of the manuscript.

Funding

The Deanship of Scientific Research (DSR) at King Abdulaziz University (KAU), Jeddah, Saudi Arabia has funded this project, under grant no. (RG-37-166-43).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the article.

Acknowledgments

The Deanship of Scientific Research (DSR) at King Abdulaziz University (KAU), Jeddah, Saudi Arabia has funded this project, under grant no. (RG-37-166-43). The authors acknowledge and give thanks to DSR for technical and financial support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Osorio, E.; Aguilera, C.; Naranjo, N.; Marín, M.; Muskus, C. Biochemical characterization of the bifunctional enzyme dihydrofolate reductase-thymidylate synthase from Leishmania (Viannia) and its evaluation as a drug target. Biomedica 2013, 33, 393–401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Tobias, A.M.; Toska, D.; Lange, K.; Eck, T.; Bhat, R.; Janson, C.A.; Rotella, D.P.; Gubler, U.; Goodey, N.M. Expression, purification, and inhibition profile of dihydrofolate reductase from the filarial nematode Wuchereria bancrofti. PLoS ONE 2018, 13, e01971732018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Raimondi, M.V.; Randazzo, O.; La Franca, M.; Barone, G.; Vignoni, E.; Rossi, D.; Collina, S. DHFR Inhibitors: Reading the Past for Discovering Novel Anticancer Agents. Molecules 2019, 24, 1140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Hawser, S.; Lociuro, S.; Islam, K. Dihydrofolate reductase inhibitors as antibacterial agents. Biochem. Pharmacol. 2006, 71, 941–948. [Google Scholar] [CrossRef] [PubMed]
  5. Rana, R.M.; Rampogu, S.; Zeb, A.; Son, M.; Park, C.; Lee, G.; Yoon, S.; Baek, A.; Parameswaran, S.; Park, S.J.; et al. In Silico Study Probes Potential Inhibitors of Human Dihydrofolate Reductase for Cancer Therapeutics. J. Clin. Med. 2019, 8, 233. [Google Scholar] [CrossRef] [Green Version]
  6. Srinivasan, B.; Tonddast-Navaei, S.; Roy, A.; Zhou, H.; Skolnick, J. Chemical space of Escherichia coli dihydrofolate reductase inhibitors: New approaches for discovering novel drugs for old bugs. Med. Res. Rev. 2019, 39, 684–705. [Google Scholar] [CrossRef]
  7. Tibon, N.S.; Ng, C.H.; Cheong, S.L. Current progress in antimalarial pharmacotherapy and multi-target drug discovery. Eur. J. Med. Chem. 2020, 188, 111983. [Google Scholar] [CrossRef]
  8. Giletti, A.; Esperon, P. Genetic markers in methotrexate treatments. Pharm. J. 2018, 18, 689–703. [Google Scholar] [CrossRef]
  9. Wróbel, A.; Arciszewska, K.; Maliszewski, D.; Drozdowska, D. Trimethoprim and other nonclassical antifolates an excellent template for searching modifications of dihydrofolate reductase enzyme inhibitors. J. Antibiot. 2020, 73, 5–27. [Google Scholar] [CrossRef] [Green Version]
  10. Chawla, P.; Teli, G.; Gill, R.K.; Narang, R.K. An Insight into Synthetic Strategies and Recent Developments of Dihydrofolate Reductase Inhibitors. ChemistrySelect 2021, 6, 12101–12145. [Google Scholar] [CrossRef]
  11. Kreutzfeld, O.; Tumwebaze, P.K.; Byaruhanga, O.; Katairo, T.; Okitwi, M.; Orena, S.; Rasmussen, S.A.; Legac, J.; Conrad, M.D.; Nsobya, S.L.; et al. Decreased Susceptibility to Dihydrofolate Reductase Inhibitors Associated With Genetic Polymorphisms in Ugandan Plasmodium falciparum Isolates. J. Infect. Dis. 2022, 225, 696–704. [Google Scholar] [CrossRef]
  12. Krucinska, J.; Lombardo, M.N.; Erlandsen, H.; Estrada, A.; Si, D.; Viswanathan, K.; Wright, D.L. Structure-guided functional studies of plasmid-encoded dihydrofolate reductases reveal a common mechanism of trimethoprim resistance in Gram-negative pathogens. Commun. Biol. 2022, 5, 459. [Google Scholar] [CrossRef]
  13. Arya, H.; Coumar, M.S. Chapter 4—Lead identification and optimization. In The Design & Development of Novel Drugs and Vaccines; Bhatt, T.K., Nimesh, S., Eds.; Academic Press: Boston, MA, USA, 2021; pp. 31–63. [Google Scholar] [CrossRef]
  14. Aguayo-Ortiz, R.; Fernández-de Gortari, E. Chapter 2—Overview of Computer-Aided Drug Design for Epigenetic Targets. In Epi-Informatics; Medina-Franco, J.L., Ed.; Academic Press: Boston, MA, USA, 2016. [Google Scholar] [CrossRef]
  15. Nastrucci, C.; Cesario, A.; Russo, P. Anticancer Drug Discovery from the Marine Environment. Recent Pat. Anti-Cancer Drug Discov. 2012, 7, 218–232. [Google Scholar] [CrossRef]
  16. Abdelmohsen, U.R.; Balasubramanian, S.; Oelschlaeger, T.A.; Grkovic, T.; Pham, N.B.; Quinn, R.J.; Hentschel, U. Potential of marine natural products against drug-resistant fungal, viral, and parasitic infections. Lancet Infect. Dis. 2017, 17, e30–e41. [Google Scholar] [CrossRef]
  17. Lu, W.-Y.; Li, H.-J.; Li, Q.-Y.; Wu, Y.-C. Application of marine natural products in drug research. Bioorganic Med. Chem. 2021, 35, 116058. [Google Scholar] [CrossRef]
  18. Pereira, F.; Aires-de-Sousa, J. Computational Methodologies in the Exploration of Marine Natural Product Leads. Mar. Drugs 2018, 16, 236. [Google Scholar] [CrossRef] [Green Version]
  19. Chen, G.; Seukep, A.J.; Guo, M. Recent Advances in Molecular Docking for the Research and Discovery of Potential Marine Drugs. Mar. Drugs 2020, 18, 545. [Google Scholar] [CrossRef]
  20. Marcos, I.S.; Conde, A.; Moro, R.F.; Basabe, P.; Diez, D.; Urones, J.G. Quinone/Hydroquinone Sesquiterpenes. Mini-Rev. Org. Chem. 2010, 7, 230–254. [Google Scholar] [CrossRef]
  21. Kumar, M.; Dagar, A.; Gupta, V.K.; Sharma, A. In silico docking studies of bioactive natural plant products as putative DHFR antagonists. Med. Chem. Res. 2014, 23, 810–817. [Google Scholar] [CrossRef]
  22. Herrmann, F.C.; Sivakumar, N.; Jose, J.; Costi, M.P.; Pozzi, C.; Schmidt, T.J. In Silico Identification and In Vitro Evaluation of Natural Inhibitors of Leishmania major Pteridine Reductase I. Molecules 2017, 22, 2166. [Google Scholar] [CrossRef] [Green Version]
  23. Jose, S.; Devi, S.S.; Al-Khafaji, K. Phytochemical constituents of Inula britannica as potential inhibitors of dihydrofolate reductase: A strategic approach against shigellosis. J. Biomol. Struct. Dyn. 2022, 40, 11932–11947. [Google Scholar] [CrossRef] [PubMed]
  24. Possart, K.; Herrmann, F.C.; Jose, J.; Costi, M.P.; Schmidt, T.J. Sesquiterpene Lactones with Dual Inhibitory Activity against the Trypanosoma brucei Pteridine Reductase 1 and Dihydrofolate Reductase. Molecules 2022, 27, 149. [Google Scholar] [CrossRef]
  25. Herrera-Acevedo, C.; Monroy-Velandia, D.; Flores-Gaspar, A.; Coy-Barrera, E. In Silico Studies to Evaluate Interactions Between Kaurane-Type Diterpenes and the Dihy-drofolate Reductase–Thymilidine Synthase of Three Leishmania Species. In Proceedings of the International Conference on Multidisciplinary Sciences (MOL2NET 2018), 4th edition UFPB, Paraiba, Brasil, 15 January 2018. [Google Scholar]
  26. Kwon, Y.-J.; Sohn, M.-J.; Kim, H.-J.; Kim, W.-G. The Lactone form of stachybotrydial: A new inhibitor of dihydrofolate reductase from stachybotrys sp. FN298. Biol. Pharm. Bull. 2014, 37, 1406–1410. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Alea, G.; Carroll, A.R.; Bowden, B.F. Coscinoquinol, a New Cytotoxic Sesterterpene From a Dictyoceratid Sponge, Coscinoderma sp. Aust. J. Chem. 1994, 47, 191–194. [Google Scholar] [CrossRef]
  28. Zubía, E.; Ortega, M.J.; Luis Carballo, J.; Salvá, J. Sesquiterpene hydroquinones from the sponge Reniera mucosa. Tetrahedron 1994, 50, 8153–8160. [Google Scholar] [CrossRef]
  29. Ibrahim, S.R.M.; Fadil, S.A.; Fadil, H.A.; Hareeri, R.H.; Alolayan, S.O.; Abdallah, H.M.; Mohamed, G.A. Dactylospongia elegans—A Promising Drug Source: Metabolites, Bioactivities, Biosynthesis, Synthesis, and Structural-Activity Relationship. Mar. Drugs 2022, 20, 221. [Google Scholar] [CrossRef]
  30. Klon, A.E.; Héroux, A.; Ross, L.J.; Pathak, V.; Johnson, C.A.; Piper, J.R.; Borhani, D.W. Atomic Structures of Human Dihydrofolate Reductase Complexed with NADPH and Two Lipophilic Antifolates at 1.09Å and 1.05Å Resolution. J. Mol. Biol. 2002, 320, 677–693. [Google Scholar] [CrossRef]
  31. Olsson, M.H.M.; Søndergaard, C.R.; Rostkowski, M.; Jensen, J.H. PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions. J. Chem. Theory Comput. 2011, 7, 525–537. [Google Scholar] [CrossRef]
  32. LigPrep; Schrödinger, LLC: New York, NY, USA, 2021.
  33. Friesner, R.A.; Banks, J.L.; Murphy, R.B.; Halgren, T.A.; Klicic, J.J.; Mainz, D.T.; Repasky, M.P.; Knoll, E.H.; Shelley, M.; Perry, J.K.; et al. Glide:  A New Approach for Rapid, Accurate Docking and Scoring. 1. Method and Assessment of Docking Accuracy. J. Med. Chem. 2004, 47, 1739–1749. [Google Scholar] [CrossRef]
  34. Friesner, R.A.; Murphy, R.B.; Repasky, M.P.; Frye, L.L.; Greenwood, J.R.; Halgren, T.A.; Sanschagrin, P.C.; Mainz, D.T. Extra Precision Glide:  Docking and Scoring Incorporating a Model of Hydrophobic Enclosure for Protein−Ligand Complexes. J. Med. Chem. 2006, 49, 6177–6196. [Google Scholar] [CrossRef] [Green Version]
  35. Glide; Schrödinger, LLC: New York, NY, USA, 2021.
  36. Hou, T.; Wang, J.; Li, Y.; Wang, W. Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics Simulations. J. Chem. Inf. Model. 2011, 51, 69–82. [Google Scholar] [CrossRef] [PubMed]
  37. Koushki, E.H.; Abolghasemi, S.; Mollica, A.; Aghaeepoor, M.; Moosavi, S.S.; Farshadfar, C.; Hasanpour, B.; Feyzi, B.; Abdi, F.; Mirzaie, S. Structure-based virtual screening, molecular docking and dynamics studies of natural product and classical inhibitors against human dihydrofolate reductase. Netw. Model. Anal. Health Inform. Bioinform. 2020, 9, 49. [Google Scholar] [CrossRef]
  38. QikProp; Schrödinger, LLC.: New York, NY, USA, 2021.
  39. Maestro; Schrödinger, LLC: New York, NY, USA, 2021.
  40. Desmond Molecular Dynamics System, Maestro-Desmond Interoperability Tools; Schrödinger: New York, NY, USA, 2021.
  41. Hollingsworth, S.A.; Dror, R.O. Molecular Dynamics Simulation for All. Neuron 2018, 99, 1129–1143. [Google Scholar] [CrossRef] [Green Version]
  42. Leelananda, S.P.; Lindert, S. Computational methods in drug discovery. Beilstein J. Org. Chem. 2016, 12, 2694–2718. [Google Scholar] [CrossRef] [PubMed]
Molecules 28 01292 g001 550
Figure 1. Structures of stachybotrydial, panicein F1, coscinoquinol, cynaropicrin, 3α-cinnamoyloxy-ent-kaur-16-en-19-oic acid, and aulocarpin.
Figure 1. Structures of stachybotrydial, panicein F1, coscinoquinol, cynaropicrin, 3α-cinnamoyloxy-ent-kaur-16-en-19-oic acid, and aulocarpin.
Molecules 28 01292 g001
Molecules 28 01292 g002 550
Figure 2. Chemical structures of compounds 1–24.
Figure 2. Chemical structures of compounds 1–24.
Molecules 28 01292 g002
Molecules 28 01292 g003 550
Figure 3. Chemical structures of compounds 25–44.
Figure 3. Chemical structures of compounds 25–44.
Molecules 28 01292 g003
Molecules 28 01292 g004 550
Figure 4. Chemical structures of compounds 45–70.
Figure 4. Chemical structures of compounds 45–70.
Molecules 28 01292 g004
Molecules 28 01292 g005 550
Figure 5. Chemical structures of compounds 71–87.
Figure 5. Chemical structures of compounds 71–87.
Molecules 28 01292 g005
Molecules 28 01292 g006 550
Figure 6. Superimposition of 3D structures with element color atoms of co-crystallized (carbon atoms customized blue) over redocked (carbon atoms customized green) reference inhibitor SRI-9662 (named Ref_1KMV in the figure).
Figure 6. Superimposition of 3D structures with element color atoms of co-crystallized (carbon atoms customized blue) over redocked (carbon atoms customized green) reference inhibitor SRI-9662 (named Ref_1KMV in the figure).
Molecules 28 01292 g006
Molecules 28 01292 g007 550
Figure 7. Molecular docking of co-crystallized inhibitor SRI-9662 in hDHFR (PDB: 1KMV). (A) 2D representation of binding interactions of SRI-9662 with amino acid residues in the active site within a 3 Å distance; (B) 3D representation of SRI-9662 in green within the hDHFR active site. The H-bond and aromatic-hydrogen interactions are in yellow and cyan dotted lines, respectively. A light blue dotted line represents the π-π stacking between aromatic rings.
Figure 7. Molecular docking of co-crystallized inhibitor SRI-9662 in hDHFR (PDB: 1KMV). (A) 2D representation of binding interactions of SRI-9662 with amino acid residues in the active site within a 3 Å distance; (B) 3D representation of SRI-9662 in green within the hDHFR active site. The H-bond and aromatic-hydrogen interactions are in yellow and cyan dotted lines, respectively. A light blue dotted line represents the π-π stacking between aromatic rings.
Molecules 28 01292 g007
Molecules 28 01292 g008 550
Figure 8. Molecular docking of co-crystallized inhibitor compound 34 in hDHFR (PDB: 1KMV). (A) 2D representation of binding interactions of 34 with amino acid residues in the active site within a 3 Å distance; (B) 3D representation of 34 in green within the hDHFR active site. The H-bonds and ionic interactions are in yellow and purple dotted lines, respectively.
Figure 8. Molecular docking of co-crystallized inhibitor compound 34 in hDHFR (PDB: 1KMV). (A) 2D representation of binding interactions of 34 with amino acid residues in the active site within a 3 Å distance; (B) 3D representation of 34 in green within the hDHFR active site. The H-bonds and ionic interactions are in yellow and purple dotted lines, respectively.
Molecules 28 01292 g008
Molecules 28 01292 g009 550
Figure 9. Molecular docking of co-crystallized inhibitor compound 28 in hDHFR (PDB: 1KMV). (A) 2D representation of binding interactions of 28 with amino acid residues in the active site within a 3 Å distance; (B) 3D representation of 28 in green within the hDHFR active site. The H-bonds and ionic interactions are in yellow and purple dotted lines, respectively.
Figure 9. Molecular docking of co-crystallized inhibitor compound 28 in hDHFR (PDB: 1KMV). (A) 2D representation of binding interactions of 28 with amino acid residues in the active site within a 3 Å distance; (B) 3D representation of 28 in green within the hDHFR active site. The H-bonds and ionic interactions are in yellow and purple dotted lines, respectively.
Molecules 28 01292 g009
Molecules 28 01292 g010 550
Figure 10. Molecular docking of co-crystallized inhibitor compound 44 in hDHFR (PDB: 1KMV). (A) 2D representation of binding interactions of 44 with amino acid residues in the active site within a 3 Å distance; (B) 3D representation of 44 in green within the hDHFR active site. The H-bonds and π-π stacking between aromatic rings are in yellow and light blue dotted lines, respectively.
Figure 10. Molecular docking of co-crystallized inhibitor compound 44 in hDHFR (PDB: 1KMV). (A) 2D representation of binding interactions of 44 with amino acid residues in the active site within a 3 Å distance; (B) 3D representation of 44 in green within the hDHFR active site. The H-bonds and π-π stacking between aromatic rings are in yellow and light blue dotted lines, respectively.
Molecules 28 01292 g010
Molecules 28 01292 g011 550
Figure 11. (A) The RMSD plot for the reference inhibitor SRI-9662 complexed with hDHFR (PDB: 1KMV) over a 100 ns simulation time; (B) Stability of the hDHFR secondary structure over 100 ns of MD simulation when complexed with SRI-9662. Protein secondary structure elements (SSE) like alpha-helices (orange color) and beta-strands (light blue color) were monitored during the simulation. The upper plot reported SSE distribution by residue index across the protein structure, and the plot at the bottom monitored each residue and its SSE assignment over the simulation time.
Figure 11. (A) The RMSD plot for the reference inhibitor SRI-9662 complexed with hDHFR (PDB: 1KMV) over a 100 ns simulation time; (B) Stability of the hDHFR secondary structure over 100 ns of MD simulation when complexed with SRI-9662. Protein secondary structure elements (SSE) like alpha-helices (orange color) and beta-strands (light blue color) were monitored during the simulation. The upper plot reported SSE distribution by residue index across the protein structure, and the plot at the bottom monitored each residue and its SSE assignment over the simulation time.
Molecules 28 01292 g011
Molecules 28 01292 g012 550
Figure 12. (A) The RMSD plot for compound 34 complexed with hDHFR (PDB: 1KMV) over a 100 ns simulation time; (B) Stability of the hDHFR secondary structure over 100 ns of MD simulation when complexed with 34. Protein secondary structure elements (SSE) like alpha-helices (orange color) and beta-strands (light blue color) were monitored during the simulation. The upper plot reported SSE distribution by residue index across the protein structure, and the plot at the bottom monitored each residue and its SSE assignment over the simulation time.
Figure 12. (A) The RMSD plot for compound 34 complexed with hDHFR (PDB: 1KMV) over a 100 ns simulation time; (B) Stability of the hDHFR secondary structure over 100 ns of MD simulation when complexed with 34. Protein secondary structure elements (SSE) like alpha-helices (orange color) and beta-strands (light blue color) were monitored during the simulation. The upper plot reported SSE distribution by residue index across the protein structure, and the plot at the bottom monitored each residue and its SSE assignment over the simulation time.
Molecules 28 01292 g012
Molecules 28 01292 g013 550
Figure 13. (A) The RMSD plot for compound 28 complexed with hDHFR (PDB: 1KMV) over a 100 ns simulation time; (B) Stability of the hDHFR secondary structure over 100 ns of MD simulation when complexed with 28. Protein secondary structure elements (SSE) like alpha-helices (orange color) and beta-strands (light blue color) were monitored during the simulation. The upper plot reported SSE distribution by residue index across the protein structure, and the plot at the bottom monitored each residue and its SSE assignment over the simulation time.
Figure 13. (A) The RMSD plot for compound 28 complexed with hDHFR (PDB: 1KMV) over a 100 ns simulation time; (B) Stability of the hDHFR secondary structure over 100 ns of MD simulation when complexed with 28. Protein secondary structure elements (SSE) like alpha-helices (orange color) and beta-strands (light blue color) were monitored during the simulation. The upper plot reported SSE distribution by residue index across the protein structure, and the plot at the bottom monitored each residue and its SSE assignment over the simulation time.
Molecules 28 01292 g013
Molecules 28 01292 g014 550
Figure 14. (A) The RMSD plot for compound 44 complexed with hDHFR (PDB: 1KMV) over a 100 ns simulation time; (B) Stability of the hDHFR secondary structure over 100 ns of MD simulation when complexed with 44. Protein secondary structure elements (SSE) like alpha-helices (orange color) and beta-strands (light blue color) were monitored during the simulation. The upper plot reported SSE distribution by residue index across the protein structure, and the plot at the bottom monitored each residue and its SSE assignment over the simulation time.
Figure 14. (A) The RMSD plot for compound 44 complexed with hDHFR (PDB: 1KMV) over a 100 ns simulation time; (B) Stability of the hDHFR secondary structure over 100 ns of MD simulation when complexed with 44. Protein secondary structure elements (SSE) like alpha-helices (orange color) and beta-strands (light blue color) were monitored during the simulation. The upper plot reported SSE distribution by residue index across the protein structure, and the plot at the bottom monitored each residue and its SSE assignment over the simulation time.
Molecules 28 01292 g014
Molecules 28 01292 g015 550
Figure 15. (A) Stacked bar graph of hDHFR (PDB: 1KMV) interactions with reference inhibitor SRI-9662 throughout the simulation. The values were normalized over the course of the trajectory. (B) A diagram illustrating the detailed connections of SRI-9662 with hDHFR. The coordination that appears at least 30% of the designated time period (0.00 to 100.00 nsec) will be presented. It is possible to observe more than 100% as some residues might be linked to the same ligand atom more than once; the pink colored arrows represent the hydrogen bond interaction while the % on the arrows represent the coordination % during the simulation time. (C) Timeline of hDHFR- SRI-9662 interactions presented in (A). The top panel presents the total number of specific interactions the protein made with the ligand over the course of the trajectory. The bottom panel presents the residues interacting with the ligand in each trajectory frame. The dark orange color indicates more than one specific interaction is made between some residues and the ligand. As the orange color becomes lighter, the number of contacts between the residue and ligand decreases. The white color indicates no interaction is formed.
Figure 15. (A) Stacked bar graph of hDHFR (PDB: 1KMV) interactions with reference inhibitor SRI-9662 throughout the simulation. The values were normalized over the course of the trajectory. (B) A diagram illustrating the detailed connections of SRI-9662 with hDHFR. The coordination that appears at least 30% of the designated time period (0.00 to 100.00 nsec) will be presented. It is possible to observe more than 100% as some residues might be linked to the same ligand atom more than once; the pink colored arrows represent the hydrogen bond interaction while the % on the arrows represent the coordination % during the simulation time. (C) Timeline of hDHFR- SRI-9662 interactions presented in (A). The top panel presents the total number of specific interactions the protein made with the ligand over the course of the trajectory. The bottom panel presents the residues interacting with the ligand in each trajectory frame. The dark orange color indicates more than one specific interaction is made between some residues and the ligand. As the orange color becomes lighter, the number of contacts between the residue and ligand decreases. The white color indicates no interaction is formed.
Molecules 28 01292 g015
Molecules 28 01292 g016 550
Figure 16. (A) Stacked bar graph of hDHFR (PDB: 1KMV) interactions with compound 34 throughout the simulation. The values were normalized over the course of the trajectory. (B) A diagram illustrating the detailed connections of 34 with hDHFR. The coordination that appears at least 30% of the designated time period (0.00 to 100.00 nsec) will be presented. It is possible to observe more than 100% as some residues might be linked to the same ligand atom more than once; the pink colored arrows represent the hydrogen bond interaction while the % on the arrows represent the coordination % during the simulation time. (C) Timeline of hDRFR-34 interactions presented in A. The top panel presents the total number of specific interactions the protein made with the ligand over the course of the trajectory. The bottom panel presents the residues interacting with the ligand in each trajectory frame. The dark orange color indicates more than one specific interaction is made between some residues and the ligand. As the orange color becomes lighter, the number of contacts between the residue and ligand decreases. The white color indicates no interaction is formed.
Figure 16. (A) Stacked bar graph of hDHFR (PDB: 1KMV) interactions with compound 34 throughout the simulation. The values were normalized over the course of the trajectory. (B) A diagram illustrating the detailed connections of 34 with hDHFR. The coordination that appears at least 30% of the designated time period (0.00 to 100.00 nsec) will be presented. It is possible to observe more than 100% as some residues might be linked to the same ligand atom more than once; the pink colored arrows represent the hydrogen bond interaction while the % on the arrows represent the coordination % during the simulation time. (C) Timeline of hDRFR-34 interactions presented in A. The top panel presents the total number of specific interactions the protein made with the ligand over the course of the trajectory. The bottom panel presents the residues interacting with the ligand in each trajectory frame. The dark orange color indicates more than one specific interaction is made between some residues and the ligand. As the orange color becomes lighter, the number of contacts between the residue and ligand decreases. The white color indicates no interaction is formed.
Molecules 28 01292 g016
Molecules 28 01292 g017 550
Figure 17. (A) Stacked bar graph of hDHFR (PDB: 1KMV) interactions with compound 28 throughout the simulation. The values were normalized over the course of the trajectory. (B) A diagram illustrating the detailed connections of 28 with hDHFR. The coordination that appears at least 30% of the designated time period (0.00 to 100.00 nsec) will be presented. It is possible to observe more than 100% as some residues might be linked to the same ligand atom more than once; the pink colored arrows represent the hydrogen bond interaction while the % on the arrows represent the coordination % during the simulation time. (C) Timeline of hDRFR-28 interactions presented in (A). The top panel presents the total number of specific interactions the protein made with the ligand over the course of the trajectory. The bottom panel presents the residues interacting with the ligand in each trajectory frame. The dark orange color indicates more than one specific interaction is made between some residues and the ligand. As the orange color becomes lighter, the number of contacts between the residue and ligand decreases. The white color indicates no interaction is formed.
Figure 17. (A) Stacked bar graph of hDHFR (PDB: 1KMV) interactions with compound 28 throughout the simulation. The values were normalized over the course of the trajectory. (B) A diagram illustrating the detailed connections of 28 with hDHFR. The coordination that appears at least 30% of the designated time period (0.00 to 100.00 nsec) will be presented. It is possible to observe more than 100% as some residues might be linked to the same ligand atom more than once; the pink colored arrows represent the hydrogen bond interaction while the % on the arrows represent the coordination % during the simulation time. (C) Timeline of hDRFR-28 interactions presented in (A). The top panel presents the total number of specific interactions the protein made with the ligand over the course of the trajectory. The bottom panel presents the residues interacting with the ligand in each trajectory frame. The dark orange color indicates more than one specific interaction is made between some residues and the ligand. As the orange color becomes lighter, the number of contacts between the residue and ligand decreases. The white color indicates no interaction is formed.
Molecules 28 01292 g017
Molecules 28 01292 g018 550
Figure 18. (A) Stacked bar graph of hDHFR (PDB: 1KMV) interactions with compound 44 throughout the simulation. The values were normalized over the course of the trajectory. (B) A diagram illustrating the detailed connections of 44 with hDHFR. The coordination that appears at least 30% of the designated time period (0.00 to 100.00 nsec) will be presented. It is possible to observe more than 100% as some residues might be linked to the same ligand atom more than once; the pink colored arrows represent the hydrogen bond interaction while the % on the arrows represent the coordination % during the simulation time. (C) Timeline of hDRFR-44 interactions is presented in (A). The top panel presents the total number of specific interactions the protein made with the ligand over the course of the trajectory. The bottom panel presents the residues interacting with the ligand in each trajectory frame. The dark orange color indicates more than one specific interaction is made between some residues and the ligand. As the orange color becomes lighter, the number of contacts between the residue and ligand decreases. The white color indicates no interaction is formed.
Figure 18. (A) Stacked bar graph of hDHFR (PDB: 1KMV) interactions with compound 44 throughout the simulation. The values were normalized over the course of the trajectory. (B) A diagram illustrating the detailed connections of 44 with hDHFR. The coordination that appears at least 30% of the designated time period (0.00 to 100.00 nsec) will be presented. It is possible to observe more than 100% as some residues might be linked to the same ligand atom more than once; the pink colored arrows represent the hydrogen bond interaction while the % on the arrows represent the coordination % during the simulation time. (C) Timeline of hDRFR-44 interactions is presented in (A). The top panel presents the total number of specific interactions the protein made with the ligand over the course of the trajectory. The bottom panel presents the residues interacting with the ligand in each trajectory frame. The dark orange color indicates more than one specific interaction is made between some residues and the ligand. As the orange color becomes lighter, the number of contacts between the residue and ligand decreases. The white color indicates no interaction is formed.
Molecules 28 01292 g018
Table
Table 1. Reported promising bio-activities of Dactylospongia elegans sesquiterpenes [29].
Table 1. Reported promising bio-activities of Dactylospongia elegans sesquiterpenes [29].
ActivityCompound Name
Cytotoxicity(−)-Ilimaquinone (1), 5-(+)-epi-ilimaquinone (2), (−)-dactyloquinone B (5), (+)-isospongiaquinone (11), mamanuthaquinone (14), hyatellaquinone (15), (+)-isohyatellaquinone (16), neomamanuthaquinone (18), 9-epi-7,8-dehydrocyclospongiaquinone-2 (20), smenospongine (25), smenospongimine (27), smenospongine b (28), smenospongine c (29), smenospongorine (30), smenospongiarine (32), 5-(+)-epi-smenospongiarine (33), smenospongidine (35), 5-(+)-epi-smenospongidine (36), isosmenospongine (39), nakijiquinone A (40), nakijiquinone B (41), nakijiquinone G (43), 5-epi-nakijiquinone Q (44), (+)-dictyoceratin A (49), (+)-19-methoxy-dictyoceratin-A (50), (+)-dictyoceratin B (51), (+)-dictyoceratin C (52), nakijinol B (65), (+)-dactylospene C (78), dactylospongenone A (81)
Antibacterial(−)-ilimaquinone (1), 5-(+)-epi-ilimaquinone (2), (−)-dactyloquinone A (4), (−)-dactyloquinone B (5), (−)-dactyloquinone C (7), (−)-dactyloquinone D (8), (+)-dactyloquinone E (9), (+)-isospongiaquinone (11), smenospongine C (29), smenospongorine (30), 5-(+)-epi-smenospongidine (36), isosmenospongine (39), nakijiquinone A (40), nakijiquinone B (41), 5-epi-nakijiquinone Q (44), (+)-dictyoceratin A (49), (+)-dictyoceratin B (51), (−)-xishaeleganin C (54), (+)-xishaeleganin d (55), (−)-xishaeleganin B (57), pelorol (64), (−)-dactyltronic acid A (71), (−)-dactyltronic acid B (72),
Anti-inflammatory(+)-dactylospene B (77), (+)-dactylospene C (78)
Antitrypanosomal(−)-ilimaquinone (1), pelorol (64)
Antimalarial(−)-ilimaquinone (1), pelorol (64)
β-Secretase 1 inhibition(−)-ilimaquinone (1), smenospongine (25)
Table
Table 2. In silico docking results of sesquiterpene metabolites with hDHFR (PDB: 1KMV) compared to the reference inhibitor SRI-9662.
Table 2. In silico docking results of sesquiterpene metabolites with hDHFR (PDB: 1KMV) compared to the reference inhibitor SRI-9662.
CompoundsPubChem CID/ChemSpider IDXP gscoreglide gscorePrime Energy
3411744241−12.431−12.431−7441.3
2850994611−11.502−11.502−7452.5
4431130045 *−10.62−10.62−7460.7
Ref_ SRI-9662 −10.432−10.432−7788.8
5121681043−9.874−9.874−7489.4
6171726095−9.71−9.71−7448.6
50-−9.598−9.598−7483.7
8410369050 *−9.341−9.341−7458.4
6550994610−8.932−8.932−7475
36-−8.884−8.884−7459.1
30101245402−8.755−8.755−7453.3
1510361056−8.705−8.705−7455.7
330770979 *−8.592−8.592−7465.9
14495021−8.494−8.494−7475.6
5221589795−8.42−8.42−7522.7
41457734−8.407−8.407−7399.1
8564063250 *−8.358−8.358−7467.9
56-−8.288−8.288−7546.5
221727418−8.256−8.256−7465.1
4210071409 −8.154−8.154−7385.6
1824691897−7.953−7.953−7447.5
82132498497−7.949−7.949−7431.7
80-−7.881−7.881−7569.4
6410067895−7.876−7.876−7478.7
6010316629 *−7.867−7.867−7472.2
1624905924−7.836−7.836−7444.7
55-−7.825−7.825−7459.6
3514286425−7.816−7.816−7435.1
1727023531−7.794−7.794−7421.4
1314526059−7.711−7.711−7477.8
31-−7.711−7.711−7426.8
4324710044 *−7.59−7.59−7373.1
57-−7.58−7.58−7511.6
48101605919−7.56−7.56−7499.8
499885835−7.56−7.56−7499.8
172291−7.5−7.5−7460.7
40457733−7.368−7.368−7374.2
63102015226−7.335−7.335−7514
3844188455−7.31−7.31−7446.5
39-−7.211−7.211−7435.4
3321727419−7.208−7.208−7445.6
23-−7.147−7.147−7460.2
62189164−7.136−7.136−7520.3
79-−7.116−7.116−7556.1
2223424798−7.071−7.071−7443
2125211413−7.062−7.062−7463.5
6-−6.942−6.942−7423.3
83132498496−6.719−6.719−7459.3
253081931−6.583−6.583−7467.6
2610617363−6.567−6.567−7469.5
77-−6.54−6.54−7545.1
78-−6.316−6.316−7553.1
7611811553−6.303−6.303−7523.8
54-−6.274−6.274−7479.5
510915278−6.194−6.194−7435.8
539977451−6.145−6.145−7486.6
912972982−6.101−6.101−7439
7065790999 *−6.055−6.055−7459.8
2950994612−5.927−5.927−7442
3210313302 *−5.704−5.704−7450.8
411035675−5.579−5.579−7423.7
2027023533 *−5.371−5.371−7440.4
24-−5.314−5.314−7422.7
7511090757−5.228−5.228−7534.9
712972980−5.134−5.134−7462.2
81102284910−5.034−5.034−7436.4
10637868−4.462−4.462−7442.6
3773930387 *−4.239−4.239−7446.9
45132606991−3.689−3.689−7417.7
1210066979−3.676−3.676−7363
46132606990−3.433−3.433−7414.9
47107805883 *−3.308−3.308−7429.9
27132578684−3.139−3.139−7406.7
812972981−3.13−3.13−7439.4
5810316627 *−2.619−2.619−7455.7
7310873154−2.199−2.199−7530
7411014966−0.852−0.852−7471.3
72-−0.122−0.122−7578.1
71547297140.840.84−7579.9
* ChemSpider ID.
Table
Table 3. ADMET prediction of the tested sesquiterpenes using QickProp.
Table 3. ADMET prediction of the tested sesquiterpenes using QickProp.
Moleculemol_MW#Stars#rtvFGCNSSASAdonorHBaccptHBQPlogPo/wQPlogHERGQPPCacoQPlogBB#metabQPlogKhsaPercent HumanOral Absorption
Recommended Range(130–725)(0.0–5.0)(0–2)(−2 Inactive) (+2 Active)(300–1000)(0–6)(2.0–20.0)(−2–6.5)Concern Below −5<25 Poor, >500 Great(−3–1.2)(1–8)(−1.5–1.5)(<25% Poor; >80% High)
1358.480.002.000.00590.811.005.503.46−3.661182.68−0.514.000.46100.00
2358.480.002.000.00587.021.005.503.40−3.591076.04−0.544.000.45100.00
3358.480.002.000.00603.351.005.503.45−3.87902.98−0.644.000.49100.00
4356.461.000.000.00573.090.005.503.29−3.501942.90−0.123.000.27100.00
5356.461.000.000.00567.600.005.503.23−3.381659.36−0.173.000.27100.00
6356.461.000.000.00575.080.005.503.26−3.561778.31−0.163.000.27100.00
7356.461.000.000.00570.260.005.503.26−3.501751.98−0.153.000.27100.00
8356.461.000.000.00582.170.005.503.30−3.661521.91−0.223.000.31100.00
9356.461.000.000.00567.440.005.503.24−3.321484.16−0.203.000.29100.00
10356.461.000.000.00587.870.005.503.29−3.781451.01−0.263.000.31100.00
11358.480.002.000.00593.801.005.503.46−3.751173.47−0.525.000.46100.00
12358.480.002.00−1.00621.901.005.503.47−4.21737.72−0.775.000.52100.00
13358.480.002.00−1.00609.271.005.503.39−3.85768.06−0.723.000.50100.00
14358.480.002.000.00604.811.005.503.53−3.781069.35−0.575.000.51100.00
15358.480.002.00−1.00620.061.005.503.45−4.10701.23−0.785.000.52100.00
16358.480.002.00−1.00607.511.005.503.43−3.91760.42−0.716.000.50100.00
17358.480.002.00−1.00607.821.005.503.42−3.93749.67−0.726.000.50100.00
18344.450.002.00−1.00613.231.005.503.30−4.02715.47−0.776.000.46100.00
19356.461.000.000.00598.490.005.503.34−3.941279.73−0.324.000.35100.00
20356.461.000.000.00598.600.005.503.38−3.901348.28−0.304.000.36100.00
21358.481.000.000.00595.100.005.503.30−3.721200.63−0.342.000.35100.00
22358.481.000.000.00596.110.005.503.33−3.741332.29−0.302.000.36100.00
23358.480.002.000.00603.691.005.503.40−3.76861.92−0.663.000.49100.00
24356.461.000.000.00572.490.005.503.24−3.411899.93−0.132.000.27100.00
25343.470.002.00−1.00569.452.505.752.50−3.62347.02−0.964.000.2787.04
26343.470.002.00−1.00561.622.505.752.43−3.46314.54−0.984.000.2685.85
27357.490.002.000.00605.482.005.753.23−3.86834.01−0.675.000.44100.00
28401.500.000.00−2.00647.793.007.752.60−2.3223.48−1.766.000.0266.70
29415.530.002.00−2.00687.222.006.753.68−2.6624.51−1.886.000.3473.34
30399.570.000.00−1.00697.892.004.504.94−4.541304.45−0.754.000.97100.00
31399.570.000.00−1.00665.462.004.504.80−4.021148.84−0.744.000.93100.00
32413.600.000.00−1.00724.462.004.505.26−4.691304.85−0.834.001.06100.00
32413.602.000.00−1.00732.982.004.505.27−4.801146.67−0.914.001.08100.00
33413.600.000.00−1.00712.202.004.505.25−4.541450.92−0.764.001.05100.00
34445.560.000.00−2.00649.733.007.203.39−1.7427.35−1.736.000.2272.53
35447.620.002.00−1.00750.391.007.004.95−5.551120.27−0.795.000.88100.00
36447.620.002.00−1.00738.691.007.004.90−5.31959.48−0.825.000.89100.00
37437.580.002.00−2.00733.503.007.253.79−5.31363.49−1.376.000.5894.97
38343.470.002.00−2.00597.782.505.752.48−4.13216.02−1.255.000.3083.25
39343.470.002.00−1.00569.702.505.752.49−3.65343.89−0.975.000.2686.93
40401.500.000.00−2.00648.043.007.752.60−2.3523.48−1.767.000.0266.69
41443.581.000.00−2.00712.863.007.753.79−2.5458.16−1.497.000.3480.71
42431.530.002.00−2.00673.743.008.452.52−2.3911.60−2.198.00−0.0260.74
43437.580.002.00−2.00728.413.007.253.60−5.16254.87−1.527.000.5691.10
44447.621.002.00−1.00760.651.007.004.89−5.82772.87−0.986.000.89100.00
45355.480.000.000.00605.151.004.503.86−4.031096.57−0.463.000.74100.00
46397.561.000.000.00695.141.004.505.06−4.591512.13−0.513.001.12100.00
47411.581.000.000.00726.021.004.505.47−4.861727.51−0.543.001.23100.00
48372.500.001.00−1.00621.472.003.504.29−3.98613.88−0.864.000.89100.00
49372.500.001.00−1.00621.472.003.504.29−3.98613.88−0.864.000.89100.00
50402.530.001.00−1.00647.582.004.254.54−3.921026.72−0.725.000.89100.00
51388.500.001.00−2.00653.042.003.254.34−4.28283.29−1.315.000.9896.23
52356.500.001.000.00612.171.002.755.02−4.101659.62−0.393.001.10100.00
53386.530.001.000.00650.861.003.505.17−4.121806.41−0.444.001.13100.00
54372.500.000.001.00600.931.003.005.10−3.725101.410.125.001.05100.00
55356.500.000.000.00617.641.002.255.33−3.983852.76−0.015.001.21100.00
56386.532.001.00−2.00824.221.003.506.38−5.911254.62−1.2211.001.43100.00
57390.520.001.00−2.00671.573.004.253.88−4.50304.62−1.334.000.7694.12
58372.500.001.00−2.00663.722.003.504.42−4.62458.92−1.096.000.97100.00
59356.501.001.000.00652.581.002.755.16−4.721254.20−0.585.001.20100.00
60388.500.001.00−2.00670.712.003.254.39−4.54239.87−1.447.001.0295.22
61388.501.001.00−1.00659.931.003.254.69−4.40510.35−0.913.001.19100.00
62314.471.000.001.00559.041.001.504.90−4.013001.200.042.001.11100.00
63314.472.000.001.00563.251.001.504.90−4.043001.810.032.001.11100.00
64372.500.001.00−1.00627.452.003.504.18−4.07621.36−0.773.000.93100.00
65355.481.000.000.00569.692.003.503.91−3.821051.25−0.535.000.67100.00
66623.924.000.000.00896.542.005.258.42−5.211691.25−0.668.002.37100.00
67623.925.000.00−1.00906.242.005.258.43−5.351545.04−0.719.002.38100.00
68624.904.000.000.00761.511.005.257.56−3.601130.21−0.657.002.07100.00
69623.924.000.000.00896.542.005.258.42−5.211691.25−0.668.002.37100.00
70623.924.000.00−1.00905.652.005.258.43−5.301511.67−0.728.002.39100.00
71362.470.001.00−1.00617.020.006.453.05−3.77832.11−0.732.000.08100.00
72362.470.001.00−1.00621.090.006.452.84−4.05724.55−0.842.000.0194.77
73404.500.002.000.00650.040.006.503.49−3.92925.43−0.603.000.35100.00
74404.500.002.000.00643.970.006.503.55−3.761033.49−0.533.000.36100.00
75404.500.002.000.00656.210.006.503.53−4.01893.31−0.633.000.37100.00
76404.500.002.000.00616.220.006.503.46−3.361469.50−0.353.000.27100.00
77400.600.001.000.00671.710.004.705.26−3.911954.56−0.516.000.93100.00
78400.600.001.000.00691.190.004.705.44−4.051954.69−0.536.001.02100.00
79432.640.001.000.00718.220.005.455.47−4.131939.73−0.664.000.91100.00
80432.640.001.000.00750.600.005.455.72−4.491950.56−0.704.001.03100.00
81390.520.000.000.00653.421.005.504.27−4.181628.90−0.494.000.71100.00
82390.520.000.000.00634.281.005.504.08−3.841090.57−0.624.000.69100.00
83390.520.000.000.00618.001.005.504.05−3.471475.16−0.474.000.66100.00
84390.520.000.000.00601.761.005.503.92−3.221276.69−0.504.000.62100.00
85404.501.001.000.00640.061.007.203.48−3.991408.26−0.593.000.32100.00
86390.520.000.000.00686.361.005.504.45−4.671597.54−0.565.000.78100.00
87250.381.000.00−1.00476.681.002.003.84−1.14441.09−0.182.000.3496.78
Abbreviations: molecular weight (mol_MW), drug-likeness (#Stars), total solvent accessible surface area (SASA), number of hydrogen bond donors and acceptors (donorHB and acceptHB), predicted octanol-water partitioning (QPlogPo/w), estimated binding to human serum albumin (QPlogKhsa), number of the possible metabolites (# metab), predicted blood-brain partitioning (QPlogBB), percentage of human oral absorption, predicted IC50 for inhibiting HERG-K+ channels (QPogHERG), predicted apparent Caco-2 cell permeability in nm/s for gut-blood barrier (QPPCaco), central nervous system activity (CNS), number of reactive functional groups present (#rtvFG), and percent human oral absorption.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Omar, A.M.; Mohammad, K.A.; Sindi, I.A.; Mohamed, G.A.; Ibrahim, S.R.M. Unveiling the Efficacy of Sesquiterpenes from Marine Sponge Dactylospongia elegans in Inhibiting Dihydrofolate Reductase Using Docking and Molecular Dynamic Studies. Molecules 2023, 28, 1292. https://doi.org/10.3390/molecules28031292

AMA Style

Omar AM, Mohammad KA, Sindi IA, Mohamed GA, Ibrahim SRM. Unveiling the Efficacy of Sesquiterpenes from Marine Sponge Dactylospongia elegans in Inhibiting Dihydrofolate Reductase Using Docking and Molecular Dynamic Studies. Molecules. 2023; 28(3):1292. https://doi.org/10.3390/molecules28031292

Chicago/Turabian Style

Omar, Abdelsattar M., Khadijah A. Mohammad, Ikhlas A. Sindi, Gamal A. Mohamed, and Sabrin R. M. Ibrahim. 2023. "Unveiling the Efficacy of Sesquiterpenes from Marine Sponge Dactylospongia elegans in Inhibiting Dihydrofolate Reductase Using Docking and Molecular Dynamic Studies" Molecules 28, no. 3: 1292. https://doi.org/10.3390/molecules28031292

APA Style

Omar, A. M., Mohammad, K. A., Sindi, I. A., Mohamed, G. A., & Ibrahim, S. R. M. (2023). Unveiling the Efficacy of Sesquiterpenes from Marine Sponge Dactylospongia elegans in Inhibiting Dihydrofolate Reductase Using Docking and Molecular Dynamic Studies. Molecules, 28(3), 1292. https://doi.org/10.3390/molecules28031292

Article Metrics

Back to TopTop