Computational basis of SARS-CoV 2 main protease inhibition: an insight from molecular dynamics simulation based findings

A computational approach to identify phytochemicals as potential inhibitor of acetylcholinesterase: Molecular docking, ADME profiling and molecular dynamics simulations

Inhibition of acetylcholinesterase (AChE) is a crucial target in the treatment of Alzheimer's disease (AD). Common anti-acetylcholinesterase drugs such as Galantamine, Rivastigmine, Donepezil, and Tacrine have significant inhibition potential. Due to side effects and safety concerns, we aimed to investigate a wide range of phytochemicals and structural analogues of these compounds. Compounds similar to the established drugs, and phytochemicals were investigated as potential inhibitors for AChE in treating AD. A total of 2,270 compound libraries were generated for further analysis. Initial virtual screening was performed using Pyrx software, resulting in 638 molecules showing higher binding affinities compared to positive controls Tacrine (-9.0 kcal/mol), Donepezil (-7.3 kcal/mol), Galantamine (-8.3 kcal/mol), and Rivastigmine (-6.4 kcal/mol). Subsequently, ADME properties were assessed, including blood-brain barrier permeability and Lipinski's rule of five violations, leading to 88 compounds passing the ADME analysis. Among the rivastigmine analogous, [3-(1-methylpiperidin-2-yl)phenyl] N,N-diethylcarbamate showed interaction with Tyr123, Tyr336, Tyr340, Phe337, Trp285 residues of AChE. Tacrine similar compounds, such as 4-amino-2-styrylquinoline, exhibited bindings with Tyr123, Phe337, Tyr336, Trp285, Trp85, Gly119, and Gly120 residues. A phytocompound (bisdemethoxycurcumin) showed interaction with Trp285, Tyr340, Trp85, Tyr71, and His446 residues of AChE with favourable binding. These findings underscore the potential of these compounds as novel inhibitors of AChE, offering insights into alternative therapeutic avenues for AD. A 100ns simulation analysis confirmed the stability of protein-ligand complex based on the RMSD, RMSF, ligand properties, PCA, DCCM and MMGBS parameters. The investigation suggested 3 ligands as a potent inhibitor of AChE which are [3-(1-methylpiperidin-2-yl)phenyl] N,N-diethylcarbamate, 4-Amino-2-styrylquinoline and bisdemethoxycurcumin. Furthermore, investigation, including in-vitro and in-vivo studies, is needed to validate the efficacy, safety profiles, and therapeutic potential of these compounds for AD treatment.

Overview of the SARS-CoV-2 nucleocapsid protein

Since the emergence of SARS-CoV in 2003, researchers worldwide have been toiling away at deciphering this virus's biological intricacies. In line with other known coronaviruses, the nucleocapsid (N) protein is an important structural component of SARS-CoV. As a result, much emphasis has been placed on characterizing this protein. Independent research conducted by a variety of laboratories has clearly demonstrated the primary function of this protein, which is to encapsidate the viral genome. Furthermore, various accounts indicate that this particular protein disrupts diverse intracellular pathways. Such observations imply its vital role in regulating the virus as well. The opening segment of this review will expound upon these distinct characteristics succinctly exhibited by the N protein. Additionally, it has been suggested that the N protein possesses diagnostic and vaccine capabilities when dealing with SARS-CoV. In light of this fact, we will be reviewing some recent headway in the use cases for N protein toward clinical purposes within this article's concluding segments. This forward movement pertains to both developments of COVID-19-oriented therapeutic targets as well as diagnostic measures. The strides made by medical researchers offer encouragement, knowing they are heading toward a brighter future combating global pandemic situations such as these.

Discovering potential inhibitors of Raf proto-oncogene serine/threonine kinase 1: a virtual screening approach towards anticancer drug development

Raf proto-oncogene serine/threonine kinase 1 (RAF1 or c-Raf) is a serine/threonine protein kinase crucial in regulating cell growth, differentiation, and survival. Any disruption or overexpression of RAF1 can result in neoplastic transformation and other disorders such as cardiomyopathy, Noonan syndrome, leopard syndrome, etc. RAF1 has been identified as a potential therapeutic target in drug development against various complex diseases, including cancer, due to its remarkable role in disease progression. Here, we carried out a multitier virtual screening study involving different in-silico approaches to discover potential inhibitors of RAF1. After applying the Lipinski rule of five, we retrieved all phytocompounds from the IMPPAT database based on their physicochemical properties. We performed a molecular docking-based virtual screening and got top hits with the best binding affinity and ligand efficiency. Then we screened out the selected hits using the PAINS filter, ADMET properties, and other druglike features. Eventually, PASS evaluation identifies two phytocompounds, Moracin C and Tectochrysin, with appreciable anti-cancerous properties. Finally, all-atom molecular dynamics simulation (MDS) followed by interaction analysis was performed on the elucidated compounds in complex with RAF1 for 200 ns to investigate their time-evolution dynamics and interaction mechanism. Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and Dynamical Cross-Correlation Matrix (DCCM) analyses then followed these results from the simulated trajectories. According to the results, the elucidated compounds stabilize the RAF1 structure and lead to fewer conformational alterations. The results of the current study indicated that Moracin C and Tectochrysin could serve as potential inhibitors of RAF1 after required validation.Communicated by Ramaswamy H. Sarma.

