Repurposing of FDA-approved drugs against active site and potential allosteric drug-binding sites of COVID-19 main protease

In Silico Identification of Putative Allosteric Pockets and Inhibitors for the KRASG13D-SOS1 Complex in Cancer Therapy

RAS mutations occur in about 30% of human cancers, leading to enhanced RAS signaling and tumor growth. KRAS is the most commonly mutated oncogene in human tumors, especially lung, pancreatic, and colorectal cancers. Direct targeting of KRAS is difficult due to its highly conserved sequence; but, its complex with the guanine nucleotide exchange factor Son of Sevenless (SOS) 1 promises an attractive target for inhibiting RAS-mediated signaling. Here, we first revealed putative allosteric binding sites of the SOS1, KRASG12C-SOS1 complex, and the ternary KRASG13D-SOS1 complex structures using two network-based models, the essential site scanning analysis and the residue interaction network model. The results enabled us to identify two new putative allosteric pockets for the ternary KRASG13D-SOS1 complex. These were then screened together with the known ligand binding site against the natural compounds in the InterBioScreen (IBS) database using the Glide software package developed by Schrödinger, Inc. The docking poses of seven hit compounds were assessed using 400 ns long molecular dynamics (MD) simulations with two independent replicas using Desmond, coupled with thermal MM-GBSA calculations for the estimation of the binding free energy values. The structural skeleton of the seven proposed compounds consists of different functional groups and heterocyclic rings that possess anti-cancer activity and exhibit persistent interactions with key residues in binding pockets throughout the MD simulations. STOCK1N-09823 was determined as the most promising hit that promoted the disruption of the interactions R73 (chain A)/N879 and R73 (chain A)/Y884, which are key for SOS1-mediated KRAS activation.

Computational Assessment of Clinical Drugs against SARS-CoV-2: Foreseeing Molecular Mechanisms and Potent Mpro Inhibitors

The emergence of new SARS-CoV-2 variants of concern (VOC) is a propulsion for accelerated potential therapeutic discovery. SARS-CoV-2's main protease (Mpro), essential for host cell viral replication, is a pre-eminent druggable protein target. Here, we perform extensive drug re-profiling of the comprehensive Excelra database, which compiles various under-trial drug candidates for COVID-19 treatment. For mechanistic understanding, the most promising screened-out molecules with targets are subjected to molecular dynamics (MD) simulations. Post-MD analyses demonstrate Darunavir, Ponatinib, and Tomivosertib forming a stable complex with Mpro, characterized by less fluctuation of Cα atoms, smooth and stable root-mean-square deviation (RMSD), and robust contact with the active site residues. Likewise, they all have lower binding free energy with Mpro, demonstrating strong affinity. In free energy landscape profiles, the distances from His41 and Cys145 exhibit a single energy minima basin, implying their preponderance in proximity to Mpro's catalytic dyad. Overall, the computational assessment earmarks promising candidates from the Excelra database, emphasizing on carrying out exhaustive biochemical experiments along with clinical trials. The work lays the foundation for potential therapeutic interventions in treating COVID-19.

Unveiling the mechanism of action of a novel natural dual inhibitor of SARS-CoV-2 Mpro and PLpro with molecular dynamics simulations

In the twenty-first century, we have witnessed multiple coronavirus pandemics. Despite declining SARS-CoV-2 cases, continued research remains vital. We report the discovery of sydowiol B, a natural product, as a dual inhibitor of SARS-CoV-2 main protease (Mpro) and papain-like protease (PLpro). Sydowiol B interacts with the nano-channel at the Mpro dimer interface and the PLpro active site. Molecular dynamics simulations suggest that sydowiol B inhibits Mpro by limiting active site expansion rather than inducing collapse. Furthermore, sydowiol B binding may amplify the fluctuation of two loops coordinating with the structural Zn2+ in PLpro, displacing Zn2+ from the zinc finger domain to the S2 helix. Sydowiol B and its analogue, violaceol I, exhibit broad-spectrum antiviral activity against homologous coronaviruses. Given the conservation of Mpro and PLpro, sydowiol B and violaceol I are promising leads for designing and developing anti-coronavirus therapies.

Toward the Design of Allosteric Effectors: Gaining Comprehensive Control of Drug Properties and Actions

While the therapeutic potential of allosteric drugs is increasingly realized, the discovery of effectors is largely incidental. The rational design of allosteric effectors requires new state-of-the-art approaches to account for the distinct characteristics of allosteric ligands and their modes of action. We present a broadly applicable computational framework for obtaining allosteric site-effector pairs, providing targeted, highly specific, and tunable regulation to any functional site. We validated the framework using the main protease from SARS-CoV-2 and the K-RasG12D oncoprotein. High-throughput per-residue quantification of the energetics of allosteric signaling and effector binding revealed known drugs capable of inducing the required modulation upon binding. Starting from fragments of known well-characterized drugs, allosteric effectors and binding sites were designed and optimized simultaneously to achieve targeted and specific signaling to distinct functional sites, such as, for example, the switch regions of K-RasG12D. The generic framework proposed in this work will be instrumental in developing allosteric therapies aligned with a precision medicine approach.

