Molecular docking and ADMET study of bioactive compounds of Glycyrrhiza glabra against main protease of SARS-CoV2.

Recent pandemic situation of COVID-19 is caused due to SARS-CoV2 and almost all the countries of the world have been affected by this highly contagious virus. Main protease (Mpro) of this virus is a highly attractive drug target among various other enzymes due to its ability to process poly-protein that is the translated product of the SARS-CoV2 RNA. The present study demonstrates molecular docking study of Glycyrrhiza glabra (Gg) active compounds such as Glycyrrhizic acid (GA), Liquiritigenin (L) and Glabridin (G) against the Mpro. Docking studies shows that these active compounds bind strongly with some of the amino acid residues in the active site of Mpro and inhibits the enzyme strongly. GA, L, and G are proposed to be strong inhibitors of the enzyme and the amino acids: His41, Gly143, Gln189, Glu 166, Cys 145, Thr25, Asn142, Met49, Cys44, Thr45 and pro168 present in the active site of Mpro were shown to make non-covalent interaction with these compounds. In silico ADMET properties prediction also shows that Gg active compounds had good solubility, absorption, permeation, non-toxic, and non- carcinogenic characteristics. Our finding concludes that all of the three active compounds of Gg have the potential to be strong inhibitors for Mpro of SARS-CoV2 but glycyrrhizic acid has a high binding affinity and a good ADMET properties than the other two.

Introduction

Wuhan, Hubei province of China was the very first location where the severe acute respiratory syndrome coronavirus2 (SARS-CoV2) outbreak recorded [1]. From Wuhan city of China this virus had spread across more than 200 countries with America and Brazil as the epicentre as of now. COVID-19 caused by SARS-CoV2 has affected 8, 242, 999 people worldwide and killed as many as 445, 535 people, global data as of 18th June 2020, 4.57 pm CEST. (WHO coronavirus disease Dashboard). World Health Organization (WHO) has stated the COVID-19, a public health emergency of international concern on January 30, 2020 [2]. SARS-CoV [3], MERS-CoV [4] and SARS-CoV2 [5] are the three categories of coronavirus that causes the deadly pneumonia in humans. Severe acute respiratory syndrome-2 (SARS-CoV2) is an enveloped, positive-sense single-strand RNA placed in a beta coronavirus family of Coronaviridae [6]. As compared to the other two, SARS-CoV2 is highly contagious and spread the infection more quickly and severely, leading to large number of deaths and thus arising a pandemic situation [7], [8].

The SARS-CoV2 has a genome ranging from 29 to 30 kb in length [9]. The recently proposed drugs targets the main protease (Mpro also known as 3CLpro) which is one of the best targets among the other enzymes in the virus [10]. This enzyme participates in the processing of polyprotein that is translated from the RNA of the same virus [11]. These polyprotein are processed and converted into mature non-structural proteins using two proteases PLpro and CLpro encoded by ORF1 (open reading frame1) [12]. A non-structural protein of coronavirus plays a significant role in replication and transcription during the infection [13]. Main protease operates at more than 10 cleavage sites on the polyprotein 1ab which is nearly 790 kDa [14]. Mpro is a homodimer having two identical chains and the recognition sequence is Leu-Gln ↓ Ser, Ala, Gly and cleaved through the Cys-His dyad protease in which -SH group of cysteine serve as a nucleophile in the proteolysis [15]. Inhibiting the main protease activity would block transcription or replication. The amino acid residues of the active site of Mpro participating in substrate binding includes Cys145 –His41 dyad ,and Phe140, Thr45, Arg188, Asp187, Met49, Asn142, Met165, Gln189, His172, Glu166 [16]. An enormous amount of data is available for describing the efficiency of an array of various active components of medicinal plants effective against the viruses. Among these medicinal plant, Glycyrrhiza glabra is the one plant, its active compounds may have the potential to act against various viruses [17]. Glycyrrhiza glabra (Gg) commonly called ‘sweet wood’ and ‘liquorice’ belongs to the family Fabaceae. It has been effective against many RNA viruses like H5N1virus, Hepatitis C virus, influenza A virus, Newcastle disease virus, H1N1 virus, Rotavirus, severe acute respiratory virus etc. and DNA viruses such as Varicella Zoster virus, Epstein-Barr virus and Herpes Simplex virus etc. [18].