Structure and dynamic simulation-based interactions of benzenoids, pyrroles and organooxygen compounds for effective targeting of GPX4 in ischemic stroke

The discovery of a novel drug for ischemic stroke is plagued by expensive and unsuccessful outcomes. FDA-approved drugs could be a viable repurposing strategy for stroke therapy. Emerging evidence suggests the regulating role of Glutathione peroxidase (GPX4) in stroke and attracts as a potential target. To overcome limited therapeutic interventions, a drug repurposing in silico investigation of FDA-approved drugs is proposed for the GPX4 receptor in distinctive species (Homo sapiens and Mus musculus). The GPX4 UniProt wild type ids, that is, P36969 (Homo sapiens), P36970 (Rattus norvegicus) and O70325 (Mus musculus) are Swiss modelled, and resultant templates are 2OBI and 6HN3 for Homo sapiens, and 5L71 for Mus musculus with a sequence identity of ∼88%. Enrichment analysis reveals high sensitivity and ranked actives with ROC and AUC values of 0.59 and 0.61, respectively. Virtual screening at extra precision resulted hit Acarbosum, is similar between 2OBI and 6HN3, demonstrating a multiple-target specificity and Iopromide, targeting 2OBI. MD simulation at 100 ns following trajectory analysis provides RMSD (∼1.2-1.8Å), RMSF (∼1.6-2.7Å), Rgyr (∼15-15.6Å) depicting stabilisation of receptor-ligand complexes. Furthermore, average B-factor value of 2OBI, 6HN3 and 5L71 is 25Å, 24Å and 60Å with a defined resolution of 1.55Å, 1.01Å and 1.80Å, respectively, depicting the thermodynamic stability of the protein structures. The dynamic cross-correlation and principal component analysis of residual fluctuations reveal more positive correlation, high atomic displacements and greater residual clustering of residues from atomic coordinates. Therefore, Acarbosum, an FDA-approved drug, could act as a potential repurposing drug with a multi-target approach translating from preclinical to clinical stages.Communicated by Ramaswamy H. Sarma.

An insight into the simulation directed understanding of the mechanism in SARS CoV-2 N-CTD, dimer integrity, and RNA-binding: Identifying potential antiviral inhibitors

Coronavirus 2019 is a transmissible disease and has caused havoc throughout the world. The present study identifies the novel potential antiviral inhibitors against the nucleocapsid C-terminal domain that aids in RNA-binding and replication. A total of 485,629 compounds were screened, and MD was performed. The trajectory analysis (DCCM & PCA), structural integrity, and degree of compaction depicted the protein-ligand complex stability (PDB-PISA and R_gyr). Results obtained from screening shortlists 13 compounds possessing high Docking score. Further, seven compounds had a permissible RMSD limit (3 Å), with robust RMSF. Post-MD analysis of the top two compounds (204 and 502), DCCM & PCA analysis show a positive atomic displacements correlation among residues of active sites-dimer (Chain A and Chain B) & residual clustering. The ΔG_int of RNA-bound (-83.5 kcal/mol) and drug-bound N-CTD-204 (-40.8 kcal/mol) and 502(-39.7 kcal/mol) as compared to Apo (-35.95 kcal/mol) suggests stabilization of protein, with less RNA-binding possibility. The R_gyr values depict the loss of compactness on RNA-binding when compared to the drug-bound N-CTD complex. Further, overlapping the protein complexes (0 ns and 100 ns) display significant changes in RMSD of the protein (204-2.07 Å and 502-1.89 Å) as compared to the Apo (1.72 Å) and RNA-bound form (1.76 Å), suggesting strong interaction for compound 204 as compared to 502. ADMET profiling indicates that these compounds can be used for further experiments (in vitro and pre-clinical). Compound 204 could be a promising candidate for targeting the N-protein-RNA assembly and viral replication.

Quantitative structure-activity relationships, molecular docking and molecular dynamics simulations reveal drug repurposing candidates as potent SARS-CoV-2 main protease inhibitors

