key: cord-0700011-akckl4u1
authors: Prajapati, Jignesh; Patel, Rohit; Rao, Priyashi; Saraf, Meenu; Rawal, Rakesh; Goswami, Dweipayan
title: Perceiving SARS-CoV-2 Mpro and PLpro dual inhibitors from pool of recognized antiviral compounds of endophytic microbes: an in silico simulation study
date: 2022-04-08
journal: Struct Chem
DOI: 10.1007/s11224-022-01932-0
sha: 5ff8b02c2ce7d415c8f03f993dda16cb777422fc
doc_id: 700011
cord_uid: akckl4u1

Coronavirus disease 2019 (COVID-19) persists and shook the global population where the endgame to this pandemic is brought on by developing vaccines in record-breaking time. Nevertheless, these vaccines are far from perfect where their efficiency ranges from 65 to 90%; therefore, vaccines are not the one only solution to overcome this situation, and apart from administration of vaccines, the scientific community is at quest for finding alternative solutions to incumber SARS-CoV-2 infection. In this study, our research group is keen on identifying a bioactive molecule that is independent in its mode of action from existing vaccines which can potentially target the SARS-CoV-2 virus replicative efficacy. Papain-like protease (PLpro) and main protease (Mpro) are the most lucrative targets of COVIDs against which the drugs can be developed, as these proteases play a vital role in the replication and development of viral particles. Researchers have modelled a compound such as GRL0617 and X77 as an inhibitor of Mpro and PLpro, respectively, but use of these compounds has several limitations on hosts like toxicity and solubility. Under the current study by deploying rigorous computational assessments, pool of microbial secondary metabolites was screened and handpicked to search a structural or functional analogue of GRL0617 and X77, with an idea to identify a compound that can serve as dual inhibitor for both PLpro and Mpro. From the manually curated database of known antiviral compounds from fungal origin, we found cytonic acids A and B to potentially serve as dual inhibitor of PLpro and Mpro.

Coronavirus disease 2019 (COVID-19) manifests a range of physiological conditions on host, ranging from being entirely asymptomatic to minor influenza to cough, cold and sore throat to manifestations like acute respiratory distress, pneumonia and lung damage to even death [1, 2] . During the first wave of the infection, all affected countries endorsed various precautionary measures, including a total or partial lockdown of residents, implementation of social distancing, wearing of face masks and frequent hand sanitization. More significantly, individuals did not follow the prevention measures precisely and the number of affected individuals changed relying upon how powerful the prevention measures were. Additionally, these recommended measures have been headed to the detriment of the economy and more crucially, these severe control measures alone are not viable enough to stop the COVID-19 pandemic. This is vividly observed in both developing and developed countries, where people were hit by a third wave of coronavirus infection due to the lack of antiviral medication or vaccination, and as a result, the scientists from around the world are prompted to find effective medications and vaccines for COVID-19 treatment and prevention [1, 3] .

In the last 1 year, COVID-19 vaccines are being developed using various strategies, executed from a typical inactivated and live-attenuated vaccines to more innovative messenger RNA (mRNA) and DNA technologies, like subunit vaccines [4] [5] [6] . For instance, BBV152 is a classical inactivated whole virion vaccine, while mRNA-1273 and mRNA-BNT162b2 are mRNA-based vaccines that target SARS-CoV-2 spike protein [5] . Herd immunity may hypothetically be achieved because of the high incidence of protection from infection after vaccination; however, there is no reliable assurance to it. Besides, there are yet various vaccination-related inquiries that are presently unaddressed, which include the length of a vaccination preventive effects and cross-protection against variants of interests (VOI) and variants of concerns (VOC) (Alpha, Beta, Gamma, Delta, Delta Plus, Epsilon, Eta, Theta, Iota, Kappa and Lambda), particularly the new VOC, Omicron (B.1.1.529) [7] . Mortality related with COVID-19 is fairly reducing as vaccination rates rise, but the list of variants is increasing as new mutations continue to evolve, challenging the whole vaccine regime [8] .

While the vaccination drive continues, another promising strategy would be to foster new drugs that are compelling against SARS-CoV-2. Virus replication necessitates a certain set of functional and structural proteins. Malfunctioning or inhibition of these proteins can slow or stop infection or can also halt viral multiplication, culminating its spread, which is desired by researchers. As a result, medicines that can regulate the accessibility of these proteins should be pursued. SARS-CoV-2 is a 29-kb RNA genome that encodes two large overlapping polyprotein precursors (pp1a and pp1ab); these two translated stretches of inactive proteins are then treated upon by PLpro and Mpro which cleave pp1a and pp1b into various small fragments where each fragment serves as active viral protein that functions in replication, transcription and assembly of the virus nuclear material into the capsid protein coat (Fig. 1) . Consequently, inhibiting these proteases may have a considerable impact on the viral machinery and the overall rate of viral infection. Scientists around the globe have well considered these proteases (Mpro and PLpro) as important drug targets for developing new Fig. 1 Representation for addressing role of viral proteases PLpro and Mpro to produce non-structural proteins participating in viral replication and transcription medications or repurposing already existing drugs [9, 10] , as inhibiting these proteases directly halts the functioning and replication of SARS-CoV-2, making the virus ineffective [11, 12] .

