In silico molecular docking analysis for repurposing approved antiviral drugs against SARS-CoV-2 main protease

Approaches Towards Better Immunosuppressive Agents

Several classes of compounds are applied in clinics due to their immunosuppressive properties in transplantology and the treatment of autoimmune diseases. Derivatives of mycophenolic acid, corticosteroids and chemotherapeutics bearing heterocyclic moieties like methotrexate, azathioprine, mizoribine, and ruxolitinib are active substances with investigated mechanisms of action. However, improved synthetic approaches of known drugs and novel derivatives are still being reported to attempt better accessibility and therapeutic properties. In this review article, we present the synthesis of the designed chemical structures based on recent literature reports concerning novel compounds as promising immunosuppressive drugs. Moreover, some of the discussed derivers revealed also other types of activities with prospective medicinal potential.

A computational simulation appraisal of banana lectin as a potential anti-SARS-CoV-2 candidate by targeting the receptor-binding domain

BACKGROUND

The ongoing concern surrounding coronavirus disease 2019 (COVID-19) primarily stems from continuous mutations in the genome of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), leading to the emergence of numerous variants. The receptor-binding domain (RBD) in the S1 subunit of the S protein of the virus plays a crucial role in recognizing the host's angiotensin-converting enzyme 2 (hACE2) receptor and facilitating cell membrane fusion processes, making it a potential target for preventing viral entrance into cells. This research aimed to determine the potential of banana lectin (BanLec) proteins to inhibit SARS-CoV-2 attachment to host cells by interacting with RBD through computational modeling.

MATERIALS AND METHODS

The BanLecs were selected through a sequence analysis process. Subsequently, the genes encoding BanLec proteins were retrieved from the Banana Genome Hub database. The FGENESH online tool was then employed to predict protein sequences, while web-based tools were utilized to assess the physicochemical properties, allergenicity, and toxicity of BanLecs. The RBDs of SARS-CoV-2 were modeled using the SWISS-MODEL in the following step. Molecular docking procedures were conducted with the aid of ClusPro 2.0 and HDOCK web servers. The three-dimensional structures of the docked complexes were visualized using PyMOL. Finally, molecular dynamics simulations were performed to investigate and validate the interactions of the complexes exhibiting the highest interactions, facilitating the simulation of their dynamic properties.

RESULTS

The BanLec proteins were successfully modeled based on the RNA sequences from two species of banana (Musa sp.). Moreover, an amino acid modification in the BanLec protein was made to reduce its mitogenicity. Theoretical allergenicity and toxicity predictions were conducted on the BanLecs, which suggested they were likely non-allergenic and contained no discernible toxic domains. Molecular docking analysis demonstrated that both altered and wild-type BanLecs exhibited strong affinity with the RBD of different SARS-CoV-2 variants. Further analysis of the molecular docking results showed that the BanLec proteins interacted with the active site of RBD, particularly the key amino acids residues responsible for RBD's binding to hACE2. Molecular dynamics simulation indicated a stable interaction between the Omicron RBD and BanLec, maintaining a root-mean-square deviation (RMSD) of approximately 0.2 nm for a duration of up to 100 ns. The individual proteins also had stable structural conformations, and the complex demonstrated a favorable binding-free energy (BFE) value.

CONCLUSIONS

These results confirm that the BanLec protein is a promising candidate for developing a potential therapeutic agent for combating COVID-19. Furthermore, the results suggest the possibility of BanLec as a broad-spectrum antiviral agent and highlight the need for further studies to examine the protein's safety and effectiveness as a potent antiviral agent.

SARS-CoV-2 main protease mutation analysis via a kinematic method

The Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2) is the virus responsible for the COVID-19 pandemic. COVID-19 continues to cause millions of deaths globally in part due to immune-evading mutations. SARS-CoV-2 main protease (Mpro) is an important enzyme for viral replication and potentially an effective drug target. Mutations affect the dynamics of enzymes and thereby their activity and ability to bind ligands. Here, we use kinematic flexibility analysis (KFA) to identify how mutations and ligand binding changes the conformational flexibility of Mpro. KFA decomposes macromolecules into regions of different flexibility near-instantly from a static structure, allowing conformational dynamics analysis at scale. Altogether, we analyzed 47 mutation sites across 69 Mpro-ligand complexes resulting in more than 3300 different structures which includes 69 mutated structures with all 47 sites mutated simultaneously and 3243 single residue mutated structures. We found that mutations generally increased the conformational flexibility of the protein. Understanding the impact of mutations on the flexibility of Mpro is essential for identifying potential drug targets in the treatment of SARS-CoV-2. Further studies in this area can offer valuable insights into the mechanisms of molecular recognition.

