Submitted : 29-08-2022 Revised : 20-10-2022 Accepted : 28-11-2022 Trad. Med. September-December 2022 Vol. , p 218-230 ISSN-p : 1410-5918 ISSN-e : 2406-9086 Antifungal Activities of Phytochemicals from Annona muricate (Sour So. : Molecular Docking and Chemoinformatics Approach Misbaudeen Abdul-Hammed1*. Ibrahim Adedotun Olaide1. Hadijat Motunrayo Adegoke1,2. Monsurat Olajide1,3. Oluwasegun Johnson Osilade1. Tolulope Irapada Afolabi1. Adelayo Idayat Abdul-Hammed1,4 1 Computational and Biophysical Chemistry Unit. Department of Pure and Applied Chemistry. Ladoke Akintola University of Technology. LAUTECH. Ogbomoso. Oyo State. Nigeria 2 Genomics Unit. Helix Biogen Institute. Ogbomoso. Oyo State. Nigeria 3 Department of Chemical Sciences. Crescent University. Abeokuta. Ogun State. Nigeria 4 Department of Computer Science. Ladoke Akintola University of Technology. LAUTECH. Ogbomoso. Oyo State. Nigeria ABSTRACT Fungal infection has become a persistent problem in humans and is sometimes life-threatening in immune-compromised individuals. This work aims to study phytochemicals from Annona muricata . our so. as probable antifungal agents against Candida albicans sterol 14-demethylase target receptor by Computer Aided-Drug Design (CADD) approach using voriconazole and fluconazole as standard drugs. modern method of drug discovery by molecular docking and chemoinformatics was used to screen 131 isolated phytochemicals with medicinal properties from Annona muricata against Candida albicans Aosterol 14-demethylase, a prominent target receptor for most anti-fungal drugs, towards the development of new anti-fungal therapeutic agents and a new approach to treat patients with fungal infections. The compounds were all subjected to analyses like ADMET, drug-likeness, bioactivity, oral-bioavailability and PASS. The results of the docking simulation and chemoinformatics analyses showed that muricin M (-7. 9 kcal/mo. , chlorogenic acid (-8. 2 kcal/mo. , roseoside (-8. 5 kcal/mo. and caffeoylquinic acid (-8. 1 kcal/mo. are potential drug candidates for treating fungal infections due to their excellent properties such as binding affinities. ADMET profile, drug-likeness, bioactivity, binding mode and interactions with the target receptor. Thus, muricin M, chlorogenic acid, roseoside and caffeoylquinic acid are recommended for further analyses towards the development of further antifungal drugs. Keywords: Chemoinformatics. Phytochemicals. Annona muricata. Fungal infection. Candida albicans. INTRODUCTION Fungal infections affect a wide diversity of populations and are caused by many different fungal pathogens from various sources (Bongomin et al. , 2. Many medical specialists come across patients with fungal infections, including general practitioners like paediatricians, dermatologists, ophthalmologists, oncologists, haematologists, intensive-care-unit AIDS otolaryngologists, which complicate the provision of holistic education about fungal infections (Tudela and Denning, 2. The rate of mortality due to fungal infections and illnesses is currently at its peak. Severe fungal infections often arise because of other health issues, including acquired immunodeficiency syndrome (AIDS), cancer, asthma, diabetes, organ transplantation, and treatment with corticosteroids (Al Aboody and Mickymaray, 2. There are different types of *Corresponding author : Misbaudeen Abdul-Hammed Email : mabdul-hammed@lautech. fungal diseases such as fungal nail infections, vaginal candidiasis, ringworm, and so on. Some serious and life-threatening fungal diseases such as aspergillosis. Candida auris infection, and invasive candidiasis (Kohler et al. , 2. also affect humans with a weakened immune system . atients with HIV, cancer, organ transplants, or certain medication. (Al Aboody and Mickymaray, 2. In an attempt to treat fungal illnesses, four groups of antifungal medications are frequently available to treat fungal infections such as azoles . lotrimazole, . aspofungin, micafungin, and anidulafungi. , flucytosine . -fluorocytosin. and polyenes . imaricin, trichomycin, candicidin, amphotericin B, nystatin, methyl partrici. (Al Aboody and Mickymaray, 2. Although, the global issue of antibiotic resistance does not exclude fungus Current antifungal drugs are imperfect, and new approaches are vital including those that support patients-host defences (Tudela and DOI: 10. 22146/mot. 77380 | Traditional Medicine Journal, 27. , 2022 This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4. 0 International License. Antifungal Activities of Phytochemicals from Annona muricate Denning, 2. Sterol 14-demethylase (CYP51. P45014DM. Erg. is a member of the cytochrome P450 superfamily, which catalyzes the oxidative removal of the 14-methyl group (C-. of 15-desaturated intermediates in ergosterol biosynthesis (Sheng et , 2. Candida albicans CYP51 consists of 528 amino acids . , including the 48-amino acidlong N-terminal membrane anchor sequence (Hargrove et al. , 2. Over the past century, the phytochemicals in plants have been a channel for pharmaceutical discovery (Moghadamtousi et al. , 2. Natural products from plants have been used to help sustain the health of mankind since the dawn of The importance of the active ingredients of plants in agriculture and medicine has stimulated significant scientific interest in the biological activities of these substances due to the increasing failure of chemotherapeutic agents, and antibiotic resistance exhibited by pathogenic organisms (Moghadamtousi et al. , 2. Researchers are increasingly turning their attention to traditional medicine, screening several medicinal plants for their potential antimicrobial activities for new leads to develop a better drug against microbial infections (Moghadamtousi et al. Plants with a lengthy account of ethnomedical usage are regarded as abundant sources of vital phytoconstituents that offer therapeutic or health advantages against a variety of illnesses and diseases. Annona muricata (Sourso. is not left out among the plants with great traditional importance (Moghadamtousi et Annona muricata contains varieties of bioactive compounds and has a very potent antimicrobial ability that can be associated with rich phytoconstituents like phenols, flavonoids, alkaloids, and acetogenins among others. All parts of Annona muricata including bark, leaves, roots, and fruits are used in natural medicine (Agu and Okolie, 2. The use of computational chemistry in drug design and discovery cannot be over-emphasized. Various studies have reported the efficacy of the computational approach in proposing potential drug leads in different phytochemicals against various diseases including coronavirus (Aanouz et , 2020. Abdul-Hammed et al. , 2021. Falade et al. , breast cancer (Abdul-Hammed et al. Alzheimer's diseases among others with higher success rates. The method is quite efficient and cost-effective in the identification of potential lead compounds in drug development. In the quest to propose a promising treatment for fungi infections with the aid of natural products, this Traditional Medicine Journal, 27. , 2022 study demonstrated the inhibitory potentials of 131 phytochemicals from Annona muricata against Candida albicans serol 14-demethylase (PDB ID: 5TZ. target receptor by using the Computer Aided-Drug Design (CADD) approach. METHODOLOGY Material The software used are PyRx virtual screening tool. BIOVIA Discovery studio visualizer. Pymol. CASTp and SpartanAo 14. Methods Preparation of the Ligands One hundred and thirty-one . from Annona muricata and two standard drugs . luconazole and voriconazol. were used for this research. The two-dimensional structures of all the ligands and standards were retrieved from the PubChem database . ttps://pubchem. and converted to 3D structures using SpartanAo 14 To obtain ligands and standards with the most stable conformation for docking simulation, a conformer search was carried out using Conformer Distribution (Molecular Mechanics/ MMFF) set on SpartanAo14 software. The most stable conformers obtained were optimized using the density functional theory method (DFT) at B3LYP and 631G* basis set, to have structures with the best equilibrium geometry before molecular docking Preparation of the Target Receptor The 3D crystal structure of Sterol 14demethylase in complex with CYP51 inhibitor (PDB ID: 5TZ. (Figure . was retrieved from Protein Data Bank (PDB) . with a resolution of 2. 00 yI (Hargrove et al. , 2. recommended for a protein of excellent quality for molecular docking. The protein was prepared by removing water molecules, ligands, and other heteroatom residues present in the structure to avoid unwanted interactions and interference with the active site, and the Ramachandran Plot (Figure . of the protein was obtained using Discovery Studio Software 2019. The binding pocket of the initial inhibitors present in 5TZ1 was used in determining the binding parameters as -28. 22, and 33. 22 for x, y, and z respectively using the grid box of the inhibitor complexed with the target receptor. Identification and Validation of the Active Sites Computed Atlas for Surface Topography of Proteins (CAST. , . ttp://sts. edu/castp/ html?2. (Tian et al. , 2. , and Biovia Misbaudeen Abdul-Hammed Figure 1. Crystal structure of Sterol 14-demethylase (PDB 1D: 5TZ. The alpha helices, beta sheets, and amino acid chain turns are in magenta, yellow and pale blue colours, respectively, while all other amino acid residues are coloured white. Figure 2. The Ramachandran plot of Sterol 14-demethylase (PDB: 5TZ. from VADAR (Volume. Area, and Dihedral Angle Reporte. web server. The colour codes in the plot correspond to conformations where atoms in the polypeptide come closer than the sum of their van der Waals radii which is sterically disallowed . , regions where there are no steric clashes . and the allowed regions with alpha-helical and beta-sheet conformations . Discovery Studio . was used to determine the binding pocket, amino acids and all the ligand interactions in the active site of 5TZ1 receptor. Result obtained was then validated with the source journal of the protease from the protein data bank . (Hargrove, et al. Molecular Docking The most stable optimized ligands and standards were docked with the target receptor using PyRx virtual screening software in duplicate (Dallakyan and Olson, 2. The mean average binding affinities and standard deviation were Traditional Medicine Journal, 27. , 2022 Antifungal Activities of Phytochemicals from Annona muricate The binding affinities obtained were used to calculate the inhibition constant (K. of the docked ligands and standards against the target receptor as shown in the equations below. iG = OeRT lnKia . Ki = exp(OeiG/RT)a. A . where R = Gas constant . 987 y 10-3 kcal/mo. T = 298. 15 K . bsolute temperatur. Ki = Inhibition constant, iG = Binding energy. The docking results were viewed and analyzed using Biovia Drug discovery studio 2019 and PyMol software. RESULT AND DISCUSSION Structural and active site analysis of sterol 14Demethylase (PDB ID: 5TZ. The active sites of sterol 14Demethylase (PDB ID: 5TZ. are located at the heme loop that contains the following residues: Tyr-64 Tyr-118 Leu-121 Thr-122 Phe-126 Ile-131 Tyr-132 Phe-228 Pro-230 Phe-233 Gly-303 Ile-304 Gly-307 Gly-308 Thr-311 Leu-376 His-377. H-bond Ser-378 Phe-380 Tyr-505 Ser-507 Met-508 (Hargrove et al. , 2. binding affinity, the lower the inhibition constantAy, and vice-versa (Ferreira et al. , 2. The 131 ligands from Annona Muricata docked against sterol 14-demethylase receptor using PyRx virtual screening tool include eighty-three . acetogenins, fourteen . alkaloids, three . cyclopeptides, ten . flavonoid triglycosides, thirteen . megastigmanes and eight . The docking scores and inhibition constants are shown in Table I. The binding affinities of the acetogenins range from -8. 3 to -6. 6 kcal/mol, with muricin L and muridienin 4 having the highest and lowest values, respectively. Thirty-one of the acetogenins have better inhibitory activities than the standard voriconazole (SD. and fluconazole (SD. The binding affinities of the alkaloids range 9 to -7. 5 kcal/mol, with anonaine and atherospemine having the highest and the lowest values respectively (Table I). Twelve alkaloids have better inhibitory activities than SD1 and SD2. The binding affinities of the cyclopeptides range 5 to -6. 6 kcal/mol, with annomuricatin C and annomuricatin B having the highest and the lowest values respectively. Two cyclopeptides have better inhibitory activities than SD1 and SD2. Flavonoid triglycosides have binding affinities ranging from -9. 4 to -5. 7 kcal/mol, with kaempferol-3-0-rutinoside and gallic acid having the highest and the lowest values, respectively, while nine of the flavonoid triglycosides have better inhibitory activities than SD1 and SD2. Also, the binding affinities of the docked megastigmanes and the phenolics range from -9. 3 to -6. 1 kcal/mol, 8 to -6. 1 kcal/mol, respectively with kaempferol-3-0-rutinoside and anoinol C having the highest and the lowest values, respectively for kaempferol-3-0rutinoside and anoinol C for the phenolics. In each case, five of the docked ligands have better inhibitory activities than SD1 and SD2. However, binding affinity is not enough to validate the reliability of the ligand as a potential drug candidate against the studied fungal infection. Therefore, all the ligands were subjected to other chemoinformatics analyses and toxicity tests. Molecular Docking Analysis Molecular docking has shown great advancement in the predictions of therapeutic interventions by predicting bound conformations and free energies of binding for small-molecule ligands to macromolecular targets . (ElHachem et al. , 2. It is used in the evaluation of the binding mode of inhibitors and receptors to guide the optimization of lead compounds (Dong et , 2. The rule of thumb is Authe higher the Assessment of the ADMET Analysis Absorption, excretion, and toxicity (ADMET) prediction can guide the selection and optimization of lead compounds in the early stage of drug development. In the parameter settings, the aqueous solubility, blood-brain barrier penetration, hepatotoxicity, human intestinal absorption (HIA), and plasma protein binding were derived from the ADMETSAR2 web server (Dong et al. , 2. ADMET and Drug-Likeness Studies The . bsorption, distribution, metabolism, excretion, and toxicity. ADMET) of the ligands were predicted admetSAR2 . ttp://lmmd. cn/admetsar2/) drug-likeness prediction of ligands was carried out by Lipinski filter using molecular inspiration server . ttps://w. com/) (Dain a et al. , 2. Prediction of Activity Spectra for Substances (PASS) and Oral Bioavailability Assessment The biological activities of the ligands from Annona muricata were analyzed for their antifungal property using the PASS software, . ttp://pharmexpert. ru/passonline/) while other properties related to the oral bioavailability properties of the ligands were obtained Swiss-ADME . ttp://w. ch/). Traditional Medicine Journal, 27. , 2022 Misbaudeen Abdul-Hammed Table Ia. Ligands from Annona muricata and Standard drugs showing docking scores and the Inhibition Constants of the Interaction of Ligands with the Crystal structure of Sterol 14 demethylase (PDB ID: 5TZ. Ligands Muricin L Annocatacin A Arianacin Cis-annonacin-10-One Gigantetronenin Muricin G Muricin E Annonacin Annopentocin A Muricin M Goniothalamicin Muricatocin C Annocatalin Annomontacin Annomuricin A Cis-Annomontacin Corepoxylone Javoricin Annonancin-10-One Annopentocin B Cis-goniothalamicin Cis-reticulatacin-10-one Cis-solamin Cis-Uvariamicin 1 Muricin K Murihexocin A Murisolin Muricatetrocin A Muricatetrocin B Muricatocin A Muricin C Anonaine Isolaureline Coreximine Xylopine Anomuricine Anomurine Annonamine Annomuricatin C Annomuricatin A Binding Affinity . G), kcal/mol Inhibition (K. AAM Ligands Acetogenins Muricin B Muricin H Xylomaticin Muricapentocin Muricin A Muricin D Muricin I Annomuricin C Annonacin A Annoreticuin-9-One Cis-Corrosolone Cis-Uvariamicin Longifolicin Annomutacin Cis-Panatellin Corossolin Corossolone Muricin J Annomuricin B Gigantetrocin B Isoannonacinone Annomuricin E Muricahexocin A Muricatocin B Acetogenins (ContAo. 4A0. Murihexocin B 4A0. Cis-Annonacin 4A0. Cohibin D 4A0. Epomuricenin A 4A0. Epomusenin B 4A0. Isoannonacin 4A0. Muricin N Alkaloids 9A0. (R)-4-O-methyl 6A0. Nornuciferine 5A0. Reticuline 4A0. Stepharine 2A0. Asimilobine 2A0. (S)-Norcorydine 2A0. Atherospemine Cyclopeptides 5A0. 25A0. 3A0. 1A0. 0A0. 0A0. 0A0. 0A0. 0A0. 9A0. 9A0. 9A0. 9A0. 9A0. 8A0. 8A0. 8A0. 8A0. 8A0. 8A0. 7A0. 7A0. 7A0. 7A0. 7A0. 7A0. Binding Affinity . G). Kcal/Mol Inhibition (K. AAM 7A0. 7A0. 7A0. 7A0. 7A0. 7A0. 7A0. 6A0. 6A0. 6A0. 6A0. 6A0. 6A0. 5A0. 5A0. 5A0. 5A0. 5A0. 5A0. 5A0. 5A0. 4A0. 4A0. 4A0. 4A0. 3A0. 3A0. 3A0. 3A0. 3A0. 3A0. 1A0. 0A0. 0A0. 9A0. 7A0. 6A0. 5A0. Traditional Medicine Journal, 27. , 2022 Antifungal Activities of Phytochemicals from Annona muricate Table Ib. Ligands from Annona muricata and Standard drugs showing docking scores and the Inhibition Constants of the Interaction of Ligands with the Crystal structure of Sterol 14 demethylase (PDB ID: 5TZ. Ligands Binding Affinity . G), kcal/mol Kaempferol-3-0Rutinoside Quercetin3-0-Rutinosid Quercetin3-ONeohesperidoside Quercertin 3-0Glucoside Quercetin 4A0. Dicaffeoylquinic acid 4-Feruloyl-5caffeoylquinic Acid Dihydrokaempferolhexoside 8A0. 3A0. Voriconazole (SD. Fluconazole (SD. Inhibition (K. AAM Ligands Flavonoid triglycosides Chlorogenic acid Binding Affinity . G). Kcal/Mol Inhibition (K. AAM 2A0. 4A0. 3A0. Epicatechine Catechine 1A0. 0A0. 5A0. Kaempferol 9A0. 2A0. 9A0. Megastigmanes Caffeoylquinic acid Feruloylglycoside 1A0. 7A0. Standard Antiviral drugs 65A0. 3A0. ADMET profile predicts that a drug molecule should have good human intestinal absorption ( HIA), ability to cross the blood-brain barrier ( . , solubility (Log S) within the range of -1 to -5, non-inhibitors of cytochrome enzyme (P. , and toxicity in terms of Ames mutagenesis should be negative or non-Ames toxic. Also, it must exhibit non-carcinogenicity, non-inhibition of hERG, and no or low level of toxicity (Tsaioun and Kates. Table II presents the ADMET analysis result. Out of the 131 ligands subjected to ADMET analysis, only 16 lead compounds with an inhibition constant estimated to be O 1. 0 were subjected to ADMET analysis across the six . cetogenins, megastigmanes, and phenolic. from Annona The sixteen ligands are non-genotoxic and non-carcinogenic. Also, all the ligands possess type i . lightly toxi. and they could be easily converted to type IV . on-toxi. during lead The human ether-a-go-go- related gene . ERG) functions as repolarization of cardiac, blockage of this could result from some molecules present in drug or inherited mutation thereby leading to QT syndrome and eventual death Traditional Medicine Journal, 27. , 2022 (Sanguinetti Tristanifirouzi. Fortunately, all 16 ligands were non-inhibitors of hERG. The Ames toxicity value reveals the potential of a drug molecule to cause mutation in DNA and could be a reason for excluding a drug molecule during the discovery process (Falade et al. , 2. As part of the major reason for screening, all the 16 ligands passed Ames mutagenicity with negative All the ligands were well absorbed in the human intestine and they possess the ability to cross the blood-brain barrier and also have aqueous solubility values within the recommended range showing that the ligands possess excellent absorption and distribution properties. Similarly, all 16 ligands are non-inhibitors of microsomal enzymes which indicates good metabolism of the drug . ttp://lmmd. cn/admetsar2/). Therefore, all 16 lead ligands have excellent ADMET properties and were subjected to further analysis for continuous validation of their potential as probable inhibitors of the receptor under study. Drug-Likeness Analysis The most common approach in drug discovery is to use some type of molecular descriptors linked with pattern recognition or interpolation technique, to distinguish between a Misbaudeen Abdul-Hammed Table II. ADMET analysis of phytochemicals from Annona muricate Absorption &Distribution b HIA ( /-) ( /-) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) (-) ( ) (-) ( ) ( ) ( ) (-) ( ) (-) ( ) (-) (-) (-) ( ) ( ) ( ) Ligands L10 L11 L12 L13 L14 L15 L16 SD1 SD2 Metabolism (CYP450 Inhibitor. 2C19 (LogS) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) (-) Extn. Toxicity 2D6 1A2 ( /-) AOT EI EC hI i i i i i i i i i i i i i i i i i i b = Blood Brain Barrier. HIA= Human Intestinal Absorption. AS= Aqueous Solubility. Extn. = Excretion. B= Biodegradation ( /-) Biodegradable ( ). Non-biodegradable (-). AM= Ames mutagenesis ( /-) Acute toxicity ( ). Non-Toxic (-). AOT= Acute Oral Toxicity. EC= Eye corrosion. EI=Eye Irritation. hI= Human either-a-gogo inhibition. C= Carcinogenicity. L1=Annocatacin A. L2=Arianacin. L3=Cis-annonacin-10-one. L4=Muricin E. L5=Annonacin. L6=Annopentocin A. L7=Muricin m. L8=Goniothalamicin. L9=Muricatocin C. L10= Coreximine. L11=Annomuricatacin A. L12= Chlorgenic acid. L13= Roseoside. L14= Dicaffeoylqunic acid. L15= Dihydrokaempferol-hexoside. L16= Caffeoylquinic acid. SD1= Voriconazole. SD2= Fluconazole data set of drugs and one of the non-drugs. Druglikeness of phytochemicals follows the pioneering Aurule of fiveAy by Lipinski (Brustle et al. , 2. The following rules must be obeyed with one . violation at most. For a molecule to obey Authe rule of fiveAy, it must exhibit molecular weight (MW) O 500 Da as an oral bioavailability criterium, hydrogen bond donor (HBD. O 5, hydrogen bond acceptor (HBA. O 10, and LogP . ctanol-water partition coefficien. O 5. These descriptors of oral Traditional Medicine Journal, 27. , 2022 Antifungal Activities of Phytochemicals from Annona muricate Table i. Drug-likeness Analysis Ligands Molecular Formula Muricin M C24H42O7 Coreximine C19H21NO4 Chlorogenic acid C16H18O9 Caffeoylquinic acid C16H18O9 Roseoside C19H30O8 Voriconazole C16H14F3N5O Fluconazole C13H12F2N6O LipinskiAos rule of five Properties Value Molecular weight(O500D. LogP(O. H-Bond donor(O. H-Bond acceptor(O. Violation Molecular weight(O500D. LogP(O. H-Bond donor(O. H-Bond acceptor(O. Violation Molecular weight(O500D. LogP(O. H-Bond donor(O. H-Bond acceptor(O. Violation Molecular weight(O500D. LogP(O. H-Bond donor(O. H-Bond acceptor(O. Violation Molecular weight(O500D. LogP(O. H-Bond donor(O. H-Bond acceptor(O. Violation Molecular weight(O500D. LogP(O. H-Bond donor(O. H-Bond acceptor(O. Violation Molecular weight(O500D. LogP(O. H-Bond donor(O. -Bond acceptor(O. Violation bioavailability are important as they predict the permeability and absorption of such drug across a biological membrane such as an epithelium cell, partition coefficient value . which is important in predicting intestinal absorption of such drug (Aucamp et al. , 2. The drug-likeness analysis of the sixteen leads showed that only five of the screened ligands and the standard drugs obeyed LipinskiAos rule of five with excellent drug-like properties (Table . Although, among the five ligands, chlorogenic acid and caffeoylquinic acid have one violation each. This indicates that the lead compounds have good Traditional Medicine Journal, 27. , 2022 oral bioavailability and permeability. So, they could be analyzed further as potential drugs. Bioactivity Analysis The bioactivity test is the base for therapeutic utilization and the potentially harmful effects of products prepared from plants (CoriaTyllez et al. , 2. The Ligand efficiency (LE). Fit Quality (FQ), and Ligand Efficiency-dependence lipophilicity (LELP) of the five . excellent lead compounds are presented in Table IV. Inhibition constant values of these compounds range from 46 to 1. 3 M, this showed that all of these Misbaudeen Abdul-Hammed Table IV. Bioactivity properties Ligands SD1 SD2 AutoDock Vina docking score cal/mo. (AAM) miLOG cal/mol/ heavy ato. E-Scale LELP L1= Coreximine. L2= Roseoside. L3= Chlorogenic. L4= Caffeoylquinic acid. L5 = Muricin M. SD1= Voriconazole. SD2= Fluconazole compounds are within the acceptable range of 0 to 0 M, except muricin M with a slightly higher value of 1. 3 M (Table IV). The heavy atoms in the structure of the ligands were used to calculate ligand efficiency (LE), with a recommended value Ou Fit Quality (FQ) was evaluated with a recommended value Ou 0. 8, and Ligand Efficiencydependence lipophilicity (LELP) was also determined with the expected recommended value in the range of -10 to 10 using equations 1 to Notably, all five ligands and the two standard drugs have excellent bioactivity properties in the right proportion as recommended, it is important to note that none of the ligands fails the bioactivity Ligand Efficiency (LE) = OeB. E y Heavy atoms (H. LE scale = 0. 873AOe0. AOe0. FQ = LE y LE scalea. LELP = LogP y LEa. AA. Binding Mode and Molecular Interactions of the Selected Hits and Standards The identification of the binding site and evaluation of the binding interaction of the ligands in the active pocket of the target receptor is important in the drug discovery process. This helps in ligand modification during the lead optimization stage of drug discovery (Falade et al. , 2. The molecular docking study predicted the residues at the interacting site of the proteins as well as their corresponding alignments (Adeoye et al. , 2. Table V and VI show the binding affinities and molecular interactions of selected lead compounds muricin M forms a conventional hydrogen bond with Met508. Tyr505. Ser507, and His377, pi-alkyl interaction with Phe380. Leu121. Pro230. Val234. Tyr132, and Leu88, unfavourable donor bond with Ser378 (Table V and VI). Upon examining the active site of the target receptor . TZ. , it was revealed that its binding pocket . ctive sit. was located in the heme loop and helix All the other residues of muricin N reported in Table VI were found at the active site of the target receptor . terol 14-demethylase . PDB ID 5TZ. except Leu88 and Val234 that was mentioned earlier in the structural and active site analysis of sterol 14-demethylase (PDB ID: 5TZ. This implies that muricin M shares the same pocket and interacts effectively with the active site of the target receptor, hence it could be a potent inhibitor of the receptor. Similarly, chlorogenic acid, roseoside, and caffeoylquinic acid showed strong interactions with 5TZ1 when compared with the standards (Table V and VI). Chlorogenic acid formed a conventional hydrogen bond with Pro230. Gly65. Tyr505. Met508. His377. Ser378, hydrogen carbon bond interactions with Ser507, pi-sulphur bond interactions with Met508, and unfavourable pibond interactions with His377 (Table V and VI). Roseoside formed a conventional hydrogen bond with Ser378. Met508. Pro230. Ser507. His377. Phe58, and pi-alkyl interactions with Tyr64. Phe380. Phe233. Leu87. Leu88. Val234, while Caffeoylquinic acid formed a conventional hydrogen bond with Ile304. His468. Arg381. Ile471. Gly308, and alkyl interactions with Lys143. Ile471. Ile131. Cys470, carbonhydrogen interactions with Gly464, and unfavourable donor interactions with Arg381. The excellent binding affinities of these ligands influence the molecular interactions (Tables V and VI). Therefore, chlorogenic acid, roseoside, and caffeoylquinic acid were found to have most of their residues at the active site of sterol 14-demethylase better than the two standard Traditional Medicine Journal, 27. , 2022 Antifungal Activities of Phytochemicals from Annona muricate Table V. Binding mode and molecular interactions of the selected hits against 5TZ1 Ligands Binding pocket Interaction L1-Muricin M L2-Chlorogenic Acid L3-Roseoside L4-Caffeolyquinic acid SD1 Voriconazole SD2- Fluconazole Traditional Medicine Journal, 27. , 2022 Misbaudeen Abdul-Hammed Table VI. Docking scores, binding sites and inhibition constants of the selected hit compounds and standards with sterol 14 -demethylase Ligands Muricin M Binding Affinity . G), kcal/mol 9A0. Chlorogenic Acid 2A0. Roseoside 5A0. Caffeoylquinic 1A0. Voriconazole 65A0. Fluconazole 3A0. CONCLUSION Receptor amino acids forming H-bond with Electrostatic/Hydrophobic ligands (H-bond Interactions Distance, yI) Met508 . 8 yI). Phe380. Ser378. Leu121. Met508 . 2 yI) Pro230. Val234. Tyr132. Leu88 Tyr505. 0 yI) Ser507. 1 yI) His377 . 9 yI) Pro230. 8 yI) Ser506. Met508. His377 Gly65. 9 yI) Pro230. 5 yI) Tyr505. 9 yI) Tyr505. 0 yI) Met508. 9 yI) His377. 7 yI) Ser378. 2 yI) Ser378. 8 yI) Ser378. 8 yI) Tyr64. Phe380. Phe233. Ser378. 2 yI) Leu87. Leu88. Val234 Met508. 0 yI) Pro230. 0 yI) Ser507. 0 yI) His377. 1 yI) Phe58. 2 yI) Ser378. 9 yI) Ala61. Tyr64. Ser378. His377. 4 yI) His377. Gly65 Tyr505. 1 yI) Met508. 8 yI) Met508. 8 yI) Ser378. 0 yI) Tyr132. 2yI) Val234. Ala61. Leu87. Arg469. 3 yI) Leu88. Trp57. Tyr505 Phe463. 3 yI) Tyr132. 4 yI) Phe463 Tyr132 Arg469 This research work evaluated one-hundred and thirty-one . ligands from Annona muricata . our so. against sterol 14demethylase using an in-silico approach . tructure-based drug desig. All these ligands were screened with the PyRx visual screening tool. ADMET SAR-2. Molinspiration web server, and many others. CASTp web tool was also used to validate their active site. The results obtained exceptionally favour muricin M (-7. 9 kcal/mo. , chlorogenic acid (-8. 2 kcal/mo. , roseoside (-8. kcal/mo. , and caffeoylquinic acid (-8. 1 kcal/mo. as probable inhibitors of sterol 14-demethylase Inhibition (K. AAM due to their outstanding binding energies. ADMET profile, drug-likeness. Bioactivity, favourable binding mode and molecular interactions with the target receptor . TZ. muricin M, chlorogenic acid, roseoside, and caffeoylquinic acid had shown favourable properties as potential inhibitors of 5TZ1 than voriconazole and fluconazole which have been in use for decades as first-line drugs for treating fungal infections (Hargrove et al. , 2. Since computational drug design is widely accepted in the world of modern drug design and development, the aforementioned ligands could be developed further and subjected to both preclinical studies and clinical trials toward the Traditional Medicine Journal, 27. , 2022 Antifungal Activities of Phytochemicals from Annona muricate development of a new anti-fungal drug. This research work has identified muricin M, chlorogenic acid, roseoside, and caffeoylquinic acid as potential inhibitors of sterol 14-demethylase (PDB: 5TZ. It is hereby recommended that structural modification of these ligands could be carried out via a lead optimization process to pharmacokinetics, and reduce their toxicity in human trials toward the development of new antifungal therapeutic agents. ACKNOWLEDGEMENT The authors acknowledge the members of the Computational and Biophysical Chemistry Research Group at the Department of Pure and Applied Chemistry. Ladoke Akintola University of Technology (LAUTECH). Ogbomoso. Oyo State. Nigeria. REFERENCES