Bioactive and natural products derived from plants, microorganisms and animals are most often recommended to treat infectious disorders [13, 14] . The term 'endophytes' signifies microorganisms living inside the tissue of plant with imparting mutualistic behaviour. Considering the published literature, endophytes could be an untapped reservoir full of diverse resources of bioactive compounds [15] . Identifying a compound from reservoir that can block PLpro and Mpro is the rationale of this research which in itself is a challenging as well as exciting avenue.

The likelihood of discovering novel inhibitors using microbial databases of natural compounds has increased as computer hardware and software technologies have progressed. Several investigations have used computer-aided drug designing (CADD) tools where in its pipeline first the compounds are screened by 'molecular docking' approach which identifies compounds that can interact with important amino acids in the active site of protein, followed by 'molecular dynamics (MD) simulations' that perform a real-time simulation of interaction between screened ligand and protein under virtual physiological environment, and a blend of advanced in silico approaches to discover natural bioactive compounds from large databases that could be used as possible lead compounds [16] [17] [18] [19] . The novelty of this piece of research work lies in the fact that antiviral metabolites from endophytic origin are very less investigated against the SARS-CoV-2 proteases (Mpro and PLpro), the reason being the scarcity of published literature and unavailability of off-the-shelf databases. To date, there are several secondary bioactive active metabolites identified from endophytes possessing ability to inhibit viral proteases, and they are emodin, ω-hydroxyemodin, ( +) -sclerotiorin, phomopsone B and phomopsone C [20] [21] [22] [23] , but such study for proteases of SARS-CoV-2 is lacking which serves as a research loophole which we are targeting in this research piece.

For this study, we have nit-picked literary evidence of antiviral compounds obtained from endophytic microorganisms and prepared a library to virtually screen them against SARS-CoV-2 Mpro and PLpro, as a dual inhibitor. This is done to ensure the optimal potential of the ligand library to act simultaneously on the activity of both proteases altogether. After screening the entire library, cytonic acids A and B were identified as PLpro and Mpro dual protease inhibitor. Under this research piece, screening hierarchy followed the use of protein-ligand docking, followed by pharmacophore hypothesis generation, which was then followed by absorption-distribution-metabolism-excretion-toxicity (ADMET) profiling and protein-ligand binding Gibbs free energy change calculation through molecular mechanics generalised Born surface area (MM-GBSA) assessment. Additionally, the dynamic nature of protein-ligand contacts was analysed using MD simulations. The simulation was undertaken for cytonic acid derivatives and control compounds for a 100-ns-long timescale. We predict that cytonic acids A and B as a dual-target therapeutic approach for COVID-19 could be more effective and would substantially reduce the use of combinational drugs as well as prevent the multi-drug dose load on the already immune-compromised and weakened host system.

PLpro protein of SARS-CoV-2 (PDB ID: 7CMD) co-crystalized form with GRL0617 and Mpro of SARS-CoV-2 (PDB ID: 6W63) co-crystalized form with X77 were retrieved from the protein databank (PDB) [24] [25] [26] . Here, GRL0617 and X77 are the known inhibitors for respective proteins and were used as reference. Proteins were imported to Schrödinger Maestro and were prepared in protein preparation wizard of Schrödinger Maestro [27] . In which, the addition of hydrogens was performed, and het state at pH 7.0 was created using Epik. Then, the protein was optimized with PROPKA keeping its parameter as pH = 7.0. On successful completion of these two steps, the minimization was performed using force field OPLS2005.

From the available literature of research articles, the data of endophytic microorganisms and endophytic microorganisms and their bioactive metabolites was manually collected. Based on experimental antiviral effects of metabolites against various viruses such as influenza, H1N1, HIV, HCV, and EV-71, the 45 antiviral compounds were retrieved from published studies. Table 1 lists the 45 compounds present in the library, along with their isolation sources and bioactivity. All the ligands were retrieved in SDF format and were imported to the Schrödinger Maestro for ligand preparation for docking, where the LipPrep wizard was used, addition of hydrogens was performed, and ligand minimization was performed using force field OPLS2005. The files were exported after ligand minimization and subsequently used for molecular docking. Generation of E-pharmacophore model E-pharmacophore hypothesis was generated separately using the crystal structure of (i) SARS-CoV-2 PLpro bound to an inhibitor GRL0617 with a resolution of 2.1 Å (PDB ID: 7CMD) and (ii) SARS-CoV-2 Mpro bound to an inhibitor X77 with a resolution of 2.1 Å (PDB ID: 6W63). During preprocessing of the complex, water molecules within 5 Å distance of the ligand were eliminated from the structure and the minimised complex was used to develop E-pharmacophore model [28, 29] . The E-pharmacophore model was developed using the Phase module (Phase, Schrödinger, LLC, New York, NY, 2021) [30] . For developing hypothesis, default pharmacophore properties like hydrogen bond acceptor (A), hydrogen bond donor (D), aromatic ring (R) and hydrophobicity (H) were mapped. On performing E-pharmacophore, two hypotheses were generated, one for each protein-ligand complex. These E-pharmacophore hypotheses were then used for screening of ligands.

A E-pharmacophore-based virtual screening was carried out on the prepared library of 45 compounds using the Phase module of the Schrodinger suite to develop a collection of compounds having requisite ligand characteristics aimed at best binding to PLpro and Mpro. This was achieved by mapping E-pharmacophore models on the library of 45 compounds.