The current outbreak of COVID-19 is leading an unprecedented scientific effort focusing on targeting SARS-CoV-2 proteins critical for its viral replication. Herein, we performed high-throughput virtual screening of more than eleven thousand FDA-approved drugs using backpropagation-based artificial neural networks (q²_LOO = 0.60, r² = 0.80 and r²_pred = 0.91), partial-least-square (PLS) regression (q²_LOO = 0.83, r² = 0.62 and r²_pred = 0.70) and sequential minimal optimization (SMO) regression (q²_LOO = 0.70, r² = 0.80 and r²_pred = 0.89). We simulated the stability of Acarbose-derived hexasaccharide, Naratriptan, Peramivir, Dihydrostreptomycin, Enviomycin, Rolitetracycline, Viomycin, Angiotensin II, Angiotensin 1-7, Angiotensinamide, Fenoterol, Zanamivir, Laninamivir and Laninamivir octanoate with 3CL_pro by 100 ns and calculated binding free energy using molecular mechanics combined with Poisson-Boltzmann surface area (MM-PBSA). Our QSAR models and molecular dynamics data suggest that seven repurposed-drug candidates such as Acarbose-derived Hexasaccharide, Angiotensinamide, Dihydrostreptomycin, Enviomycin, Fenoterol, Naratriptan and Viomycin are potential SARS-CoV-2 main protease inhibitors. In addition, our QSAR models and molecular dynamics simulations revealed that His41, Asn142, Cys145, Glu166 and Gln189 are potential pharmacophoric centers for 3CL_pro inhibitors. Glu166 is a potential pharmacophore for drug design and inhibitors that interact with this residue may be critical to avoid dimerization of 3CL_pro. Our results will contribute to future investigations of novel chemical scaffolds and the discovery of novel hits in high-throughput screening as potential anti-SARS-CoV-2 properties.Communicated by Ramaswamy H. Sarma.

Nelfinavir: An Old Ally in the COVID-19 Fight?

After almost three years of the pandemic, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is still spreading around the world, causing notable sanitary and social issues. New antiviral therapies are constantly under investigation. However, few options have been approved for the treatment of COVID-19. Clinical trials are currently ongoing to evaluate the efficacy of nelfinavir on mild−moderate COVID-19. This study aims to investigate the activity of this compound on SARS-CoV-2 “Variants of Concern” (VOCs), comparing its effectiveness with the approved drugs remdesivir and molnupiravir. The experiments were conducted in a biosafety level 3 facility. In this study, we used a Vero-E6-cell-based infection assay to investigate the in vitro activity of nelfinavir, molnupiravir, and remdesivir. Four strains of SARS-CoV-2 were tested: 20A.EU1, B.1.1.7, P.1, and B.1.617.2. All compounds reached micromolar/submicromolar EC50, EC90, and EC99. Furthermore, the Cmax/EC50 and Cmax/EC90 ratios were >1 for all compounds and all variants tested. Our study demonstrated that nelfinavir, as molnupiravir, and remdesivir are effective in vitro on SARS-CoV-2 variants.

Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2

Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.

In-silico investigation of some recent natural compounds for their potential use against SARS-CoV-2: a DFT, molecular docking and molecular dynamics study

Since its first appearance in December 2019, SARS-CoV-2 has infected many people all over the world, causing serious health problems in many people and causing many deaths, but no specific drug has been developed yet. SARS-CoV-2 main protease (SARS-CoV-2 M^pro) has an important role in viral replication and transcription, so inhibition of this enzyme is proposed to be an attractive route for the treatment of COVID-19. Natural compounds have been used in the treatment of many diseases throughout the history. In this study, it was aimed to investigate SARS-CoV-2 M^pro inhibition abilities, thus the therapeutic potentials of some novel phytochemicals which have recently been entered the literature. For this purpose, eleven novel phytochemicals obtained from various natural resources have been investigated for their potential antiviral activity against SARS-CoV-2 with the use of in silico methods. In the first part of the study, DFT (density functional theory) calculations were performed on the investigated compounds. In this part, geometry optimizations, vibrational analyses, and MEP (molecular electrostatic potential) map calculations were performed. In the second part, molecular docking calculations and then molecular dynamics (MD) simulations were performed to investigate how these natural compounds interact with SARS-CoV-2 M^pro which is a promising target for COVID-19 treatments. In this part, MM-PBSA (molecular mechanics with Poisson-Boltzmann surface area) calculations were also performed to determine the binding free energies of the investigated compounds. Results showed that most of the investigated compounds interacted with SARS-CoV-2 M^pro effectively and can be promising structures for drug development studies for COVID-19.Communicated by Ramaswamy H. Sarma.

