Computational Analysis Reveals Monomethylated Triazolopyrimidine as a Novel Inhibitor of SARS-CoV-2 RNA-Dependent RNA Polymerase (RdRp)

The human population is still facing appalling conditions due to several outbreaks of Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) virus. The absence of specific drugs, appropriate vaccines for mutants, and knowledge of potential therapeutic agents makes this situation more difficult. Several 1, 2, 4-triazolo [1, 5-a] pyrimidine (TP)-derivative compounds were comprehensively studied for antiviral activities against RNA polymerase of HIV, HCV, and influenza viruses, and showed immense pharmacological interest. Therefore, TP-derivative compounds can be repurposed against the RNA-dependent RNA polymerase (RdRp) protein of SARS-CoV-2. In this study, a meta-analysis was performed to ensure the genomic variability and stability of the SARS-CoV-2 RdRp protein. The molecular docking of natural and synthetic TP compounds to RdRp and molecular dynamic (MD) simulations were performed to analyse the dynamic behaviour of TP compounds at the active site of the RdRp protein. TP compounds were also docked against other non-structural proteins (NSP1, NSP2, NSP3, NSP5, NSP8, NSP13, and NSP15) of SARS-CoV-2. Furthermore, the inhibition potential of TP compounds was compared with Remdesivir and Favipiravir drugs as a positive control. Additionally, TP compounds were analysed for inhibitory activity against SARS-CoV RdRp protein. This study demonstrates that TP analogues (monomethylated triazolopyrimidine and essramycin) represent potential lead molecules for designing an effective inhibitor to control viral replication. Furthermore, in vitro and in vivo studies will strengthen the use of these inhibitors as suitable drug candidates against SARS-CoV-2.     

Targeting the RdRp of Emerging RNA Viruses: The Structure-Based Drug Design Challenge

The RNA-dependent RNA polymerase (RdRp) is an essential enzyme for the viral replication process, catalyzing the viral RNA synthesis using a metal ion-dependent mechanism. In recent years, RdRp has emerged as an optimal target for the development of antiviral drugs, as demonstrated by recent approvals of sofosbuvir and remdesivir against Hepatitis C virus (HCV) and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), respectively. In this work, we overview the main sequence and structural features of the RdRp of emerging RNA viruses such as Coronaviruses, Flaviviruses, and HCV, as well as inhibition strategies implemented so far. While analyzing the structural information available on the RdRp of emerging RNA viruses, we provide examples of success stories such as for HCV and SARS-CoV-2. In contrast, Flaviviruses’ story has raised attention about how the lack of structural details on catalytically-competent or ligand-bound RdRp strongly hampers the application of structure-based drug design, either in repurposing and conventional approaches.     

CRISPR-Mediated Profiling of Viral RNA at Single-Nucleotide Resolution

Mass pathogen screening is critical to preventing the outbreaks and spread of infectious diseases. The large-scale epidemic of COVID-19 and the rapid mutation of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus have put forward new requirements for virus detection and identification techniques. Here, we report a CRISPR-based Amplification-free Viral RNA Electrical Detection platform (CAVRED) for the rapid detection and identification of SARS-CoV-2 variants. A series of CRISPR RNA assays were designed to amplify the CRISPR-Cas system‘s ability to discriminate between mutant and wild RNA genomes with a single-nucleotide difference. The identified viral RNA information was converted into readable electrical signals through field-effect transistor biosensors for the achievement of highly sensitive detection of single-base mutations. CAVRED can detect the SARS-CoV-2 virus genome as low as 1 cp μL⁻¹ within 20 mins without amplification, and this value is comparable to the detection limit of real-time quantitative polymerase chain reaction. Based on the excellent RNA mutation detection ability, an 8-in-1 CAVRED array was constructed and realized the rapid identification of 40 simulated throat swab samples of SARS-CoV-2 variants with a 95.0 % accuracy. The advantages of accuracy, sensitivity, and fast speed of CAVRED promise its application in rapid and large-scale epidemic screening.     