The prepared proteins as shown in earlier step were used for molecular docking. 

MM-GBSA determines the Gibbs free energy change in kcal/mol that predicts how spontaneous the interaction reaction can be between selected screened ligand and protein.

The greater negativity in the values shows more spontaneous reaction. This ΔG value can be compared for with that of reference native ligand. If the values shown by screened selected ligand are more negative than by the reference ligand, then the binding spontaneity of the screened ligand is stronger than of the reference native ligand [31] [32] [33] [34] . MM-GBSA assessment was performed using the Prime module of Schrodinger Maestro making use of OPLS2005 force field. The ΔGBind is the overall Gibbs free energy change value of the system, which is also bifurcated further into segments as ΔGCoulomb-Coulomb energy, ΔGHbondhydrogen-bonding correction, ΔGLipo-lipophilic energy, ΔGPacking-pi-pi packing correction, ΔGvdW-Van der Waals energy, etc.

The ADMET properties of the GRL0617 and X77 (as controls) and top screened ligands (test endophytic molecules) NA not available were predicted using the pkCSM-pharmacokinetics server [35] which provides the theoretical values that are obtained by proven animal model turned in to finely tuned algorithms. Values for each of the parameters like water solubility, CaCO 2 permeability, intestinal absorption (human), skin permeability, P-glycoprotein substrate, P-glycoprotein I and II inhibitor, blood-brain barrier (BBB) permeability, central nervous system (CNS) permeability, cytochrome p450 (CYP)2D6 and CYP3A4 substrate, CYP1A2-CYP2C19-CYP2C9-CYP2D6-CYP3A4 inhibitor, renal clearance, AMES toxicity, max. tolerated dose, hERG I and II inhibitor, oral rat acute and chronic toxicity, hepatotoxicity, skin sensitisation, T. pyriformis toxicity and minnow toxicity are obtained for each ligand under study by this assessment. This helps to understand and compare drugs, like properties of screened selected ligands in comparison with the reference molecules.

The Desmond module of Schrödinger Maestro was used to perform simulations of SARS-CoV-2 PLpro-ligand and Mproligand complexes. Cytonic acid A and cytonic acid B were having efficient binding energy and proper ADMET profile and both were able to mimic interactions of native ligand with PLpro as well as Mpro; MD simulations were carried out for the following complexes: (I) PLpro-GRL0617, (II) PLprocytonic acid A, (III) PLpro-cytonic acid B, (IV) Mpro-X77, (V) Mpro-cytonic acid A and (VI) Mpro-cytonic acid B. For the simulation assessments, the complexes PLpro-GRL0617 and Mpro-X77 were taken as reference controls against which the simulation results of cytonic acid A and cytonic acid B with both the proteins were compared with. For MD simulations, the virtual simulation box was prepared into which complex was placed, and then, the simulation box was filled with water molecules with TIP3P model, after which the ion Cl − or Na + was placed in the simulation box for neutralisation. MD simulations were performed with NPT ensemble at 300 K and 1.01 bar for 100 ns. Post MD simulation assessments were performed by assessing ligand receptor interaction profiles with respect to time, types and nature; the root mean square deviation (RMSD) and root mean square fluctuation (RMSF) assessments of trajectories were also analysed to evaluate the stability of protein-ligand complexes during the simulations run.

The prepared library for this in silico analysis included specialized nit-picked literary evidenced endophytic compounds having in vitro anti-viral activity observed against influenza, H1N1, HIV, HCV and EV-71 viruses. Forty-five antiviral compounds obtained from endophytic microorganisms were used to virtually screen them against SARS-CoV-2 Mpro and PLpro, as a dual inhibitor. Table 1 contains information on the names of the compounds, their endophyte sources and their antiviral efficacy against particular viruses.

Receptor-pharmacophore based on the 3D structure of a target protein can provide detailed and accurate information on ligand interaction attributes. In totality, there are about 8 descriptors used for pharmacophore modelling; amongst which, hydrogen bond donors, acceptors, lipophilic regions and negative and positive ionizable groups along with aromatic rings are the most commonly used descriptors in pharmacophore modelling. For this study, 3D structure-based e-pharmacophores were prepared with the 'Phase' module which utilises the receptor-ligand pharmacophore generation procedure and discovered amino acids involved in ligand binding for a co-crystal GRL0617 inhibitor within the binding site of PLpro (Fig. 2a) and for a co-crystal X77 inhibitor within the binding site of Mpro (Fig. 2c) . Hydrogen bond acceptor, hydrogen bond donor and pi-pi stacking of aromatic ring were amongst the five key 3D attributes of the generated e-pharmacophore for the evaluated protein. Figure 2b depicts the five 3D pharmacophore characteristics for the PLpro-GRL0617 complex that includes two donor hydrogen bonds, one acceptor hydrogen bond and one aromatic ring sphere, while Fig. 2d depicts the five 3D pharmacophore characteristics for the Mpro-X77 complex, including two acceptor hydrogen bonds and three aromatic ring spheres.

