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.

Mechanistic insights into the Japanese encephalitis virus RNA dependent RNA polymerase protein inhibition by bioflavonoids from Azadirachta indica.

Japanese encephalitis (JE) virus is a flavivirus causing encephalitis causing neurological damage. RNA-dependent-RNA-polymerase (RdRp) is responsible for genome replication making it excellent anti-viral target. In this study, the crystal structure of JE RdRp (jRdRp) and bioflavonoids reported in Azadirachta indica were retrieved from specific databases. Structure-based virtual screening was employed using MTiOpenScreen server and top four compounds selected with the most negative docking scores. Conformations were redocked using AutoDock Vina; these complexes showed mechanistic interactions with Arg474, Gly605, Asp668, and Trp800 residues in the active site of jRdRp, i.e., guanosine-5'-triphosphate. Furthermore, 100 ns classical molecular dynamics simulation and binding free energy calculation showed stability of docked bioflavonoids in the active jRdRp pocket and significant contribution of van-der-Waals interactions for docked complex stability during simulation. Therefore, this study predicted the anti-viral activity of Gedunin, Nimbolide, Ohchinin acetate, and Kulactone against jRdRp and can be considered for further antiviral drug development.

Natural plant products as potential inhibitors of RNA dependent RNA polymerase of Severe Acute Respiratory Syndrome Coronavirus-2.

Drug repurposing studies targeting inhibition of RNA dependent RNA polymerase (RdRP) of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) have exhibited the potential effect of small molecules. In the present work a detailed interaction study between the phytochemicals from Indian medicinal plants and the RdRP of SARS-CoV-2 has been performed. The top four phytochemicals obtained through molecular docking were, swertiapuniside, cordifolide A, sitoindoside IX, and amarogentin belonging to Swertia chirayita, Tinospora cordifolia and Withania somnifera. These ligands bound to the RdRP were further studied using molecular dynamics simulations. The principal component analysis of these systems showed significant conformational changes in the finger and thumb subdomain of the RdRP. Hydrogen bonding, salt-bridge and water mediated interactions supported by MM-GBSA free energy of binding revealed strong binding of cordifolide A and sitoindoside IX to RdRP. The ligand-interacting residues belonged to either of the seven conserved motifs of the RdRP. These residues were polar and charged amino acids, namely, ARG 553, ARG 555, ASP 618, ASP 760, ASP 761, GLU 811, and SER 814. The glycosidic moieties of the phytochemicals were observed to form favourable interactions with these residues. Hence, these phytochemicals may hold the potential to act as RdRP inhibitors owing to their stability in binding to the druggable site.

In Silico Screening of the DrugBank Database to Search for Possible Drugs against SARS-CoV-2.

Coronavirus desease 2019 (COVID-19) is responsible for more than 1.80 M deaths worldwide. A Quantitative Structure-Activity Relationships (QSAR) model is developed based on experimental pIC50 values reported for a structurally diverse dataset. A robust model with only five descriptors is found, with values of R2 = 0.897, Q2LOO = 0.854, and Q2ext = 0.876 and complying with all the parameters established in the validation Tropsha's test. The analysis of the applicability domain (AD) reveals coverage of about 90% for the external test set. Docking and molecular dynamic analysis are performed on the three most relevant biological targets for SARS-CoV-2: main protease, papain-like protease, and RNA-dependent RNA polymerase. A screening of the DrugBank database is executed, predicting the pIC50 value of 6664 drugs, which are IN the AD of the model (coverage = 79%). Fifty-seven possible potent anti-COVID-19 candidates with pIC50 values > 6.6 are identified, and based on a pharmacophore modelling analysis, four compounds of this set can be suggested as potent candidates to be potential inhibitors of SARS-CoV-2. Finally, the biological activity of the compounds was related to the frontier molecular orbitals shapes.

A novel virtual screening procedure identifies Pralatrexate as inhibitor of SARS-CoV-2 RdRp and it reduces viral replication in vitro.

The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus poses serious threats to the global public health and leads to worldwide crisis. No effective drug or vaccine is readily available. The viral RNA-dependent RNA polymerase (RdRp) is a promising therapeutic target. A hybrid drug screening procedure was proposed and applied to identify potential drug candidates targeting RdRp from 1906 approved drugs. Among the four selected market available drug candidates, Pralatrexate and Azithromycin were confirmed to effectively inhibit SARS-CoV-2 replication in vitro with EC50 values of 0.008µM and 9.453 µM, respectively. For the first time, our study discovered that Pralatrexate is able to potently inhibit SARS-CoV-2 replication with a stronger inhibitory activity than Remdesivir within the same experimental conditions. The paper demonstrates the feasibility of fast and accurate anti-viral drug screening for inhibitors of SARS-CoV-2 and provides potential therapeutic agents against COVID-19.