Structure-based virtual screening and molecular dynamics of phytochemicals derived from Saudi medicinal plants to identify potential COVID-19 therapeutics

Coronavirus disease 2019 (COVID-19) has affected almost every country in the world by causing a global pandemic with a high mortality rate. Lack of an effective vaccine and/or antiviral drugs against SARS-CoV-2, the causative agent, has severely hampered the response to this novel coronavirus. Natural products have long been used in traditional medicines to treat various diseases, and purified phytochemicals from medicinal plants provide a valuable scaffold for the discovery of new drug leads. In the present study, we performed a computational screening of an in-house database composed of ~1000 phytochemicals derived from traditional Saudi medicinal plants with recognised antiviral activity. Structure-based virtual screening was carried out against three druggable SARS-CoV-2 targets, viral RNA-dependent RNA polymerase (RdRp), 3-chymotrypsin-like cysteine protease (3CL^pro) and papain like protease (PL^pro) to identify putative inhibitors that could facilitate the development of potential anti-COVID-19 drug candidates. Computational analyses identified three compounds inhibiting each target, with binding affinity scores ranging from −9.9 to −6.5 kcal/mol. Among these, luteolin 7-rutinoside, chrysophanol 8-(6-galloylglucoside) and kaempferol 7-(6″-galloylglucoside) bound efficiently to RdRp, while chrysophanol 8-(6-galloylglucoside), 3,4,5-tri-O-galloylquinic acid and mulberrofuran G interacted strongly with 3CL^pro, and withanolide A, isocodonocarpine and calonysterone bound tightly to PL^pro. These potential drug candidates will be subjected to further in vitro and in vivo studies and may assist the development of effective anti-COVID-19 drugs.     

Nucleotide Analogues Bearing a C2′ or C3′-Stereogenic All-Carbon Quaternary Center as SARS-CoV-2 RdRp Inhibitors

The design of novel nucleoside triphosphate (NTP) analogues bearing an all-carbon quaternary center at C2′ or C3′ is described. The construction of this all-carbon stereogenic center involves the use of an intramoleculer photoredox-catalyzed reaction. The nucleoside analogues (NA) hydroxyl functional group at C2′ was generated by diastereoselective epoxidation. In addition, highly enantioselective and diastereoselective Mukaiyama aldol reactions, diastereoselective N-glycosylations and regioselective triphosphorylation reactions were employed to synthesize the novel NTPs. Two of these compounds are inhibitors of the RNA-dependent RNA polymerase (RdRp) of SARS-CoV-2, the causal virus of COVID-19.     

Molecular docking revealed the binding of nucleotide/side inhibitors to Zika viral polymerase solved structures

A new Zika virus (ZIKV) outbreak started in 2015. According to the World Health Organization, 84 countries confirmed ZIKV infection. RNA-dependent RNA polymerase (RdRp) was an appealing target for drug designers during the last two decades. Through molecular docking, we screened 16 nucleotide/side inhibitors against ZIKV RdRp. While the mode of interaction with ZIKV is different from that in the hepatitis C virus (HCV), nucleotide/side inhibitors in this study (mostly anti-HCV) showed promising binding affinities (?6.2 to ?9.7 kcal/mol calculated by AutoDock Vina) to ZIKV RdRp. Setrobuvir, YAK and, to a lesser extent, IDX-184 reveal promising results compared to other inhibitors in terms of binding ZIKV RdRp. These candidates would be powerful anti-ZIKV drugs.     