The 45 potential anti-viral compounds of endophytes origin (Table 1) were screened using the E-pharmacophore hypothesis created from the GRL0617-PLpro complex and X77-Mpro. The 'Phase' screen score and corresponding binding site variables were used to filter these compounds. The Phase screen score is utilised to score and rank the matching compounds. This scoring function assesses both the quantity (if partial matching is permitted) and quality of ligand feature matching, where quality is defined by site, vector and volume scoring components. After the execution of the proposed pharmacophore hypothesis, a total of 15 molecules managed to pass this filter for the PLpro-GRL0167 complex (Table 2 ) and 18 molecules for the Mpro-X77 complex ( Table 3) .

After using Glide's Receptor Grid Generation tool in Maestro to define the grid box, the minimised 3D molecular structure of the co-crystallised inhibitor GRL0617 was docked into the active site of viral PLpro. The molecular interaction of the co-crystallised inhibitor (GRL0617) with the binding site of the viral protease SARS-CoV-2-PLpro was discovered. In general, hydrophobic interactions induce GRL0617's binding to SARS-PLpro, which results in protein inhibition. The aromatic rings' 1-naphthyl moiety frames the primary contact as a pi-pi interaction with Tyr264 and Tyr268, is less exposed to solvent and fits into the cavity at the locus where the leucine at the P4 position is required. Furthermore, the 1-naphthyl moiety interacts with Pro247's and Pro248's side chains. The (R)methyl group at GRL0167's stereocenter is located inside the protein, between Tyr264 and Thr301 residues, in a small polar region. GRL0617 features a single aromatic ring with NH2 at R3 position, which is positioned at the cavity's aperture, in addition to the 1-naphthyl moiety. The cavity is polar because of the presence of several polar groups, such as Gln269's side chain oxygen and Tyr268's hydroxyl group, which function as hydrogen bond acceptors in the interaction. After pharmacophore screening, a total of 15 microbial metabolites were investigated for molecular docking utilising Schrödinger's Glide XP lead optimization methodology. Four fungal metabolites of endophytic origin, aspergilline E, cytonic acid A, cytonic acid B and 4-dehydroxyaltersolanol A, produced significant binding energies along with multiple conserved interactions with Tyr268 of SARS-CoV-2-PLpro. Each of these molecules produced the binding energies of − 7.7, − 7.655, − 7.292 and − 6.859 kcal/mol, respectively ( Table 4 ). The interactions between these antiviral compounds and the protein (PLpro) are depicted in Fig. 3 , and the structural features of GRL0617 and the top four antiviral metabolites are listed in Table 5 . Consequently, after successively establishing the grid box with Glide's Receptor Grid Generation tool in Maestro, the minimised 3D molecular structure of the co-crystallised inhibitor X77 was docked into the active site of viral Mpro. The molecular interaction of the co-crystallised inhibitor (X77) with the binding site of the viral protease SARS-CoV-2-Mpro was observed. The co-crystallised binding site appears to be engaged in a number of hydrogen bonding, Pi stacking and hydrophobic interactions. The non-covalent inhibitor X77 forms a hydrogen bond with Asn142, Gly143, Glu166 and His41. Furthermore, the ligand is sandwiched between the benzene ring of Phe140 and the imidazole ring of His172, forming a pi-pi connection with Cys145 and His41. Also hydrophobic interactions with different amino acids were observed including Thr25, Leu27, Ser144, Phe140, Leu141, His172, Leu167, Tyr54, His164, Gln189, Met165 and Pro168 within the active site.

The top four molecules from endophytes that showed to effectively bind with Mpro are cytonic acid B, cytonic acid A, asperphenalenone D and asperphenalenone A, respectively, with the binding energy of − 9.998 kcal/ mol, − 9.798 kcal/mol, − 8.969 kcal/mol and − 6.2417 kcal/ mol ( Table 4 ). The compound with the best docking score, cytonic acid B, makes six hydrogen bonds with the amino acids, namely, Thr24, Thr26, Asn142, Glu166, Arg188 and Thr190; along with it forms Pi-anion interaction with Cys145. Moreover, hydrophobic interactions in the form of alkyl/pi-alkyl interaction with His41, Met49 and Met165 are being formed. Apart from these, Van der Waals interaction with this ligand is formed by 13 different amino acids as represented in Fig. 4 . Cytonic acid A tends to form hydrogen bonds with Thr24, Asn142, Glu166 and Arg188, while alkyl/pi-alkyl interaction involves residues Cys44, Thr26, Met49 and Met165. Van der Waals interactions are also formed with several amino acid residues of the active Fig. 4 . Moreover, the ligand properties of the reference inhibitors and top screened ligands of endophytic origins for both PLpro and Mpro are represented in Table 5 .

The interaction of a receptor and a ligand modifies the energies of both the free receptor and the ligand (ΔGBind); these energies also have a substantial influence on the stability of the receptor-ligand complex. In general, negative energy reflects a system's greater stability. Table 4 shows the various energies produced by MM-GBSA for 7CMD with native ligand GRL0617, as well as the top four docking hits. In this case, the binding energies of GRL0617 with both 7CMD (− 67.56 kcal/mol) are negative, indicating that the interactions are more spontaneous. Similarly, cytonic acid A is the most stable of the studied inhibitors, with a binding energy of − 58.96 kcal/mol, while cytonic acid B is the second most stable, with a binding energy of − 50.9 kcal/mol. The MM-GBSA energies of the remaining compounds examined were likewise negative, although the values were significantly higher, indicating a low occurrence rate for the complex. The compound aspergilline E had the greatest docking score with SARS-CoV-2-PLpro; however, its ΔGBind energy (− 50.81 kcal/mol) was determined to be greater than that of cytonic acid derivatives. In a similar manner, for Mpro, considering the ΔGBind scores of reference compounds and ligands, the reference ligand X77 has the best ΔGBind score of − 79.89 kcal/mol, while of all the other test ligands of endophytic origins, cytonic acid A has the best ΔGBind score of − 75.30 kcal/ mol, followed by cytonic acid B, asperphenalenone D and Table 4 Docking scores and the contributing binding residues of known PLpro inhibitor GRL-0617, Mpro inhibitor X77 and selected top antiviral metabolites generated using XP docking