In silico studies evidenced the role of structurally diverse plant secondary metabolites in reducing SARS-CoV-2 pathogenesis.

Plants are endowed with a large pool of structurally diverse small molecules known as secondary metabolites. The present study aims to virtually screen these plant secondary metabolites (PSM) for their possible anti-SARS-CoV-2 properties targeting four proteins/ enzymes which govern viral pathogenesis. Results of molecular docking with 4,704 ligands against four target proteins, and data analysis revealed a unique pattern of structurally similar PSM interacting with the target proteins. Among the top-ranked PSM which recorded lower binding energy (BE), > 50% were triterpenoids which interacted strongly with viral spike protein-receptor binding domain, > 32% molecules which showed better interaction with the active site of human transmembrane serine protease were belongs to flavonoids and their glycosides, > 16% of flavonol glycosides and > 16% anthocyanidins recorded lower BE against active site of viral main protease and > 13% flavonol glycoside strongly interacted with active site of viral RNA-dependent RNA polymerase. The primary concern about these PSM is their bioavailability. However, several PSM recorded higher bioavailability score and found fulfilling most of the drug-likeness characters as per Lipinski's rule (Coagulin K, Kamalachalcone C, Ginkgetin, Isoginkgetin, 3,3'-Biplumbagin, Chrysophanein, Aromoline, etc.). Natural occurrence, bio-transformation, bioavailability of selected PSM and their interaction with the target site of selected proteins were discussed in detail. Present study provides a platform for researchers to explore the possible use of selected PSM to prevent/ cure the COVID-19 by subjecting them for thorough in vitro and in vivo evaluation for the capabilities to interfering with the process of viral host cell recognition, entry and replication.

COSMO-RS-Based Descriptors for the Machine Learning-Enabled Screening of Nucleotide Analogue Drugs against SARS-CoV-2.

Chemical similarity-based approaches employed to repurpose or develop new treatments for emerging diseases, such as COVID-19, correlates molecular structure-based descriptors of drugs with those of a physiological counterpart or clinical phenotype. We propose novel descriptors based on a COSMO-RS (short for conductor-like screening model for real solvents) σ-profiles for enhanced drug screening enabled by machine learning (ML). The descriptors' performance is hereby illustrated for nucleotide analogue drugs that inhibit the ribonucleic acid-dependent ribonucleic acid polymerase, key to viral transcription and genome replication. The COSMO-RS-based descriptors account for both chemical reactivity and structure, and are more effective for ML-based screening than fingerprints based on molecular structure and simple physical/chemical properties. The descriptors are evaluated using principal component analysis, an unsupervised ML technique. Our results correlate with the active monophosphate forms of the leading drug remdesivir and the prospective drug EIDD-2801 with nucleotides, followed by other promising drugs, and are superior to those from molecular structure-based descriptors and molecular docking. The COSMO-RS-based descriptors could help accelerate drug discovery for the treatment of emerging diseases.

Unraveling the anti-influenza effect of flavonoids: Experimental validation of luteolin and its congeners as potent influenza endonuclease inhibitors.

The biological effects of flavonoids on mammal cells are diverse, ranging from scavenging free radicals and anti-cancer activity to anti-influenza activity. Despite appreciable effort to understand the anti-influenza activity of flavonoids, there is no clear consensus about their precise mode-of-action at a cellular level. Here, we report the development and validation of a screening assay based on AlphaScreen technology and illustrate its application for determination of the inhibitory potency of a large set of polyols against PA N-terminal domain (PA-Nter) of influenza RNA-dependent RNA polymerase featuring endonuclease activity. The most potent inhibitors we identified were luteolin with an IC50 of 72 ± 2 nM and its 8-C-glucoside orientin with an IC50 of 43 ± 2 nM. Submicromolar inhibitors were also evaluated by an in vitro endonuclease activity assay using single-stranded DNA, and the results were in full agreement with data from the competitive AlphaScreen assay. Using X-ray crystallography, we analyzed structures of the PA-Nter in complex with luteolin at 2.0 Å resolution and quambalarine B at 2.5 Å resolution, which clearly revealed the binding pose of these polyols coordinated to two manganese ions in the endonuclease active site. Using two distinct assays along with the structural work, we have presumably identified and characterized the molecular mode-of-action of flavonoids in influenza-infected cells.

