key: cord-0822696-mrgdw07v authors: Allam, Ahmed E.; Assaf, Hamdy K.; Hassan, Heba Ali; Shimizu, Kuniyoshi; Elshaier, Yaseen A. M. M. title: An in silico perception for newly isolated flavonoids from peach fruit as privileged avenue for a countermeasure outbreak of COVID-19 date: 2020-08-13 journal: RSC advances DOI: 10.1039/d0ra05265e sha: 30ce4e0c0ca6a82cf6a15ee6ca4de16fda09999e doc_id: 822696 cord_uid: mrgdw07v 3′-Hydroxy-4′-methoxy-chroman-7-O-β-d-glucopyranoside 4 was first isolated from a natural source, together with three known compounds, the ferulic acid heptyl ester 1, naringenin 2, and 4,2′,4′-trihydroxy-6′-methoxychalcone-4′-O-β-d-glucopyranoside 3, which were isolated from peach [Prunus persica (L.) Batsch] fruits. These compounds were subjected to different virtual screening strategies in order to examine their activity to combat the COVID-19 outbreak. The study design composed of some major aspects: (a) docking with main protease (M(pro)), (b) docking with spike protein, (c) 3D shape similarity study (Rapid Overlay Chemical Similarity-ROCS) to the clinically used drugs in COVID-19 patients, and finally, (d) the rule of five and the estimated pre-ADMT properties of the separated flavonoids. Docking study with M(pro) of SARS-CoV-2 (PDB ID:6LU7, and 6Y2F) showed that compound 3, its aglycone part, and compound 4 have a strong binding mode to a protease receptor with key amino acids, especially Gln:166AA, and having a similar docking pose to co-crystalized ligands. Docking with the spike protein of SARS-CoV-2 illustrated that compounds 3 and 4 have a good binding affinity to PDB ID:6VSB through the formation of HBs with Asp:467A and Asn:422A. According to ROCS analysis, compounds 1, 3, and 4 displayed high similarities to drugs that prevent SARS-Co2 entry to the lung cells or block the inflammatory storm causing lung injury. Compounds 3 and 4 are good candidates for drug development especially because they showed predicted activity against SARS-CoV-2 through different mechanisms either by preventing genome replication or by blocking inflammatory storm that trigger lung injury. These compounds were isolated from peach fruit, and the study supports data and continues with the recommendation of peach fruits in controlling and managing COVID-19 cases. Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV-2) also known as 2019-nCoV is the causative microbe for the current pandemic COVID-19. 1 Coronaviruses (CoV) are a large family of viruses that lead to infection and sickness starting from common cold to severe diseases. 2 Middle East Respiratory Syndrome (MERS-CoV) 2012 and Severe Acute Syndrome (SARS-CoV) 2003 are well-known diseases that had arisen from this family. COVID-19 is the third recognized spillover of coronavirus to humans in the last two decades. The International Committee on Taxonomy of Viruses has considered the placement of the human pathogen, tentatively called 2019-nCoV, within the Coronaviridae family, 3 Fig. 1 . Although studies on the mechanism of replication as well as the pathogenesis of several coronaviruses have been very active, this family of coronaviruses received much attention because of the new human coronavirus that was responsible for severe acute respiratory syndrome (SARS), a contagious and fatal illness. 4 According to the World Health Organization (WHO), from December 2019 to now, about 10 million conrmed cases and about 0.5 million deaths have been reported. The spread of respiratory droplets through sneezing, coughing, or close contact between individuals is the main cause for the transmission of SARS-CoV-2. The symptoms of COVID-19 range from fever, headache, dry coughing, dyspnea, diarrhea, loss of taste, and fatigue to viral pneumonia, acute respiratory distress, and hypoxia. In severe cases, COVID-19 patients require mechanical ventilation. 5, 6 To date, there is no benecial therapeutics and limited effective treatment opportunities persist. targeting SARS-CoV-2 and these drugs, as shown in Fig. 2 , are prescribed under what is known as "Drug repurposing approaches or Drug repositioning approaches". 7 1.2. According to their molecular mechanisms, these drugs could be categorized as follows (1) Drugs pointing at SARS-CoV-2 directly, these drugs are either; (a) The inhibitors of crucial viral enzymes responsible for genome replication through inhibition of RNA polymerase or viral protease. (b) Hindering viral entrance to human cells through the inhibition of spike proteins. (2) Drugs as a human immune system modulator, these drugs are either; (a) Boosting the inherent response, which has a signicant role mainly in controlling viruses' replication. (b) Preventing the inammatory storm that causes lung injury. (3) Symptomatic regulator using VEGF inhibitors as Bevacizumab clinical trial (NCT04275414), and Intravenous Immunoglobulin. Here in Table 1 , the ongoing therapeutic choices are represented briey that may lead to combat this novel pandemic outbreak (Fig. 2) . Other therapies, involving mesenchymal stem cells, monoclonal antibody, methylprednisolone, recombinant interferon, This journal is © The Royal Society of Chemistry 2020 RSC Adv., 2020, 10, 29983-29998 | 29985 Table 1 Category of drugs recommended in the treatment of COVID-19 Category Notes Inhibiting the RNAdependent RNA polymerase It was originally developed by Gilead Sciences (USA) against the Ebola virus as a prodrug because its structure resembles endogenous adenosine 8 It has been previously shown to display antiviral activities against coronaviruses, especially SARS-CoV 9 Remdesivir is now being tested in different countries, two randomized phases, III trials 10, 11 Favipiravir It was developed by Toyama Chemical (Japan) for the treatment of inuenza It's a prodrug structurally resembling the endogenous guanine 12 Comparing to remdesivir, less preclinical studies have been established for favipiravir to kill SARS-CoV-2 Currently, it was approved by the National Medical Products Administration of China as the rst anti-COVID-19 drug in China, as the clinical trial had demonstrated efficacy with minimal side effects Ivermectin Inhibiting the viral Protease Anti-parasitic drug and was reported to inhibit both human immunodeciency virus (HIV) and dengue virus 13, 14 A recent in vitro study has proven its capability to reduce the viral RNA of SARS-CoV-2 (ref. 1) Currently, establishing a safety prole is the next step to verify Ivermectin's ability for curing patients diagnosed with COVID-19 (ref. 15 ) Lopinavir/ritonavir Both are used in a combination protocol for the treatment of HIV Although coronaviruses encode a different enzymatic class of protease, a number of clinical and in vitro model studies have been completed using this combination on SARS and MERS viruses. 16, 17 Currently, both drugs are prescribed in a clinical trial in early-stage COVID-19; however, no benet was observed beyond standard care 18 Hydroxychloroquine (HCQ) Blocking virus-cell membrane fusion A well-known antimalarial and anti-autoimmune agent with basic character It blocks the viral infection of SARS-CoV-2 by altering the endosomal pH required for the membrane fusion between the virus and the host cell 19 or by interfering in the glycosylation of its cellular receptor, ACE2 (ref. 20) It was approved in different countries to control the spread of COVID-19. 21 However, larger randomized controlled trials are required for further evaluation 22 On 17 June 2020, WHO announced that the use of HCQ in COVID-19 treatment was ceased Umifenovir It is a fully-functionalized indole with an antiviral activity against inuenza infection. It directs hemagglutinin glycoprotein on the surface of the inuenza virus and subsequently prevents its fusion with endosome aer endocytosis Currently, it is undergoing trials for COVID-19 (ref. It is a globular glycoprotein found in mammalian milk It acquires a wide spectrum of biological activities especially immunological properties In some countries e.g. Egypt, it has been prescribed in the protocol for the treatment of COVID-19 patients LF prevents coronavirus to get attached to Heparan Sulfate Proteoglycans (HSPGs) present at the surface of host cells as these viruses are considered to bind to the host cell by binding rst to HSPGs 24 Recombinant human Angiotensin-converting enzyme 2-APN01(rhACE2) Preventing the inammatory storm It blocks S protein from interacting with cellular ACE2 and so inhibits SARS-CoV-2 replication The administration of rhACE2 can decrease the serum level of angiotensin II and so there was no further activation of ACE2 receptor Currently, a pilot study is now evaluating the role of rhACE2 in COVID-19 Interleukin (IL)-6 inhibitors e.g. sarilumab, siltuximab, and tocilizumab Interleukin (IL)-6, IL-1, and TNF-a are the most important pro-inammatory cytokines in the human body Interleukin (IL) inhibitors may ameliorate severe damage to lung tissue caused by cytokine release in patients with serious SARS-CoV-2 infections. Several studies have indicated a "cytokine storm" with the release of IL-6, IL-1, IL-12, and IL-18, along with tumor necrosis factor-alpha (TNFa) and other inammatory mediators. The increased pulmonary inammatory response may result in increased alveolar-capillary gas exchange, making oxygenation difficult in patients with severe illness 25, 26 Fingolimod It is an immunomodulation drug, mostly used for treating multiple sclerosis 27 It is a highly potent functional antagonist of S1P1 receptors in lymph node T cells Currently, it is under clinical trials for the treatment of SARS-CoV-2, NCT04280588, MRCTA, and ECFAH of FMU [2020]027 Antibiotics as azithromycin Known macrolide antibiotics Reduce viremia to zero when used in combination with other drugs and natural killer cells have different functions to boost inherent response. Currently, corticosteroids have been approved in patients with COVID-19 with a hope of preventing lung brosis in patients with unresolved acute respiratory distress syndrome. Drug discovery approaches like virtual screening (VS), drug repositioning, quantitative structure-activity relationship (QSAR), and articial intelligence (also, called machine learning), are required more as the world is facing unwelcomed and uncontrolled scenario of current pandemic caused by SARS-CoV-2. 31, 32 Presently, virtual screening libraries, including ligand-based or target-based approaches, are a focal point for medicinal chemist to trigger time and money efficiency. 7 Because the discovery of new drugs requires a long time and expense, searching for new compounds from natural sources known for their high safety and applicability will be a good avenue to treat SARS-CoV-2. In order to quickly discovere lead compounds especially from food sources for clinical trials, a virtual screening study was initiated to identify new drugs targeting SARS-CoV-2. Among these food sources, peach [Prunus persica (L.) Batsch] is one of the nutritionally edible important fruits in the world. In terms of its biological activity, the peach fruit has been reported to have plenty of phenolic acids, avanones, and chalcones, which have a wide range of activity against viral infections. 33 The fruit is also involved in different biological activities, such as growth-inhibition activity against different breast cell lines, 34 antioxidant, anti-lipase, and anti-dementia activities. 35 Recently, Traditional Chinese Medicine (TCM) has recommended peach, under the Chinese name "Taoren", as a contributor in the prevention and treatment of COVID-19, 36 that evoked us to carry out this study. The total ethanolic extract of Prunus persica L. fruits afforded 3 0hydroxy-4 0 -methoxy-chroman-7-O-b-D-glucopyranoside 4, and it is the rst time it has been isolated from a natural source, together with three known compounds, ferulic acid heptyl ester 1, 37 naringenin 2, 38 and 4,2 0 ,4 0 -trihydroxy-6 0 -methoxychalcone- 35 as shown in Fig. 3 . The structures of the known compounds were identied by the comparison of their spectroscopic data with those reported in the literature. Compound 4 was obtained as a yellow amorphous powder, positive FAB-MS analysis showed [M + Na] + at m/z 471, corresponding to molecular weight 448 m/z and molecular formula The IR spectrum of compound 4 indicated the presence of hydroxyl (3400 cm À1 ) and carbonyl (1690 cm À1 ) groups. The UV spectrum of compound 4 exhibited a maximum at 286 and 330 (sh) nm and indicated that 4 is a avanone. 39 The bathochromic shis induced by AlCl 3 (70 nm) and NaOMe (74 nm) are typical of avanone with two hydroxyl groups at C-3 0 and C-4 0 . No UV bathochromic shi was observed with NaOAc suggesting the absence of free 7-hydroxyl group. 39 The 1 H-NMR spectrum of 4 reveals the presence of two sets of ABX-spin systems, one reminiscent to the avanone heterocyclic C ring at d H 5.2 (1H, dd, J ¼ 3, 12.6 Hz, H-2), d H 2.7 (1H, dd, J 41 Also, the 1 H-NMR spectrum showed an anomeric proton signal at d H 4.7 (1H, d, J ¼ 7.8 Hz, H-1 00 ), which indicated the b-conguration of the glucose moiety. Additionally, the 1 H-NMR spectrum exhibited a singlet at d H 3.80 (3H), indicating the presence of a methoxy group in the B ring. 42 The 13 C-NMR spectrum showed the presence of a quaternary carbon C-4 at d C 192.77, an oxymethine carbon C-2 at d C 80.00, an oxymethylene carbon C-3 at d C 46.43, and other carbons of two benzene rings, which indicated the presence of a avanone moiety. 43 Also, the 13 C-NMR spectrum showed an anomeric Approved for hospitalized cases Nitazoxanide 28 Orally active broad-spectrum antiparasitic Prodrug transformed rapidly to active metabolites, such as, tizoxanide and tizoxanide Known to potentiate interferon-alfa and interferon-beta assembly Previously revealed an in vitro activity against MERS-CoV and other coronaviruses Recommended in combination with azithromycin Thalidomide Thalidomide has been repurposed with different pharmacological effects It was reported as an anti-inammatory drug as it prevents the synthesis of TNF-a, 29 treating H1N1-infected mice by reducing pro-inammatory cytokines 30 Present studies are based on its immunomodulatory properties to treat COVID-19 (NCT04273529, NCT04273581), phase 2 This journal is © The Royal Society of Chemistry 2020 RSC Adv., 2020, 10, 29983-29998 | 29987 carbon at d C 105.12, which conrmed the presence of glucose moiety. 44 The HMBC spectrum showed a signicant correlation between the anomeric proton (H-1 00 ) at d H 4.7 of glucose moiety and the carbon (C-7) at d C 162.52. 43, 45 The elucidation of the absolute conguration at C-2 was based on the values of coupling constants with methylenic protons H-3 a,b (Jax-ax ¼ 12.6 and Jax-eq ¼ 3 Hz). The close similarity of H-2 chemical shis, as can be seen in Table 2 , with those of the literature, 46,47 thus conrmed the S-conguration of C-2. From the above-mentioned spectroscopic data, compound 4 was identied as (À)-(2S)-3 0 -hydroxy-4 0 -methoxy-chroman-7-O-b-D-glucopyranoside, and it is the rst time it has been isolated from a natural source. 2.2.1. Molecular docking study. SARS-CoV-2 main protease (M pro ) and spike glycoprotein (S) are crucial elements in the infectious route of the virus, and they have been recognized as important targets for therapeutic strategies. 48 SARS-CoV-2 is a single-stranded positive-sense RNA genome. The virus envelop contains spike protein (S) which regulates viral pass into the host cells. Then, two polyproteins i.e. pp1a and pp1ab are promptly translated upon entry into the host cells, and nally, these are separated by two viral protease enzymes to make viral replication. Thus, these molecular targets can be is considered as druggable targets. 49 Numerous molecular dynamic simulation studies were addressed to gure out the active sites of the SARS-CoV-2 main protease (M pro ) or spike protein 50 respectively, Fig. 4 . Furthermore, a docking study with the spike protein of SARS-CoV-2 (PDB ID:6vsb 53 ) was represented. By this tactic, we will indicate the possible mechanism of these naturally occurring avonoids to ght COVID-19. The standard ligands N3 and a-ketoamide in addition to clinically used protease inhibitor drugs, such as remdisever, lopinavir, and ritonavir were docked with the available co-crystallized main protease (M pro ). The experimental docking of these standards presented similarity in their binding mode inside the receptor especially HBs with the key peptides and The literature survey displayed that avonoids were previously reported to inhibit the main protease of SARS-CoV-1. 4, 54 Recently, there are studies for the activity of phenolic compounds and avonoids against SARS-CoV-2. 26 To well understand the expected activity of the isolated avonoids from peach fruits, docking analyses were performed using the experimental crystal structure of SARS-CoV-2, which was determined recently. Consensus score for the isolated compounds and some selected antiviral drugs are illustrated in Table 3 . Among the isolated compounds and based on the consensus score value, the avanone 4, chalcone 3, and its aglycone bind strongly with M pro receptor. To validate our study, the reported ligand was re-docked with the active pocket. N3 showed interactions with receptor having HB formations with Gln 192: A and Glu 166: A, Fig. 5a . Compound 2, the (S) isomer, displayed HBs by its carbonyl function which is arranged close to the Glu 166: A. The hydroxyl group at C-8 participated in HB interaction with Ser 144: A. The hydroxyl group attached to C2-phenyl moiety at P-position occupied the amino acid cles through hydrophobic-hydrophobic interactions, Fig. 5b . In order to ensure the effect of compound's absolute conguration on docking mode (enantioselective docking relationship) and subsequently explore their structure activity relationship (SAR), the (R) isomer demonstrated different binding pose and mode with the receptor. Docking results clarify that both the enantiomers of compound 2 (Naringenin) displayed dissimilarity in their binding pose inside the active site. Taxifolin is a dihydro quercetin avonoid with an interesting spectrum of biological activities. 55 Recently, in a virtual screening study, taxifolin showed the lowest binding free energy against M pro of SARS-CoV-2. 56 Our attention was directed to examine the effect of different pharmacophoric features of taxifolin and compound 2 that will support the structureactivity relationship and drug design strategy. Taxifolin presented a binding with the receptor by hydroxyl functionality at C-3, which is arranged close to Glu 166: A with the formation of strong HB. Whereas the phenyl moiety and chroman ring were well-inserted through hydrophobic-hydrophobic interaction in Regarding the ferulic acid heptyl ester 1, as indicated from its structure, it formed hydrophobic-hydrophobic interactions without the formation of HB, Fig. 5c . This result emphasizes the importance of physicochemical parameters, such as log P, which will be elucidated later. Likewise, compound 3 disclosed strong hydrophobic-hydrophobic interactions, Fig. 5d . To perceive the effect of glycosidation, the 3D shape of compound 3 inside the receptor represented dissimilarity with its aglycone part. The aglycone moiety disclosed strong bonding interaction (consensus score 80). It formed HBs with Glu 166: A, Hist 163: A, and Ser 144: A through the carbonyl at C-2 and OH at C-6 (glycoside linkage), respectively, Fig. 5e . Concerning compound 4, it facilitated HB formation with Ser. 144: A through the OH group of phenyl at C-2, and the sugar moiety also participated in HB formation with Gln 192: A, Fig. 5f Fig. 7a . The Umifernovir drug exhibited HB with Asn:422A (consensus score 174), ESI. † Among the separated avonoids, compounds 3 and 4 were observed to have a good binding affinity to PDB ID:6vsb. Compound 3 (consensus score 464) showed HB interaction with Asp:467A and Asn:422A, Fig. 7b ; while the corresponding aglycone illustrated a docking outside the receptor (consensus score 626), ESI. † This feature indicates that glycosylation is important to provide the compound some polarity essential for binding interactions. Compound 4 formed HBs with Ilu:468A and hydrophobic interactions (consensus score 396), Fig. 7c . Although compound 2 formed HBs with Ser:494A, Val:350A (consensus score 606) but it interacts with the S protein outside the created receptor, Fig. 7d . However, the analog R isomer appeared without docking and was visualized outside the receptor (ESI †), and ferulic ester (1) occupied the receptor with hydrophobic interactions only (consensus score 880), Fig. 7e . Despite all the prospective approaches of molecular docking, the chemistry of ligands (tautomerism and ionization), the exibility of the receptor, the availability of receptor crystal structure, and scoring function persisted the challenge. In this regard, although molecular docking is an attractive technique to recognize drug biomolecular interactions for drug design and discovery, many challenges for docking study could occur. 2.2.2. Shape similarity and lead hopping TC scores. Rapid Overlay Chemical Structure (ROCS) is a virtual relative screening tool used to observe the resemblance between chemical entities based on their three-dimensional shapes. [57] [58] [59] High match in shape reects a high match in biology while the high match in biology does not reect the relationship between biological activity and the 3D shape structure. Different applications of ROCS can be applied, such as virtual screening, lead hopping, molecular alignment, pose generation, and structural predictions. 34, 35, 37 The approach of this study is to identify molecules that can adopt shapes extraordinarily similar to reported drugs recommended in COVID-19 disease. ROCS was applied to identify the similarity of our compounds with different drugs or cover other proteins related to SARS-CoV-2 not co-crystallized yet, Table 4 . ROCS application involves two les, 38 which must be in most stable conformer generated by omega algorism: (a) database le, which is a compound collection for this study; (b) query le here is the standard drugs. In this work, we decided to examine our separated compounds with different queries in order to gure out the most possible mechanism in targeting SARS-CoV-2, Table 5 . The study depicts a match between the database le and query drug in a separate run. This match is based on the volume overlap of optimally aligned molecules, which are virtually independent of the atom types and bonding patterns of the query. Two main outputs generated from ROCS analysis: (1) shape similarity which includes shape counter, shape atoms, and color atom labels for database set (query compounds) and was visualized by vROCS and VIDA applications. We added this feature which has been written in Fig. 8a-e; and (2) a set of scores expressed in Tanimoto scores. The most imperative score is Tanimoto Combo (TC) that contains both shape ts and color. This has a value between 0 and 2 and the score is used for ranking the hit list, Table 4 . Among the separated compounds, compound 4 has the highest similarity to remdesivir, lopinavir, and hydroxychloroquine (HCQ) with TC ¼ 0.71, 0.58, and 0.76 respectively. Table 4 signies that this compound has high shape similarity to these drugs that prevent the genome of the virus to replicate. Fig. 8a and b illustrated color, volume and shape of remdesivir and lopinavir respectively. rand their overlay with compound 4 with remdesivir. However, the two benzyl moieties of lopinavir occupied the outside volume of compound 4. Moreover, chalcone 3 has high similarity to protease inhibitor a-ketomide (TC ¼ 0.64). The aglycone part of compound 3 (considered as a metabolite, Fig. 8c ) exhibited high similarity to Hydroxychloroquine HCQ (TC 0.85). Additionally, it has also a high similarity to umifenovir (TC 0.87). These results emphasize the hypothesis that this compound (3) has the ability for blocking the entry of virus. Ferulic ester 1 has the highest similarity to favipiravr, umifenover, and ngolamide, demonstrating that this compound has a high likeliness towards different mechanisms for SARS-CoV-2 treatment, Fig. 8d and e. In response to their high privilege character but with limited expected drug-likeness properties, these combinations of characters encouraged us to add predicted physicochemical parameter values as will be discussed later to study the drug ability and low toxicity of compounds. 2.2.3. Predicted pharmacokinetic and pharmacodynamic parameters. 60 The drug-likeness score of oral drugs is determined based on their physicochemical properties. The drug candidates with high drug-likeness scores illustrated higher absorption and bioavailability with lower doses and have fewer drug-drug interaction warnings. 61 Absorption, distribution, metabolism, excretion, and toxicity (ADMET) calculations for drug candidates are required for drug development. These values contribute to determining the failure of approximately 60% of all drugs in clinical trial phases. In this regard, ADMET is determined at the beginning of drug discovery phases to eliminate molecules with poor ADMET properties from an earlier drug discovery pipeline with the aim to save research costs. In this respect, chemoinformatic tools were used to predict ADMET parameters. As illustrated in Table 5 , compounds 3 and 4 were considered as isomers as both of them have the same molecular weight and are similar in most ADMET and physicochemical parameters (drug-likeness score ¼ 0.56 and 0.34, respectively). The aglycone part of compounds 3 and 4 showed higher acceptable physicochemical and Pr-ADMT parameters. According to their 2D structure (compounds 3 and 4), these compounds originated from two different architecture systems. Compound 1 has a higher lipophilicity value than compound 2, then compound 3 and nally compound 4. Their lipophilicity scores were represented by log P values and BBB. Compound 1 exhibited high log P value of 5.08 and therefore low hydrophilicity causes poor absorption or permeation with a low drug-likeness score of À0.51, highlighting that lipophilicity is an important parameter in the drug development process. However, compound 2, aglycones of 3, and compound 4 well-exhibited good log P values (2.38-2.81), indicating that the hydrolysis of the glycoside linkage is valuable for compound activity. Silica gel and RP-18 silica gel were purchased from Wako (Osaka, Japan). Thin-layer chromatography (TLC) silica gel G 60 F 254 was purchased from Merck (Darmstadt, Germany). All other ingredients used were of the highest grade available. The fresh materials were freeze-dried and milled below 1 mm. Then, their extracts were prepared by shaking (200 rpm) with solvent (ethanol) at room temperature for 48 h. The extracted solution was ltered and evaporated using ADVANTEC no. 2 lter paper (Toyo Roshi Kaisha, Ltd, Tokyo). The ethanol extract was suspended in distilled water and divided between n-hexane, ethyl acetate, and the remaining water to give the n-hexane fraction (12 g), ethyl acetate fraction (5.1 g), and the remaining aqueous fraction (45.5 g). We attempted to isolate the active compounds from the ethanol extract derived from peach fruits. The ethyl acetate fraction from the ethanol extract was sub-fractionated on a silica gel column using chloroform-methanol gradient elution (25%, 50%, 75%, and 100%; 2 L each). The fraction eluted by 25% methanol (1.4 g) was further separated by chromatography on an ODS column (80 Â 200 mm; Cosmosil 140 C-18 PREP, Nacalai Tesque, Tokyo, Japan) using six mobile phase systems of methanol-water (10%, 25%, 40%, 50%, 70%, and 90% v/v; elution volume: 1.5 L of each) to give six corresponding fractions. The fraction eluted with 10% methanol (310 mg) was further chromatographed by column chromatography on silica gel and eluted by a stepwise gradient of chloroform-methanol This journal is © The Royal Society of Chemistry 2020 RSC Adv., 2020, 10, 29983-29998 | 29995 (ratios of 9 : 1, 85 : 15, 8 : 2, 7 : 3, and 6 : 4; v/v elution volume: 200 ml each) to give ve corresponding fractions. The fraction eluted with 85 : 15 chloroform-methanol resulted in the elution of compounds 1 and 2 (55 and 40 mg, respectively). The fraction eluted with 8 : 2 resulted in the elution of compound 3 (35 mg) and that eluted with 7 : 3 chloroform-methanol resulted in the elution of compound 4 (28 mg). The structures of all compounds were elucidated by 1D and 2D NMR spectroscopy, including 1H, 13C, HSQC, and HMBC experiments, Fig. 3 . The 1 H, 13 C-NMR, and 2D spectra of the isolated compounds were recorded using a Bruker DRX 600 NMR spectrometer (Bruker Daltonics, Billerica, MA). The molecular docking and shape similarity known as ROCS studies were carried out using the OpenEye Modeling soware [Fast Rigid Exhaustive Docking (FRED) Receptor, version 2.2.5; OpenEye Scientic Soware, SantaFe, NM (USA); http:// www.eyesopen.com]. A virtual library of target compounds was used and their energies were minimized using the MMFF94 force eld, followed by the generation of multi-conformers using the OMEGA application. Data were visualized by vida commands. Lipinski's rule (rule of ve) and molecular property prediction was calculated by the free accesses to website https:// www.molso.com/servers.html. PreADMET estimation was determined by utilizing the free access of website https://preadmet.bmdrc.kr/. The isolated compounds from peach fruits showed a high resemblance to drugs that prevent the multiplication of SARS-CoV-2. Although compound 1 is very similar to some of these drugs, such as ngolamide, however, it showed less number of physicochemical parameters and low drug-likeness scores. The aglycones of compounds 3 and 4 illustrated strong binding interactions with M pro and spike protein of SARS-CoV-2. Compounds originated from avanone core systems of 2, 4, and those from a,b-unsaturated acetophenone core architecture 3 have high similarity in their 3D structures and their predicted biological activities. This study pointed out the importance of taking peach fruits as a source of naturally occurring compounds isolated from food sources in alleviating SARS-CoV-2, and hence the advice of eating peach fruits to COVID-19 patients. Additionally, it will help us to afford semisynthetic new scaffolds with potential activities against this disease. Our outcomes prove the value of this screening study, which can lead to the rapid discovery of drugs from natural food sources and prioritize them as potential multitargets acting on pandemic COVID-19. We declare that we have no conicts of interest. Techniques of avonoid identication