Exploring Druggable Binding Sites on the Class A GPCRs Using the Residue Interaction Network and Site Identification by Ligand Competitive Saturation

G protein-coupled receptors (GPCRs) play a central role in cellular signaling and are linked to many diseases. Accordingly, computational methods to explore potential allosteric sites for this class of proteins to facilitate the identification of potential modulators are needed. Importantly, the availability of rich structural data providing the locations of the orthosteric ligands and allosteric modulators targeting different GPCRs allows for the validation of approaches to identify new allosteric binding sites. Here, we validate the combination of two computational techniques, the residue interaction network (RIN) model and the site identification by ligand competitive saturation (SILCS) method, to predict putative allosteric binding sites of class A GPCRs. RIN analysis identifies hub residues that mediate allosteric signaling within a receptor and have a high capacity to alter receptor dynamics upon ligand binding. The known orthosteric (and allosteric) binding sites of 18 distinct class A GPCRs were successfully predicted by RIN through a dataset of 105 crystal structures (91 ligand-bound, 14 unbound) with up to 77.8% (76.9%) sensitivity, 92.5% (95.3%) specificity, 51.9% (50%) precision, and 86.2% (92.4%) accuracy based on the experimental and theoretical binding site data. Moreover, graph spectral analysis of the residue networks revealed that the proposed sites were located at the interfaces of highly interconnected residue clusters with a high ability to coordinate the functional dynamics. Then, we employed the SILCS-Hotspots method to assess the druggability of the novel sites predicted for 7 distinct class A GPCRs that are critical for a variety of diseases. While the known orthosteric and allosteric binding sites are successfully explored by our approach, numerous putative allosteric sites with the potential to bind drug-like molecules are proposed. The computational approach presented here promises to be a highly effective tool to predict putative allosteric sites of GPCRs to facilitate the design of effective modulators.

A potential allosteric inhibitor of SARS-CoV-2 main protease (Mpro) identified through metastable state analysis

Anti-COVID19 drugs, such as nirmatrelvir, have been developed targeting the SARS-CoV-2 main protease, Mpro, based on the critical requirement of its proteolytic processing of the viral polyproteins into functional proteins essential for viral replication. However, the emergence of SARS-CoV-2 variants with Mpro mutations has raised the possibility of developing resistance against these drugs, likely due to therapeutic targeting of the Mpro catalytic site. An alternative to these drugs is the development of drugs that target an allosteric site distant from the catalytic site in the protein that may reduce the chance of the emergence of resistant mutants. Here, we combine computational analysis with in vitro assay and report the discovery of a potential allosteric site and an allosteric inhibitor of SARS-CoV-2 Mpro. Specifically, we identified an Mpro metastable state with a deformed catalytic site harboring potential allosteric sites, raising the possibility that stabilization of this metastable state through ligand binding can lead to the inhibition of Mpro activity. We then performed a computational screening of a library (∼4.2 million) of drug-like compounds from the ZINC database and identified several candidate molecules with high predicted binding affinity. MD simulations showed stable binding of the three top-ranking compounds to the putative allosteric sites in the protein. Finally, we tested the three compounds in vitro using a BRET-based Mpro biosensor and found that one of the compounds (ZINC4497834) inhibited the Mpro activity. We envisage that the identification of a potential allosteric inhibitor of Mpro will aid in developing improved anti-COVID-19 therapy.

A computational workflow to determine drug candidates alternative to aminoglycosides targeting the decoding center of E. coli ribosome

The global antibiotic resistance problem necessitates fast and effective approaches to finding novel inhibitors to treat bacterial infections. In this study, we propose a computational workflow to identify plausible high-affinity compounds from FDA-approved, investigational, and experimental libraries for the decoding center on the small subunit 30S of the E. coli ribosome. The workflow basically consists of two molecular docking calculations on the intact 30S, followed by molecular dynamics (MD) simulations coupled with MM-GBSA calculations on a truncated ribosome structure. The parameters used in the molecular docking suits, Glide and AutoDock Vina, as well as in the MD simulations with Desmond were carefully adjusted to obtain expected interactions for the ligand-rRNA complexes. A filtering procedure was followed, considering a fingerprint based on aminoglycoside's binding site on the 30S to obtain seven hit compounds either with different clinical usages or aminoglycoside derivatives under investigation, suggested for in vitro studies. The detailed workflow developed in this study promises an effective and fast approach for the estimation of binding free energies of large protein-RNA and ligand complexes.

Identifying and Assessing Putative Allosteric Sites and Modulators for CXCR4 Predicted through Network Modeling and Site Identification by Ligand Competitive Saturation