Repurposing of neprilysin inhibitor 'sacubitrilat' as an anti-cancer drug by modulating epigenetic and apoptotic regulators

Modifications in the epigenetic landscape have been considered a hallmark of cancer. Histone deacetylation is one of the crucial epigenetic modulations associated with the aggressive progression of various cancer subtypes. Herein, we have repurposed the neprilysin inhibitor sacubitrilat as a potent anticancer agent using in-silico protein-ligand interaction profiler (PLIP) analysis, molecular docking, and in vitro studies. The screening of PLIP profiles between vorinostat/panobinostat and HDACs/LTA4H followed by molecular docking resulted in five (Sacubitrilat, B65, BDS, BIR, and NPV) FDA-approved, experimental and investigational drugs. Sacubitrilat has demonstrated promising anticancer activity against colorectal cancer (SW-480) and triple-negative breast cancer (MDA-MB-231) cells, with IC₅₀ values of 14.07 μg/mL and 23.02 μg/mL, respectively. FACS analysis revealed that sacubitrilat arrests the cell cycle at the G0/G1 phase and induces apoptotic-mediated cell death in SW-480 cells. In addition, sacubitrilat inhibited HDAC isoforms at the transcriptomic level by 0.7-0.9 fold and at the proteomic level by 0.5-0.6 fold as compared to the control. Sacubitrilat increased the protein expression of tumor-suppressor (p53) and pro-apoptotic makers (Bax and Bid) by 0.2-2.5 fold while decreasing the expression of anti-apoptotic Bcl2 and Nrf2 proteins by 0.2-0.5 fold with respect to control. The observed cleaved PARP product indicates that sacubitrilat induces apoptotic-mediated cell death. This study may pave the way to identify the anticancer potential of sacubitrilat and can be explored in human clinical trials.

An investigation for the interaction of gamma oryzanol with the Mpro of SARS-CoV-2 to combat COVID-19: DFT, molecular docking, ADME and molecular dynamics simulations

COVID-19 has affected more or less every nation across the world and affected the economy very badly. Infection of this virus in human took the life of millions. We have already faced the first and the second waves of COVID-19 and recently, the nations or humanity is afraid of new strain, that is, OMICRON. Considered to highly infectious than the previous strains. Therefore, the researchers are working to find a promising molecule with no or permissible toxicity. In the present work, authors have chosen 10 molecules including the molecules used in curing the infection from nCoV. All the molecules were docked against Mpro of nCoV using iGemdock, a reliable computational tool. Based on the binding energy obtained, it can be seen that only latermovir; remdesivir; zanamivir showed better binding affinity than the gamma oryzanol, the molecule of interest in this work. These three molecules are already in use to cure the patients siffering from the infection of nCoV. But, we need a cost effective and easily available molecule to fight against this viral infection. The binding energy obtained for the formation of complex of gamma oryzanol with Mpro of nCoV through molecular docking is -118.787 kcal/mol. It forms conventional hydrogen bonds with the CYS145 (2.51 Å), LEU141 (3.01 Å) and SER144 (3.09 Å); forms C-H bonds with PHE140 (3.37 Å) and HIS163 (2.91 Å), forms alkyl interactions with ALA191 (3.59 and 4.74 Å), CYS145 (4.90 Å). One interesting information is obtained that the value of log Kp of gamma oryzanol is least means more permeable to skin in comparison of other molecules used in the work. Gamma oryzanol in known for to its biological potency like it can modulate the oxidative stress as well as inflammation. DFT calculations of gamma oryzanol (GO) was made at different temperature and no change in the delocalization of electron density as well no change in free energy is observed. Molecular dynamics (MD) simulations of gamma oryzanol with the Mpro of nCoV at different temperatures was performed. The formation of the complex between GO and Mpro of CoV at 290 K, 300 K, 310 K and 320 K for 100 ns was investigated. It has been observed that the effective binding is observed at 290 K, therefore, it can be said that the inhibition of the Mpro of nCoV with GO is maximum at 290 K.Communicated by Ramaswamy H. Sarma.

Potential antiviral peptides targeting the SARS-CoV-2 spike protein

BACKGROUND

