Journal of Food Bioactives, ISSN 2637-8752 print, 2637-8779 online
Journal website www.isnff-jfb.com

Original Research

Volume 25, March 2024, pages 81-94


In silico ADME and molecular simulation studies of pharmacological activities of phytoconstituents of Annona muricata (L.) Fruit

Figures

Figure 1.
Figure 1. Protein-protein interaction of A. muricata molecular targets.
Figure 2.
Figure 2. Overall molecular target genes network.
Figure 3.
Figure 3. Average rank of kinases and transcription factors across all the library.
Figure 4.
Figure 4. Interaction of the binding poses of PDCD4 with: (a) Annonaine; (b) Asimilobine; (c) Dicaffeoylquinic acid; (d) Dihydrokaempferol-hexoside; (e) fisetin; (f) Kaempferol; (g) Kaempferol-3-O-rutinoside; (h) Luteolin-3,7-di-O-glucoside; (i) Morin; (j) Myricetin; (k) Nornuciferine; (l) Dexamethasone.
Figure 5.
Figure 5. Protein-ligand complex simulation results (a) RMSD of Annonaine and PDCD4 (b) RMSF of PDCD4 on binding to Orientin. (c) Interaction profile of the contact between Annonaine and PDCD4 (d) RMSD of Dexamethasone and PDCD4. (e) RMSF of PDCD4 on binding to Quercetin; (f) Interaction profile of the contact between Dexamethasone and PDCD4.
Figure 6.
Figure 6. A schematic of detailed ligand-protein interactions (a) Annonaine and PDCD4 (b) Dexamethasone and PDCD4.

Tables

Table 1. Predicted pharmacokinetics properties of selected ligands
 
SNLigandsPredicted ADME Parameter from SWISSADME
MWMRTPSA (Å2)Log PESOL Log SESOL ClassGIABBB permeantP-gpCYPs inhibitorLipBSSA
Physicochemical properties: Molecular weight (MW), Molar Refractivity (MR), Total polar surface area (TPSA). Lipophilicity: Consensus Log P. Water Solubility: ESOL Log S, ESOL Class. Pharmacokinetics: Gastrointestinal absorption (GIA), Blood-brain barrier (BBB), P-glycoprotein (P-gp) substrate, Inhibition of Cytochrome P450 (CYPs) type CYP1A2, CYP2C19, CYP2C9, CYP2D6, and CYP3A4. Druglikeness: Lipinski (Lip), Bioavailability Score (BS), Medicinal Chemistry: Synthetic accessibility (SA).
1Annonacin596.88172.67116.456.92−7.13Poorly solubleLowNoNoCYP3A410.557.36
2Annonacin-10-one594.86171.71113.296.83−6.8Poorly solubleLowNoYesCYP3A410.557.18
3Annonaine265.3180.1330.492.88−3.71SolubleHighYesYesCYP1A2, CYP2D6, CYP3A400.553.36
4Cis-Annoreticuin596.88172.67116.456.94−7.13Poorly solubleLowNoNoCYP3A410.557.36
5Asimilobine267.3282.5941.492.65−3.54SolubleHighYesYesCYP1A2, CYP2D6, CYP3A400.553.24
6Cinnamic acid148.1643.1137.31.79−2.37SolubleHighYesNoNONE00.851.67
7Corossolone578.86170.5593.067.75−7.66Poorly solubleLowNoYesNONE20.176.9
8Coumaric acid164.1645.1357.531.26−2.02SolubleHighYesNoNONE00.851.61
9Dicaffeoylquinic acid516.45126.9211.280.91−3.65SolubleLowNoYesNONE30.114.83
10Dihydrokaempferol-hexoside450.39105.12186.37−0.5−2.8SolubleLowNoYesNONE20.175.24
11Epoxymurin-A530.86167.5538.8310.35−10.03InsolubleLowNoYesNONE20.176.47
12Epoxymurin-B530.86167.5538.8310.37−10.03InsolubleLowNoYesNONE20.176.47
13Epomusenin-A558.92177.1738.8311.09−10.75InsolubleLowNoYesNONE20.176.73
14Epomusenin-B558.92177.1738.8311.08−10.75InsolubleLowNoYesNONE20.176.73
15Fisetin286.2476.01111.131.55−3.35SolubleHighNoNoCYP1A2, CYP2D6, CYP3A400.553.16
16Kaempferol286.2476.01111.131.58−3.31SolubleHighNoNoCYP1A2, CYP2D6, CYP3A400.553.14
17Kaempferol 3-O-rutinoside594.52139.36249.2−0.73−3.42SolubleLowNoYesNONE30.176.48
18Luteolin 3,7-di-O-glucoside610.52140.26269.43−1.49−3.22SolubleLowNoYesNONE30.176.39
19Montecristin574.92180.0566.769.92−9.64Poorly solubleLowNoYesNONE20.176.82
20Morin302.2478.03131.361.2−3.16SolubleHighNoNoCYP1A2, CYP2D6, CYP3A400.553.25
21Muricatetrocin B470.64129.41116.453.8−4.1Moderately solubleHighNoNoCYP3A400.556.21
22Muricatocin A612.88173.84136.686.17−6.61Poorly solubleLowNoNoCYP3A410.557.5
23Myricetin318.2480.06151.590.79−3.01SolubleLowNoNoCYP1A2, CYP3A410.553.27
24N-methylcoclaurine299.3690.5252.932.59−3.82SolubleHighYesYesCYP2D600.552.92
25Nornuciferine281.3587.0630.493.01−3.74SolubleHighYesYesCYP1A2, CYP2D6, CYP3A400.553.35
26Reticuline329.3997.0162.162.6−3.88SolubleHighYesYesCYP2D600.553.07
27Sabadelin530.86167.5538.8310.35−10.03InsolubleLowNoYesNONE20.176.47
28Xylomatenin622.92181.81116.457.47−7.43Poorly solubleLowNoNoCYP3A410.557.58

 