Several active compounds from Gg have been reported that inhibit the virus’s growth. Some studies showed that Glycyrrhizin (glycyrrhizinate; glycyrrhizic acid) inhibit the binding of the virus to the target cells, and also beneficial in controlling the viral replication and found to have an important antiviral activity. Studies on clinically isolates of severe acute respiratory syndrome Virus [FFM-1 and FFM-2] gave important information about the anti-viral property of glycyrrhizic acid [19], [20]. Another active compound, Glabridin (G) isolated from Gg roots. It has been used as a traditional medicine in Asia and has anti-inflammatory, anti-osteoporotic, anti-nephritic, anti-antheroenic, estrogenic, and anti-oxidant, regulation of metabolism, neuro-protective and skin-whitening properties. This is also reported to have vasodilation property and used for the treatment of many diseases like bacterial, bronchitis infection [27], [28]. Liquiritigenin (L) has anti-cancer, anti-oxidative, antimicrobial activities, hepato-protective, immune regulatory, and cardio-protective effects and also has anti-platelet aggregation, anti-tumorigenic, anti-angiogenic, anti-allergic properties [29].

For the virtual screening and structure-based drug designing, molecular docking is commonly used in bioinformatics based research and is accustomed to predict non-covalent interaction (mainly hydrogen bonding) of macromolecule (receptor) with the ligand (drug molecule). AutoDock is a suited program for docking and virtual screening which calculates the grid internally (automatically), for the atom types that are required. A more improved version, AutoDock Vina achieves within the average accuracy of the binding pose prediction and bring order of magnitude faster than previously developed AutoDock4 [30].

OSIRIS Property Explorer is employed to determine drug Likeness. It provides drug relevant information whenever a structure is valid and underlines various properties with a high risk of inadmissible effects such as poor intestinal absorption, mutagenicity or drug conform behaviour through different customization characteristics. This also determines the clog P (Logarithm of compounds partition coefficient between water and n-Octanol) which determines the hydrophilicity of the compounds. Larger the clogP value, lower the hydrophilicity and thus poor the permeation and absorption. Log S suggest solubility; lower the log S value, higher the solubility which enhances the absorption. The Topological polar surface area (TPSA) suggest the surfaces that belong to the polar atoms and molecules in the compound. Larger TPSA value is related with the least permeability of the membrane. The compounds that have larger TPSA value will be a better substrate for p-glycoprotein which is liable for efflux of drug from the cell and thus reduced TPSA was beneficial for drug-like property. Some studies also anticipated that a molecule having a better penetration property through CNS should have lesser value of TPSA [26].

A good drug candidate can be consumed in the desired time as well as distributed through the entire body for its efficient action and metabolism. Another factor is toxicity which is often very vital and dominates the Absorption, Distribution, Metabolism, Excretion behaviour of the drugs. Some studies have shown that drug failed at clinical trials due to its adverse side effect and their toxicity which has been proven to be very detrimental and expensive in the drug development. In silico, ADMET and drug-likeness prediction help in the discovery of new targets and compounds with anticipated biological activities [31]. Human intestinal absorption [HIA] is an important ADMET properties. The exploitation of drugs in the body is such a complex process that can hardly analysed by statistical models. HIA is one of the important steps during compounds (drug) transporting to their desired target [31]. Blood-brain barrier [BBB] is used to predict various features of the CNS vasculature. CNS vessel does not have pores in the endothelial cell and have various properties that regulate the transport of ions, molecules, and cells which make CNS vessel very prohibitive in nature and provides an interference for the delivery of compounds into the central nervous system [32]. Carcinogenicity means the ability of the substance to cause cancer. P-glycoprotein substrate is a compound that utilizes the P-glycoprotein transporter for several exercises like excretion and absorption of drugs and other vital exercises [33]. These properties are often used for the study of ADMET behaviour of the drugs.