The coronavirus disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection became an international pandemic and created a public health crisis. The binding of the viral Spike glycoprotein to the human cell receptor angiotensin-converting enzyme 2 (ACE2) initiates viral infection. The development of efficient treatments to combat coronavirus disease is considered essential.

METHODS

An in silico approach was employed to design amino acid peptide inhibitor against the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein. The designed inhibitor (SARS-CoV-2 PEP 49) consists of amino acids with the α1 helix and the β4 - β5 sheets of ACE2. The PEP-FOLD3 web tool was used to create the 3D structures of the peptide amino acids. Analyzing the interaction between ACE2 and the RBD of the Spike protein for three protein data bank entries (6M0J, 7C8D, and 7A95) indicated that the interacting amino acids were contained inside two regions of ACE2: the α1 helical protease domain (PD) and the β4 - β5 sheets.

RESULTS

Molecular docking analysis of the designed inhibitor demonstrated that SARS-CoV-2 PEP 49 attaches directly to the ACE2 binding site of the Spike protein with a binding affinity greater than the ACE2, implying that the SARS-CoV-2 PEP 49 model may be useful as a potential RBD binding blocker.

Screening of Phytoconstituents from Traditional Plants against SARSCoV- 2 using Molecular Docking Approach

Background:

The emergence of COVID-19 as a fatal viral disease encourages researchers to develop effective and efficient therapeutic agents. The intervention of in silico studies has led to revolutionary changes in the conventional method of testing the bioactivity of plant constituents.

Objective:

The current study deals with the investigation of some traditional immunomodulators of plant origin to combat this ailment.

Materials and Methods:

A total of 151 phytomolecules of 12 immunomodulatory plants were evaluated for their inhibitory action against the main protease (PDB ID: 7D1M) and NSP15 endoribonuclease (PDB ID: 6WLC) by structure-based virtual screening. In addition, the promising molecules with ligand efficiency of more than -0.3(kcal/mol)/heavy atoms were further predicted for pharmacokinetic properties and druggability using the SwissADME web server, and their toxicity was also evaluated using Protox-II.

Result:

Myricetin-3-O-arabinofuranoside of cranberry plant was found to be the most potential candidate against both enzymes: main protease (–14.2 kcal/mol) and NSP15 endoribonuclease (–12.2 kcal/mol).

Conclusion:

The promising outcomes of the current study may be implemented in future drug development against coronavirus. The findings also help in the development of lead candidates of plant origin with a better ADMET profile in the future.

A comparative study of 5- fluorouracil, doxorubicin, methotrexate, paclitaxel for their inhibition ability for Mpro of nCoV: Molecular docking and molecular dynamics simulations

A new corona virus (nCoV) is aetiological agent responsible for the viral pneumonia epidemic. Three is no specific therapeutic medicines available for the treatment of this condition and also effective treatment choices are few. In this work author tried to investigate some repurposing drug such as 5- fluorouracil, doxorubicin, methotrexate and paclitaxel against the main protease (Mpro) of nCoV by the computational model. Molecular docking was performed to screen out the best compound and doxorubicin was found to have minimum binding energy −121.89 kcal/mol. To further study, MD simulations were performed at 300 K and the result successfully corroborate the energy obtained by molecular docking. Temperature dependent MD simulation of the best molecule that is doxorubicin obtained from docking result was performed to check the variation in structural changes in Mpro of nCoV at 290 K, 310 K, 320 K and 325 K. It is sound that doxorubicin binds effectively with Mpro of nCoV at 290 K. Further ADME properties of the 5- fluorouracil, doxorubicin, methotrexate and paclitaxel were also evaluated to understand the bioavailability.

Exploring the anti-SARS-CoV-2 main protease potential of FDA approved marine drugs using integrated machine learning templates as predictive tools

Since the inception of COVID-19 pandemic in December 2019, socio-economic crisis begins to rise globally and SARS-CoV-2 was responsible for this outbreak. With this outbreak, currently, world is in need of effective and safe eradication of COVID-19. Hence, in this study anti-SAR-Co-2 potential of FDA approved marine drugs (Biological macromolecules) data set is explored computationally using machine learning algorithm of Flare by Cresset Group, Field template, 3D-QSAR and activity Atlas model was generated against FDA approved M-pro SARS-CoV-2 repurposed drugs including Nafamostat, Hydroxyprogesterone caporate, and Camostat mesylate. Data sets were categorized into active and inactive molecules on the basis of their structural and biological resemblance with repurposed COVID-19 drugs. Then these active compounds were docked against the five different M-pro proteins co-crystal structures. Highest LF VS score of Holichondrin B against all main protease co-crystal structures ranked it as lead drug. Finally, this new technique of drug repurposing remained efficient to explore the anti-SARS-CoV-2 potential of FDA approved marine drugs.