In silico study of novel niclosamide derivatives,SARS-CoV-2 nonstructural proteins catalytic residue-targeting small molecules drug candidates

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2)-mediated coronavirus disease 2019 (COVID-19) infection remains a global pandemic and health emergency with overwhelming social and economic impacts throughout the world. Therapeutics for COVID-19 are limited to only remdesivir; therefore, there is a need for combined, multidisciplinary efforts to develop new therapeutic molecules and explore the effectiveness of existing drugs against SARS-CoV-2. In the present study, we reported eight (SCOV-L-02, SCOV-L-09, SCOV-L-10, SCOV-L-11, SCOV-L-15, SCOV-L-18, SCOV-L-22, and SCOV-L-23) novel structurally related small-molecule derivatives of niclosamide (SCOV-L series) for their targeting potential against angiotensin-converting enzyme-2 (ACE2), type II transmembrane serine protease (TMPRSS2), and SARS-COV-2 nonstructural proteins (NSPs) including NSP5 (3CLpro), NSP3 (PLpro), and RdRp. Our correlation analysis suggested that ACE2 and TMPRSS2 modulate host immune response via regulation of immune-infiltrating cells at the site of tissue/organs entries. In addition, we identified some TMPRSS2 and ACE2 microRNAs target regulatory networks in SARS-CoV-2 infection and thus open up a new window for microRNAs-based therapy for the treatment of SARS-CoV-2 infection. Our in vitro study revealed that with the exception of SCOV-L-11 and SCOV-L-23 which were non-active, the SCOV-L series exhibited strict antiproliferative activities and non-cytotoxic effects against ACE2- and TMPRSS2-expressing cells. Our molecular docking for the analysis of receptor-ligand interactions revealed that SCOV-L series demonstrated high ligand binding efficacies (at higher levels than clinical drugs) against the ACE2, TMPRSS2, and SARS-COV-2 NSPs. SCOV-L-18, SCOV-L-15, and SCOV-L-09 were particularly found to exhibit strong binding affinities with three key SARS-CoV-2’s proteins: 3CLpro, PLpro, and RdRp. These compounds bind to the several catalytic residues of the proteins, and satisfied the criteria of drug-like candidates, having good adsorption, distribution, metabolism, excretion, and toxicity (ADMET) pharmacokinetic profile. Altogether, the present study suggests the therapeutic potential of SCOV-L series for preventing and managing SARs-COV-2 infection and are currently under detailed investigation in our lab.     

Molecular dynamics simulation study of PTP1B with allosteric inhibitor and its application in receptor based pharmacophore modeling

Allosteric inhibition of protein tyrosine phosphatase 1B (PTP1B), has paved a new path to design specific inhibitors for PTP1B, which is an important drug target for the treatment of type II diabetes and obesity. The PTP1B_1–282-allosteric inhibitor complex crystal structure lacks α7 (287–298) and moreover there is no available 3D structure of PTP1B_1–298 in open form. As the interaction between α7 and α6–α3 helices plays a crucial role in allosteric inhibition, α7 was modeled to the PTP1B_1–282 in open form complexed with an allosteric inhibitor (compound-2) and a 5 ns MD simulation was performed to investigate the relative orientation of the α7–α6–α3 helices. The simulation conformational space was statistically sampled by clustering analyses. This approach was helpful to reveal certain clues on PTP1B allosteric inhibition. The simulation was also utilized in the generation of receptor based pharmacophore models to include the conformational flexibility of the protein-inhibitor complex. Three cluster representative structures of the highly populated clusters were selected for pharmacophore model generation. The three pharmacophore models were subsequently utilized for screening databases to retrieve molecules containing the features that complement the allosteric site. The retrieved hits were filtered based on certain drug-like properties and molecular docking simulations were performed in two different conformations of protein. Thus, performing MD simulation with α7 to investigate the changes at the allosteric site, then developing receptor based pharmacophore models and finally docking the retrieved hits into two distinct conformations will be a reliable methodology in identifying PTP1B allosteric inhibitors. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Potential Achilles heels of SARS-CoV-2 are best displayed by the base order-dependent component of RNA folding energy