Table 2. Target prediction results
 
SNCompoundTarget gene codeTarget descriptionp-valueMTC
1AnnonacinPDCD4Programmed cell death protein 44.236e-740.36
PTGER2Prostaglandin E2 receptor EP2 subtype1.586e-060.32
2Annonacin-10-onePDCD4Programmed cell death protein 46.804e-700.33
IL6STInterleukin-6 receptor subunit beta4.887e-230.30
PPM1AProtein phosphatase 1A7.318e-230.28
SLCO2A1Solute carrier organic anion transporter family member 2A18.104e-180.28
3AnnonaineTAS1R1Taste receptor type 1 member 13.462e-060.29
4Cis-AnnoreticuinPDCD4Programmed cell death protein 44.236e-740.36
5AsimilobineTHTyrosine 3-monooxygenase1.001e-390.32
DRD1D(1A) dopamine receptor1.11e-160.39
MMP26Matrix metalloproteinase-261.384e-070.32
6Cinnamic acidHCAR2Hydroxycarboxylic acid receptor 28.423e-101.00
NR0B2Nuclear receptor subfamily 0 group B member 29.432e-460.31
GPR183G-protein coupled receptor 1832.941e-420.39
SENP2Sentrin-specific protease 23.568e-410.33
CYP1B1Cytochrome P450 1B16.792e-400.48
RCOR3REST corepressor 31.157e-390.31
PAMPeptidyl-glycine alpha-amidating monooxygenase1.372e-380.52
CHRNA10Neuronal acetylcholine receptor subunit alpha-101.097e-350.28
7CorossolonePDCD4Programmed cell death protein 43.7e-330.30
PPM1AProtein phosphatase 1A2.063e-230.29
8Coumaric acidCA3Carbonic anhydrase 31.966e-291.00
CA6Carbonic anhydrase 64.201e-291.00
CA5BCarbonic anhydrase 5B, mitochondrial5.523e-271.00
CA5ACarbonic anhydrase 5A, mitochondrial2.466e-261.00
CA14Carbonic anhydrase 147.845e-211.00
CA7Carbonic anhydrase 79.975e-201.00
AKR1B1Aldo-keto reductase family 1 member B15.551e-161.00
9Dicaffeoylquinic acidNSD2Histone-lysine N-methyltransferase NSD20.00089821.00
CXCL12Stromal cell-derived factor 11.437e-760.38
NR0B2Nuclear receptor subfamily 0 group B member 23.814e-470.33
MYOCMyocilin1.277e-310.39
10Dihydrokaempferol-hexosideTOP1DNA topoisomerase 16.893e-610.36
SLC5A2Sodium/glucose cotransporter 28.932e-580.43
SLC28A3Solute carrier family 28 member 39.656e-530.42
TYRTyrosinase5.641e-360.42
IL2Interleukin-25.276e-300.40
CBSCystathionine beta-synthase2.038e-270.48
11Epoxymurin-APDCD4Programmed cell death protein 41.016e-450.36
LYPLA2Acyl-protein thioesterase 21.4e-220.29
12Epoxymurin-BPDCD4Programmed cell death protein 41.016e-450.36
LYPLA2Acyl-protein thioesterase 21.4e-220.29
13Epomusenin-APDCD4Programmed cell death protein 41.016e-450.36
LYPLA2Acyl-protein thioesterase 21.4e-220.29
14Epomusenin-BPDCD4Programmed cell death protein 41.016e-450.36
LYPLA2Acyl-protein thioesterase 21.4e-220.29
15FisetinELAVL3ELAV-like protein 32.048e-650.66
HSD17B3Testosterone 17-beta-dehydrogenase 31.047e-520.43
ALDH2Aldehyde dehydrogenase, mitochondrial4.45e-500.41
CYP1B1Cytochrome P450 1B16.994e-390.66
16KaempferolELAVL3ELAV-like protein 31.407e-810.78
PTPRSReceptor-type tyrosine-protein phosphatase S5.492e-590.75
CBSCystathionine beta-synthase9.017e-510.45
CREB1Cyclic AMP-responsive element-binding protein 15.513e-450.30
P4HBProtein disulfide-isomerase6.101e-410.52
17Kaempferol 3-O-rutinosideP4HBProtein disulfide-isomerase3.535e-711.00
SLC28A3Solute carrier family 28 member 31.353e-430.35
ELAVL3ELAV-like protein 32.773e-430.39
IL2Interleukin-21.262e-310.42
TYRTyrosinase9.763e-300.34
XDHXanthine dehydrogenase/oxidase1.157e-250.65
18Luteolin 3,7-di-O-glucosideCREB1Cyclic AMP-responsive element-binding protein 15.266e-690.35
SLC5A2Sodium/glucose cotransporter 23.962e-600.41
SLC28A3Solute carrier family 28 member 33.077e-500.40
IL2Interleukin-25.836e-400.89
19MontecristinPDCD4Programmed cell death protein 41.128e-770.41
SLCO2A1Solute carrier organic anion transporter family member 2A12.172e-380.32
POLHDNA polymerase eta3.331e-160.31
20MorinPTPRSReceptor-type tyrosine-protein phosphatase S6.131e-561.00
MPGDNA-3-methyladenine glycosylase1.089e-501.00
DAPK1Death-associated protein kinase 11.344e-211.00
ELAVL3ELAV-like protein 38.402e-820.68
P4HBProtein disulfide-isomerase3.742e-360.44
21Muricatetrocin BPDCD4Programmed cell death protein 42.688e-590.33
22Muricatocin APDCD4Programmed cell death protein 41.009e-810.40
23MyricetinELAVL3ELAV-like protein 32.196e-1011.00
XDHXanthine dehydrogenase/oxidase5.844e-361.00
24N-methylcoclaurineDRD1D(1A) dopamine receptor1.982e-950.58
DHCR77-dehydrocholesterol reductase2.187e-430.50
TUBB1Tubulin beta-1 chain5.991e-260.35
ABCB1ATP-dependent translocase ABCB19.037e-260.45
SLC18A2Synaptic vesicular amine transporter4.29e-240.35
25NornuciferineTUBBTubulin beta chain4.386e-250.31
26ReticulineDRD1D(1A) dopamine receptor9.256e-980.56
DHCR77-dehydrocholesterol reductase3.993e-440.51
TUBB1Tubulin beta-1 chain1.432e-430.40
FSHRFollicle-stimulating hormone receptor1.151e-290.30
NDUFA1NADH dehydrogenase [ubiquinone] 1 alpha subcomplex subunit 15.531e-170.32
27SabadelinPDCD4Programmed cell death protein 41.016e-450.36
LYPLA2Acyl-protein thioesterase 21.4e-220.29
28XylomateninPDCD4Programmed cell death protein 49.378e-700.41
SLCO2A1Solute carrier organic anion transporter family member 2A11.494e-360.31

 