The chemokine receptor CXCR4 is a critical target for the treatment of several cancer types and HIV-1 infections. While orthosteric and allosteric modulators have been developed targeting its extracellular or transmembrane regions, the intramembrane region of CXCR4 may also include allosteric binding sites suitable for the development of allosteric drugs. To investigate this, we apply the Gaussian Network Model (GNM) to the monomeric and dimeric forms of CXCR4 to identify residues essential for its local and global motions located in the hinge regions of the protein. Residue interaction network (RIN) analysis suggests hub residues that participate in allosteric communication throughout the receptor. Mutual residues from the network models reside in regions with a high capacity to alter receptor dynamics upon ligand binding. We then investigate the druggability of these potential allosteric regions using the site identification by ligand competitive saturation (SILCS) approach, revealing two putative allosteric sites on the monomer and three on the homodimer. Two screening campaigns with Glide and SILCS-Monte Carlo docking using FDA-approved drugs suggest 20 putative hit compounds including antifungal drugs, anticancer agents, HIV protease inhibitors, and antimalarial drugs. In vitro assays considering mAB 12G5 and CXCL12 demonstrate both positive and negative allosteric activities of these compounds, supporting our computational approach. However, in vivo functional assays based on the recruitment of β-arrestin to CXCR4 do not show significant agonism and antagonism at a single compound concentration. The present computational pipeline brings a new perspective to computer-aided drug design by combining conformational dynamics based on network analysis and cosolvent analysis based on the SILCS technology to identify putative allosteric binding sites using CXCR4 as a showcase.

Evaluations of FDA-approved Drugs Targeting 3CLP of SARS-CoV-2 Employing a Repurposing Strategy

BACKGROUND

The SARS-CoV-2 coronavirus (COVID-19) has raised innumerable global concerns, and few effective treatment strategies have yet been permitted by the FDA to lighten the disease burden. SARS-CoV-2 3C-like proteinase (3CLP) is a crucial protease and plays a key role in the viral life cycle, as it controls replication, and thus, it is viewed as a target for drug design.

METHODS

In this study, we performed structure-based virtual screening of FDA drugs approved during 2015-2019 (a total of 220 drugs) for interaction with the active site of 3CLP (PDB ID 6LU7) using AutoDock 4.2. We report the top ten drugs that outperform the reported drugs against 3CLP (Elbasvir and Nelfinavir), particularly Cefiderocol, having the highest affinity among the compounds tested, with a binding energy of -9.97 kcal/mol. H-bond (LYS102:HZ2-ligand: O49), hydrophobic (ligand-VAL104), and electrostatic (LYS102:NZ-ligand: O50) interactions were observed in the cefiderocol-3CLP complex. The docked complex was subjected to a 50 ns molecular dynamics study to check its stability, and stable RMSD and RMSF graphs were observed.

RESULTS

Accordingly, we suggest cefiderocol might be effective against SARS-CoV-2 and urge that experimental validation be performed to determine the antiviral efficacy of cefiderocol against SARS-CoV-2.

DISCUSSION

Along with these, cefiderocol is effective for treating respiratory tract pathogens and a wide range of gram-negative bacteria for whom there are limited therapeutic alternatives.

CONCLUSION

This article aimed to explore the FDA-approved drugs as a repurposing study against 3CLP for COVID-19 management.

Antiviral Effect of 5'-Arylchalcogeno-3-aminothymidine Derivatives in SARS-CoV-2 Infection

The understanding that zidovudine (ZDV or azidothymidine, AZT) inhibits the RNA-dependent RNA polymerase (RdRp) of SARS-CoV-2 and that chalcogen atoms can increase the bioactivity and reduce the toxicity of AZT has directed our search for the discovery of novel potential anti-coronavirus compounds. Here, the antiviral activity of selenium and tellurium containing AZT derivatives in human type II pneumocytes cell model (Calu-3) and monkey kidney cells (Vero E6) infected with SARS-CoV-2, and their toxic effects on these cells, was evaluated. Cell viability analysis revealed that organoselenium (R3a-R3e) showed lower cytotoxicity than organotellurium (R3f, R3n-R3q), with CC50 ≥ 100 µM. The R3b and R3e were particularly noteworthy for inhibiting viral replication in both cell models and showed better selectivity index. In Vero E6, the EC50 values for R3b and R3e were 2.97 ± 0.62 µM and 1.99 ± 0.42 µM, respectively, while in Calu-3, concentrations of 3.82 ± 1.42 µM and 1.92 ± 0.43 µM (24 h treatment) and 1.33 ± 0.35 µM and 2.31 ± 0.54 µM (48 h) were observed, respectively. The molecular docking calculations were carried out to main protease (Mpro), papain-like protease (PLpro), and RdRp following non-competitive, competitive, and allosteric inhibitory approaches. The in silico results suggested that the organoselenium is a potential non-competitive inhibitor of RdRp, interacting in the allosteric cavity located in the palm region. Overall, the cell-based results indicated that the chalcogen-zidovudine derivatives were more potent than AZT in inhibiting SARS-CoV-2 replication and that the compounds R3b and R3e play an important inhibitory role, expanding the knowledge about the promising therapeutic capacity of organoselenium against COVID-19.

Repurposing clinically available drugs and therapies for pathogenic targets to combat SARS-CoV-2