Targeting SARS-CoV-2 main protease by teicoplanin: A mechanistic insight by docking, MM/GBSA and molecular dynamics simulation

First emerged in late December 2019, the outbreak of novel severe acute respiratory syndrome corona virus-2 (SARS-CoV-2) pandemic has instigated public-health emergency around the globe. Till date there is no specific therapeutic agent for this disease and hence, the world is craving to identify potential antiviral agents against SARS-CoV-2. The main protease (M^Pro) is considered as an attractive drug target for rational drug design against SARS-CoV-2 as it is known to play a crucial role in the viral replication and transcription. Teicoplanin is a glycopeptide class of antibiotic which is regularly used for treating Gram-positive bacterial infections, has shown potential therapeutic efficacy against SARS-CoV-2 in vitro. Therefore, in this study, a mechanistic insight of intermolecular interactions between teicoplanin and SARS-CoV-2 M^Pro has been scrutinized by molecular docking. Both monomeric and dimeric forms of M^Pro was used in docking involving blind as well as defined binding site based on the known inhibitor. Binding energies of teicoplanin-M^Pro complexes were estimated by Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) computations from docking and simulated trajectories. The dynamic and thermodynamics constraints of docked drug in complex with target proteins under specific physiological conditions was ascertained by all-atom molecular dynamics simulation of 100 ns trajectory. Root mean square deviation and fluctuation of carbon α chain justified the stability of the bound complex in biological environments. The outcomes of current study are supposed to be fruitful in rational design of antiviral drugs against SARS-CoV-2.

Exploring biological efficacy of novel benzothiazole linked 2,5-disubstituted-1,3,4-oxadiazole hybrids as efficient α-amylase inhibitors: Synthesis, characterization, inhibition, molecular docking, molecular dynamics and Monte Carlo based QSAR studies

In an effort to explore a class of novel antidiabetic agents, we have made an effort to synergize the α-amylase inhibitory potential of 1,3-benzothiazole and 1,3,4-oxadiazole scaffolds by combining the two into a single structure via an ether linkage. The structure of synthesized benzothiazole clubbed oxadiazole derivatives are established by different spectral techniques. The synthesized hybrids are evaluated for their in vitro inhibitory potential against α-amylase. Compound 8f is found to be the most potent with a significant inhibition (87.5 ± 0.74% at 50 μg/mL, 82.27 ± 1.85% at 25 μg/mL and 79.94 ± 1.88% at 12.5 μg/mL) when compared to positive control acarbose (77.96 ± 2.06%, 71.17 ± 0.60%, 67.24 ± 1.16% at 50 μg/mL, 25 μg/mL and 12.5 μg/mL concentration). Molecular docking of the most potent enzyme inhibitor, 8f, shows promising interaction with the binding site of biological macromolecule Aspergillus oryzae α-amylase (PDB ID: 7TAA) and human pancreatic α-amylase (PDB ID: 3BAJ). To a step further, in-depth QSAR studies show a significant correlation between the experimental and the predicted inhibitory activities with the best R_validation²= 0.8701. The developed QSAR model can provide ample information about the structural features responsible for the increase and decrease of inhibitory activity. The mechanistic interpretation of the structure-activity relationship (SAR) is done with the help of combined computational calculations i.e. molecular docking and QSAR. Finally, molecular dynamic simulations are performed to get an insight into the binding mode of the most potent derivative with α-amylase from A. oryzae (PDB ID: 7TAA) and human pancreas (PDB ID: 3BAJ).

Identification of homologous human miRNAs as antivirals towards COVID-19 genome

The COVID-19 fatality rate is ~57% worldwide. The investigation of possible antiviral therapy using host microRNA (miRNA) to inhibit viral replication and transmission is the need of the hour. Computational techniques were used to predict the hairpin precursor miRNA (pre-miRNAs) of COVID-19 genome with high homology towards human (host) miRNA. Top 21 host miRNAs with >80% homology towards 18 viral pre miRNAs were identified. The Gibbs free energy (ΔG) between host miRNAs and viral pre-miRNAs hybridization resulted in the best 5 host miRNAs having the highest base-pair complementarity. miR-4476 had the strongest binding with viral pre-miRNA (ΔG = -21.8 kcal/mol) due to maximum base pairing in the seed sequence. Pre-miR-651 secondary structure was most stable due to the (1) least minimum free energy (ΔG = -24.4 kcal/mol), energy frequency, and noncanonical base pairing and (2) maximum number of stem base pairing and small loop size. Host miRNAs-viral mRNAs interaction can effectively inhibit viral transmission and replication. Furthermore, miRNAs gene network and gene-ontology studies indicate top 5 host miRNAs interaction with host genes involved in transmembrane-receptor signaling, cell migration, RNA splicing, nervous system formation, and tumor necrosis factor-mediated signaling in respiratory diseases. This study identifies host miRNA/virus pre-miRNAs strong interaction, structural stability, and their gene-network analysis provides strong evidence of host miRNAs as antiviral COVID-19 agents.