Table 3. Molecular docking results
 
S.NPhytochemicalsPubChem CIDPDCD4 (AlphaFold ID: AF-Q53EL6) Binding Affinity ΔG (kcal.mol−1)
Docking parameter: PDCD4 [spacing: 0.525, box size: 126 × 126 × 126, center: −9.248 × 0.161 × 1.314].
1Annonacin354398−4.741
2Annonacin-10-one180161−4.591
3Annonaine160597−7.294
4Cis-Annoreticuin72778911−4.869
5Asimilobine160875−6.409
6Cinnamic acid444539−3.96
7Corossolone4366126−4.615
8Coumaric acid637542−5.183
9Dicaffeoylquinic acid12358846−6.775
10Dihydrokaempferol-hexoside10478918−7.012
11Epoxymurin-A5281161−4.416
12Epoxymurin-B131752983−3.789
13Epomusenin-A10507050−3.842
14Epomusenin-B10698082−3.614
15Fisetin5281614−6.979
16Kaempferol5280863−6.749
17Kaempferol 3-O-rutinoside5318767−6.629
18Luteolin 3,7-di-O-glucoside5490298−7.65
19Montecristin102083640−3.861
20Morin5281670−6.884
21Muricatetrocin B393472−5.362
22Muricatocin A133072−5.191
23Myricetin5281672−6.538
24N-methylcoclaurine440595−5.946
25Nornuciferine41169−6.226
26Reticuline439653−5.907
27Sabadelin101006011−3.426
28Xylomatenin10484035−4.5
STDDexamethasone5743−6.682

 

Table 4. Prime MMGBSA binding energy of interaction of PDCD4 with Annonaine and Dexamethasone respectively, before and after molecular dynamics simulation.
 
ComplexSimulation Time (ns)MMGBSA ΔGbind (kcal.mol−1)
CoulombCovalentHbondLipoPackingSolv_GBvdWGbind (Total)
Covalent: Covalent binding energy; Coulomb: Coulomb energy; Lipo: Lipophilic energy; Hbond: Hydrogen bonding energy; Packing: Pi-pi packing correction; Solv GB: Generalized Born electrostatic solvation energy; vdW: Van der Waals energy; Total: Total energy (Prime energy).
Annonaine and PDCD40−64.4282.791−0.60−14.323−0.54169.403−28.149−35.851
100−52.4711.676−0.797−15.184−1.51952.998−23.720−39.019
Dexamethasone and PDCD400.9130.092−0.420−7.605011.380−32.850−28.489
100−4.693−0.012−0.027−8.555015.696−30.694−28.284