The coronavirus disease 2019 (COVID-19) pandemic has affected a large portion of the global population, both physically and mentally. Current evidence suggests that the rapidly evolving coronavirus subvariants risk rendering vaccines and antibodies ineffective due to their potential to evade existing immunity, with enhanced transmission activity and higher reinfection rates that could lead to new outbreaks across the globe. The goal of viral management is to disrupt the viral life cycle as well as to relieve severe symptoms such as lung damage, cytokine storm, and organ failure. In the fight against viruses, the combination of viral genome sequencing, elucidation of the structure of viral proteins, and identifying proteins that are highly conserved across multiple coronaviruses has revealed many potential molecular targets. In addition, the time- and cost-effective repurposing of preexisting antiviral drugs or approved/clinical drugs for these targets offers considerable clinical advantages for COVID-19 patients. This review provides a comprehensive overview of various identified pathogenic targets and pathways as well as corresponding repurposed approved/clinical drugs and their potential against COVID-19. These findings provide new insight into the discovery of novel therapeutic strategies that could be applied to the control of disease symptoms emanating from evolving SARS-CoV-2 variants.

Exploring the Potential of Black Soldier Fly Larval Proteins as Bioactive Peptide Sources through in Silico Gastrointestinal Proteolysis: A Cheminformatic Investigation

Despite their potential as a protein source for human consumption, the health benefits of black soldier fly larvae (BSFL) proteins following human gastrointestinal (GI) digestion are poorly understood. This computational study explored the potential of BSFL proteins to release health-promoting peptides after human GI digestion. Twenty-six proteins were virtually proteolyzed with GI proteases. The resultant peptides were screened for high GI absorption and non-toxicity. Shortlisted peptides were searched against the BIOPEP-UWM and Scopus databases to identify their bioactivities. The potential of the peptides as inhibitors of myeloperoxidase (MPO), NADPH oxidase (NOX), and xanthine oxidase (XO), as well as a disruptor of Keap1–Nrf2 protein–protein interaction, were predicted using molecular docking and dynamics simulation. Our results revealed that about 95% of the 5218 fragments generated from the proteolysis of BSFL proteins came from muscle proteins. Dipeptides comprised the largest group (about 25%) of fragments arising from each muscular protein. Screening of 1994 di- and tripeptides using SwissADME and STopTox tools revealed 65 unique sequences with high GI absorption and non-toxicity. A search of the databases identified 16 antioxidant peptides, 14 anti-angiotensin-converting enzyme peptides, and 17 anti-dipeptidyl peptidase IV peptides among these sequences. Results from molecular docking and dynamic simulation suggest that the dipeptide DF has the potential to inhibit Keap1–Nrf2 interaction and interact with MPO within a short time frame, whereas the dipeptide TF shows promise as an XO inhibitor. BSFL peptides were likely weak NOX inhibitors. Our in silico results suggest that upon GI digestion, BSFL proteins may yield high-GI-absorbed and non-toxic peptides with potential health benefits. This study is the first to investigate the bioactivity of peptides liberated from BSFL proteins following human GI digestion. Our findings provide a basis for further investigations into the potential use of BSFL proteins as a functional food ingredient with significant health benefits.

Virtual Screening of FDA-Approved Drugs for Enhanced Binding with Mitochondrial Aldehyde Dehydrogenase

Mitochondrial aldehyde dehydrogenase (ALDH2) is a potential target for the treatment of substance use disorders such as alcohol addiction. Here, we adopted computational methods of molecular dynamics (MD) simulation, docking, and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) analysis to perform a virtual screening of FDA-approved drugs, hitting potent inhibitors against ALDH2. Using MD-derived conformations as receptors, butenafine (net charge q = +1 e) and olaparib (q = 0) were selected as promising compounds with a low toxicity and a binding strength equal to or stronger than previously reported potent inhibitors of daidzin and CVT-10216. A few negatively charged compounds were also hit from the docking with the Autodock Vina software, while the MM-PBSA analysis yielded positive binding energies (unfavorable binding) for these compounds, mainly owing to electrostatic repulsion in association with a negatively charged receptor (q = -6 e for ALDH2 plus the cofactor NAD⁺). This revealed a deficiency of the Vina scoring in dealing with strong charge-charge interactions between binding partners, due to its built-in protocol of not using atomic charges for electrostatic interactions. These observations indicated a requirement of further verification using MD and/or MM-PBSA after docking prediction. The identification of key residues for the binding implied that the receptor residues at the bottom and entrance of the substrate-binding hydrophobic tunnel were able to offer additional interactions with different inhibitors such as π-π, π-alkyl, van der Waals contacts, and polar interactions, and that the rational use of these interactions is beneficial to the design of potent inhibitors against ALDH2.

Agathisflavone, a natural biflavonoid that inhibits SARS-CoV-2 replication by targeting its proteases