The base order-dependent component of folding energy has revealed a highly conserved region in HIV-1 genomes that associates with RNA structure. This corresponds to a packaging signal that is recognized by the nucleocapsid domain of the Gag polyprotein. Long viewed as a potential HIV-1 "Achilles heel," the signal can be targeted by a new antiviral compound. Although SARS-CoV-2 differs in many respects from HIV-1, the same technology displays regions with a high base order-dependent folding energy component, which are also highly conserved. This indicates structural invariance (SI) sustained by natural selection. While the regions are often also protein-encoding (e. g. NSP3, ORF3a), we suggest that their nucleic acid level functions can be considered potential "Achilles heels" for SARS-CoV-2, perhaps susceptible to therapies like those envisaged for AIDS. The ribosomal frameshifting element scored well, but higher SI scores were obtained in other regions, including those encoding NSP13 and the nucleocapsid (N) protein.     

In silico identification of noncompetitive inhibitors targeting an uncharacterized allosteric site of falcipain-2

Falcipain-2 (FP-2) is a Plasmodium falciparum hemoglobinase widely targeted in the search for antimalarials. FP-2 can be allosterically modulated by various noncompetitive inhibitors that have been serendipitously identified. Moreover, the crystal structures of two inhibitors bound to an allosteric site, termed site 6, of the homolog enzyme human cathepsin K (hCatK) suggest that the equivalent region in FP-2 might play a similar role. Here, we conduct the rational identification of FP-2 inhibitors through virtual screenings (VS) of compounds into several pocket-like conformations of site 6, sampled during molecular dynamics (MD) simulations of the free enzyme. Two noncompetitive inhibitors, ZINC03225317 and ZINC72290660, were confirmed using in vitro enzymatic assays and their poses into site 6 led to calculated binding free energies matching the experimental ones. Our results provide strong evidence about the allosteric inhibition of FP-2 through binding of small molecules to site 6, thus opening the way toward the discovery of new inhibitors against this enzyme.

     

Pharmacophore model-aided virtual screening combined with comparative molecular docking and molecular dynamics for identification of marine natural products as SARS-CoV-2 papain-like protease inhibitors

Targeting SARS-CoV-2 papain-like protease using inhibitors is a suitable approach for inhibition of virus replication and dysregulation of host anti-viral immunity. Engaging all five binding sites far from the catalytic site of PLpro is essential for developing a potent inhibitor. We developed and validated a structure-based pharmacophore model with 9 features of a potent PLpro inhibitor. The pharmacophore model-aided virtual screening of the comprehensive marine natural product database predicted 66 initial hits. This hit library was downsized by filtration through a molecular weight filter of ≤ 500 g/mol. The 50 resultant hits were screened by comparative molecular docking using AutoDock and AutoDock Vina. Comparative molecular docking enables benchmarking docking and relieves the disparities in the search and scoring functions of docking engines. Both docking engines retrieved 3 same compounds at different positions in the top 1 % rank, hence consensus scoring was applied, through which CMNPD28766, aspergillipeptide F emerged as the best PLpro inhibitor. Aspergillipeptide F topped the 50-hit library with a pharmacophore-fit score of 75.916. Favorable binding interactions were predicted between aspergillipeptide F and PLpro similar to the native ligand XR8-24. Aspergillipeptide F was able to engage all the 5 binding sites including the newly discovered BL2 groove, site V. Molecular dynamics for quantification of Cα-atom movements of PLpro after ligand binding indicated that it exhibits highly correlated domain movements contributing to the low free energy of binding and a stable conformation. Thus, aspergillipeptide F is a promising candidate for pharmaceutical and clinical development as a potent SARS-CoV-2 PLpro inhibitor.     

Coronavirus Infection-Associated Cell Death Signaling and Potential Therapeutic Targets

