Ligand and Structure-Based In Silico Determination of the Most Promising SARS-CoV-2 nsp16-nsp10 2′-o-Methyltransferase Complex Inhibitors among 3009 FDA Approved Drugs

As a continuation of our earlier work against SARS-CoV-2, seven FDA-approved drugs were designated as the best SARS-CoV-2 nsp16-nsp10 2′-o-methyltransferase (2′OMTase) inhibitors through 3009 compounds. The in silico inhibitory potential of the examined compounds against SARS-CoV-2 nsp16-nsp10 2′-o-methyltransferase (PDB ID: (6W4H) was conducted through a multi-step screening approach. At the beginning, molecular fingerprints experiment with SAM (S-Adenosylmethionine), the co-crystallized ligand of the targeted enzyme, unveiled the resemblance of 147 drugs. Then, a structural similarity experiment recommended 26 compounds. Therefore, the 26 compounds were docked against 2′OMTase to reveal the potential inhibitory effect of seven promising compounds (Protirelin, (1187), Calcium folinate (1913), Raltegravir (1995), Regadenoson (2176), Ertapenem (2396), Methylergometrine (2532), and Thiamine pyrophosphate hydrochloride (2612)). Out of the docked ligands, Ertapenem (2396) showed an ideal binding mode like that of the co-crystallized ligand (SAM). It occupied all sub-pockets of the active site and bound the crucial amino acids. Accordingly, some MD simulation experiments (RMSD, RMSF, R_g, SASA, and H-bonding) have been conducted for the 2′OMTase—Ertapenem complex over 100 ns. The performed MD experiments verified the correct binding mode of Ertapenem against 2′OMTase exhibiting low energy and optimal dynamics. Finally, MM-PBSA studies indicated that Ertapenem bonded advantageously to the targeted protein with a free energy value of −43 KJ/mol. Furthermore, the binding free energy analysis revealed the essential amino acids of 2′OMTase that served positively to the binding. The achieved results bring hope to find a treatment for COVID-19 via in vitro and in vivo studies for the pointed compounds.     

In Silico Exploration of Potential Natural Inhibitors against SARS-Cov-2 nsp10

In continuation of our previous effort, different in silico selection methods were applied to 310 naturally isolated metabolites that exhibited antiviral potentialities before. The applied selection methods aimed to pick the most relevant inhibitor of SARS-CoV-2 nsp10. At first, a structural similarity study against the co-crystallized ligand, S-Adenosyl Methionine (SAM), of SARS-CoV-2 nonstructural protein (nsp10) (PDB ID: 6W4H) was carried out. The similarity analysis culled 30 candidates. Secondly, a fingerprint study against SAM preferred compounds 44, 48, 85, 102, 105, 182, 220, 221, 282, 284, 285, 301, and 302. The docking studies picked 48, 182, 220, 221, and 284. While the ADMET analysis expected the likeness of the five candidates to be drugs, the toxicity study preferred compounds 48 and 182. Finally, a density-functional theory (DFT) study suggested vidarabine (182) to be the most relevant SARS-Cov-2 nsp10 inhibitor.     

Structural Homology-Based Drug Repurposing Approach for Targeting NSP12 SARS-CoV-2

The severe acute respiratory syndrome coronavirus 2, also known as SARS-CoV-2, is the causative agent of the COVID-19 global pandemic. SARS-CoV-2 has a highly conserved non-structural protein 12 (NSP-12) involved in RNA-dependent RNA polymerase (RdRp) activity. For the identification of potential inhibitors for NSP-12, computational approaches such as the identification of homologous proteins that have been previously targeted by FDA-approved antivirals can be employed. Herein, homologous proteins of NSP-12 were retrieved from Protein DataBank (PDB) and the evolutionary conserved sequence and structure similarity of the active site of the RdRp domain of NSP-12 was characterized. The identified homologous structures of NSP-12 belonged to four viral families: Coronaviridae, Flaviviridae, Picornaviridae, and Caliciviridae, and shared evolutionary conserved relationships. The multiple sequences and structural alignment of homologous structures showed highly conserved amino acid residues that were located at the active site of the RdRp domain of NSP-12. The conserved active site of the RdRp domain of NSP-12 was evaluated for binding affinity with the FDA-approved antivirals, i.e., Sofosbuvir and Dasabuvir in a molecular docking study. The molecular docking of Sofosbuvir and Dasabuvir with the active site that contains conserved motifs (motif A-G) of the RdRp domain of NSP-12 revealed significant binding affinity. Furthermore, MD simulation also inferred the potency of Sofosbuvir and Dasabuvir. In conclusion, targeting the active site of the RdRp domain of NSP-12 with Dasabuvir and Sofosbuvir might reduce viral replication and pathogenicity and could be further studied for the treatment of SARS-CoV-2.     

Marine Brominated Tyrosine Alkaloids as Promising Inhibitors of SARS-CoV-2

There have been more than 150 million confirmed cases of SARS-CoV-2 since the beginning of the pandemic in 2019. By June 2021, the mortality from such infections approached 3.9 million people. Despite the availability of a number of vaccines which provide protection against this virus, the evolution of new viral variants, inconsistent availability of the vaccine around the world, and vaccine hesitancy, in some countries, makes it unreasonable to rely on mass vaccination alone to combat this pandemic. Consequently, much effort is directed to identifying potential antiviral treatments. Marine brominated tyrosine alkaloids are recognized to have antiviral potential. We test here the antiviral capacity of fourteen marine brominated tyrosine alkaloids against five different target proteins from SARS-CoV-2, including main protease (M^pro) (PDB ID: 6lu7), spike glycoprotein (PDB ID: 6VYB), nucleocapsid phosphoprotein (PDB ID: 6VYO), membrane glycoprotein (PDB ID: 6M17), and non-structural protein 10 (nsp10) (PDB ID: 6W4H). These marine alkaloids, particularly the hexabrominated compound, fistularin-3, shows promising docking interactions with predicted binding affinities (S-score = −7.78, −7.65, −6.39, −6.28, −8.84 Kcal/mol) for the main protease (M^pro) (PDB ID: 6lu7), spike glycoprotein (PDB ID: 6VYB), nucleocapsid phosphoprotein (PDB ID: 6VYO), membrane glycoprotein (PDB ID: 6M17), and non-structural protein 10 (nsp10) (PDB ID: 6W4H), respectively, where it forms better interactions with the protein pockets than the native interaction. It also shows promising molecular dynamics, pharmacokinetics, and toxicity profiles. As such, further exploration of the antiviral properties of fistularin-3 against SARS-CoV-2 is merited.     

Synthesis and Biochemical Evaluation of 8H-Indeno[1,2-d]thiazole Derivatives as Novel SARS-CoV-2 3CL Protease Inhibitors

The COVID-19 pandemic caused by SARS-CoV-2 is a global burden on human health and economy. The 3-Chymotrypsin-like cysteine protease (3CL^pro) becomes an attractive target for SARS-CoV-2 due to its important role in viral replication. We synthesized a series of 8H-indeno[1,2-d]thiazole derivatives and evaluated their biochemical activities against SARS-CoV-2 3CL^pro. Among them, the representative compound 7a displayed inhibitory activity with an IC₅₀ of 1.28 ± 0.17 μM against SARS-CoV-2 3CL^pro. Molecular docking of 7a against 3CL^pro was performed and the binding mode was rationalized. These preliminary results provide a unique prototype for the development of novel inhibitors against SARS-CoV-2 3CL^pro.     

Ensemble Docking Coupled to Linear Interaction Energy Calculations for Identification of Coronavirus Main Protease (3CLpro) Non-Covalent Small-Molecule Inhibitors

SARS-CoV-2, or severe acute respiratory syndrome coronavirus 2, represents a new strain of Coronaviridae. In the closing 2019 to early 2020 months, the virus caused a global pandemic of COVID-19 disease. We performed a virtual screening study in order to identify potential inhibitors of the SARS-CoV-2 main viral protease (3CL^pro or M^pro). For this purpose, we developed a novel approach using ensemble docking high-throughput virtual screening directly coupled with subsequent Linear Interaction Energy (LIE) calculations to maximize the conformational space sampling and to assess the binding affinity of identified inhibitors. A large database of small commercial compounds was prepared, and top-scoring hits were identified with two compounds singled out, namely 1-[(R)-2-(1,3-benzimidazol-2-yl)-1-pyrrolidinyl]-2-(4-methyl-1,4-diazepan-1-yl)-1-ethanone and [({(S)-1-[(1H-indol-2-yl)methyl]-3-pyrrolidinyl}methyl)amino](5-methyl-2H-pyrazol-3-yl)formaldehyde. Moreover, we obtained a favorable binding free energy of the identified compounds, and using contact analysis we confirmed their stable binding modes in the 3CL^pro active site. These compounds will facilitate further 3CL^pro inhibitor design.     

Establishing an Analogue Based In Silico Pipeline in the Pursuit of Novel Inhibitory Scaffolds against the SARS Coronavirus 2 Papain-Like Protease

The ongoing coronavirus pandemic has been a burden on the worldwide population, with mass fatalities and devastating socioeconomic consequences. It has particularly drawn attention to the lack of approved small-molecule drugs to inhibit SARS coronaviruses. Importantly, lessons learned from the SARS outbreak of 2002–2004, caused by severe acute respiratory syndrome coronavirus 1 (SARS-CoV-1), can be applied to current drug discovery ventures. SARS-CoV-1 and SARS-CoV-2 both possess two cysteine proteases, the main protease (M^pro) and the papain-like protease (PL^pro), which play a significant role in facilitating viral replication, and are important drug targets. The non-covalent inhibitor, GRL-0617, which was found to inhibit replication of SARS-CoV-1, and more recently SARS-CoV-2, is the only PL^pro inhibitor co-crystallised with the recently solved SARS-CoV-2 PL^pro crystal structure. Therefore, the GRL-0617 structural template and pharmacophore features are instrumental in the design and development of more potent PL^pro inhibitors. In this work, we conducted scaffold hopping using GRL-0617 as a reference to screen over 339,000 ligands in the chemical space using the ChemDiv, MayBridge, and Enamine screening libraries. Twenty-four distinct scaffolds with structural and electrostatic similarity to GRL-0617 were obtained. These proceeded to molecular docking against PL^pro using the AutoDock tools. Of two compounds that showed the most favourable predicted binding affinities to the target site, as well as comparable protein-ligand interactions to GRL-0617, one was chosen for further analogue-based work. Twenty-seven analogues of this compound were further docked against the PL^pro, which resulted in two additional hits with promising docking profiles. Our in silico pipeline consisted of an integrative four-step approach: (1) ligand-based virtual screening (scaffold-hopping), (2) molecular docking, (3) an analogue search, and, (4) evaluation of scaffold drug-likeness, to identify promising scaffolds and eliminate those with undesirable properties. Overall, we present four novel, and lipophilic, scaffolds obtained from an exhaustive search of diverse and uncharted regions of chemical space, which may be further explored in vitro through structure-activity relationship (SAR) studies in the search for more potent inhibitors. Furthermore, these scaffolds were predicted to have fewer off-target interactions than GRL-0617. Lastly, to our knowledge, this work contains the largest ligand-based virtual screen performed against GRL-0617.     

Identification of Persuasive Antiviral Natural Compounds for COVID-19 by Targeting Endoribonuclease NSP15: A Structural-Bioinformatics Approach

SARS-CoV-2 is a positive-stranded RNA virus that bundles its genomic material as messenger-sense RNA in infectious virions and replicates these genomes through RNA intermediates. Several virus-encoded nonstructural proteins play a key role during the viral life cycle. Endoribonuclease NSP15 is vital for the replication and life cycle of the virus, and is thus considered a compelling druggable target. Here, we performed a combination of multiscoring virtual screening and molecular docking of a library of 1624 natural compounds (Nuclei of Bioassays, Ecophysiology and Biosynthesis of Natural Products (NuBBE) database) on the active sites of NSP15 (PDB:6VWW). After sequential high-throughput screening by LibDock and GOLD, docking optimization by CDOCKER, and final scoring by calculating binding energies, top-ranked compounds NuBBE-1970 and NuBBE-242 were further investigated via an indepth molecular-docking and molecular-dynamics simulation of 60 ns, which revealed that the binding of these two compounds with active site residues of NSP15 was sufficiently strong and stable. The findings strongly suggest that further optimization and clinical investigations of these potent compounds may lead to effective SARS-CoV-2 treatment.     

Identification of SARS-CoV-2 Receptor Binding Inhibitors by In Vitro Screening of Drug Libraries

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the coronavirus disease 2019 (COVID-19) global pandemic. The first step of viral infection is cell attachment, which is mediated by the binding of the SARS-CoV-2 receptor binding domain (RBD), part of the virus spike protein, to human angiotensin-converting enzyme 2 (ACE2). Therefore, drug repurposing to discover RBD-ACE2 binding inhibitors may provide a rapid and safe approach for COVID-19 therapy. Here, we describe the development of an in vitro RBD-ACE2 binding assay and its application to identify inhibitors of the interaction of the SARS-CoV-2 RBD to ACE2 by the high-throughput screening of two compound libraries (LOPAC^®1280 and DiscoveryProbeTM). Three compounds, heparin sodium, aurintricarboxylic acid (ATA), and ellagic acid, were found to exert an effective binding inhibition, with IC50 values ranging from 0.6 to 5.5 µg/mL. A plaque reduction assay in Vero E6 cells infected with a SARS-CoV-2 surrogate virus confirmed the inhibition efficacy of heparin sodium and ATA. Molecular docking analysis located potential binding sites of these compounds in the RBD. In light of these findings, the screening system described herein can be applied to other drug libraries to discover potent SARS-CoV-2 inhibitors.     

Repositioning of Ligands That Target the Spike Glycoprotein as Potential Drugs for SARS-CoV-2 in an In Silico Study

The worldwide health emergency of the SARS-CoV-2 pandemic and the absence of a specific treatment for this new coronavirus have led to the use of computational strategies (drug repositioning) to search for treatments. The aim of this work is to identify FDA (Food and Drug Administration)-approved drugs with the potential for binding to the spike structural glycoprotein at the hinge site, receptor binding motif (RBM), and fusion peptide (FP) using molecular docking simulations. Drugs that bind to amino acids are crucial for conformational changes, receptor recognition, and fusion of the viral membrane with the cell membrane. The results revealed some drugs that bind to hinge site amino acids (varenicline, or steroids such as betamethasone while other drugs bind to crucial amino acids in the RBM (naldemedine, atovaquone, cefotetan) or FP (azilsartan, maraviroc, and difluprednate); saquinavir binds both the RBM and the FP. Therefore, these drugs could inhibit spike glycoprotein and prevent viral entry as possible anti-COVID-19 drugs. Several drugs are in clinical studies; by focusing on other pharmacological agents (candesartan, atovaquone, losartan, maviroc and ritonavir) in this work we propose an additional target: the spike glycoprotein. These results can impact the proposed use of treatments that inhibit the first steps of the virus replication cycle.     

Virtual Screening and Molecular Dynamics Simulation Study of Influenza Polymerase PB2 Inhibitors

Influenza A virus is the main cause of worldwide epidemics and annual influenza outbreaks in humans. In this study, a virtual screen was performed to identify compounds that interact with the PB2 cap-binding domain (CBD) of influenza A polymerase. A virtual screening workflow based on Glide docking was used to screen an internal database containing 8417 molecules, and then the output compounds were selected based on solubility, absorbance, and structural fingerprints. Of the 16 compounds selected for biological evaluation, six compounds were identified that rescued cells from H1N1 virus-mediated death at non-cytotoxic concentrations, with EC50 values ranging from 2.5–55.43 μM, and that could bind to the PB2 CBD of H1N1, with Kd values ranging from 0.081–1.53 μM. Molecular dynamics (MD) simulations of the docking complexes of our active compounds revealed that each compound had its own binding characteristics that differed from those of VX-787. Our active compounds have novel structures and unique binding modes with PB2 proteins, and are suitable to serve as lead compounds for the development of PB2 inhibitors. An analysis of the MD simulation also helped us to identify the dominant amino acid residues that play a key role in binding the ligand to PB2, suggesting that we should focus on increasing and enhancing the interaction between inhibitors and these major amino acids during lead compound optimization to obtain more active PB2 inhibitors.     

Olive-Derived Triterpenes Suppress SARS COV-2 Main Protease: A Promising Scaffold for Future Therapeutics

SARS CoV-2 pandemic is still considered a global health disaster, and newly emerged variants keep growing. A number of promising vaccines have been recently developed as a protective measure; however, cost-effective treatments are also of great importance to support this critical situation. Previously, betulinic acid has shown promising antiviral activity against SARS CoV via targeting its main protease. Herein, we investigated the inhibitory potential of this compound together with three other triterpene congeners (i.e., ursolic acid, maslinic acid, and betulin) derived from olive leaves against the viral main protease (M^pro) of the currently widespread SARS CoV-2. Interestingly, betulinic, ursolic, and maslinic acids showed significant inhibitory activity (IC₅₀ = 3.22–14.55 µM), while betulin was far less active (IC₅₀ = 89.67 µM). A comprehensive in-silico analysis (i.e., ensemble docking, molecular dynamic simulation, and binding-free energy calculation) was then performed to describe the binding mode of these compounds with the enzyme catalytic active site and determine the main essential structural features required for their inhibitory activity. Results presented in this communication indicated that this class of compounds could be considered as a promising lead scaffold for developing cost-effective anti-SARS CoV-2 therapeutics.     

Genomic variation and point mutations analysis of Indian COVID-19 patient samples submitted in GISAID database

Corona virus disease 2019 (COVID-19) endemic has havoc on the world; the causative virus of the pandemic is SARS CoV-2. Pharmaceutical companies and academic institutes are in continuous efforts to identify anti-viral therapy or vaccines, but the most significant challenge faced is the highly evolving genome of SARS CoV-2, which is imparting evolutionary selective benefits to the virus. To understand the viral mutations, we have retrieved nine hundred and thirty-four samples from different states of India via the GISAID database and analyzed the frequency of all types of point mutation in all structural, non-structural proteins, and accessory factors of SARS CoV-2. Spike glycol protein, nsp3, nsp6, nsp12, N and NS3 were the most evolving proteins. High frequency point mutations were Q496P (nsp2), A380V (nsp4), A994D (nsp3), L37F (nsp6), P323L & A97V (nsp12), Q57H (ns3), D614G (S), P13L (N), R203K (N), G204R (N) and S194L (N).     

In Silico Studies of Some Isoflavonoids as Potential Candidates against COVID-19 Targeting Human ACE2 (hACE2) and Viral Main Protease (Mpro)

The Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) caused the “COVID-19” disease that has been declared by WHO as a global emergency. The pandemic, which emerged in China and widespread all over the world, has no specific treatment till now. The reported antiviral activities of isoflavonoids encouraged us to find out its in silico anti-SARS-CoV-2 activity. In this work, molecular docking studies were carried out to investigate the interaction of fifty-nine isoflavonoids against hACE2 and viral M^pro. Several other in silico studies including physicochemical properties, ADMET and toxicity have been preceded. The results revealed that the examined isoflavonoids bound perfectly the hACE-2 with free binding energies ranging from −24.02 to −39.33 kcal mol⁻¹, compared to the co-crystallized ligand (−21.39 kcal mol^–1). Furthermore, such compounds bound the M^pro with unique binding modes showing free binding energies ranging from −32.19 to −50.79 kcal mol^–1, comparing to the co-crystallized ligand (binding energy = −62.84 kcal mol^–1). Compounds 33 and 56 showed the most acceptable affinities against hACE2. Compounds 30 and 53 showed the best docking results against M^pro. In silico ADMET studies suggest that most compounds possess drug-likeness properties.     

Synthesis and Antimicrobial Activity Screening of Piperazines Bearing N,N′-Bis(1,3,4-thiadiazole) Moiety as Probable Enoyl-ACP Reductase Inhibitors

A new N,N′-disubstituted piperazine conjugated with 1,3,4-thiadiazole and 1,2,4-triazole was prepared and the chemical structures were identified by IR, NMR and elemental analysis. All the prepared compounds were tested for their antimicrobial activity. The antimicrobial results indicated that the tested compounds showed significant antibacterial activity against gram-negative strains, especially E. coli, relative to gram-positive bacteria. Docking analysis was performed to support the biological results; binding modes with the active site of enoyl reductase amino acids from E. coli showed very good scores, ranging from −6.1090 to −9.6184 kcal/mol. Correlation analysis was performed for the inhibition zone (nm) and the docking score.     

Synthesis and molecular docking study of novel COVID-19 inhibitors

In 2020, the world tried to combat the corona virus (COVID-19) pandemic. A proven treatment method specific to Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) is still not found. In this study, seven new antiviral compounds were designed for COVID-19 treatment. The ability of these compounds to inhibit COVID-19’s RNA processing was calculated by the molecular docking study. It has been observed that the compounds can have high binding affinities especially against NSP12 (between -9.06 and -8.00 kcal/mol). The molecular dynamics simulation of NSP12-ZG 7 complex proved the stability of interaction. The synthesis of two most active molecules was performed by one-pot reaction and characterized by FT-IR, ¹H-NMR, ¹³C-NMR, and mass spectroscopy. The compounds presented with their synthesis are inhibitory core structures against SARS-CoV-2 infection.     

In Vitro Evaluation of In Silico Screening Approaches in Search for Selective ACE2 Binding Chemical Probes

New models for ACE2 receptor binding, based on QSAR and docking algorithms were developed, using XRD structural data and ChEMBL 26 database hits as training sets. The selectivity of the potential ACE2-binding ligands towards Neprilysin (NEP) and ACE was evaluated. The Enamine screening collection (3.2 million compounds) was virtually screened according to the above models, in order to find possible ACE2-chemical probes, useful for the study of SARS-CoV2-induced neurological disorders. An enzymology inhibition assay for ACE2 was optimized, and the combined diversified set of predicted selective ACE2-binding molecules from QSAR modeling, docking, and ultrafast docking was screened in vitro. The in vitro hits included two novel chemotypes suitable for further optimization.