Despite the fast development of vaccines, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) still circulates through variants of concern (VoC) and escape the humoral immune response. SARS-CoV-2 has provoked over 200,000 deaths/months since its emergence and only a few antiviral drugs showed clinical benefit up to this moment. Thus, chemical structures endowed with anti-SARS-CoV-2 activity are important for continuous antiviral development and natural products represent a fruitful source of substances with biological activity. In the present study, agathisflavone (AGT), a biflavonoid from Anacardium occidentale was investigated as a candidate anti-SARS-CoV-2 compound. In silico and enzymatic analysis indicated that AGT may target mainly the viral main protease (M^pro) and not the papain-like protease (PL^pro) in a non-competitive way. Cell-based assays in type II pneumocytes cell lineage (Calu-3) showed that SARS-CoV-2 is more susceptible to AGT than to apigenin (APG, monomer of AGT), in a dose-dependent manner, with an EC₅₀ of 4.23 ± 0.21 μM and CC₅₀ of 61.3 ± 0.1 μM and with a capacity to inhibit the level of pro-inflammatory mediator tumor necrosis factor-alpha (TNF-α). These results configure AGT as an interesting chemical scaffold for the development of novel semisynthetic antivirals against SARS-CoV-2.

Ivermectin: A Controversial Focal Point during the COVID-19 Pandemic

The SARS-CoV-2 pandemic has confirmed the apocalyptic predictions that virologists have been making for several decades. The challenge the world is facing is that of trying to find a possible treatment, and a viable and expedient option for addressing this challenge is the repurposing of drugs. However, in some cases, although these drugs are approved for use in humans, the mechanisms of action involved are unknown. In this sense, to justify its therapeutic application to a new disease, it is ideal, but not necessary, to know the basic mechanisms of action involved in a drug's biological effects. This review compiled the available information regarding the various effects attributed to Ivermectin. The controversy over its use for the treatment of COVID-19 is demonstrated by this report that considers the proposal unfeasible because the therapeutic doses proposed to achieve this effect cannot be achieved. However, due to the urgent need to find a treatment, an exhaustive and impartial review is necessary in order to integrate the knowledge that exists, to date, of the possible mechanisms through which the treatment may be helpful in defining safe doses and schedules of Ivermectin.

Integrated bioinformatics-cheminformatics approach toward locating pseudo-potential antiviral marine alkaloids against SARS-CoV-2-Mpro

The emergence of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) with the most contagious variants, alpha (B.1.1.7), beta (B.1.351), delta (B.1.617.2), and Omicron (B.1.1.529) has continuously added a higher number of morbidity and mortality, globally. The present integrated bioinformatics-cheminformatics approach was employed to locate potent antiviral marine alkaloids that could be used against SARS-CoV-2. Initially, 57 antiviral marine alkaloids and two repurposing drugs were selected from an extensive literature review. Then, the putative target enzyme SARS-CoV-2 main protease (SARS-CoV-2-Mpro) was retrieved from the protein data bank and carried out a virtual screening-cum-molecular docking study with all candidates using PyRx 0.8 and AutoDock 4.2 software. Further, the molecular dynamics (MD) simulation of the two most potential alkaloids and a drug docking complex at 100 ns (with two ligand topology files from PRODRG and ATB server, separately), the molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) free energy, and contributions of entropy were investigated. Then, the physicochemical-toxicity-pharmacokinetics-drug-likeness profiles, the frontier molecular orbitals energies (highest occupied molecular orbital, lowest unoccupied molecular orbital, and ΔE), and structural-activity relationship were assessed and analyzed. Based on binding energy, 8-hydroxymanzamine (-10.5 kcal/mol) and manzamine A (-10.1 kcal/mol) from all alkaloids with darunavir (-7.9 kcal/mol) and lopinavir (-7.4 kcal/mol) against SARS-CoV-2-Mpro were recorded. The MD simulation (RMSD, RMSF, Rg, H-bond, MM/PBSA binding energy) illustrated that the 8-hydroxymanzamine exhibits a static thermodynamic feature than the other two complexes. The predicted physicochemical, toxicity, pharmacokinetics, and drug-likeness profiles also revealed that the 8-hydroxymanzamine could be used as a potential lead candidate individually and/or synergistically with darunavir or lopinavir to combat SARS-CoV-2 infection after some pharmacological validation.

Molecular dynamics and docking studies on potentially active natural phytochemicals for targeting SARS-CoV-2 main protease

In the present study, we screened eighty seven novel phytochemical compounds from four popular herbs, such as, Aegle Marmelos, Coleus Amboinicus, Aerva Lanata and Biophytum Sensitivum and identified the best three for targeting the main protease (M^pro) receptor of SARS-CoV-2. After categorizing all the phytochemicals based upon LibDock scores and hydrogen bonding interactions, the top ranked 26 compounds were further subjected for detailed Molecular dynamics (MD) study. From these screening we identified that Aegelinosides B leads the list with a high LibDock value of 142.50 (binding energy: -8.54 kcal/mol), which is better than several popular reference compounds namely, Tipranavir (LibDock score, 141.50), Saquinavir (125.34), Zopicole (122.9), Pirenepine (122.70), (115.37), Metixene (109.18), Oxiconazole Pimozide (138.00) and Rimonabant (91.88). Detailed analysis for structural stability (RMSD), Cα fluctuations (RMSF), intermolecular hydrogen bond interactions, effect of solvent accessibility (SASA) and compactness (Rg) factors were performed for the best six compounds and it is found that they are very stable and exhibit folding behavior. Apart from the docking and MD tests, through further drug-likeness and toxicity tests, three compounds, such as, Aegelinosides B, Epicatechin, and Feruloyltyramine (LibDock scores, respectively, 142.50, 124.33 and 129.06) can be suggested for fighting SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Scope of repurposed drugs against the potential targets of the latest variants of SARS-CoV-2