COVID-19 is the name of the disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection that occurred in 2019. The virus–host-specific interactions, molecular targets on host cell deaths, and the involved signaling are crucial issues, which become potential targets for treatment. Spike protein, angiotensin-converting enzyme 2 (ACE2), cathepsin L-cysteine peptidase, transmembrane protease serine 2 (TMPRSS2), nonstructural protein 1 (Nsp1), open reading frame 7a (ORF7a), viral main protease (3C-like protease (3CLpro) or Mpro), RNA dependent RNA polymerase (RdRp) (Nsp12), non-structural protein 13 (Nsp13) helicase, and papain-like proteinase (PLpro) are molecules associated with SARS-CoV infection and propagation. SARS-CoV-2 can induce host cell death via five kinds of regulated cell death, i.e., apoptosis, necroptosis, pyroptosis, autophagy, and PANoptosis. The mechanisms of these cell deaths are well established and can be disrupted by synthetic small molecules or natural products. There are a variety of compounds proven to play roles in the cell death inhibition, such as pan-caspase inhibitor (z-VAD-fmk) for apoptosis, necrostatin-1 for necroptosis, MCC950, a potent and specific inhibitor of the NLRP3 inflammasome in pyroptosis, and chloroquine/hydroxychloroquine, which can mitigate the corresponding cell death pathways. However, NF-κB signaling is another critical anti-apoptotic or survival route mediated by SARS-CoV-2. Such signaling promotes viral survival, proliferation, and inflammation by inducing the expression of apoptosis inhibitors such as Bcl-2 and XIAP, as well as cytokines, e.g., TNF. As a result, tiny natural compounds functioning as proteasome inhibitors such as celastrol and curcumin can be used to modify NF-κB signaling, providing a responsible method for treating SARS-CoV-2-infected patients. The natural constituents that aid in inhibiting viral infection, progression, and amplification of coronaviruses are also emphasized, which are in the groups of alkaloids, flavonoids, terpenoids, diarylheptanoids, and anthraquinones. Natural constituents derived from medicinal herbs have anti-inflammatory and antiviral properties, as well as inhibitory effects, on the viral life cycle, including viral entry, replication, assembly, and release of COVID-19 virions. The phytochemicals contain a high potential for COVID-19 treatment. As a result, SARS-CoV-2-infected cell death processes and signaling might be of high efficacy for therapeutic targeting effects and yielding encouraging outcomes.     

Chemical Proteomic Discovery of Isotype-Selective Covalent Inhibitors of the RNA Methyltransferase NSUN2

5-Methylcytosine (m⁵C) is an RNA modification prevalent on tRNAs, where it can protect tRNAs from endonucleolytic cleavage to maintain protein synthesis. The NSUN family (NSUN1-7 in humans) of RNA methyltransferases are capable of installing the methyl group onto the C⁵ position of cytosines in RNA. NSUNs are implicated in a wide range of (patho)physiological processes, but selective and cell-active inhibitors of these enzymes are lacking. Here, we use cysteine-directed activity-based protein profiling (ABPP) to discover azetidine acrylamides that act as stereoselective covalent inhibitors of human NSUN2. Despite targeting a conserved catalytic cysteine in the NSUN family, the NSUN2 inhibitors show negligible cross-reactivity with other human NSUNs and exhibit good proteome-wide selectivity. We verify that the azetidine acrylamides inhibit the catalytic activity of recombinant NSUN2, but not NSUN6, and demonstrate that these compounds stereoselectively disrupt NSUN2-tRNA interactions in cancer cells, leading to a global reduction in tRNA m⁵C content. Our findings thus highlight the potential to create isotype-selective and cell-active inhibitors of NSUN2 with covalent chemistry targeting a conserved catalytic cysteine.     

Site mapping and small molecule blind docking reveal a possible target site on the SARS-CoV-2 main protease dimer interface