The main objective of the present study is to carry out molecular docking analysis of Glycyrrhiza glabra active compounds, Glycyrrhizic acid, Liquiritigenin and Glabridin against the main protease (Mpro) one by one followed by molecular interaction study (hydrogen bond prediction between target and drugs), drug-likeness behaviour and ADMET prediction to confirm the efficiency and efficacy of these active compound against SARS-CoV2.

Methods and materials

Selection of ligand

Glycyrrhizic acid (PubChem CID: 14982), Liquiritigenin (PubChem CID: 114829), and Glabridin (PubChem CID: 124052) were recognized as potential SARS-CoV inhibitors from literature survey and maybe also effective against SARS-CoV2. The structure of Liquiritigenin and Glabridin, and Glycyrrhizic acid were downloaded from PubChem (http://pubchem.ncbi.nlm.nih.gov) two control drug that is lopinavir (PubChem CID: 92727) and rotinavir (PubChem CID: 39622) also downloaded and converted their file format from SDF file to PDB file [Supplementary Fig. 1 (a) – (c)].

Selection of target molecule (Receptor/Mpro): Active site and sequence analysis

The crystal structure of the main protease (PDB ID- 6Y2E) downloaded from PDB database https://www.rcsb.org/structure/6Y2E, DOI: https://doi.org//10.2210/pdb6Y2E/pdb [25] and saved as PDB format. Its functional domain and the number and name of amino acid residues present in the pocket of the active site of 6Y2E were verified using ScanProsite (ExPASy) (http://prosite.expas.org) (supplementary figure S4) [21] and CDD webserver was used to determine the conserved sequence of 6Y2E [22], Sequence analysis of the main protease of SARS-CoV2 was retrieved using ProtParam, ProtScale (http://web.expasy.org) [23].

Preparation of ligands (Glycyrrhiza glabra active compounds) and main protease (Mpro)

The Ligand was prepared by adjusting ionization, torsion, degree of freedom and stereo-chemical variation using AutoDock MG Tool. Energy of the main protease was minimized based on GROMAS6 43B1 algorithm using Swiss Model Viewer, water was removed and Polar hydrogens were added to the main protease of SARS-CoV2 followed by calculation of gasteiger Charges using the AutoDock MG Tool and saved in PDBQT file format. The receptor grid was developed using AutoDock Grid tool with grid dimensions of 30A0 × 30A0 × 30 A0 with 0.375 A0 spacing and the grid box was set at −15.692, –32.457, 2.498 in X, Y and Z dimension [24].

Molecular docking of Glycyrrhiza glabra (Gg) active compounds with Mpro

The crystal structure main protease (Mpro) downloaded from RCSB Protein Data Bank (PDB ID: 6Y2E) with a resolution of 1.7 A0 and R-value free of 0.222 and R-value work of 0.171. Molecular docking was performed using in silico docking software, AutoDock Vina [30] with the three active compounds of Gg, glycyrrhizic acid, liquiritigenin and glabridin one by one against 6Y2E. The docking was done using command Prompt and the other prerequisite conditions before the docking like the ligand and enzyme preparation has been established. PyMol and BIOVIA Discovery Studio Visualizer was used to visualize the binding interaction between the active compounds of Glycyrrhiza glabra with 3D structure of M pro (6Y2E) of SARS-CoV2 [25] and Simulation was done with CHARMM36 force field in Discovery Studio visualisation tool.

Drug likeness prediction

OSIRIS Property Explorer (http://www.organic-chemistry.org/prog/peo/) [26] is employed to determine likeness of drugs. Properties like Log S calculation, TPSA, Clog P calculation, molecular mass, and drug-likeness based on fragment and drug score of all the three active compounds have been established.

ADMET properties prediction

ADMET refers to Absorption, Distribution, Metabolism, Excretion, and Toxicity. It contains the pharmacokinetic profile of a compound (drug molecule) and plays a significant role in determining its pharmacodynamics activities. Taking all the three active compounds into account, properties like bioavailability, brain penetration, oral absorption, carcinogenicity, and other human intestinal absorption properties of the active compounds have been determined using admetSAR webserver (http://Immd.ecust.edu.cn:8000/) [31], [32], [33], [34].

Results and discussion

Sequence analysis and active site identification

Sequence analysis of SARS-CoV2 Mpro revealed a theoretical PI value of 5.95 with an extinction coefficient of 33640 M−1 cm−1 and absorbance (1 g/L) computed 0.995. The Aliphatic index of the Mpro was 82.12 and shows a stability index of 27.65. All these result shows that SARS-CoV2 Mpro is a very stable macromolecule. The secondary structure of SARS-CoV2 main protease revealed that the alpha helix and the beta strand were maximally present (supplementary table 1). The ExPASy Prosite and CDD webserver result shows that His41 and Cys145 has a vital role in the catalytic activity of the main protease of SARS-CoV2 and also has been conserved in all the different types of coronaviruses as depicted in Fig. 1 .

Fig 1

Sequence alignment using CDD webserver represent His41 and Cys145 as highly conserved amino acids (denoted by #) in all coronaviruses and actively participated in cleavage of polyprotein of SARS-CoV2 into structural proteins.

Molecular docking analysis of control drugs with the Mpro of the SARS-CoV2

The molecular docked pose of the minimum binding affinity conformer of the Lopinavir and Ritonavir (control drugs) demonstrated that they bind to the active site of the Mpro protein of the SARS-CoV2. Fig. 2 showed that Lopinavir binds through two conventional H-bond with Asn142, two carbon H-bond with Thr25 and Gln189, and two alkyl bond with His41 respectively. However, Ritonavir also bind firmly to the active site of the main protease and it stabilized the Mpro active site through four hydrogen Bond in which two are conventional hydrogen bond and rest two are carbon hydrogen bond. Asn142 and Thr25 amino acid residue formed conventional and carbon hydrogen bond respectively. Whereas Cys44 and Cys145 formed alkyl bond and His41 and Pro168 shows Pi-alkyl interaction, besides other residues also found to be bound to control drugs through Vander Waals interaction depicted in 2D (two dimensional) plot (Fig. 3 ). The binding affinity energy of the control drugs are mentioned in table 1 . If we compare the binding affinity energy of control drugs, Lopinavir and Ritonavir, we observed that Lopinavir have highest negative binding affinity with −7.6 Kcal/mol followed by Ritonavir with −7.3 Kcal/mol [35].

Fig. 2

3-D (three-dimensional) diagram based on hydrogen bond donor and acceptor characteristics of amino acid residues (A) and 2-D (two dimensional) (B) diagram showing several interaction between Mpro and Lopinavir. CYS: Cysteine, THR: Threonine, ASN: Asparagine.

Fig. 3

3-D (three dimensional) diagram based on hydrogen bond donor and acceptor characteristics of amino acid residues (A) and 2-D (two dimensional) (B) diagram showing several interaction between Mpro and Ritonavir. PRO: Proline, THR: Threonine, LEU: Leucine, HIS: Histidine, ASN: Asparagine.

Table 1

Molecular Docking analysis: Binding energy and the type of molecular interaction of control drugs against Mpro. UNK: Glycyrrhiza glabra active compounds.

Control Drugs	Binding Affinity (Kcal/mol)	Donor Atom	Distance	Acceptor Atom	Category	Types of interaction
Lopinavir	−7.6	ASN142:HD22	2.39	UNK: O	Hydrogen Bond	Conventional H-bond
		UNK:H	1.75	ASN142:OD1	Hydrogen Bond	Conventional H-bond
		UNK:H	2.63	THR25:OG1	Hydrogen Bond	Conventional H-bond
		UNK:H	1.72	CYS44	Hydrogen Bond	Conventional H-bond
		UNK	3.98	UNK	Hydrophobic	Pi-Pi Stacked
		UNK	4.98	CYS145	Hydrophobic	Pi-Alkyl
Ritonavir	−7.3	ASN142:HD22	2.63	UNK: O	Hydrogen Bond	Conventional H-bond
		UNK:H1	2.43	ASN142:OD1	Hydrogen Bond	Conventional H-bond
		UNK:H13	2.61	THR25:OG1	Hydrogen Bond	Carbon H-bond
		UNK:H34	2.23	GLN189:OE1	Hydrogen Bond	Carbon H-bond
		UNK:C12	5.14	LEU27	Hydrophobic	Alkyl
		UNK:C12	3.87	CYS145	Hydrophobic	Alkyl
		HIS41	4.55	UNK:C12	Hydrophobic	Pi-Alkyl
		UNK	5.33	PRO168	Hydrophobic	Pi-Alkyl

Molecular docking analysis of Gg active compounds with the Mpro of the SARS-CoV2

The molecular docked pose of the minimum binding affinity conformer of the Gg active compounds that is glycyrrhizic acid, Glabridin and Liquiritigenin demonstrated that they also firmly goes and bind to the active site of the Mpro protein of the SARS-CoV-2 (Fig. 4 and table 2 ). Glycyrrhizic acid (GA) binds within the Mpro active site through several non-covalent interactions. GA binds active site through one conventional H-bond with Glu166 and three carbon H-bond with Ser46, Pro168 and Asn142. GA also found to be bound to active site through alkyl and Pi-alkyl bond with Cys145, Leu27 and His41 respectively (Fig. 7). Glabridin stabilizes the Mpro active site through hydrogen bond with His41 and Glu166 and through alkyl and Pi-alkyl with Cys145 and Pro168 respectively and play an important role in the stabilisation of the active site of the main protease (Fig. 6 ). Likewise Liquiritigenin also binds to the active site and formed several non-covalent interaction in which Thr25 and Thr45 are found to bind with the Liquiritigenin through hydrogen bond and, Met49, and Cys145 formed Pi-sulfur interaction to Liquiritigenin L also firmly binds to the His41 through Pi-Alkyl that stabilize the active site (Fig. 5). On comparing the binding affinity of the control drugs with the binding affinity of Gg active compounds, we found that GA have the highest negative binding affinity energy than Lopinavir and ritonavir with −8.0 Kcal/mol, G have lower binding affinity than Lopinavir but have similar binding affinity with ritonavir with −7.3 Kcal/mol whereas L have lowest binding affinity among all ligands [36].

Fig 4

Three dimensional structure of main Protease based on the hydrophobicity of amino acid residues and docked compounds are shown in black.

Table 2

Molecular Docking analysis: Binding energy and the type of hydrogen bond interaction of all the three active compounds of Gg against Mpro. UNK: Glycyrrhiza glabra active compounds.

Gg active compounds	Binding Affinity (Kcal/mol)	Donor Atom	Distance	Acceptor Atom	Category	Types of interaction
G	−7.3	HIS41:HE2	2.37	UNK:O1	Hydrogen Bond	Conventional H-bond
		UNK:H	1.84	GLU166: O	Hydrogen Bond	Conventional H-bond
		MET49:SD	5.46	UNK	Other	Pi-Sulfur
		CYS145	4.43	UNK	Hydrophobic	Alkyl
		UNK:C12	3.99	CYS145	Hydrophobic	Alkyl
		HIS41	5.34	UNK	Hydrophobic	Pi-Alkyl
		HIS41	5.32	UNK:C12	Hydrophobic	Pi-Alkyl
		UNK	5.2	PRO168	Hydrophobic	Pi-Alkyl
L	−7.0	UNK:H	2.17	THR25: O	Hydrogen Bond	Conventional H-bond
		THR45:H	2.84	UNK:O2	Hydrogen Bond	Carbon H-bond
		MET49:S	5.6	UNK	Other	Pi-Sulfur
		CYS145:S	3.64	UNK	Other	Pi-Sulfur
		HIS41	5.82	UNK	Hydrophobic	Pi-Pi Stacked
GA	−8.0	UNK:H7	2.51	GLU166: O	Hydrogen Bond	Conventional H-Bond
		SER46:HB2	2.47	UNK:O12	Hydrogen Bond	Carbon H-Bond
		PRO168:HD2	2.8	UNK:O14	Hydrogen Bond	Carbon H-Bond
		UNK:H54	2.46	ASN142:OD1	Hydrogen Bond	Carbon H-Bond
		CYS145	5.21	UNK	Hydrophobic	Alkyl
		UNK:C28	5.19	LEU27	Hydrophobic	Alkyl
		HIS41	4.58	UNK:C28	Hydrophobic	Pi-Alkyl

Fig. 7

3-D (three dimensional) diagram based on hydrogen bond donor and acceptor characteristics of amino acid residues (a) and 2-D (two dimensional) (b) diagram showing several interaction between Mpro and Glycyrrhizic acid (GA). HIS: Histidine, GLU: Glutamic acid, LEU: Leucine, ASN: Asparagine, CYS: Cysteine, PRO: Proline.

Fig. 6

3-D (three dimensional) diagram based on hydrogen bond donor and acceptor characteristics of amino acid residues (A) and 2-D (two dimensional) (B) diagram showing several interaction between Mpro and Glabridin (G). PRO: Proline, MET: Methionine, GLU: Glutamic acid, CYS: Cysteine, HIS: Histidine.

Fig. 5

3-D (three dimensional) diagram based on hydrogen bond donor and acceptor characteristics of amino acid residues (a) and 2-D (two dimensional) (b) diagram showing several interaction between Mpro and Liquiritigenin (L).

The result of the OSIRIS Property Explorer (http://www.organic-chemistry.org/prog/peo/) shown in table 3 . Out of the three active compounds, GA shows a very low Clog P than G and L and thus GA has high hydrophilicity and can be absorbed easily by the cell of the SARS-CoV2 thus showing drug efficacy but in terms of solubility, GA is least soluble in blood than the other two. L shows maximum solubility in the blood. TPSA of GA is also higher as compared with the other two compounds which tell the efficiency of the compounds in terms of membrane permeability. Here again GA has more potential to cross the lipid bilayer of the virus and thus can strongly bind and inhibits the enzyme (Mpro). Our predicted values lie well within the drug-likeness.

Table 3

Drug likeness prediction using OSIRIS Property Explorer, clog P (O/W): Logarithm of partition coefficient between n-octanol and water, logS: aqueous solubility, TPSA: Topological polar surface area.

Gg active compounds	clog P(≤5)	Solubility: log S(-6.5 to 0.5)	Molecular mass	TPSA	Overall drug score
GA	0.99	−5.14	822	267	0.19
G	4.31	−4.15	324	58.92	0.24
L	2.5	−2.94	256	66.76	0.42

ADMET prediction

ADMET (Absorption, Distribution, Metabolism, Excretion and Toxicity) properties were retrieved using admetSAR server, which shows that all Glycyrrhiza glabra active compounds i.e. Glycyrrhizic acid (GA) (Accession Number: DB13751), Liquiritigenin (L) (PubChem CID: 114829) and Glabridin (G) (PubChem CID: 124052) are having better human intestinal absorption score (table 4 ). Higher human intestinal absorption indicates that the compound of interest may be better consumed from the gastrointestinal tract upon oral administration. The penetration through the Blood-Brain Barrier is computed to be the best for G and comparatively higher than L and almost similar to GA. In terms of predicting the efflux by P-glycoprotein from the cell, GA comes out to be inhibitor and non-substrate of P-glycoprotein, G comes out to be non-inhibitor and substrate and L comes out to be non-inhibitor and non-substrate. An inhibitor of P-glycoprotein means that drug will inhibit the efflux process from the cell and enhances the bioavailability and non-inhibitor of P-glycoprotein mean that drug will efflux from the cell by P-glycoprotein and limits the bioavailability by pumping back into the lumen and may promote the elimination of that drug into the bile and urine [34]. In the case of metabolism, GA and G found to be non-inhibitor and substrate while L found to be inhibitor and non-substrate. A non-inhibitor of cytochrome P450 means that the molecule will not hinder the biotransformation of the compound (drug) metabolized by cytochrome P450 [31]. The Carcinogenic profile also shows that GA, G, and L are non-carcinogenic so it can be applied as drugs for treating COVID-19 as there would not be any bioaccumulation of compounds in the human body and these compounds would not lead to cancer in future if the patient is treated for a long duration. Acute oral toxicity was found to be higher in G followed by GA and L. It shows that L has less oral toxicity than GA which further suggests that the drugs L are not metabolized in the gastrointestinal tract before reaching the desired target. All these results have also been depicted in the form of a graph (Fig. 8 A, B and table 4).

Table 4

ADMET profile of Gg active compounds, HIA: Human intestinal absorption, BBB: Blood Brain Barrier, CYP450: Cytochrome P450.

S.NO	Gg active compounds	HIA (probability)	BBB (probability)	Acute oral Toxicity (Kg/mol)	CYP450 3A4 inhibitor/substrate	Carcinogenicity
1.	GA	0.8125	0.8793	2.13283	Non-Inhibitor/substrate	Non carcinogenic
2.	G	0.9948	0.9312	3.11555	Non-Inhibitor/substrate	Non carcinogenic
3.	L	0.9947	0.688	1.188737	Inhibitor/Non-substrate	Non carcinogenic

Fig. 8

A. Comparative BBB, HIA and B. Acute oral toxicity of Gg active compounds.

Conclusion

The present studies conclude the presence of His and Cys as the conserved amino acid of Mpro of SARS-CoV2. Molecular docking indicated that the three active compounds of Glycyrrhiza glabra namely glycyrrhizic acid, Liquiritigenin, and Glabridin successfully docked with the amino acid molecule at the catalytic site of the Mpro with a high negative binding affinity and formed several molecular interaction with the main protease of SARS-CoV2. They may form a complex with the Mpro of SARS-CoV2 causing inhibition of the catalytic activity of the main protease. His41 and Cys145 are conserved amino acids found in different types of coronavirus family and also plays a crucial role in enzyme catalysis and our docking result also showed that His41 and Cys145 are participated in the non-covalent interaction with all the three active compounds of Glycyrrhiza glabra. So from this study it is concluded that these active compounds exactly occupy the same location and binds these amino acid residues of the enzyme which otherwise would occupied by the natural substrates for the enzyme (Mpro), so these compounds can be replaced as a proposed drug for the main protease of SARS-COV2. admetSAR results show that glycyrrhizic acid would be more suitable drug candidate among the other two compounds as it is more selective, potent, non-carcinogenic, and non-tumorigenic. OSIRIS property explorer had shown Glabridin to have toxicity to the reproductive system and Liquiritigenin to be mutagenic and hence less likely be a good drug candidate than glycyrrhizic acid. The outbreak of SARS-CoV2 in December 2019 having no antiviral drug attainable for instant relief in the initial stage, This computational study may provide some promise for the designing of new antiviral drugs with least side effect and more target oriented. Our future aim would be to confirm our in-silico results by in-vitro and in-vivo studies.

CRediT authorship contribution statement

Vivek Srivastava: Conceptualization, Supervision. Ankush Yadav: Methodology, Software, Data curation, Writing - original draft. Paratpar Sarkar: Visualization, Investigation, Writing - review & editing.

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.