Synthesis, crystal structure, DFT and molecular docking studies of N-acetyl-2,4-[diaryl-3-azabicyclo[3.3.1]nonan-9-yl]-9-spiro-4'-acetyl-2'-(acetylamino)-4',9-dihydro-[1',3',4']-thiadiazoles: A potential SARS-nCoV-2 Mpro (COVID-19) inhibitor

In this paper, we describe the synthesis and crystal structure analysis of N-acetyl-2,4-[diphenyl-3-azabicyclo[3.3.1]nonan-9-yl]-9-spiro-4′-acetyl-2′-(acetylamino)-4′,9-dihydro-[1′,3′,4′]-thiadiazole (3a) and N-acetyl- 2,4-[bis(p-methoxyphenyl)-3-azabicyclo[3.3.1]nonan-9-yl]-9-spiro-4′-acetyl-2′-(acetylamino)-4′,9-dihydro-[1′,3′,4′]-thiadiazole (3b). The title compounds 3a and 3b are characterized by 1D NMR and single crystal x-ray diffraction analysis. Non-covalent interactions in a molecule were identified by Hirshfeld surface (d_norm contacts and 2D fingerprint plot) analysis. In addition, the existence of chalcogen bond (S•••O bond) in the molecular structures (3a and 3b) are described by NCI-RDG and QTAIM analysis. NBO analysis is employed to describe the orbital interactions and electron transfer between sulfur and oxygen atoms. Molecular docking is carried out for compounds 3a and 3b with COVID-19 viral protein SARS-nCoV-2 M^pro (PDB ID: 6LU7).

Development of a database of RNA helicase inhibitors (VHIMDB) of pathogenic viruses and in silico screening for the potential drug molecules

The pathogenic RNA virus that infects human beings contains the RNA helicase enzyme, responsible for the replication of the viral genome. The enzyme is used as a suitable target against which the drug molecule acts. Therefore, the identification and proposal the novel compounds that can be targeted toward the helicase enzymes to stop the functioning of the enzyme is desirable. Although many viral helicase inhibitor molecules have been identified, still yet no unique database is available for these compounds. This research work envisages developing a curated database of RNA helicase inhibitors. The database contains in total of 353 entries that are computationally predicted and experimentally verified RNA helicase inhibitors. The database contains information like compound name, chemical properties, chemical format, and name of the target virus to which it acts against it with a user-friendly menu-driven search engine. Presently, the database is freely available at: https://vhimdb.rsatpathy.in/. Further, in silico screening of the whole database by drug-likeness and toxicity resulted in 14 potential drug molecules. The selected molecules were analyzed for their effectiveness in binding by using molecular docking score and interaction with the helicase enzymes of three categories of pathogenic viruses (SARS-CoV-2, SARS-CoV, and MERS-CoV).

Seeking antiviral drugs to inhibit SARS-CoV-2 RNA dependent RNA polymerase: A molecular docking analysis

COVID-19 outbreak associated with the severe acute respiratory syndrome coronavirus (SARS-CoV-2) raised health concerns across the globe and has been considered highly transmissible between people. In attempts for finding therapeutic treatment for the new disease, this work has focused on examining the polymerase inhibitors against the SARS-CoV-2 nsp12 and co-factors nsp8 and nsp7. Several polymerase inhibitors were examined against PDB ID: 6M71 using computational analysis evaluating the ligand's binding affinity to replicating groove to the active site. The findings of this analysis showed Cytarabine of -5.65 Kcal/mol with the highest binding probability (70%) to replicating groove of 6M71. The complex stability was then examined over 19 ns molecular dynamics simulation suggesting that Cytarabine might be possible potent inhibitor for the SARS-CoV-2 RNA Dependent RNA Polymerase.

Active Compounds in Zingiber officinale as Possible Redox Inhibitors of 5-Lipoxygenase Using an In Silico Approach