In silico ADMET and molecular docking study on searching potential inhibitors from limonoids and triterpenoids for COVID-19.

Virtual screening of phytochemicals was performed through molecular docking, simulations, in silico ADMET and drug-likeness prediction to identify the potential hits that can inhibit the effects of SARS-CoV-2. Considering the published literature on medicinal importance, 154 phytochemicals with analogous structure from limonoids and triterpenoids were selected to search potential inhibitors for the five therapeutic protein targets of SARS-CoV-2, i.e., 3CLpro (main protease), PLpro (papain-like protease), SGp-RBD (spike glycoprotein-receptor binding domain), RdRp (RNA dependent RNA polymerase) and ACE2 (angiotensin-converting enzyme 2). The in silico computational results revealed that the phytochemicals such as glycyrrhizic acid, limonin, 7-deacetyl-7-benzoylgedunin, maslinic acid, corosolic acid, obacunone and ursolic acid were found to be effective against the target proteins of SARS-CoV-2. The protein-ligand interaction study revealed that these phytochemicals bind with the amino acid residues at the active site of the target proteins. Therefore, the core structure of these potential hits can be used for further lead optimization to design drugs for SARS-CoV-2. Also, the medicinal plants containing these phytochemicals like licorice, neem, tulsi, citrus and olives can be used to formulate suitable therapeutic approaches in traditional medicines.

Analysis of SARS-CoV-2 RNA-dependent RNA polymerase as a potential therapeutic drug target using a computational approach.

BACKGROUND: The Severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) outbreak originating in Wuhan, China, has raised global health concerns and the pandemic has now been reported on all inhabited continents. Hitherto, no antiviral drug is available to combat this viral outbreak. METHODS: Keeping in mind the urgency of the situation, the current study was designed to devise new strategies for drug discovery and/or repositioning against SARS-CoV-2. In the current study, RNA-dependent RNA polymerase (RdRp), which regulates viral replication, is proposed as a potential therapeutic target to inhibit viral infection. RESULTS: Evolutionary studies of whole-genome sequences of SARS-CoV-2 represent high similarity (> 90%) with other SARS viruses. Targeting the RdRp active sites, ASP760 and ASP761, by antiviral drugs could be a potential therapeutic option for inhibition of coronavirus RdRp, and thus viral replication. Target-based virtual screening and molecular docking results show that the antiviral Galidesivir and its structurally similar compounds have shown promise against SARS-CoV-2. CONCLUSIONS: The anti-polymerase drugs predicted here-CID123624208 and CID11687749-may be considered for in vitro and in vivo clinical trials.

Phylogenetic Analysis and Structural Perspectives of RNA-Dependent RNA-Polymerase Inhibition from SARs-CoV-2 with Natural Products.

Most recently, an outbreak of severe pneumonia caused by the infection of SARS-CoV-2, a novel coronavirus first identified in Wuhan, China, imposes serious threats to public health. Upon infecting host cells, coronaviruses assemble a multi-subunit RNA-synthesis complex of viral non-structural proteins (nsp) responsible for the replication and transcription of the viral genome. Therefore, the role and inhibition of nsp12 are indispensable. A cryo-EM structure of RdRp from SARs-CoV-2 was used to identify novel drugs from Northern South African medicinal compounds database (NANPDB) by using computational virtual screening and molecular docking approaches. Considering Remdesivir as the control, 42 compounds were shortlisted to have docking score better than Remdesivir. The top 5 hits were validated by using molecular dynamics simulation approach and free energy calculations possess strong inhibitory properties than the Remdesivir. Thus, this study paved a way for designing novel drugs by decoding the architecture of an important enzyme and its inhibition with compounds from natural resources. This disclosing of necessary knowledge regarding the screening and the identification of top hits could help to design effective therapeutic candidates against the coronaviruses and design robust preventive measurements.

The potential chemical structure of anti-SARS-CoV-2 RNA-dependent RNA polymerase.