The unprecedented outbreak of the severe acute respiratory syndrome (SARS) Coronavirus-2, across the globe, triggered a worldwide uproar in the search for immediate treatment strategies. With no specific drug and not much data available, alternative approaches such as drug repurposing came to the limelight. To date, extensive research on the repositioning of drugs has led to the identification of numerous drugs against various important protein targets of the coronavirus strains, with hopes of the drugs working against the major variants of concerns (alpha, beta, gamma, delta, omicron) of the virus. Advancements in computational sciences have led to improved scope of repurposing via techniques such as structure-based approaches including molecular docking, molecular dynamic simulations and quantitative structure activity relationships, network-based approaches, and artificial intelligence-based approaches with other core machine and deep learning algorithms. This review highlights the various approaches to repurposing drugs from a computational biological perspective, with various mechanisms of action of the drugs against some of the major protein targets of SARS-CoV-2. Additionally, clinical trials data on potential COVID-19 repurposed drugs are also highlighted with stress on the major SARS-CoV-2 targets and the structural effect of variants on these targets. The interaction modelling of some important repurposed drugs has also been elucidated. Furthermore, the merits and demerits of drug repurposing are also discussed, with a focus on the scope and applications of the latest advancements in repurposing.

Allosteric Binding Sites of the SARS-CoV-2 Main Protease: Potential Targets for Broad-Spectrum Anti-Coronavirus Agents

The current pandemic caused by the COVID-19 disease has reached everywhere in the world and has affected every aspect of our lives. As of the current data, the World Health Organization (WHO) has reported more than 300 million confirmed COVID-19 cases worldwide and more than 5 million deaths. M^pro is an enzyme that plays a key role in the life cycle of the SARS-CoV-2 virus, and it is vital for the disease progression. The M^pro enzyme seems to have several allosteric sites that can hinder the enzyme catalytic activity. Furthermore, some of these allosteric sites are located at or nearby the dimerization interface which is essential for the overall M^pro activity. In this review paper, we investigate the potential of the M^pro allosteric site to act as a drug target, especially since they interestingly appear to be resistant to mutation. The work is illustrated through three subsequent sections: First, the two main categories of M^pro allosteric sites have been explained and discussed. Second, a total of six pockets have been studied and evaluated for their druggability and cavity characteristics. Third, the experimental and computational attempts for the discovery of new allosteric inhibitors have been illustrated and discussed. To sum up, this review paper gives a detailed insight into the feasibility of developing new M^pro inhibitors to act as a potential treatment for the COVID-19 disease.

In Silico Screening of Novel TMPRSS2 Inhibitors for Treatment of COVID-19

COVID-19, a pandemic caused by the virus SARS-CoV-2, has spread globally, necessitating the search for antiviral compounds. Transmembrane protease serine 2 (TMPRSS2) is a cell surface protease that plays an essential role in SARS-CoV-2 infection. Therefore, researchers are searching for TMPRSS2 inhibitors that can be used for the treatment of COVID-19. As such, in this study, based on the crystal structure, we targeted the active site of TMPRSS2 for virtual screening of compounds in the FDA database. Then, we screened lumacaftor and ergotamine, which showed strong binding ability, using 100 ns molecular dynamics simulations to study the stability of the protein–ligand binding process, the flexibility of amino acid residues, and the formation of hydrogen bonds. Subsequently, we calculated the binding free energy of the protein–ligand complex by the MM-PBSA method. The results show that lumacaftor and ergotamine interact with residues around the TMPRSS2 active site, and reached equilibrium in the 100 ns molecular dynamics simulations. We think that lumacaftor and ergotamine, which we screened through in silico studies, can effectively inhibit the activity of TMPRSS2. Our findings provide a basis for subsequent in vitro experiments, having important implications for the development of effective anti-COVID-19 drugs.

Commercially Available Flavonols Are Better SARS-CoV-2 Inhibitors Than Isoflavone and Flavones

Despite the fast development of vaccines, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is still circulating and generating variants of concern (VoC) that escape the humoral immune response. In this context, the search for anti-SARS-CoV-2 compounds is still essential. A class of natural polyphenols known as flavonoids, frequently available in fruits and vegetables, is widely explored in the treatment of different diseases and used as a scaffold for the design of novel drugs. Therefore, herein we evaluate seven flavonoids divided into three subclasses, isoflavone (genistein), flavone (apigenin and luteolin) and flavonol (fisetin, kaempferol, myricetin, and quercetin), for COVID-19 treatment using cell-based assays and in silico calculations validated with experimental enzymatic data. The flavonols were better SARS-CoV-2 inhibitors than isoflavone and flavones. The increasing number of hydroxyl groups in ring B of the flavonols kaempferol, quercetin, and myricetin decreased the 50% effective concentration (EC₅₀) value due to their impact on the orientation of the compounds inside the target. Myricetin and fisetin appear to be preferred candidates; they are both anti-inflammatory (decreasing TNF-α levels) and inhibit SARS-CoV-2 mainly by targeting the processability of the main protease (M^pro) in a non-competitive manner, with a potency comparable to the repurposed drug atazanavir. However, fisetin and myricetin might also be considered hits that are amenable to synthetic modification to improve their anti-SARS-CoV-2 profile by inhibiting not only M^pro, but also the 3′–5′ exonuclease (ExoN).