Glide score (kcal/mol) Contributing binding residues

GRL-0617 (control) − 6.441  ASP164, GLN269, TYR268, TYR264, PRO247, PRO248, LEU162, GLY163, ARG166,  MET208, TYR273, GLY271, CYS270, THR301  Aspergilline E  − 7.77  PRO247, TYR268, ASP164, GLY163, GLY271, TYR264, LEU162, TYR273, GLN269 Table 6 .

Seven distinct models were used to measure the absorption property. CaCO 2 permeability and intestinal absorption were assessed; except for cytonic acid A and cytonic acid B, all of the substances have greater CaCO 2 permeability values.

Furthermore, all of the tested compounds including cytonic acid A and cytonic acid B as well as control substances had the poorest skin permeability, implying low absorption. For the P-glycoprotein assay, the lead compounds were found to be the compatible substrates as they could potentially be able to pass across the membrane using ATP-binding cassette (ABC) transporter. All substances tested as inhibitors of P-glycoproteins I and II were negative, except aspergilline E. The VDss assay is used to optimise the measurement of the total quantity of medications required for uniform distribution of pharmaceuticals throughout the blood. The assay values for all the drugs were less than 0.45 log L/kg, indicating that less medication volume will be required, but the values for the reference inhibitors GRL0617 and X77 were positive integer values. All the chemicals studied had negative CNS permeability values, indicating that they are projected to have decreased CNS permeability. Metabolism of the test drugs within the body was evaluated and except for asperphenalenone D and asperphenalenone A, which are anticipated to interact with CYP3A4, all other test chemicals were determined to be negative for all seven distinct models of metabolism. The total clearance and renal OCT2 assays were used to determine the drugs' excretion from the body. All the test chemicals, including GRL0167, tested negative for the renal OCT2 assay, indicating that none of them can be eliminated by organic cation transporter 2, apart from X77. Out of all ten distinct models that were used to measure toxicity, AMES test revealed that all four test molecules were negative, indicating that they are not carcinogenic or mutagenic, although GRL0617 was reported positive for this test. For cytonic acid A and cytonic acid B, maximum recommended tolerated dose values were found to be highest, that is 0.44 and 0.467 log(mg/kg/day), respectively. Table 7 contains the values of the ADMET analysis and properties of the test compounds along with reference compounds GRL0617 and X77.

Two isomeric chemicals, cytonic acids A and B, previously reported as inhibitors of human cytomegalovirus protease (hCMV), were observed to produce docking scores and MM-GBSA energies for both PLpro and Mpro. As a result, both compounds were subjected to an MD simulations experiment to determine the stability of the protein-ligand combination for 100 ns. In addition to the test substance (cytonic acids A and B), simulations of GRL0617 with SARS-CoV-2-PLpro and X77 with SARS-CoV-2-Mpro were run as a control set for each protein.