An outbreak of coronavirus disease 2019 (COVID-19) occurred in Wuhan and it has rapidly spread to almost all parts of the world. For coronaviruses, RNA-dependent RNA polymerase (RdRp) is an important polymerase that catalyzes the replication of RNA from RNA template and is an attractive therapeutic target. In this study, we screened these chemical structures from traditional Chinese medicinal compounds proven to show antiviral activity in severe acute respiratory syndrome coronavirus (SARS-CoV) and the similar chemical structures through a molecular docking study to target RdRp of SARS-CoV-2, SARS-CoV, and Middle East respiratory syndrome coronavirus (MERS-CoV). We found that theaflavin has a lower idock score in the catalytic pocket of RdRp in SARS-CoV-2 (-9.11 kcal/mol), SARS-CoV (-8.03 kcal/mol), and MERS-CoV (-8.26 kcal/mol) from idock. To confirm the result, we discovered that theaflavin has lower binding energy of -8.8 kcal/mol when it docks in the catalytic pocket of SARS-CoV-2 RdRp by using the Blind Docking server. Regarding contact modes, hydrophobic interactions contribute significantly in binding and additional hydrogen bonds were found between theaflavin and RdRp. Moreover, one π-cation interaction was formed between theaflavin and Arg553 from the Blind Docking server. Our results suggest that theaflavin could be a potential SARS-CoV-2 RdRp inhibitor for further study.

Discovery of a small molecule inhibitor targeting dengue virus NS5 RNA-dependent RNA polymerase.

Dengue is a mosquito-borne viral infection that has spread globally in recent years. Around half of the world's population, especially in the tropics and subtropics, is at risk of infection. Every year, 50-100 million clinical cases are reported, and more than 500,000 patients develop the symptoms of severe dengue infection: dengue haemorrhagic fever and dengue shock syndrome, which threaten life in Asia and Latin America. No antiviral drug for dengue is available. The dengue virus (DENV) non-structural protein 5 (NS5), which possesses the RNA-dependent RNA polymerase (RdRp) activity and is responsible for viral replication and transcription, is an attractive target for anti-dengue drug development. In the present study, 16,240 small-molecule compounds in a fragment library were screened for their capabilities to inhibit the DENV type 2 (DENV2) RdRp activities in vitro. Based on in cellulo antiviral and cytotoxity assays, we selected the compound RK-0404678 with the EC50 value of 6.0 µM for DENV2. Crystallographic analyses revealed two unique binding sites for RK-0404678 within the RdRp, which are conserved in flavivirus NS5 proteins. No resistant viruses emerged after nine rounds of serial passage of DENV2 in the presence of RK-0404678, suggesting the high genetic barrier of this compound to the emergence of a resistant virus. Collectively, RK-0404678 and its binding sites provide a new framework for antiviral drug development.

Nucleoside Analogs with Selective Antiviral Activity against Dengue Fever and Japanese Encephalitis Viruses.

Dengue virus (DENV) and Japanese encephalitis virus (JEV) are important arthropod-borne viruses from the Flaviviridae family. DENV is a global public health problem with significant social and economic impacts, especially in tropical and subtropical areas. JEV is a neurotropic arbovirus endemic to east and southeast Asia. There are no U.S. FDA-approved antiviral drugs available to treat or to prevent DENV and JEV infections, leaving nearly one-third of the world's population at risk for infection. Therefore, it is crucial to discover potent antiviral agents against these viruses. Nucleoside analogs, as a class, are widely used for the treatment of viral infections. In this study, we discovered nucleoside analogs that possess potent and selective anti-JEV and anti-DENV activities across all serotypes in cell-based assay systems. Both viruses were susceptible to sugar-substituted 2'-C-methyl analogs with either cytosine or 7-deaza-7-fluoro-adenine nucleobases. Mouse studies confirmed the anti-DENV activity of these nucleoside analogs. Molecular models were assembled for DENV serotype 2 (DENV-2) and JEV RNA-dependent RNA polymerase replication complexes bound to nucleotide inhibitors. These models show similarities between JEV and DENV-2, which recognize the same nucleotide inhibitors. Collectively, our findings provide promising compounds and a structural rationale for the development of direct-acting antiviral agents with dual activity against JEV and DENV infections.

Structure(s), function(s), and inhibition of the RNA-dependent RNA polymerase of noroviruses.

Noroviruses belong to the Caliciviridae family of single-stranded positive-sense RNA viruses. The genus Norovirus includes seven genogroups (designated GI-GVII), of which GI, GII and GIV infect humans. Human noroviruses are responsible for widespread outbreaks of acute gastroenteritis and represent one of the most common causes of foodborne illness. No vaccine or antiviral treatment options are available for norovirus infection. The RNA-dependent RNA polymerase (RdRp) of noroviruses is a key enzyme responsible for transcription and replication of the viral genome. Here, we review the progress made in understanding the structures and functions of norovirus RdRp and its use as a target for small molecule inhibitors. Crystal structures of the RdRp at different stages of substrate interaction have been determined, which shed light on its multi-step catalytic cycle. The in vitro assays and in vivo animal models that have been developed to identify and characterize inhibitors of norovirus RdRp are also summarized, followed by an update on the current antiviral research targeting different regions of norovirus RdRp. In the future, structure-based drug design and rational optimization of known nucleoside and non-nucleoside inhibitors of norovirus RdRp may pave the way towards the next generation of direct-acting antivirals.