Drug repurposing against main protease and RNA-dependent RNA polymerase of SARS-CoV-2 using molecular docking, MM-GBSA calculations and molecular dynamics

A virus called severe acute respiratory distress syndrome coronavirus type 2 (SARS‐CoV‐2) is the causing organism of coronavirus disease 2019 (COVID-19), which has severely affected human life and threatened public health. The pandemic took millions of lives worldwide and caused serious negative effects on human society and the economy. SARS-CoV-2 main protease (Mpro) and RNA-dependent RNA polymerase (RdRp) are interesting targets due to their crucial role in viral replication and growth. Since there is only one approved therapy for COVID-19, drug repurposing is a promising approach to finding molecules with potential activity against COVID-19 in a short time and at minimal cost. In this study, virtual screening was performed on the ChEMBL library containing 9923 FDA-approved drugs, using various docking filters with different accuracy. The best drugs with the highest docking scores were further examined for molecular dynamics (MD) studies and MM-GBSA calculations. The results of this study suggest that nadide, cangrelor and denufosol are promising potential candidates against COVID-19. Further in vitro, preclinical and clinical studies of these candidates would help to discover safe and effective anti-COVID-19 drugs.

Molecular Docking as a Potential Approach in Repurposing Drugs Against COVID-19: a Systematic Review and Novel Pharmacophore Models

Purpose of Review

This article provides a review of the recent literature related to the FDA-approved drugs that had been repurposed as potential drug candidates against COVID-19. Moreover, we performed a quality pharmacophore study for frequently studied targets, namely, the main protease, RNA-dependent RNA polymerase, and spike protein.

Recent Findings

Ever since the COVID-19 pandemic, the whole spectrum of scientific community is still unable to invent an absolute therapeutic agent for COVID-19. Considering such a fact, drug repurposing strategies seem a truly viable approach to develop novel therapeutic interventions.

Summery

Drug repurposing explores previously approved drugs of known safety and pharmacokinetics profile for possible new effects, reducing the cost, time, and predicting prospective side effects and drug interactions. COVID-19 virulent machinery appeared similar to other viruses, making antiviral agents widely repurposed in pursuit for curative candidates. Our main protease pharmacophoric study revealed multiple features and could be a probable starting point for upcoming research.

Prospective mode of action of Ivermectin: SARS-CoV-2

The well-known anti-helminthic drug ivermectin (IVM) has been established as an example of drug repurposing for the management of SARS-CoV-2 infection. Various study has been done to understand the inhibitory mechanism of IVM against SARS-CoV-2 targets. Broadly, IVM has been categorized as a host-directed agent and the proposed mechanism involves inhibition of the IMPα/ß1-mediated nuclear import of viral proteins. In addition, in vitro/in vivo and molecular docking/dynamic simulation studies suggested multitargets mechanism of IVM against SARS-CoV-2. Present manuscript attempts to provide an overview of the detailed mechanism of action based on experimental and computational studies. The knowledge of binding interaction of IVM and SARS-CoV-2 targets will give the direction to developed new and potential anti-COVID agents.

Exploring species-specific inhibitors with multiple target sites on S. aureus pyruvate kinase using a computational workflow

