Saurov Mahanta1†
Tufan Naiya2†
Kunal Biswas3†
Liza Changkakoti1†
Yugal Kishore Mohanta4
Bhaben Tanti5
Awdhesh Kumar Mishra6*
Tapan Kumar Mohanta7*
Nanaocha Sharma8*- 1National Institute of Electronics and Information Technology (NIELIT), Guwahati, India
- 2Department of Biotechnology, Maulana Abul Kalam Azad University of Technology, West Bengal, India
- 3Centre for Nanoscience and Nanotechnology, Sathyabama Institute of Science and Technology, Chennai, India
- 4Department of Applied Biology, School of Biological Sciences, University of Science and Technology Meghalaya (USTM), Baridua, India
- 5Department of Botany, Gauhati University, Guwahati, India
- 6Department of Biotechnology, Yeungnam University, Gyeongsan, South Korea
- 7Natural and Medical Sciences Research Centre, University of Nizwa, Nizwa, Oman
- 8Institute of Bioresources and Sustainable Development, Imphal, India
SARS-CoV-2 is the virus responsible for causing COVID-19 disease in humans, creating the recent pandemic across the world, where lower production of Type I Interferon (IFN-I) is associated with the deadly form of the disease. Membrane protein or SARS-CoV-2 M proteins are known to be the major reason behind the lower production of human IFN-I by suppressing the expression of IFNβ and Interferon Stimulated Genes. In this study, 7,832 compounds from 32 medicinal plants of India possessing traditional knowledge linkage with pneumonia-like disease treatment, were screened against the Homology-Modelled structure of SARS-CoV-2 M protein with the objective of identifying some active phytochemicals as inhibitors. The entire study was carried out using different modules of Schrodinger Suite 2020-3. During the docking of the phytochemicals against the SARS-CoV-2 M protein, a compound, ZIN1722 from Zingiber officinale showed the best binding affinity with the receptor with a Glide Docking Score of −5.752 and Glide gscore of −5.789. In order to study the binding stability, the complex between the SARS-CoV-2 M protein and ZIN1722 was subjected to 50 ns Molecular Dynamics simulation using Desmond module of Schrodinger suite 2020-3, during which the receptor-ligand complex showed substantial stability after 32 ns of MD Simulation. The molecule ZIN1722 also showed promising results during ADME-Tox analysis performed using Swiss ADME and pkCSM. With all the findings of this extensive computational study, the compound ZIN1722 is proposed as a potential inhibitor to the SARS-CoV-2 M protein, which may subsequently prevent the immunosuppression mechanism in the human body during the SARS-CoV-2 virus infection. Further studies based on this work would pave the way towards the identification of an effective therapeutic regime for the treatment and management of SARS-CoV-2 infection in a precise and sustainable manner.
Introduction
Research Question
Owning to the massive global threat posed by Covid-19, discovery of novel and effective therapeutic strategies against the disease becomes the need of the hour (Nand et al., 2020). The present study was designed to computationally identify phytocompound(s) against Membrane Glycoprotein (M protein) of SARS-CoV-2 with a view to targeting the antagonistic effect M protein exerts on Type I Interferon induced anti-viral immunity in the host cell.
The Disease and Its Causative
Since the first recorded case in November 2019 in Wuhan, China, Covid-19 has gone on to infect citizens in all the continents (Puttaswamy et al., 2020; Elekofehinti et al., 2021). The National Health Commission of China went on to confirm the human-to-human transmission of the disease in late January 2020, following which on 31 January, 2020, WHO pronounced COVID-19 as a ‘public health emergency of international concern (11 March 2020)’, and subsequently as a pandemic in March 2020 (Aouidate et al., 2021; Elfiky, 2021; Sridhar et al., 2021). As of 21 October, 2021, there have been a total of 241,886,635 recorded Covid-19 cases and 4,919,755 reported deaths (WHO Coronavirus (COVID-19) Dashboard | WHO Coronavirus (COVID-19) Dashboard with Vaccination Data). Owing to its virulence and high rate of infection, the Covid-19 virus, SARS-CoV-2, has caused large-scale damage, leading to global economic losses (Andrianov et al., 2020; Jiménez-Alberto et al., 2020; Saxena, 2020; Elekofehinti et al., 2021).
Belonging to the RNA virus family, Coronaviridae, SARS-CoV-2 is a Betacoronavirus (genus) like the member coronaviruses, SARS and MERS (Aouidate et al., 2021; Elfiky, 2021). Covid-19, the recently identified Coronavirus disease caused by SARS-CoV-2, is way more fatal to humans compared to the previously identified strains (Mohanta, 2020; Puttaswamy et al., 2020).
SARS-CoV-2 enters the host cells with the help of ACE2, Angiotensin-Converting Enzyme 2, which is found attached to the cell membrane of several organs in the human body. SARS-CoV-2 endocytoses its way into a host cell via means of interacting with the ACE2 receptor using its unique surface glycoprotein, spike protein. The virus then goes on to hijack the replication machinery of the host cell and produces a new set of viral particles ready to infect a new set of cells (Elekofehinti et al., 2021).
SARS-CoV-2 primarily affects the lungs. Additionally, the virus affects other organs such as intestine, kidneys, blood vessels among several others (Palmeira et al., 2020; Elekofehinti et al., 2021).
Infected people display a wide range of symptoms varying from person-to-person including fever, dyspnea, dry cough, loss of appetite, vomiting, loss of sense of taste and smell, which in certain cases progresses to pneumonia leading to acute respiratory distress syndrome (ARDS), multi-organ failure, usually followed by death (Gimeno et al., 2020; Puttaswamy et al., 2020; Elekofehinti et al., 2021). Aggressive inflammatory response to SARS-CoV-2 infection may result in Cytokine Storm which mediates widespread lung infection and is often followed by septic shock and multi-organ failure, making Cytokine Storm one of the major causes of fatality among the patients (Mohanta T. K. et al., 2020; Tay et al., 2020).
Coronaviruses classify as an extremely diverse group of enveloped viruses having positive-sense, single-stranded RNA (ssRNA), possessing a genome bearing 27–32 kilobases (Khan et al., 2020; Shamsi et al., 2020). The genome contains at least 6 ORFs which code for a minimum of 4 structural proteins and 16 non-structural proteins (Jiménez-Alberto et al., 2020).
The four structural proteins include Nucleocapsid (N) protein, Envelope (E) protein, Membrane (M) protein, and Spike (S) protein (Sridhar et al., 2021). These four proteins together form the envelope of SARS-CoV-2 and are responsible primarily for the assembly of the virus and subsequently budding. These proteins also modulate the host’s cellular response to the infection caused by the virus (Elekofehinti et al., 2021). Specifically, the N protein is responsible for the processing of the viral genome, and the S proteins are the components that polymerize to give coronaviruses its popular crown-like shape in addition to being the primary component facilitating the virus entry into the host cell by the process of receptor-mediated endocytosis (Nand et al., 2020; Elekofehinti et al., 2021). Further, the M protein, through its transmembrane domain, binds to the nucleocapsid of the virus and provides a shape to the virion. The E protein plays a crucial role in viral pathogenesis and the virus assembly and budding. Whereas, the 16 non-structural proteins essentially play crucial roles in the virus replication machinery (Khan et al., 2020).
The genome of SARS-CoV-2 produces a ∼800 kDa polypeptide as its transcription product which is proteolytically cleaved by enzymes, 3-Chymotrypsin-like protease (3CLpro) and Papain-like protease (PLpro). Cleavage by 3CLpro essentially leads to the formation of the various non-structural proteins mentioned prior. Hence, 3CLpro was identified as a crucial protein (Nand et al., 2020). Owing to their roles, these proteins were identified as drug targets, however, among these, S protein and 3CLpro, identify as the most popular drug targets (Nand et al., 2020).
Therapeutics and Its Status
The dramatic unfolding of the pandemic has drawn global concern towards the development of therapeutic measures against the virus (Morgon et al., 2021). Various researchers worldwide have adopted various strategies towards finding an effective drug molecule. Most approaches include the use of data from preliminary in vitro studies and clinical trials, computational techniques, high-throughput screening and drug-repurposing (Puttaswamy et al., 2020; Wu et al., 2020). Using these approaches, several studies had suggested potential drug candidates, which include now popular names like Remdesivir, Oseltamivir, Tocilizumab, Ribavirin, Chloroquine, Hydroxychloroquine, among others (Gimeno et al., 2020; Al-Sehemi et al., 2021; Tazikeh-Lemeski et al., 2021). These are mostly antivirals, corticosteroids, antibiotics and immunoglobulins (Morgon et al., 2021). There are, however, several issues around these compounds. The efficacy of most of these remains controversial (Gimeno et al., 2020). Several of the repurposed drugs displayed potent results in in vitro studies but failed to show encouraging results in clinical trials (Gogoi et al., 2021). Various others have failed to specifically interact with the viral protease (Elekofehinti et al., 2021). All in all, the absence of a safe, effective, potent and definitive drug molecule so far greatly emphasizes need to find newer potential candidates (Palmeira et al., 2020; Gogoi et al., 2021; Morgon et al., 2021; Tazikeh-Lemeski et al., 2021).
Even with vaccines out now, it’ll be 2 years before citizens of the poorest countries gain access to their doses (Padma, 2021). Additionally, the dramatically evolving variants of the virus resulting from mutations and ineffectiveness of the proposed molecules towards these variants (Andrianov et al., 2020), point to the need for identifying novel therapeutic approaches.
Membrane Protein
M protein of SARS-CoV-2 plays a crucial role in the virus assembly (Bhowmik et al., 2020). M protein is inherently a glycoprotein that is generally conserved across the genus of Betacoronaviruses, especially across SARS CoVs (Bianchi et al., 2020). Despite the crucial role played by membrane proteins, few drugs are being developed targeting them (Ita, 2021). To mediate membrane fusion, SARS-CoV-2 spike glycoprotein must be sequentially cleaved at its S1/S2 and S2’ cleavage sites. Furin can cleave the polybasic insertion (PRRAR) at the S1/S2 cleavage site of SARS-CoV-2. (Peacock et al., 2021).
Innate Immunity, which is the first line of defence in the host against viral organisms, is triggered via PAMPs (Pathogen-Associated Molecular Patterns), which also includes ssRNAs, leading to the production of Interferons (IFNs): Type I (α, β) and Type III. These Type I and Type III IFNs, through a series of signalling cascades, go on to facilitate the expression of Interferon Stimulated Genes (ISGs) which confer antiviral immunity to the host cell (Zheng et al., 2020). In such a scenario, encoding of IFN antagonists by the viruses appears to be a common strategy adopted by the viruses to evade the host’s antiviral immunity. In recent screening studies, it was found that the M protein of SARS-CoV-2, inhibited SeV (Sendai virus)-induced IFN-β gene activation. The SARS-CoV-2 infection causes declined production of Type I and Type III IFNs, coupled with lowered ISG response. Whereas, in severe cases of infection, this expression is strongly triggered throughout the course of the infection (Zheng et al., 2020). Nonetheless, it has been established that SARS-CoV-2 M protein is a Type I and Type III Interferon-antagonist. It is to note that this antagonism is achieved by negatively influencing the formation of a multiprotein complex, the RIG-I/MDA-5–MAVS–TRAF3–TBK1 signalosome and by the impairing of MAVS (Mitochondrial Antiviral Signaling Protein) activation (Zheng et al., 2020; Fu et al., 2021).
The Type I IFNs (IFN-Is) confer strong antiviral immunity to the host cells, hence the antagonistic strategies adopted by SARS-CoV-2 to evade the host immunity has been an intense matter of research. Additionally, it was found that IFN-Is facilitates the severity of the Cytokine Storms, which is a notably major complication observed in Covid-19, often leading to ARDS, followed by death (Mohanta TK. et al., 2020) (Schreiber, 2020).
Significance of Computational Studies in Covid-19 Research
After identification of a potential drug target, the development of a drug molecule right from identification of potential ligands to procurement of a marketable drug, needless to say, is an extremely time-consuming and expensive process (Al-Sehemi et al., 2020; Tazikeh-Lemeski et al., 2021) which in case of the deadly SARS-CoV-2, is not the most feasible approach (Sridhar et al., 2021). In silico methods like Structure-Based Drug Design (SBDD) or Computer-Aided Drug Discovery (CADD) techniques prove to be extremely rapid, cost-effective, robust and reliable alternatives (Al-Sehemi et al., 2020, 2021; Gimeno et al., 2020; Tatar et al., 2021). CADD techniques have been vigorously used by researchers for the identification of lead molecules for various diseases so far, including Covid-19, and presented encouraging results (Gentile et al., 2020; Gimeno et al., 2020; Naik et al., 2020; Aouidate et al., 2021).
Significance of Compounds From Medicinal Plants
Owing to the numerous side effects, including drug resistance, posed by synthetic drugs, recent drug discovery has seen a shift towards natural products as sources of drug molecules (Gogoi et al., 2021). Plants are rich sources of pharmacologically active or bioactive compounds possessing huge structural diversity (Naik et al., 2020; Puttaswamy et al., 2020; Gogoi et al., 2021). These bioactive small molecules are the plant’s secondary metabolites (Puttaswamy et al., 2020). Extensive in silico research has been conducted since the outbreak towards identifying phytocompounds possessing anti-Covid-19 properties. Some of the proposed candidates include curcumin, epi-catechin-gallate, kaempferol, catechin, quercetin, garlic essential oil, compounds from Houttuynia cordata, among several others (Puttaswamy et al., 2020; Thuy et al., 2020; Das et al., 2021).
Summary of the Current Study
In the present study, 7,832 compounds from 32 medicinal plants of India possessing traditionally established linkage with pneumonia-like disease treatment, were curated. These phytocompounds were subsequently virtually screened, using CADD protocols, against a Homology Modelled structure of SARS-CoV-2 M protein. The ligands were subjected to Molecular Docking and subsequently to Molecular Dynamics simulation of 50 nanoseconds (ns), both of which were carried out in Schrödinger Suite 2020-3. The finally shortlisted ligands were evaluated for their ADMET properties. Following the experiments, the ligand, ZIN1722, from Zingiber officinale, emerged as the best binding ligand, indicating its potential anti-Covid-19 property.
Materials and Methods
Prediction of the 3D Structure of SARS-CoV-2 M Protein
Due to the absence of the structure of the SARS-CoV-2 M protein in Protein Data Bank, the 3D-structure of the protein was determined via Homology Modelling (Venselaar et al., 2010).
The sequence for the protein, encoded by gene M, was retrieved from UniProt SwissProt database (SwissProt accession number P0DTC5.1) (Boeckmann et al., 2003). The 222-amino acid-long protein was subjected to Homology Modelling using homology/ab initio hybrid online server BhageerathH+ (Jayaram et al., 2014).
The analysis of the quality of the predicted structure of M protein was performed using ProSA analysis (Wiederstein and Sippl, 2007) and subsequently, the Ramachandran plot was generated for the predicted structure (Ho et al., 2009) through ProCheck (Roman Laskowski et al., 1983).
Preparation of the Receptor Target: SARS-CoV-2 M Protein
The predicted 3D structure of the SARS-CoV-2 M protein was processed and prepared for Molecular Docking in Protein Preparation Wizard of Epik module of Schrödinger Suite, 2020-3 (Paul et al., 2021). The similar binding sites and excess water molecules in the protein structure were removed, and the bond orders were refined. Using the Prime module of Schrodinger Suite, 2020-3, the missing chain atoms were added (Suganya et al., 2014). Further, the 3D structure of the protein was minimized using Optimized Potentials For Liquid Simulations-3 (OPLS3e) molecular force field with RMSD (Root Mean Square Deviation) of crystallographic heavy atoms which were kept at 0.3Å (Paul et al., 2021).
The structure was further subjected to active site analysis on AADS which is an Automated Active Site Identification server developed by the Supercomputing Facility for Bioinformatics and Computational Biology, IIT Delhi (Singh et al., 2011).
Selection and Preparation of the Ligands
A ligand library of 7,832 compounds from 32 traditionally used medicinal plants of India was curated. To prepare the ligands for Molecular Docking and further studies, a series of activities were executed. The two-dimensional structures of the ligands were converted to their respective 3D counterparts using the LigPrep module of the Schrödinger Suite (Schrödinger, 2017). Further, to screen out the unfavourable molecules from the library, Lipinski’s rule of Drug Discovery was used as the criteria (Lipinski, 2004). Next, energy minimization and geometry optimization were performed for the 3D structures of the ligands, followed by desaltation and correction of chirality. The tautomeric and ionization states were generated between pH 6.8–7.2 using Epik module of Schrödinger suite 2020-3. The ligand structures were minimized using OPLS3e force field in Schrödinger Suite 2020-3 until a Root Mean Square Deviation (RMSD) value of 2.0Å was attained. The optimized ligands were docked against the selected receptor proteins.
Molecular Docking
Molecular Docking was performed between the receptor and the optimized ligands in Glide module of Schrödinger Suite 2020-3 using the High Throughput Virtual Screening (HTVS) mode. Taking into consideration the Glide Docking Score, fifty ligands were shortlisted and further subjected to Extra Precision (XP) docking option of Glide (Wilmes et al., 2004; Schrödinger, 2017). Prior to performing molecular docking, grid box of the receptor protein structure was defined on the basis of the information of the predicted active site.
Molecular Dynamics Simulation
The stability of the interactions within the docked complex constituted by SARS-CoV-2 M protein and the best-binding ligand was analyzed by performing Molecular Dynamic simulation of 50 ns (Schrödinger, 2017). The SARS-CoV-2 M Protein-ligand complex within the explicit solvent system with the OPLS3e force field was studied using the Desmond module of Schrödinger suite 2020-3. The atomic framework had been solvated with TIP3P crystallographic water particles with orthorhombic intermittent limit conditions for a 10 Å buffer region. The overlapping water molecules were eliminated and Na+ were added as counter ions to neutralize the entire framework of atoms. An ensemble (NPT) of Nose-Hoover thermostat and barostat was applied to maintain the constant temperature of 300 K and pressure of 1 bar in the system (Schrödinger, 2017). A hybrid energy minimization algorithm with 1,000 steps of steepest descent followed by conjugate gradient algorithms was used.
ADME-Tox Analysis
To assess the safety and pharmacological competency of the shortlisted ligand, that is, to assess the absorption, distribution, metabolism, excretion and toxicity of the shortlisted ligand, ADME-Tox test of the ligand was performed. The ADME-Tox analysis was performed using two online tools, SwissADME and pkCSM, to encompass a pool of pharmacokinetic and medicinal chemistry properties (Pires et al., 2015; Daina et al., 2017).
• Swiss ADME
SwissADME is a freely-available, user-friendly online tool which hosts numerous fast and robust predictive models for the evaluation of pharmacokinetics, physicochemical properties, drug-likeliness and medicinal chemistry friendliness.
• pkCSM
pkCSM is an integrated, freely-available, user-friendly web server that utilises graph-based signatures to perform ADMET studies. These signatures encode distance patterns between the atoms which are in turn used to represent the small molecule as well as train the predictive models.
Results
Molecular Docking
The homology modelled-3D-structure of SARS-CoV-2 M protein was docked against the selected ligands using the High Throughput Virtual Screening (HTVS) and subsequently, Extra Precision (XP) modes of the Glide docking module of the Schrödinger Suite, version 2020-3. The docking score and other relevant information regarding the best twenty (20) docked complexes have been presented in Table 1. The results of the docking experiment presented a list of the best-binding ligands with respect to the receptor. With a Glide gscore of −5.789, ZIN1722, a compound from the root of Zingiber officinale (Ginger), emerged as the best-binding ligand.
TABLE 1. Docking score and ranking of the compounds evaluated through Glide Docking module.
Various interactions were observed between the compound, ZIN1722, and the structure of M protein. A total of four hydrogen bonds were noted. Glu115 and His154 of both the chains of the receptor have formed hydrogen bonds with the ligand. The interactions have been pictorially represented in Figures 1A,B. The compound, ZIN1722, hence, was considered for further study based on its best in silico affinity towards the receptor among all the ligands.
FIGURE 1. (A) Surface view of the interactions between SARS-CoV-2 M Protein bound to the ligand ZIN1722; (B) Detailed interaction diagram of SARS-CoV-2 M Protein with ZIN 1722. [before MD simulation].
Molecular Dynamics Simulation
To assess the stability of the interactions between ZIN1722 and the M protein structure, the docked complex was studied under Molecular Dynamics Simulation of 50 ns in Desmond software, Schrödinger Suite version 2020-3. The results of the Molecular Dynamics simulation study have been presented in Figures 2–4. In Figure 2A, the left Y-axis represents the RMSD evolution of the proteins, that is, the structural conformation of the proteins throughout the simulation and the right Y-axis represents the stability of the ligand with respect to the protein and the protein’s binding pocket. The Root Mean Square Deviation (RMSD) of the ligand, ZIN1722, bound to SARS-CoV-2 M Protein, was found to be stable after 32 ns of the simulation. Whereas, the RMSD of the receptor protein, bound to the ligand, was almost stable throughout the entire 50 ns of the simulation which signifies that the structural conformation of the receptor protein in the ligand-bound state was quite stable. The Ligand RMSD, Intramolecular Hydrogen Bonds, Radius of Gyration, Molecular Surface Area, Polar Surface Area and Solvent Accessible Surface Area of ZIN1722 in complex with the receptor protein was studied further during the 50 ns-long simulation and the resultant images are presented in Figure 3.
FIGURE 2. (A) The RMSD plot for SARS-CoV-2 M protein and ZIN1722 complex, where the RMSDs of Protein backbone and the Ligands were calculated separately throughout the MD trajectory of 50 ns; (B) RMSF plot of the Protein chain in Ligand bound state.
FIGURE 3. Ligand RMSD, Radius of Gyration (rGyr), Intramolecular Hydrogen Bond (intraHB), Molecular Surface Area (MolSA), Solvent Accessible Surface Area (SASA), Polar Surface Area (PSA) as calculated during the 50 ns of MD Simulation.
In the RMSF (Root Mean Square Fluctuations) local charges throughout the protein chain were characterized. The RMSF plots of the receptors bound to the selected ligands are placed in Figure 2B. The peaks in Figure 2B indicate the protein residues that fluctuated the most during the entire 50 ns-long simulation; the green-coloured vertical bars indicate the residues of the protein that interacted with the ligand. The fluctuations around amino acids no 25, 50, 100, 125 and 175 were recorded along with the N and C terminal amino acids as expected. As a rule, the N- and C-terminal regions of proteins fluctuate more than the rest of the residues, similarly, the structured regions like alpha helices and beta strands fluctuate less than the loop regions. Throughout the trajectory, mostly water bridges, Hydrogen bonds and Hydrophobic interactions were recorded between the receptor and the ligand. The amino acids Glu115, His154, Leu135 and Asn 113 had shown Hydrogen bonding with different atoms of the ligand ZIN1722. A hydrogen bond with Glu115 was formed and maintained over more than 80% of the total duration of the simulation.
A summary of the interactions between the receptor and the ligand throughout the MD simulation has been presented in Figures 4A,B.
FIGURE 4. (A) Interactions between the SARS-CoV-2 M protein and the Ligand throughout the Simulation is depicted in bars. Different color signifies different interactions and contacts formed between the amino acids and the ligand atoms. (B) Graphical representation of the interactions and contacts throughout the Simulation.
ADME-Tox Analysis
The pharmacokinetic properties and toxicity status of ZIN1722 were assessed using pkCSM and SWISS ADME web servers. The compound passed the Lipinski rule of five and Ghose filter which indicates that it may be a good lead compound. The absorption value of ZIN1722 predicted by pkCSM server is 65.4% whereas, the volume of distribution is 1.232 (log L/kg). Total clearance, which is an indicator of excretion value, is predicted to 0.526 (log ml/min/kg) whereas, oral rat acute toxicity (Rat LD50 for toxicity) is predicted to be 2.437 (mol/kg), which is in a quite acceptable range. The detailed report of the ADME-Tox study performed on both the aforementioned servers has been presented in Table 2 and Table 3. All in all, both the in silico ADME-Tox tests present encouraging results. Hence, the molecule, ZIN1722, may be considered for in vitro validation and subsequent studies.
TABLE 2. Summary of the ADME-Tox results performed using SWISS-ADME server.
TABLE 3. Summary of the ADME-Tox study performed using pkCSM server.
The phytocompound, ZIN1722, bears the common name Myricetin. Its IUPAC name is 5,7-Dihydroxy-4-Oxo-2-(3,4,5-Trihydroxyphenyl)-4H-Chromen-3-Olate. It is obtained from the root of Zingiber officinale, commonly known as Ginger. The PubChem Compound ID (CID) is 25201643 (Kim et al., 2021).
Discussion
Ever since its outbreak in 2019, Covid-19 has caused massive global damage in all spheres. Needless to say, in the absence of a single safe, effective and definite treatment measure and the rapidly mutating strains of the virus, SARS-CoV-2, the discovery of new and potent drug candidates as well as novel therapeutic approaches, become the need of the hour (Andrianov et al., 2020; Nand et al., 2020; Palmeira et al., 2020; Gogoi et al., 2021; Morgon et al., 2021; Padma, 2021; Tazikeh-Lemeski et al., 2021). The present study was designed with a view to targeting the less popular drug target of SARS-CoV-2, Membrane Glycoprotein, also known as M Protein. M protein plays a critical role in various functions of the virus, mainly in virus assembly (Bhowmik et al., 2020). However, this study was designed targeting another functionally significant role of M protein in infecting a host. M protein facilitates evasion of the host’s immune system by means of rendering a declined production of Type I and Type III Interferon by suppressing the expression of Interferon Stimulated Genes (ISGs) (Schreiber, 2020; Zheng et al., 2020; Fu et al., 2021). Type I IFN confers strong anti-viral immunity to the host, hence, inhibiting M proteins which inhibit the production of Type I IFN, may hinder disease progression in the infected individuals. In the present study, 7,832 phytocompounds from 32 medicinal plants of India were virtually screened, resulting in the identification of compound ZIN1722 as the compound with the best binding affinity towards the receptor protein, M protein, whose structure was predicted through Homology Modelling. To study the stability of the interactions formed between ZIN1722 and the M protein, the complex was subjected to a 50 ns-long Molecular Dynamics Simulation. Throughout the trajectory, it was found that the receptor-ligand complex displayed considerable stability and hence, it can be concluded that the ligand, ZIN1722, can be considered as a compound possessing encouraging affinity towards SARS-CoV-2 M protein. Further, the safety and pharmacological competency of the phytocompound was assessed through extensive pharmacokinetic profiling performed on SwissADME and pkCSM webservers. The compound passed the Lipinski rule of five and Ghose filter, and moderate-to-excellent values of most of the pharmacokinetic parameters were assessed. Overall, results of the in silico ADME-Tox analysis indicate encouraging scope of the compound, ZIN1722, as a drug candidate. Based on these findings, ZIN1722 may be considered for further in vitro and in vivo validation as a novel anti-SARS-CoV-2 compound (Daina et al., 2017).
The phytocompound, ZIN1722, also known as Myricetin [PubChem Compound ID (CID) is 25201643; IUPAC name: 5,7-Dihydroxy-4-Oxo-2-(3,4,5-Trihydroxyphenyl)-4H-Chromen-3-Olate] is a flavonoid found mainly in the roots of Ginger (Zingiber officinale) (Yu et al., 2012; Song et al., 2021). Myricetin is also found in plants of many families and distributed in many foods like berries, fruits, vegetables, tea, etc.
From the current study, it was found that ZIN1722 or Myricetin has a considerable affinity towards SARS-CoV-2 M protein which is sufficient for inhibiting the activity of M protein which may disrupt the process of evasion of host antiviral immunity by means of inhibiting the declined production of Type I IFN and subsequently, lowered expression of ISGs. Inhibition of M protein is also expected to disrupt the process of viral assembly. Conclusively, inhibition of SARS-CoV-2 M protein by Myricetin is expected to mitigate viral infection in the infected individual, thereby, opening up novel therapeutic approaches to tackling the Covid-19 virus.
Myricetin has been previously reported to possess anti-viral, anti-bacterial, anti-tumour, anti-inflammatory activities. The compound is also known to be protective against cardiovascular diseases, neurological damage and potential liver injury (Semwal et al., 2016; Song et al., 2021).
In previous studies, Myricetin has also been reported to show high binding affinity towards ACE2 receptor, indicating the compound’s potential inhibitory activity towards viral entry into host cell during SARS-CoV-2 infection. Additionally, Myricetin has also been reported to possess inhibitory activity towards SARS-CoV helicase which supports the potential of this molecule in combating diseases caused by coronaviruses, including Covid-19 (Yu et al., 2012).
The findings of the current study are encouraging and align with previous research findings. The compound, Myricetin, identified through this in silico study as a potent anti-SARS-CoV-2 agent, possesses hopeful drug-like properties which can be further validated through in vitro and in vivo assays, and subsequent clinical studies, thereby, paving way for the identification of a potential therapeutic agent against Covid-19.
Conclusion
Two years into the Covid-19 pandemic with countless fatalities, damage of varied kinds and degrees worldwide, no definitive drug molecule or treatment strategy, quest for an effective, potent and novel therapeutic approach against the disease becomes the need of the hour. The current study is a computational approach towards identifying a novel molecule against the causative of Covid-19, SARS-CoV-2. Digressing from the popularly studied receptor targets, the current study focuses on the immunological perspective of the infection, involving the SARS-CoV-2 Membrane Glycoprotein, M protein. In addition to playing a role in virus assembly, SARS-CoV-2 M protein has been reported to cause a declined production of Type I Interferon and a downregulated expression of the Interferon Stimulated Genes (ISGs), thereby facilitating evasion of the virus by the host innate immunity and subsequently, rendering a decreased anti-viral response of the host immune system towards the infection. Keeping in view this recent, valuable finding, this study was designed to identify inhibitory molecule(s) towards the SARS-CoV-2 M protein. In this study, 7,832 compounds from 32 medicinal plants of India traditionally known to be used in the treatment of pneumonia-like diseases, were virtually screened against a Homology Modelled structure of SARS-CoV-2 M protein. The best-binding ligand from the Molecular Docking experiment was further studied for its stability with the receptor through Molecular Dynamics simulation study, and for its pharmacokinetic profile through ADME-Tox study. The phytocompound, ZIN1722, that is, Myricetin was identified from this study as a potent SARS-CoV-2 M protein inhibitor. The findings of this study align with the findings of previously reported studies. The findings of this study are preliminary but highly encouraging, hence, further evaluation of the potency and efficacy of the phytocompound through in vitro and in vivo validation assays, and subsequent clinical studies, present hopeful scope towards the discovery of a novel therapeutic molecule against SARS-CoV-2 in the future (WHO Coronavirus Dashboard).
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
SM: conceived the idea and run the experiments. LC, SM, and TN drafted the original manuscript. AM, BT and KB: revised the manuscript. YM, TM, and NS: analyzed the data.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors would like to extend their sincere thanks to their respective host institutes for providing the necessary support to conduct the work.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2022.805344/full#supplementary-material
References
Al-Sehemi, A. G., Pannipara, M., Parulekar, R. S., Patil, O., Choudhari, P. B., Bhatia, M. S., et al. (2020). Potential of NO Donor Furoxan as SARS-CoV-2 Main Protease (Mpro) Inhibitors: In Silico Analysis. J. Biomol. Struct. Dyn. 1, 1–15. doi:10.1080/07391102.2020.1790038
Al-Sehemi, A. G., Parulekar, R. S., Pannipara, M., Manzur Ali, P P, Zubaidha, P. K., Bhatia, M. S., et al. (2021). In Silico evaluation of NO Donor Heterocyclic Vasodilators as SARS-CoV-2 Mpro Protein Inhibitor. J. Biomol. Struct. Dyn. 1, 1–18. doi:10.1080/07391102.2021.2005682
Andrianov, A. M., Kornoushenko, Y. V., Karpenko, A. D., Bosko, I. P., and Tuzikov, A. V. (2020). Computational Discovery of Small Drug-like Compounds as Potential Inhibitors of SARS-CoV-2 Main Protease. J. Biomol. Struct. Dyn. 39, 1–13. doi:10.1080/07391102.2020.1792989
Aouidate, A., Ghaleb, A., Chtita, S., Aarjane, M., Ousaa, A., Maghat, H., et al. (2021). Identification of a Novel Dual-Target Scaffold for 3CLpro and RdRp Proteins of SARS-CoV-2 Using 3D-Similarity Search, Molecular Docking, Molecular Dynamics and ADMET Evaluation. J. Biomol. Struct. Dyn. 39, 4522–4535. doi:10.1080/07391102.2020.1779130
Bhowmik, D., Nandi, R., Jagadeesan, R., Kumar, N., Prakash, A., and Kumar, D. (2020). Identification of Potential Inhibitors against SARS-CoV-2 by Targeting Proteins Responsible for Envelope Formation and Virion Assembly Using Docking Based Virtual Screening, and Pharmacokinetics Approaches. Infect. Genet. Evol. 84, 104451. doi:10.1016/j.meegid.2020.104451
Bianchi, M., Benvenuto, D., Giovanetti, M., Angeletti, S., Ciccozzi, M., and Pascarella, S. (2020). “Sars-CoV-2 Envelope and Membrane Proteins: Differences from Closely Related Proteins Linked to Cross-Species Transmission?,” in Molecular Surveillance of Emerging Pathogens Using Bioinformatic and Bayesian Phylogeny View Project. doi:10.20944/preprints202004.0089.v1
Boeckmann, B., Bairoch, A., Apweiler, R., Blatter, M. C., Estreicher, A., Gasteiger, E., et al. (2003). The SWISS-PROT Protein Knowledgebase and its Supplement TrEMBL in 2003. Nucleic Acids Res. 31, 365–370. doi:10.1093/nar/gkg095
Daina, A., Michielin, O., and Zoete, V. (2017). SwissADME: A Free Web Tool to Evaluate Pharmacokinetics, Drug-Likeness and Medicinal Chemistry Friendliness of Small Molecules. Sci. Rep. 7, 42717. doi:10.1038/srep42717
Das, S. K., Mahanta, S., Tanti, B., Tag, H., and Hui, P. K. (2021). Identification of Phytocompounds from Houttuynia Cordata Thunb. As Potential Inhibitors for SARS-CoV-2 Replication Proteins through GC-MS/LC-MS Characterization, Molecular Docking and Molecular Dynamics Simulation. Mol. Divers. 26, 365–388. doi:10.1007/s11030-021-10226-2
Elekofehinti, O. O., Iwaloye, O., Josiah, S. S., Lawal, A. O., Akinjiyan, M. O., and Ariyo, E. O. (2021). Molecular Docking Studies, Molecular Dynamics and ADME/tox Reveal Therapeutic Potentials of STOCK1N-69160 against Papain-like Protease of SARS-CoV-2. Mol. Divers. 25, 1761–1773. doi:10.1007/s11030-020-10151-w
Elfiky, A. A. (2021). SARS-CoV-2 RNA Dependent RNA Polymerase (RdRp) Targeting: an In Silico Perspective. J. Biomol. Struct. Dyn. 39, 3204–3212. doi:10.1080/07391102.2020.1761882
Fu, Y. Z., Wang, S. Y., Zheng, Z. Q., Yi Huang, Yi., Li, W. W., Xu, Z. S., et al. (2021). SARS-CoV-2 Membrane Glycoprotein M Antagonizes the MAVS-Mediated Innate Antiviral Response. Cell Mol Immunol 18, 613–620. doi:10.1038/s41423-020-00571-x
Gentile, D., Patamia, V., Scala, A., Sciortino, M. T., Piperno, A., and Rescifina, A. (2020). Putative Inhibitors of SARS-COV-2 Main Protease from a Library of marine Natural Products: A Virtual Screening and Molecular Modeling Study. Mar. Drugs 18, 225. doi:10.3390/md18040225
Gimeno, A., Mestres-Truyol, J., Ojeda-Montes, M. J., Macip, G., Saldivar-Espinoza, B., Cereto-Massagué, A., et al. (2020). Prediction of Novel Inhibitors of the Main Protease (M-Pro) of SARS-CoV-2 through Consensus Docking and Drug Reposition. Int. J. Mol. Sci. 21, 3793. doi:10.3390/ijms21113793
Gogoi, M., Borkotoky, M., Borchetia, S., Chowdhury, P., Mahanta, S., and Barooah, A. K. (2021). Black tea Bioactives as Inhibitors of Multiple Targets of SARS-CoV-2 (3CLpro, PLpro and RdRp): a Virtual Screening and Molecular Dynamic Simulation Study. J. Biomol. Struct. Dyn., 1–24. doi:10.1080/07391102.2021.1897679
Ho, B. K., Thomas, A., and Brasseur, R. (2009). Revisiting the Ramachandran Plot: Hard-Sphere Repulsion, Electrostatics, and H-Bonding in the Alpha-helix. Protein Sci. 12, 2508–2522. doi:10.1110/ps.03235203
Ita, K. (2021). Coronavirus Disease (COVID-19): Current Status and Prospects for Drug and Vaccine Development. Arch. Med. Res. 52, 15–24. doi:10.1016/j.arcmed.2020.09.010
Jayaram, B., Dhingra, P., Mishra, A., Kaushik, R., Mukherjee, G., Singh, A., et al. (2014). Bhageerath-H: A Homology/Ab Initio Hybrid Server for Predicting Tertiary Structures of Monomeric Soluble Proteins. BMC Bioinformatics 15, S7. doi:10.1186/1471-2105-15-S16-S7
Jiménez-Alberto, A., Ribas-Aparicio, R. M., Aparicio-Ozores, G., and Castelán-Vega, J. A. (2020). Virtual Screening of Approved Drugs as Potential SARS-CoV-2 Main Protease Inhibitors. Comput. Biol. Chem. 88, 107325. doi:10.1016/j.compbiolchem.2020.107325
Khan, A., Khan, M., Saleem, S., Babar, Z., Ali, A., Khan, A. A., et al. (2020). Phylogenetic Analysis and Structural Perspectives of RNA-dependent RNA-Polymerase Inhibition from SARs-CoV-2 with Natural Products. Interdiscip. Sci. 12, 335–348. doi:10.1007/s12539-020-00381-9
Kim, S., Chen, J., Cheng, T., Gindulyte, A., He, J., He, S., et al. (2021). PubChem in 2021: New Data Content and Improved Web Interfaces. Nucleic Acids Res. 49, D1388–D1395. doi:10.1093/nar/gkaa971
Lipinski, C. A. (2004). Lead- and Drug-like Compounds: The Rule-Of-Five Revolution. Drug Discov. Today Technol. 1, 337–341. doi:10.1016/j.ddtec.2004.11.007
Mohanta, T. K. (2020). Corona Virus (CoVid19) Genome: Genomic and Biochemical Analysis Revealed its Possible Synthetic Origin. J. Appl. Biotechnol. Bioeng. 7, 200–213. doi:10.15406/jabb.2020.07.00235
Mohanta, T. K., Sharma, N., Arina, P., and Defilippi, P. (2020a). Molecular Insights into the MAPK Cascade during Viral Infection: Potential Crosstalk between HCQ and HCQ Analogues. Biomed. Res. Int. 2020, 1–9. doi:10.1155/2020/8827752
Mohanta, T. K., Sharma, N., Arina, P., and Defilippi, P. (2020b). Molecular Insights into the MAPK Cascade during Viral Infection: Potential Crosstalk between HCQ and HCQ Analogues. Biomed. Res. Int. 2020, 8827752–8827759. doi:10.1155/2020/8827752
Morgon, N. H., Grandini, G. S., Yoguim, M. I., Porto, C. M., Santana, L. C., Biswas, S., et al. (2021). Potential Activity of Linezolid against SARS-CoV-2 Using Electronic and Molecular Docking Study. J. Mol. Model. 27, 222. doi:10.1007/s00894-021-04828-8
Naik, B., Gupta, N., Ojha, R., Singh, S., Prajapati, V. K., and Prusty, D. (2020). High Throughput Virtual Screening Reveals SARS-CoV-2 Multi-Target Binding Natural Compounds to lead Instant Therapy for COVID-19 Treatment. Int. J. Biol. Macromol 160, 1–17. doi:10.1016/j.ijbiomac.2020.05.184
Nand, M., Maiti, P., Joshi, T., Chandra, S., Pande, V., Kuniyal, J. C., et al. (2020). Virtual Screening of Anti-HIV1 Compounds against SARS-CoV-2: Machine Learning Modeling, Chemoinformatics and Molecular Dynamics Simulation Based Analysis. Sci. Rep. 10, 20397. doi:10.1038/s41598-020-77524-x
Padma, T. V. (2021). COVID Vaccines to Reach Poorest Countries in 2023 - Despite Recent Pledges. Nature 595, 342–343. doi:10.1038/d41586-021-01762-wAvailable from: https://www.nature.com/articles/d41586-021-01762-w (Accessed October 24, 2021)
Palmeira, A., Sousa, E., Köseler, A., Sabirli, R., Gören, T., Türkçüer, İ., et al. (2020). Preliminary Virtual Screening Studies to Identify Grp78 Inhibitors Which May Interfere with Sars-Cov-2 Infection. Pharmaceuticals (Basel) 13, 1–13. doi:10.3390/ph13060132
Paul, D., Mahanta, S., Tag, H., Das, S. K., das Gupta, D., Tanti, B., et al. (2021). Identification of Tyrosine Kinase Inhibitors from Panax Bipinnatifidus and Panax Pseudoginseng for RTK-HER2 and VEGFR2 Receptors, by In Silico Approach. Mol. Divers. doi:10.1007/s11030-021-10304-5
Peacock, T. P., Goldhill, D. H., Zhou, J., Baillon, L., Frise, R., Swann, O. C., et al. (2021). The Furin Cleavage Site in the SARS-CoV-2 Spike Protein Is Required for Transmission in Ferrets. Nat. Microbiol. 6, 899–909. doi:10.1038/s41564-021-00908-w
Pires, D. E., Blundell, T. L., and Ascher, D. B. (2015). pkCSM: Predicting Small-Molecule Pharmacokinetic and Toxicity Properties Using Graph-Based Signatures. J. Med. Chem. 58, 4066–4072. doi:10.1021/acs.jmedchem.5b00104
Puttaswamy, H., Gowtham, H. G., Ojha, M. D., Yadav, A., Choudhir, G., Raguraman, V., et al. (2020). In Silico studies Evidenced the Role of Structurally Diverse Plant Secondary Metabolites in Reducing SARS-CoV-2 Pathogenesis. Sci. Rep. 10, 20584. doi:10.1038/s41598-020-77602-0
Roman Laskowski, B. A., Macarthur, M. W., and Thornton, J. M. (1983). Computer Programs PROCHECK: A Program to Check the Stereochemicai Quality of Protein Structures.
Saxena, A. (2020). Drug Targets for COVID-19 Therapeutics: Ongoing Global Efforts. J. Biosci. 45, 87. doi:10.1007/s12038-020-00067-w
Schreiber, G. (2020). The Role of Type I Interferons in the Pathogenesis and Treatment of COVID-19. Front. Immunol. 11, 595739. doi:10.3389/fimmu.2020.595739
Semwal, D. K., Semwal, R. B., Combrinck, S., and Viljoen, A. (2016). Myricetin: A Dietary Molecule with Diverse Biological Activities. Nutrients 8, 90. doi:10.3390/nu8020090
Shamsi, A., Mohammad, T., Anwar, S., AlAjmi, M. F., Hussain, A., Rehman, M. T., et al. (2020). Glecaprevir and Maraviroc Are High-Affinity Inhibitors of SARS-CoV-2 Main Protease: Possible Implication in COVID-19 Therapy. Biosci. Rep. 40, BSR20201256. doi:10.1042/BSR20201256
Singh, T., Biswas, D., and Jayaram, B. (2011). AADS--an Automated Active Site Identification, Docking, and Scoring Protocol for Protein Targets Based on Physicochemical Descriptors. J. Chem. Inf. Model. 51, 2515–2527. doi:10.1021/ci200193z
Song, X., Tan, L., Wang, M., Ren, C., Guo, C., Yang, B., et al. (2021). Myricetin: A Review of the Most Recent Research. Biomed. Pharmacother. 134, 111017. doi:10.1016/j.biopha.2020.111017
Sridhar, S., Saha, G., Lata, S., and Mehrotra, R. (2021). Molecular Docking Studies of Indian Variants of Pathophysiological Proteins of SARS-CoV-2 with Selected Drug Candidates. J. Genet. 100, 64. doi:10.1007/s12041-021-01313-2
Suganya, R., Kalva, S., Saleena, L. M., Rathi Suganya, P., and Sudevan, K. (2014). Homology Modeling for Human ADAM12 Using Prime, I-Tasser and Easymodeller Molecular Modelling of Synthetic Ion Channels View Project Structure Based Pharmacophore Modelling and Docking with Adam 33 for Identification of Natural lead Compounds View Project HOMOLOGY MODELING for HUMAN ADAM12 USING PRIME, I-TASSER and EASYMODELLER. Available from: https://www.researchgate.net/publication/286503309.
Tatar, G., Ozyurt, E., and Turhan, K. (2021). Computational Drug Repurposing Study of the RNA Binding Domain of SARS-CoV-2 Nucleocapsid Protein with Antiviral Agents. Biotechnol. Prog. 37, e3110. doi:10.1002/btpr.3110
Tay, M. Z., Poh, C. M., Rénia, L., MacAry, P. A., and Ng, L. F. P. (2020). The trinity of COVID-19: Immunity, Inflammation and Intervention. Nat. Rev. Immunol. 20, 363–374. doi:10.1038/s41577-020-0311-8
Tazikeh-Lemeski, E., Moradi, S., Raoufi, R., Shahlaei, M., Janlou, M. A. M., and Zolghadri, S. (2021). Targeting SARS-COV-2 Non-structural Protein 16: a Virtual Drug Repurposing Study. J. Biomol. Struct. Dyn. 39, 4633–4646. doi:10.1080/07391102.2020.1779133
Thuy, B. T. P., My, T. T. A., Hai, N. T. T., Hieu, L. T., Hoa, T. T., Thi Phuong Loan, H., et al. (2020). Investigation into SARS-CoV-2 Resistance of Compounds in Garlic Essential Oil. ACS Omega 5, 8312–8320. doi:10.1021/acsomega.0c00772
Venselaar, H., Joosten, R. P., Vroling, B., Baakman, C. A., Hekkelman, M. L., Krieger, E., et al. (2010). Homology Modelling and Spectroscopy, a Never-Ending Love story. Eur. Biophys. J. 39, 551–563. doi:10.1007/s00249-009-0531-0
WHO Coronavirus (COVID-19) Dashboard | WHO Coronavirus (COVID-19) Dashboard with Vaccination Data Available from: https://covid19.who.int/ (Accessed October 22, 2021).
Wiederstein, M., and Sippl, M. J. (2007). ProSA-web: Interactive Web Service for the Recognition of Errors in Three-Dimensional Structures of Proteins. Nucleic Acids Res. 35, W407–W410. doi:10.1093/nar/gkm290
Wilmes, G. M., Archambault, V., Austin, R. J., Jacobson, M. D., Bell, S. P., and Cross, F. R. (2004). Interaction of the S-phase Cyclin Clb5 with an "RXL" Docking Sequence in the Initiator Protein Orc6 Provides an Origin-Localized Replication Control Switch. Genes Dev. 18, 981–991. doi:10.1101/gad.1202304
Wu, C., Liu, Y., Yang, Y., Zhang, P., Zhong, W., Wang, Y., et al. (2020). Analysis of Therapeutic Targets for SARS-CoV-2 and Discovery of Potential Drugs by Computational Methods. Acta Pharm. Sin B 10, 766–788. doi:10.1016/j.apsb.2020.02.008
Yu, M. S., Lee, J., Lee, J. M., Kim, Y., Chin, Y. W., Jee, J. G., et al. (2012). Identification of Myricetin and Scutellarein as Novel Chemical Inhibitors of the SARS Coronavirus Helicase, nsP13. Bioorg. Med. Chem. Lett. 22, 4049–4054. doi:10.1016/j.bmcl.2012.04.081
Zheng, Y., Zhuang, L., Han, M. L., Kang, D., Gaou, C., Wang, P. H., et al. (2020). Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Membrane (M) Protein Inhibits Type I and III Interferon Production by Targeting RIG-I/MDA-5 Signaling. Signal. Transduct Target. Thery 5, 299. doi:10.1038/s41392-020-00438-7
Keywords: homology modelling, phyto-compound, molecular dynamics, SARS-CoV2, in silico, membrane glycoprotein
Citation: Mahanta S, Naiya T, Biswas K, Changkakoti L, Mohanta YK, Tanti B, Mishra AK, Mohanta TK and Sharma N (2022) Plant Source Derived Compound Exhibited In Silico Inhibition of Membrane Glycoprotein In SARS-CoV-2: Paving the Way to Discover a New Class of Compound For Treatment of COVID-19. Front. Pharmacol. 13:805344. doi: 10.3389/fphar.2022.805344
Received: 01 November 2021; Accepted: 04 March 2022;
Published: 07 April 2022.
Edited by:
Balakumar Chandrasekaran, ITM University, IndiaReviewed by:
Sailendra Kumar Mahanta, Assam Kaziranga University, IndiaDnyaneshwar Rathod, Sant Gadge Baba Amravati University, India
Copyright © 2022 Mahanta, Naiya, Biswas, Changkakoti, Mohanta, Tanti, Mishra, Mohanta and Sharma. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Awdhesh Kumar Mishra, YXdhZGhlc2guYmlvdGVjaDA3QGdtYWlsLmNvbQ==; Tapan Kumar Mohanta, bm9zdG9jLnRhcGFuQGdtYWlsLmNvbQ==, dGFwYW4ubW9oYW50YUB1bml6d2EuZWR1Lm9t; Nanaocha Sharma, c2hhcm1hLm5hbmFvY2hhQGdtYWlsLmNvbQ==
†These authors have contributed equally to this work