Structural and conformational analysis of SARS CoV 2 N-CTD revealing monomeric and dimeric active sites during the RNA-binding and stabilization: Insights towards potential inhibitors for N-CTD

The advent of SARS-CoV-2 has become a universal health issue with no appropriate cure available to date. The coronavirus nucleocapsid (N) protein combines viral genomic RNA into a ribonucleoprotein and protects the viral genome from the host's nucleases. Structurally, the N protein comprises two independent domains: the N-terminal domain (NTD) for RNA-binding and C-terminal domain (CTD) involved in RNA-binding, protein dimerization, and nucleocapsid stabilization. The present study explains the structural aspects associated with the involvement of nucleocapsid C-terminal domain in the subunit assembly that helps the RNA binding and further stabilizing the virus assembly by protecting RNA from the hosts exonucleases degradation. The molecular dynamics (MD) simulations of the N-CTD and RNA complex suggests two active sites (site I: a monomer) and (site II: a dimer) with structural stability (RMSD: ~2 Å), Cα fluctuations (RMSF: ~3 Å) and strong protein-ligand interactions were estimated through the SiteMap module of Schrodinger. Virtual screening of 2456 FDA-approved drugs using structure-based docking identified top two leads distinctively against Site-I (monomer): Ceftaroline fosamil (MM-GBSA = -47.12 kcal/mol) and Cefoperazone (-45.84 kcal/mol); and against Site-II (dimer): Boceprevir, (an antiviral protease inhibitor, -106.78 kcal/mol) and Ceftaroline fosamil (-99.55 kcal/mol). The DCCM and PCA of drugs Ceftaroline fosamil (PC1+PC2 = 71.9%) and Boceprevir (PC1 +PC2 = 61.6%) show significant correlated residue motions which suggests highly induced conformational changes in the N-CTD dimer. Therefore, we propose N-CTD as a druggable target with two active binding sites (monomer and dimer) involved in specific RNA binding and stability. The RNA binding site with Ceftaroline fosamil binding can prevent viral assembly and can act as an antiviral for coronavirus.

COVID-19–associated rhino-orbital-cerebral mucormycosis: A systematic review, meta-analysis, and meta-regression analysis

BACKGROUND:

Till now, no meta-analysis is available to address the clinical profile, risk factors, different interventions, and outcomes among COVID-19–associated rhino-orbito-cerebral mucormycosis (C-ROCM) cases.

MATERIALS AND METHODS:

Eight literature databases were screened using appropriate keywords from November 1, 2019, to June 30, 2021. The objectives were to analyze the clinical and microbiological profile, risk factor/comorbidity, intervention, and outcome. “R-metafor package” was used for analysis.

RESULTS:

A total of 23 studies were included. The mean age of presentation of C-ROCM was 54.6 years. The most common presentation was ptosis (72.7%), lid edema (60.6%), proptosis (60.6%), ophthalmoplegia (57.3%), loss of vision (53.7%), facial edema (34.7%), and nasal-blockage (11.8%). Evidence of intracranial spread was seen in 42.8% of cases. Rhizopus was the most common fungus (57.1%) isolated in fungal culture. Among C-ROCM patients, diabetes was the commonest comorbid condition, and the use of corticosteroids related to COVID-19 treatment was the most common risk factor (85.75%). Compared to controlled diabetics, C-ROCM was significantly higher among uncontrolled diabetics (odds ratio [OR] 0.15, 95% confidence interval [C.I.] 0.041–0.544, P = 0.0010). However, no significant association was seen between C-ROCM and COVID-19 severity (OR 0.930, 95% C.I. 0.212–4.087, P = 0.923). For treatment, amphotericin-B was the most common antifungal drug used which was followed by surgical options. However, mortality was high (prevalence 0.344, 95% C.I. 0.205–0.403) despite treatment.

CONCLUSION:

Although local rhino-orbito symptoms were the first to appear, rapid intracranial extension was seen in a significant number of C-ROCM cases. Uncontrolled diabetes and excessive use of corticosteroid were the most common risk factors present among the C-ROCM cases. High index clinical suspicion is imperative (specifically among COVID-19 patients with diabetes), and routine screening may be helpful.