The RMSD values for all the frames in the trajectory were determined once the MD simulation was completed. The RMSD value of any simulation, in general, indicates the changes during the simulation of specific atoms with respect to their original state. The docked pose of the proteinligand complex was used as a baseline state or reference pose to compute the RMSD values. The RMSD value of protein atoms is represented by the Y-axis (left) in Figs. 5 and 6. The equilibration of the protein backbone for PLpro-GRL0617 (Fig. 5 ) occurs at an RMSD value of 2.3 ± 0.5 Å, whereas the equilibration of the protein backbone for Mpro-X77 (Fig. 6 ) occurs at an RMSD value of 2.2 ± 0.6 Å. For proteins that are bigger and non-globular in shape, the value might exceed beyond 4.0 Å. For the smaller and most globular proteins, RMSD value is ideally supposed to be in the limit of 1.0-4.0 Å. Here, the RMSD value after equilibration is not exceeding 4 Å for all the complexes under study at the peak end which is an advantageous parameter to screen the lead against the known ligand. For the complex of PLpro-cytonic acid A (Fig. 5) , RMSD value remains within the accepted limits of 2.0 ± 0.5 Å considering the size of protein, whereas for PLpro-cytonic acid B, RMSD value is equilibrated at 2.3 ± 0.5 Å. Similarly for complexes Mprocytonic acid A and Mpro-cytonic acid B (Fig. 6) , the protein backbone RMSD never exceeds 3.5 Å. In the graphs of Fig. 5 right Y-axis represents the ligand RMSD value. 'Lig fit Prot' signifies the RMSD values of ligand with reference to protein backbone. Throughout the whole simulation, 'Lig fit Prot' value of the PLpro-GRL0617 (Fig. 5) complex is in the range of 1.5 ± 0.5 Å, whereas for the PLpro-cytonic acid A complex this value is found to be 5.5 ± 1 Å. For, PLprocytonic acid B this value is 4.5 ± 1 Å which is slightly lower than that of cytonic acid A which suggests the relatively higher stability of the complex. On the other hand, 'Lig fit Lig' value of GRL0617 is lower than that of both cytonic acid isomers which suggests changes in the binding pose of ligand, but the difference is miniscule, hence negligible. Such results show that GRL0617 is having relatively higher stability compared to the tested ligand in binding pocket at the given docking pose. The 'Lig fit Prot' values for Mpro-cytonic acid A and Mpro-cytonic acid B the RMSD Fig. 3 Interaction profile of GRL0617 and best docked four antiviral compounds with SARS-CoV-2-PLpro (amino acids without any bond interactions are interacting with Van der Waals forces) ◂ values peak at 6.5 Å and 7.3 Å, respectively (Fig. 6) . For all the ligands, the 'Lig fit Prot' implies that the interactions between protein and ligand are quite stable. Additionally, Figs. 7, 8, 9, 10, 11 and 12 provide pictorial insights onto the entire simulation results based on collected trajectories and 10,000 captured frameshifts. Figure 7 represents the interactions made by reference ligand (GRL0617) and test compounds (cytonic acid A and cytonic acid B) with PLpro where the interaction types with percent interaction during 100-ns MD simulation run are represented. Similarly, for Mpro the interactions made by reference ligand (X77) and test compounds (cytonic acid A and cytonic acid B) are represented in Fig. 8 . Figure 9 represents the interaction fraction of reference ligand (GRL0617) and test compounds (cytonic acid A and cytonic acid B) with PLpro while Fig. 10 represents the interaction fraction of reference ligand (X77) and test compounds (cytonic acid A and cytonic acid B) with Mpro. Hydrogen bonds, hydrophobic interactions such as Pi cation, pi-pi stacking, water bridges and ionic interactions made by ligand with 

Coronaviruses (CoVs) are diverse group of viruses that have potentials to infect various mammals including humans. These viruses have positive single strand of RNA as its genetic material. These CoVs have their structure architecture composed of proteins, namely, S-, E-, M-and N-protein. α, β, γ and δ are four genera into which these CoVs are distributed, whereby α and β are the ones that have nag of infecting humans as their S-protein has structural complementary with human's ACE-2 receptor, enabling these viruses to attach human cell surface ACE-2 which is a vital interaction to cause infection [36] [37] [38] . Similarly, γ and δ CoVs have the tendency to infect Aves. As a reflection on the past, Middle East respiratory syndrome (MERS)-CoV and severe acute respiratory syndrome (SARS)-CoV deadly viruses belong to β genera of CoVs; also, the current pandemic causing SARS-CoV-2 belongs to β genera of CoVs. The recent variant Omicron that is ought to have very high infectivity is having high degree of mutations in S-protein which enables it to interact more efficiently with ACE-2; however, the target proteins used under present study, i.e. Mpro and PLpro, are identical with the previous variants of SARS-CoV-2; therefore, the research presented in this study is applicable to Omicron variant as well.