5-Lipoxygenase (5-LOX) converts arachidonic acid to lipidic inflammatory mediators such as leukotrienes (LTs). In diseases such as asthma, LTs contribute to a physiopathology that could be reverted by blocking 5-LOX. Natural products with anti-inflammatory potential such as ginger have been used as nutraceuticals since ancient times. 6-Gingerol and 6-shogaol are the most abundant compounds in the ginger rhizome; they possess anti-inflammatory, antioxidant, and chemopreventive properties. In the present study, 6-gingerol and 6-shogaol structures were analyzed and compared with two commercial 5-LOX inhibitors (zileuton and atreleuton) and with other inhibitor candidates (3f, NDGA, CP 209, caffeic acid, and caffeic acid phenethyl ester (CAPE)). The pharmacokinetics and toxicological properties of 6-gingerol, 6-shogaol, and the other compounds were evaluated. Targeted molecular coupling was performed to identify the optimal catalytic pocket for 5-LOX inhibition. The results showed that 6-gingerol and 6-shogaol follow all of the recommended pharmacokinetic parameters. These compounds could be inhibitors of 5-LOX because they present specific interactions with the residues involved in molecular inhibition. The current study demonstrated the potential of 6-gingerol and 6-shogaol as anti-inflammatory agents that inhibit 5-LOX, as they present a high level of performance in the toxicological analysis and could be catabolized by the cytochrome p450 enzymatic complex; however, 6-gingerol was superior in safety compared to 6-shogaol.

Identification of bioactive molecules from Triphala (Ayurvedic herbal formulation) as potential inhibitors of SARS-CoV-2 main protease (Mpro) through computational investigations

Severe acute respiratory syndrome coronavirus disease (SARS-CoV-2) induced coronavirus disease 2019 (COVID-19) pandemic is the present worldwide health emergency. The global scientific community faces a significant challenge in developing targeted therapies to combat the SARS-CoV-2 infection. Computational approaches have been critical for identifying potential SARS-CoV-2 inhibitors in the face of limited resources and in this time of crisis. Main protease (M^pro) is an intriguing drug target because it processes the polyproteins required for SARS-CoV-2 replication. The application of Ayurvedic knowledge from traditional Indian systems of medicine may be a promising strategy to develop potential inhibitor for different target proteins of SARS-CoV-2. With this endeavor, we docked bioactive molecules from Triphala, an Ayurvedic formulation, against M^pro followed by molecular dynamics (MD) simulation (100 ns) to investigate their inhibitory potential against SARS-CoV-2. The top four best docked molecules (terflavin A, chebulagic acid, chebulinic acid, and corilagin) were selected for MD simulation study and the results obtained were compared to native ligand X77. From docking and MD simulation studies, the selected molecules showed promising binding affinity with the formation of stable complexes at the active binding pocket of M^pro and exhibited negative binding energy during MM-PBSA calculations, indication their strong binding affinity with the target protein. The identified bioactive molecules were further analyzed for drug-likeness by Lipinski's filter, ADMET and toxicity studies. Computational (in silico) investigations identified terflavin A, chebulagic acid, chebulinic acid, and corilagin from Triphala formulation as promising inhibitors of SARS-CoV-2 M^pro, suggesting experimental (in vitro/in vivo) studies to further explore their inhibitory mechanisms.

Artecanin of Laurus nobilis is a novel inhibitor of SARS-CoV-2 main protease with highly desirable druglikeness

Main protease (M^pro) is a critical enzyme in the life cycle of severe acute respiratory syndrome Coronavirus -2 (SARS-CoV-2). Due to its essential role in the maturation of the polyproteins, the necessity to inhibit M^pro is one of the essential means to prevent the outbreak of COVID-19. In this context, this study was conducted on the natural compounds of medicinal plants that are commonly available in the Middle East to find out the most potent one to inhibit M^pro with the best bioavailability and druglikeness properties. A total of 3392 compounds of sixty-six medicinal plants were retrieved from PubChem database and docked against M^pro. Thirty compounds with the highest docking scores with M^pro were chosen for further virtual screening. Variable druglikeness and toxicity potentials of these compounds were evaluated using SwissADME and Protox servers respectively. Out of these virtually screened compounds, artecanin was predicted to exhibit the most favourable druglikeness potentials, accompanied by no predicted hepatoxicity, carcinogenicity, mutagenicity, and cytotoxicity. Molecular dynamics (MD) simulations showed that M^pro-artecanin complex exhibited comparable stability with that observed in the ligand-free M^pro. This study revealed for the first time that artecanin from Laurus nobilis provided a novel static and dynamic inhibition for M^pro with excellent safety, oral bioavailability, and pharmacokinetic profile. This study suggested the ability of artecanin to be used as a potential natural inhibitor that can be used to block or at least counteract the SARS-CoV-2 invasion.Communicated by Ramaswamy H. Sarma.