Identification, molecular cloning and expression analysis of five RNA-dependent RNA polymerase genes in Salvia miltiorrhiza.

RNA-dependent RNA polymerases (RDRs) act as key components of the small RNA biogenesis pathways and play significant roles in post-transcriptional gene silencing (PTGS) and antiviral defense. However, there is no information about the RDR gene family in Salvia miltiorrhiza, an emerging model medicinal plant with great economic value. Through genome-wide predication and subsequent molecular cloning, five full-length S. miltiorrhiza RDR genes, termed SmRDR1-SmRDR5, were identified. The length of SmRDR cDNAs varies between 3,262 (SmRDR5) and 4,130 bp (SmRDR3). The intron number of SmRDR genes varies from 3 (SmRDR1, SmRDR3 and SmRDR4) to 17 (SmRDR5). All of the deduced SmRDR protein sequences contain the conserved RdRp domain. Moreover, SmRDR2 and SmRDR4 have an additional RRM domain. Based on the phylogenetic tree constructed with sixteen RDRs from Arabidopsis, rice and S. miltiorrhiza, plant RDRs may be divided into four groups (RDR1-RDR4). The RDR1 group contains an AtRDR and an OsRDR, while includes two SmRDRs. On the contrary, the RDR3 group contains three AtRDRs and two OsRDRs, but has only one SmRDR. SmRDRs were differentially expressed in flowers, leaves, stems and roots of S. miltiorrhiza and responsive to methyl jasmonate treatment and cucumber mosaic virus infection. The results suggest the involvement of RDRs in S. miltiorrhiza development and response to abiotic and biotic stresses. It provides a foundation for further studying the regulation and biological functions of SmRDRs and the biogenesis pathways of small RNAs in S. miltiorrhiza.

Multiple virtual screening approaches for finding new hepatitis C virus RNA-dependent RNA polymerase inhibitors: structure-based screens and molecular dynamics for the pursue of new poly pharmacological inhibitors.

The RNA polymerase NS5B of Hepatitis C virus (HCV) is a well-characterised drug target with an active site and four allosteric binding sites. This work presents a workflow for virtual screening and its application to Drug Bank screening targeting the Hepatitis C Virus (HCV) RNA polymerase non-nucleoside binding sites. Potential polypharmacological drugs are sought with predicted active inhibition on viral replication, and with proven positive pharmaco-clinical profiles. The approach adopted was receptor-based. Docking screens, guided with contact pharmacophores and neural-network activity prediction models on all allosteric binding sites and MD simulations, constituted our analysis workflow for identification of potential hits. Steps included: 1) using a two-phase docking screen with Surflex and Glide Xp. 2) Ranking based on scores, and important H interactions. 3) a machine-learning target-trained artificial neural network PIC prediction model used for ranking. This provided a better correlation of IC50 values of the training sets for each site with different docking scores and sub-scores. 4) interaction pharmacophores-through retrospective analysis of protein-inhibitor complex X-ray structures for the interaction pharmacophore (common interaction modes) of inhibitors for the five non-nucleoside binding sites were constructed. These were used for filtering the hits according to the critical binding feature of formerly reported inhibitors. This filtration process resulted in identification of potential new inhibitors as well as formerly reported ones for the thumb II and Palm I sites (HCV-81) NS5B binding sites. Eventually molecular dynamics simulations were carried out, confirming the binding hypothesis and resulting in 4 hits.

Identification of influenza virus inhibitors which disrupt of viral polymerase protein-protein interactions.

Due to their ability to rapidly mutate, influenza viruses quickly develop resistance against many antiviral substances, leading to an urgent need for new compounds. The trimeric viral polymerase complex, a major target for the development of new inhibitors, must be assembled from the PB1, PB2, and PA subunits for successful infection. Here, we describe ELISA-based assays which allow the identification of peptides which impair polymerase complex formation. Since the protein-protein interaction domains of the viral polymerase are highly conserved, these inhibitors are also predicted to be active against a broad range of influenza strains. Using this method, identification of small molecules and lead compounds against influenza A and B viruses should be feasible.