Aside from immunisation improvement with vaccine research, mainstream researchers have additionally given their energies and assets to distinguish intensifies that can counter the SARS-CoV-2 infection with alternative strategies [16, 39, 40] . Computationally driven molecular docking and MD recreations have been oftentimes used to identify compounds that possess potentials to interact with natural biochemistry of virus and halt its function [2, [41] [42] [43] . Assortment of metabolites as potential drugs from bacteria, moulds, yeast and plants has been investigated to restrain the urgent viral target proteins as their natural availability makes them unfathomably accessible. For example, pyranonigrin A, flaviolin and sterenin M are some secondary metabolites screened from moulds which by molecular docking and MD recreation studies with Mpro are recommended that they can restrain the viral replication by means of associating with this protein [44] [45] [46] . Besides, there are plethora of examinations successfully completed and published with in silico workflow utilizing sub-atomic powerful MD reenactment and docking assays that are very helpful in various different everyday issues of biological sciences [47] [48] [49] [50] [51] . Furthermore, these in silico methodologies have additionally been used for strength appraisal of hydroxychloroquine against different expected target proteins of SARS-CoV-2 [52] . In these computational examinations, the reasoning is extremely fundamental that in the event that any of these metabolites serving as potential drug leads can intrude with the viral proteins, they can thwart the typical life pattern of the viral particle [53] [54] [55] . Endophytic microbes are symbiotically associated with plants and have been proven to produce novel or analogues of host bioactive metabolites exhibiting a variety of biological activities including antiviral activity [56, 57] . However, compounds of specifically endophytic origin are not explored against SARS-CoV-2. In the current study, 45 compounds that have been reported active against viruses such as influenza, H1N1, HIV, HCV, hCMV and EV-71 were retrieved through literature and screened through rigorous computation workflow. Top compounds obtained after docking and MM-GBSA analysis for each protease (PLpro and Mpro) are known to have antiviral activity against different viruses. Top ligands for both the proteins in the present study, cytonic acids A and B, are known to possess in vitro inhibitory activities towards hCMV protease with the IC 50 values of 43 µmol and 11 µmol, respectively. Both the cytonic acids were extracted from the fermented culture of endophytic Fig. 7 Protein-ligand interaction diagram for PLpro showing percent of total time involved for particular interaction during molecular dynamics fungus Cytonaema sp., which was isolated from the plant Quercus sp. (European oak) [58] . Other compounds interacting with Mpro are asperphenalenone A and asperphenalenone D, which were obtained from the ethyl acetate fraction of fermented culture of Aspergillus sp., an endophytic fungus isolated from the plant Kadsura longipedunculata. Asperphenalenone A and D displayed anti-HIV activity with IC 50 s of 4.5 and 2.4 µM, respectively [59] . Aspergilline E was extracted along with aspergillines A-D from Aspergillus versicolor, an endophyte isolated from the plant Paris polyphylla var. yunnanensis, collected from China. Aspergilline E showed weak antiviral activity towards tobacco mosaic virus with IC 50 of 33.6 ± 3.0 μM [60] . In the present study, it interacts with the PLpro with multiple interaction types. Another compound having strong interaction with PLpro is 4-dehydroxyaltersolanol A, which was obtained from endophytic fungus Nigrospora sp. YE3033, harboured in Aconitum carmichaelii. 4-Dehydroxyaltersolanol A displayed in vitro antiviral activity towards H1N1 influenza A with IC 50 of 8.35 ± 1.41 µg/mL [61] . With all these evidences, it is quiet insightful that cytonic acids A and B, asperphenalenone A, asperphenalenone D and aspergilline E have antiviral activities against different forms of viruses and could be of potential against SARS-CoV-2; of course, in vitro assays are Fig. 9 Protein-ligand interaction diagram showing interaction fraction of crucial interacting amino acids of PLpro needed to validate the findings, but such in silico study helps us to identify potential hits out of large pools of secondary metabolites and it is analogous of finding a purple wooden stick from the Amazon forest with ease.

The best docked compounds are cytonic acid A and cytonic acid B, which have significant interactions with both PLpro and Mpro proteases, according to the MD simulation studies. Considering the lack of treatment options for coronavirus infection, the current study proposes cytonic acids A and B as a dual inhibitor that interacts with critical amino acids in the binding site of PLpro and Mpro to inhibit their function, which could be useful in future in vitro and in vivo studies for COVID-19 therapeutics. Because we are dealing with infectious agents, in vitro validation of the suggested drugs mandates the use of Biosafety Level 4 (BSL4). Only an in silico investigation has been conducted. The information offered in this study will be beneficial to the scientific community. Even so, this work could potentially untap many more such exploratory work with endophytic compounds.

Curse of La Corona: unravelling the scientific and psychological conundrums of the 21st century pandemic

ABBV-744 as a potential inhibitor of SARS-CoV-2 main protease enzyme against COVID-19

Rally around your fellows: information and social trust in a real-world experiment during the corona crisis

COVID-19 vaccines: where we stand and challenges ahead

SARS-CoV-2 vaccines strategies: a comprehensive review of phase 3 candidates

A review on COVID-19 vaccines: stages of clinical trials, mode of actions and efficacy

Evolutionary analysis of the delta and Delta Plus variants of the SARS-CoV-2 viruses

Antibody resistance of SARS-CoV-2 variants B.1.351 and B.1.1.7

In search of RdRp and Mpro inhibitors against SARS CoV-2: molecular docking, molecular dynamic simulations and ADMET analysis

Druggable targets from coronaviruses for designing new antiviral drugs

In silico studies on the interaction between Mpro and PLpro from SARS-CoV-2 and ebselen, its metabolites and derivatives

An integrative in silico drug repurposing approach for identification of potential inhibitors of SARS-CoV-2 main protease

Natural products for drug discovery in the 21st century: innovations for novel drug discovery

Natural products' role against COVID-19

Isolation and characterization of bioactive metabolites from Xylaria psidii, an endophytic fungus of the medicinal plant Aegle marmelos and their role in mitochondrial dependent apoptosis against pancreatic cancer cells

Berries anthocyanins as potential SARS-CoV-2 inhibitors targeting the viral attachment and replication; molecular docking simulation

Identifying structural -functional analogue of GRL0617, the only well -established inhibitor for papain -like protease ( PLpro ) of SARS -CoV2 from the pool of fungal metabolites using docking and molecular dynamics simulation

Systematic virtual screening in search of SARS CoV-2 inhibitors against spike glycoprotein: pharmacophore screening, molecular docking, ADMET analysis and MD simulations

In search of SARS CoV-2 replication inhibitors: virtual screening, molecular dynamics simulations and ADMET analysis

Bioactive anthraquinones from endophytic fungus Aspergillus versicolor isolated from red sea algae

New azaphilones, phomopsones A-C with biological activities from an endophytic fungus Phomopsis sp. CGMCC No.5416

Different culture metabolites of the Red Sea fungus Fusarium equiseti optimize the inhibition of hepatitis C virus NS3/4A protease (HCV PR)

Azaphilone and isocoumarin derivatives from the endophytic fungus Penicillium sclerotiorum PSU-A13

Pharmacophore based virtual screening, molecular docking, molecular dynamics and MM-GBSA approach for identification of prospective SARS-CoV-2 inhibitor from natural product databases

Identifying structural-functional analogue of GRL0617, the only well-established inhibitor for papain-like protease (PLpro) of SARS-CoV2 from the pool of fungal metabolites using docking and molecular dynamics simulation

The complex structure of GRL0617 and SARS-CoV-2 PLpro reveals a hot spot for antiviral drug discovery

Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments

Three-dimensional pharmacophore methods in drug discovery

Protein and ligand preparation : parameters, protocols, and influence on virtual screening enrichments

PHASE: a novel approach to pharmacophore modeling and 3D database searching

Biomolecular simulations: recent developments in force fields, simulations of enzyme catalysis, protein-ligand, protein-protein, and proteinnucleic acid noncovalent interactions

Recent advances in free energy calculations with a combination of molecular mechanics and continuum models

Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models

Combined molecular mechanical and continuum solvent approach (MM-PBSA/GBSA) to predict ligand binding

pkCSM: Predicting small-molecule pharmacokinetic and toxicity properties using graph-based signatures

The newly emerged SARS-Like coronavirus HCoV-EMC also has an "Achilles' heel": current effective inhibitor targeting a 3C-like protease

Computational and in-vitro validation of natural molecules as potential acetylcholinesterase inhibitors and neuroprotective agents

Drug design targeting the main protease, the Achilles heel of coronaviruses

Evaluation of potential anti-RNA-dependent RNA polymerase (RdRP) drugs against the newly emerged model of COVID-19 RdRP using computational methods

Identification of potential inhibitors of coronavirus hemagglutinin-esterase using molecular docking, molecular dynamics simulation and binding free energy calculation

Drug design and repurposing with DockThor-VS web server focusing on SARS-CoV-2 therapeutic targets and their non-synonym variants

Molecular docking study of potential phytochemicals and their effects on the complex of SARS-CoV2 spike protein and human ACE2

Searching for target-specific and multi-targeting organics for Covid-19 in the Drugbank database with a double scoring approach

Proposing a fungal metabolite-flaviolin as a potential inhibitor of 3CLpro of novel coronavirus SARS-CoV-2 identified using docking and molecular dynamics

Sterenin M as a potential inhibitor of SARS-CoV-2 main protease identified from MeFSAT database using molecular docking, molecular dynamics simulation and binding free energy calculation

Reckoning a fungal metabolite, pyranonigrin A as a potential Main protease (Mpro) inhibitor of novel SARS-CoV-2 virus identified using docking and molecular dynamics simulation

Rapid identification of potential inhibitors of SARS-CoV-2 main protease by deep docking of 1.3 billion compounds

QSAR modeling of SARS-CoV Mpro inhibitors identifies sufugolix, cenicriviroc, proglumetacin, and other drugs as candidates for repurposing against SARS-CoV-2

In silico drug repurposing for COVID-19: targeting SARS-CoV-2 proteins through docking and consensus ranking

Breaking bad: better call gingerol for improving antibiotic susceptibility of Pseudomonas aeruginosa by inhibiting multiple quorum sensing pathways

Revealing the molecular interplay of curcumin as Culex pipiens acetylcholine esterase 1 (AChE1) inhibitor

Structural insight to hydroxychloroquine-3C-like proteinase complexation from SARS-CoV-2: inhibitor modelling study through molecular docking and MD-simulation study

The rise of gingerol as anti-QS molecule: darkest episode in the LuxR-mediated bioluminescence saga

Identification of promising compounds from curry tree with cyclooxygenase inhibitory potential using a combination of machine learning, molecular docking, dynamics simulations and binding free energy calculations

Comparative assessment of RNA-dependent RNA polymerase (RdRp) inhibitors under clinical trials to control SARS-CoV2 using rigorous computational workflow

Endophytic fungi: a treasure trove of novel anticancer compounds. Current Research in Pharmacology and Drug Discovery 100050

Untapped potentials of endophytic fungi: a review of novel bioactive compounds with biological applications

Cytonic acids A and B: novel tridepside inhibitors of hCMV protease from the endophytic fungus Cytonaema species

Metabolites from the plant endophytic fungus Aspergillus sp

CPCC 400735 and their anti-HIV activities

Aspergillines A−E, highly oxygenated hexacyclic indole−tetrahydrofuran−tetramic acid derivatives from Aspergillus versicolor

Antiviral anthraquinones and azaphilones produced by an endophytic fungus Nigrospora sp. from Aconitum carmichaeli

Sesquiterpenoids and xanthones derivatives produced by sponge-derived fungus Stachybotry sp. HH1 ZSDS1F1-2. The Journal of Antibiotics advance online publication

Antiviral isoindolone derivatives from an endophytic fungus Emericella sp. associated with Aegiceras corniculatum

Chloropupukeananin, the first chlorinated pupukeanane derivative, and its precursors from Pestalotiopsis fici

/5 tricyclic ring system from the endophytic fungus Periconia sp

Altertoxins with potent anti-HIV activity from Alternaria tenuissima QUE1Se, a fungal endophyte of Quercus emoryi

Xiamycin, a pentacyclic indolosesquiterpene with selective anti-HIV activity from a bacterial mangrove endophyte

Neosartoryadins A and B, fumiquinazoline alkaloids from a mangrovederived fungus Neosartorya udagawae HDN13-313

Microbial natural product alternariol 5-O-methyl ether inhibits HIV-1 integration by blocking nuclear import of the pre-integration complex