Experimental evidence indicated that bacterial pyruvate kinase of glycolysis can be evaluated as an alternative target to eliminate infections, while antibiotic resistance poses a global threat. Here, we use a computational workflow to reveal and investigate the potential allosteric sites of methicillin-resistant S. aureus PK, which can help in designing species-specific drugs to inhibit activity of this organism. Residue interaction networks point to a known allosteric site at the small C-C interface, a potential allosteric site near the small interface (site #1), and a second potential allosteric site at the large interface (site #2). 2 µs-long molecular dynamics (MD) simulations with AMBER16 generate different conformations of one narrow target site. Known and potential allosteric sites on the selected conformers are investigated using ensemble docking with AutoDock Vina and a library of 2447 FDA-approved drugs. We determine 18 hits, comprising ergot-alkaloids, anti-cancer-agents, antivirals, analgesics, cardiac glycosides, all with a high docking z-score for three sites. 5 selected compounds with high, average and low z-scores are subjected to 50 ns-long MD simulations for MM-GBSA calculations. ΔG_bind values up to -49.3 kcal/mol at the C-C interface, up to -32.7 kcal/mol at site #1, and up to -53.3 kcal/mol at site #2 support the docking calculations. We investigate mitapivat and TT-232 as reference compounds under clinical trial, targeting human PK isomers. We suggest 18 FDA-approved hits from the docking calculations and TT-232 as potential inhibitors with multiple target sites on S. aureus PK. This study also proposes pharmacophores models for de novo drug design.Communicated by Ramaswamy H. Sarma.

Computational Identification of Possible Allosteric Sites and Modulators of the SARS-CoV-2 Main Protease

In this study, we target the main protease (M^pro) of the SARS-CoV-2 virus as it is a crucial enzyme for viral replication. Herein, we report three plausible allosteric sites on M^pro that can expand structure-based drug discovery efforts for new M^pro inhibitors. To find these sites, we used mixed-solvent molecular dynamics (MixMD) simulations, an efficient computational protocol that finds binding hotspots through mapping the surface of unbound proteins with 5% cosolvents in water. We have used normal mode analysis to support our claim of allosteric control for these sites. Further, we have performed virtual screening against the sites with 361 hits from M^pro screenings available through the National Center for Advancing Translational Sciences (NCATS). We have identified the NCATS inhibitors that bind to the remote sites better than the active site of M^pro, and we propose these molecules may be allosteric regulators of the system. After identifying our sites, new X-ray crystal structures were released that show fragment molecules in the sites we found, supporting the notion that these sites are accurate and druggable.

Small molecule interactions with the SARS-CoV-2 main protease: In silico all-atom microsecond MD simulations, PELE Monte Carlo simulations, and determination of in vitro activity inhibition

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused the ongoing COVID-19 pandemic. With some notable exceptions, safe and effective vaccines, which are now being widely distributed globally, have largely begun to stabilise the situation. However, emerging variants of concern and vaccine hesitancy are apparent obstacles to eradication. Therefore, the need for the development of potent antivirals is still of importance. In this context, the SARS-CoV-2 main protease (M^pro) is a critical target and numerous clinical trials, predominantly in the private domain, are currently in progress. Here, our aim was to extend our previous studies, with hypericin and cyanidin-3-O-glucoside, as potential inhibitors of the SARS-CoV-2 M^pro. Firstly, we performed all-atom microsecond molecular dynamics simulations, which highlight the stability of the ligands in the M^pro active site over the duration of the trajectories. We also invoked PELE Monte Carlo simulations which indicate that both hypericin and cyanidin-3-O-glucoside preferentially interact with the M^pro active site and known allosteric sites. For further validation, we performed an in vitro enzymatic activity assay that demonstrated that hypericin and cyanidin-3-O-glucoside inhibit M^pro activity in a dose-dependent manner at biologically relevant (μM) concentrations. However, both ligands are much less potent than the well-known covalent antiviral GC376, which was used as a positive control in our experiments. Nevertheless, the biologically relevant activity of hypericin and cyanidin-3-O-glucoside is encouraging. In particular, a synthetic version of hypericin has FDA orphan drug designation, which could simplify potential clinical evaluation in the context of COVID-19.

Current status of structure-based drug repurposing against COVID-19 by targeting SARS-CoV-2 proteins

More than one and half years have passed, as of August 2021, since the COVID-19 caused by the novel coronavirus named SARS-CoV-2 emerged in 2019. While the recent success of vaccine developments likely reduces the severe cases, there is still a strong requirement of safety and effective therapeutic drugs for overcoming the unprecedented situation. Here we review the recent progress and the status of the drug discovery against COVID-19 with emphasizing a structure-based perspective. Structural data regarding the SARS-CoV-2 proteome has been rapidly accumulated in the Protein Data Bank, and up to 68% of the total amino acid residues encoded in the genome were covered by the structural data. Despite a global effort of in silico and in vitro screenings for drug repurposing, there is only a limited number of drugs had been successfully authorized by drug regulation organizations. Although many approved drugs and natural compounds, which exhibited antiviral activity in vitro, were considered potential drugs against COVID-19, a further multidisciplinary investigation is required for understanding the mechanisms underlying the antiviral effects of the drugs.

Effectiveness of Natural Antioxidants against SARS-CoV-2? Insights from the In-Silico World

The SARS CoV-2 pandemic has affected millions of people around the globe. Despite many efforts to find some effective medicines against SARS CoV-2, no established therapeutics are available yet. The use of phytochemicals as antiviral agents provides hope against the proliferation of SARS-CoV-2. Several natural compounds were analyzed by virtual screening against six SARS CoV-2 protein targets using molecular docking simulations in the present study. More than a hundred plant-derived secondary metabolites have been docked, including alkaloids, flavonoids, coumarins, and steroids. SARS CoV-2 protein targets include Main protease (MPro), Papain-like protease (PLpro), RNA-dependent RNA polymerase (RdRp), Spike glycoprotein (S), Helicase (Nsp13), and E-Channel protein. Phytochemicals were evaluated by molecular docking, and MD simulations were performed using the YASARA structure using a modified genetic algorithm and AMBER03 force field. Binding energies and dissociation constants allowed the identification of potentially active compounds. Ligand-protein interactions provide an insight into the mechanism and potential of identified compounds. Glycyrrhizin and its metabolite 18-β-glycyrrhetinic acid have shown a strong binding affinity for MPro, helicase, RdRp, spike, and E-channel proteins, while a flavonoid Baicalin also strongly binds against PLpro and RdRp. The use of identified phytochemicals may help to speed up the drug development and provide natural protection against SARS-CoV-2.