The SARS-CoV-2 virus is causing COVID-19 resulting in an ongoing pandemic with serious health, social, and economic implications. Much research is focused in repurposing or identifying new small molecules which may interact with viral or host-cell molecular targets. An important SARS-CoV-2 target is the main protease (M^pro), and the peptidomimetic α-ketoamides represent prototypical experimental inhibitors. The protease is characterised by the dimerization of two monomers each which contains the catalytic dyad defined by Cys¹⁴⁵ and His⁴¹ residues (active site). Dimerization yields the functional homodimer. Here, our aim was to investigate small molecules, including lopinavir and ritonavir, α-ketoamide 13b, and ebselen, for their ability to interact with the M^pro. The sirtuin 1 agonist SRT1720 was also used in our analyses. Blind docking to each monomer individually indicated preferential binding of the ligands in the active site. Site-mapping of the dimeric protease indicated a highly reactive pocket in the dimerization region at the domain III apex. Blind docking consistently indicated a strong preference of ligand binding in domain III, away from the active site. Molecular dynamics simulations indicated that ligands docked both to the active site and in the dimerization region at the apex, formed relatively stable interactions. Overall, our findings do not obviate the superior potency with respect to inhibition of protease activity of covalently-linked inhibitors such as α-ketoamide 13b in the M^pro active site. Nevertheless, along with those from others, our findings highlight the importance of further characterisation of the M^pro active site and any potential allosteric sites.     

Field-Template,QSAR, Ensemble Molecular Docking,and 3D-RISM Solvation Studies Expose Potential of FDA-Approved Marine Drugs as SARS-CoVID-2 Main Protease Inhibitors

Currently, SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) has infected people among all countries and is a pandemic as declared by the World Health Organization (WHO). SARS-CoVID-2 main protease is one of the therapeutic drug targets that has been shown to reduce virus replication, and its high-resolution 3D structures in complex with inhibitors have been solved. Previously, we had demonstrated the potential of natural compounds such as serine protease inhibitors eventually leading us to hypothesize that FDA-approved marine drugs have the potential to inhibit the biological activity of SARS-CoV-2 main protease. Initially, field-template and structure–activity atlas models were constructed to understand and explain the molecular features responsible for SARS-CoVID-2 main protease inhibitors, which revealed that Eribulin Mesylate, Plitidepsin, and Trabectedin possess similar characteristics related to SARS-CoVID-2 main protease inhibitors. Later, protein–ligand interactions are studied using ensemble molecular-docking simulations that revealed that marine drugs bind at the active site of the main protease. The three-dimensional reference interaction site model (3D-RISM) studies show that marine drugs displace water molecules at the active site, and interactions observed are favorable. These computational studies eventually paved an interest in further in vitro studies. Finally, these findings are new and indeed provide insights into the role of FDA-approved marine drugs, which are already in clinical use for cancer treatment as a potential alternative to prevent and treat infected people with SARS-CoV-2.     

Histidine affinity chromatography-based methodology for the simultaneous isolation of Escherichia coli small and ribosomal RNA

Research on RNA has led to many important biological discoveries and the improvement of therapeutic technologies. In particular, there is a great focus on small RNA and ribosomal RNA owing to their key functions in the cell, which make them excellent therapeutic targets. Although the study of these RNA classes is progressing, some limitations have been found regarding the use of suitable techniques that are able to produce and isolate biologically competent and chemically stable RNA. To address this, we have developed a novel histidine affinity chromatography-based isolation methodology for small and ribosomal RNA molecules. The new procedure involves three main steps: (1) cell lysis with guanidinium buffer, (2) RNA primary isolation with ammonium sulfate precipitation and (3) histidine affinity chromatography to specifically purify small RNA and ribosomal RNA from other Escherichia coli impurities (genomic DNA and proteins). The RNA quality assessment revealed that both RNA species were obtained with a high recovery, integrity and purity. The potential of this method to achieve a reproducible RNA isolation with appropriate quality has been demonstrated and it should have broad application in the structural, biophysical and biomedical investigation of systems involving RNA components.     

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.