Identification and characterization of novel RdRp and Nsp15 inhibitors for SARS-COV2 using computational approach

Targeting human thymidylate synthase: Ensemble-based virtual screening for drug repositioning and the role of water

A drug repositioning computational approach was carried to search inhibitors for human thymidylate synthase. An ensemble-based virtual screening of FDA-approved drugs showed the drugs Imatinib, Lumacaftor and Naldemedine to be likely candidates for repurposing. The role of water in the drug-receptor interactions was revealed by the application of an extended AutoDock scoring function that included the water forcefield. The binding affinity scores when hydrated ligands were docked were improved in the drugs considered. Further binding free energy calculations based on the Molecular Mechanics Poisson-Boltzmann Surface Area method revealed that Imatinib, Lumacaftor and Naldemedine scored -130.7 ± 28.1, -210.6 ± 29.9 and -238.0 ± 25.4 kJ/mol, respectively, showing good binding affinity for the candidates considered. Overall, the analysis of the molecular dynamics trajectory of the receptor-drug complexes revealed stable structures for Imatinib, Lumacaftor and Naldemedine, for the entire simulation time.

Discovering new potential inhibitors to SARS-CoV-2 RNA dependent RNA polymerase (RdRp) using high throughput virtual screening and molecular dynamics simulations

RNA dependent RNA polymerase (RdRp), is an essential in the RNA replication within the life cycle of the severely acute respiratory coronavirus-2 (SARS-CoV-2), causing the deadly respiratory induced sickness COVID-19. Remdesivir is a prodrug that has seen some success in inhibiting this enzyme, however there is still the pressing need for effective alternatives. In this study, we present the discovery of four non-nucleoside small molecules that bind favorably to SARS-CoV-2 RdRp over the active form of the popular drug remdesivir (RTP) and adenosine triphosphate (ATP) by utilizing high-throughput virtual screening (HTVS) against the vast ZINC compound database coupled with extensive molecular dynamics (MD) simulations. After post-trajectory analysis, we found that the simulations of complexes containing both ATP and RTP remained stable for the duration of their trajectories. Additionally, it was revealed that the phosphate tail of RTP was stabilized by both the positive amino acid pocket and magnesium ions near the entry channel of RdRp which includes residues K551, R553, R555 and K621. It was also found that residues D623, D760, and N691 further stabilized the ribose portion of RTP with U10 on the template RNA strand forming hydrogen pairs with the adenosine motif. Using these models of RdRp, we employed them to screen the ZINC database of ~ 17 million molecules. Using docking and drug properties scoring, we narrowed down our selection to fourteen candidates. These were subjected to 200 ns simulations each underwent free energy calculations. We identified four hit compounds from the ZINC database that have similar binding poses to RTP while possessing lower overall binding free energies, with ZINC097971592 having a binding free energy two times lower than RTP.

Recent insights into the structure and function of coronavirus ribonucleases

Coronaviruses use approximately two-thirds of their 30-kb genomes to encode nonstructural proteins (nsps) with diverse functions that assist in viral replication and transcription, and evasion of the host immune response. The SARS-CoV-2 pandemic has led to renewed interest in the molecular mechanisms used by coronaviruses to infect cells and replicate. Among the 16 Nsps involved in replication and transcription, coronaviruses encode two ribonucleases that process the viral RNA-an exonuclease (Nsp14) and an endonuclease (Nsp15). In this review, we discuss recent structural and biochemical studies of these nucleases and the implications for drug discovery.

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.

Genome-wide lone strand adenine methylation in Deinococcus radiodurans R1: Regulation of gene expression through DR0643-dependent adenine methylation

DNA methylation is a covalent modification of adenine or cytosine in the genome of an organism and is found in diverse microbes including the radiation resistant bacterium Deinococcus radiodurans R1. Although earlier findings have confirmed repression or de-repression of certain genes in adenine methyltransferase (DR_0643/Dam1_DR) deficient D. radiodurans mutant however, the overall regulatory aspects of Dam1_DR-mediated adenine methylation remain mostly unexplored. In the present study, we compared the genome-wide methylome and the corresponding transcriptome of D. radiodurans WT and Δdam1 mutant to explore the correlation between methylation and gene expression. In D. radiodurans, deletion of DR_0643 ORF (Δdam1) led to hypomethylation of 512 genes resulting in differential expression of 168 genes (99 genes are upregulated and 69 genes are downregulated). The modification patterns deduced for Dam1_DR (DR_0643) and Dam2_DR (DR_2267) were non-palindromic and atypical. Moreover, we observed methylation at opportunistic sites that show adenine methylation only in D. radiodurans Δdam1 and not in D. radiodurans WT. Correlation between the methylome and transcriptome suggests that hypomethylation at Dam1_DR specific sites had both negative as well as a positive effects on gene expression. Pathways such as amino acid metabolism, transport, oxidative phosphorylation, quorum sensing, signal transduction, two-component system, glycolysis/gluconeogenesis, TCA cycle, glyoxylate and dicarboxylate metabolism were modulated by Dam1_DR-mediated adenine methylation in D. radiodurans. Processes such as DNA repair, recombination, ATPase and transmembrane transporter activity were enriched when Dam1_DR mutant was subjected to radiation stress. We further evaluated the molecular interactions and mode of binding between Dam1_DR protein and S-adenosyl methionine using molecular docking followed by MD simulation. To get a better insight into the methylation mechanism, the Dam1_DR-SAM complex was also docked with a DNA molecule to elucidate DNA-Dam1_DR structural interaction during methyl-group transfer reaction. In summary, our work presents comprehensive and integrative approaches to investigate both functional and structural aspects of DNA adenine methyltransferase (Dam1_DR) in D. radiodurans biology.

In silico drug repositioning against human NRP1 to block SARS-CoV-2 host entry

Despite COVID-19 turned into a pandemic, no approved drug for the treatment or globally available vaccine is out yet. In such a global emergency, drug repurposing approach that bypasses a costly and long-time demanding drug discovery process is an effective way in search of finding drugs for the COVID-19 treatment. Recent studies showed that SARS-CoV-2 uses neuropilin-1 (NRP1) for host entry. Here we took advantage of structural information of the NRP1 in complex with C-terminal of spike (S) protein of SARS-CoV-2 to identify drugs that may inhibit NRP1 and S protein interaction. U.S. Food and Drug Administration (FDA) approved drugs were screened using docking simulations. Among top drugs, well-tolerated drugs were selected for further analysis. Molecular dynamics (MD) simulations of drugs-NRP1 complexes were run for 100 ns to assess the persistency of binding. MM/GBSA calculations from MD simulations showed that eltrombopag, glimepiride, sitagliptin, dutasteride, and ergotamine stably and strongly bind to NRP1. In silico Alanine scanning analysis revealed that Tyr²⁹⁷, Trp³⁰¹, and Tyr³⁵³ amino acids of NRP1 are critical for drug binding. Validating the effect of drugs analyzed in this paper by experimental studies and clinical trials will expedite the drug discovery process for COVID-19.

Computational methods directed towards drug repurposing for COVID-19: advantages and limitations

Novel coronavirus disease 2019 (COVID-19) has significantly altered the socio-economic status of countries. Although vaccines are now available against the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a causative agent for COVID-19, it continues to transmit and newer variants of concern have been consistently emerging world-wide. Computational strategies involving drug repurposing offer a viable opportunity to choose a medication from a rundown of affirmed drugs against distinct diseases including COVID-19. While pandemics impede the healthcare systems, drug repurposing or repositioning represents a hopeful approach in which existing drugs can be remodeled and employed to treat newer diseases. In this review, we summarize the diverse computational approaches attempted for developing drugs through drug repurposing or repositioning against COVID-19 and discuss their advantages and limitations. To this end, we have outlined studies that utilized computational techniques such as molecular docking, molecular dynamic simulation, disease-disease association, drug-drug interaction, integrated biological network, artificial intelligence, machine learning and network medicine to accelerate creation of smart and safe drugs against COVID-19.

Exploring the Druggability of Conserved RNA Regulatory Elements in the SARS-CoV-2 Genome

SARS-CoV-2 contains a positive single-stranded RNA genome of approximately 30 000 nucleotides. Within this genome, 15 RNA elements were identified as conserved between SARS-CoV and SARS-CoV-2. By nuclear magnetic resonance (NMR) spectroscopy, we previously determined that these elements fold independently, in line with data from in vivo and ex-vivo structural probing experiments. These elements contain non-base-paired regions that potentially harbor ligand-binding pockets. Here, we performed an NMR-based screening of a poised fragment library of 768 compounds for binding to these RNAs, employing three different 1 H-based 1D NMR binding assays. The screening identified common as well as RNA-element specific hits. The results allow selection of the most promising of the 15 RNA elements as putative drug targets. Based on the identified hits, we derive key functional units and groups in ligands for effective targeting of the RNA of SARS-CoV-2.

Identifying SARS-CoV-2 antiviral compounds by screening for small molecule inhibitors of nsp15 endoribonuclease

SARS-CoV-2 is responsible for COVID-19, a human disease that has caused over 2 million deaths, stretched health systems to near-breaking point and endangered economies of countries and families around the world. Antiviral treatments to combat COVID-19 are currently lacking. Remdesivir, the only antiviral drug approved for the treatment of COVID-19, can affect disease severity, but better treatments are needed. SARS-CoV-2 encodes 16 non-structural proteins (nsp) that possess different enzymatic activities with important roles in viral genome replication, transcription and host immune evasion. One key aspect of host immune evasion is performed by the uridine-directed endoribonuclease activity of nsp15. Here we describe the expression and purification of nsp15 recombinant protein. We have developed biochemical assays to follow its activity, and we have found evidence for allosteric behaviour. We screened a custom chemical library of over 5000 compounds to identify nsp15 endoribonuclease inhibitors, and we identified and validated NSC95397 as an inhibitor of nsp15 endoribonuclease in vitro. Although NSC95397 did not inhibit SARS-CoV-2 growth in VERO E6 cells, further studies will be required to determine the effect of nsp15 inhibition on host immune evasion.

Current trends in computer aided drug design and a highlight of drugs discovered via computational techniques: A review

Computer-aided drug design (CADD) is one of the pivotal approaches to contemporary pre-clinical drug discovery, and various computational techniques and software programs are typically used in combination, in a bid to achieve the desired outcome. Several approved drugs have been developed with the aid of CADD. On SciFinder®, we evaluated more than 600 publications through systematic searching and refining, using the terms, virtual screening; software methods; computational studies and publication year, in order to obtain data concerning particular aspects of CADD. The primary focus of this review was on the databases screened, virtual screening and/or molecular docking software program used. Furthermore, we evaluated the studies that subsequently performed molecular dynamics (MD) simulations and we reviewed the software programs applied, the application of density functional theory (DFT) calculations and experimental assays. To represent the latest trends, the most recent data obtained was between 2015 and 2020, consequently the most frequently employed techniques and software programs were recorded. Among these, the ZINC database was the most widely preferred with an average use of 31.2%. Structure-based virtual screening (SBVS) was the most prominently used type of virtual screening and it accounted for an average of 57.6%, with AutoDock being the preferred virtual screening/molecular docking program with 41.8% usage. Following the screening process, 38.5% of the studies performed MD simulations to complement the virtual screening and GROMACS with 39.3% usage, was the popular MD software program. Among the computational techniques, DFT was the least applied whereby it only accounts for 0.02% average use. An average of 36.5% of the studies included reports on experimental evaluations following virtual screening. Ultimately, since the inception and application of CADD in pre-clinical drug discovery, more than 70 approved drugs have been discovered, and this number is steadily increasing over time.

Evidence of Omics, Immune Infiltration, and Pharmacogenomic for SENP1 in the Pan-Cancer Cohort

Sentrin specific-protease 1 (SENP1) is a protein involved in deSUMOylation that is almost overexpressed in cancer. SENP1 has a determinative role in the activation of transcription programs in the innate immune responses and the development B of and C lymphocytes. We found, SENP1 possibly plays a critical role in immune infiltration and acts as an expression marker in PAAD, ESCA, and THYM. CD4+ T cells, CD8+ T cells, and macrophages were more key-related immune cells, indicating that SENP1 might be introduced as a potential target for cancer immunotherapy. We further showed that dysregulation of SENP1 is powerfully associated with decreased patient survival and clinical stage. Total SENP1 protein also increases in cancer. SENP1 is also controlled by transcription factors (TFs) CREB1, KDM5A, REST, and YY1 that regulates apoptosis, cell cycle, cell proliferation, invasion, tumorigenesis, and metastasis. These TFs were in a positive correlation with SENP1. MiR-138-5p, miR-129-1-3p, and miR-129-2-3p also inhibit tumorigenesis through targeting of SENP1. The SENP1 expression level positively correlated with the expression levels of UBN1, SP3, SAP130, NUP98, NUP153 in 32 tumor types. SENP1 and correlated and binding genes: SAP130, NUP98, and NUP153 activated cell cycle. Consistent with this finding, drug analysis was indicated SENP1 is sensitive to cell cycle, apoptosis, and RTK signaling regulators. In the end, SENP1 and its expression-correlated and functional binding genes were enriched in cell cycle, apoptosis, cellular response to DNA damage stimulus. We found that the cell cycle is the main way for tumorigenesis by SENP1. SENP1 attenuates the effect of inhibitory drugs on the cell cycle. We also introduced effective FDA-Approved drugs that can inhibit SENP1. Therefore in the treatments in which these drugs are used, SENP1 inhibition is a suitable approach. This study supplies a wide analysis of the SENP1 across The Cancer Genome Atlas (CGA) cancer types. These results suggest the potential roles of SENP1 as a biomarker for cancer. Since these drugs and the drugs that cause to resistance are applied to cancer treatment, then these two class drugs can use to inhibition of SENP1.

Computational Simulations Identified Marine-Derived Natural Bioactive Compounds as Replication Inhibitors of SARS-CoV-2

The rapid spread of COVID-19, caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a worldwide health emergency. Unfortunately, to date, a very small number of remedies have been to be found effective against SARS-CoV-2 infection. Therefore, further research is required to achieve a lasting solution against this deadly disease. Repurposing available drugs and evaluating natural product inhibitors against target proteins of SARS-CoV-2 could be an effective approach to accelerate drug discovery and development. With this strategy in mind, we derived Marine Natural Products (MNP)-based drug-like small molecules and evaluated them against three major target proteins of the SARS-CoV-2 virus replication cycle. A drug-like database from MNP library was generated using Lipinski's rule of five and ADMET descriptors. A total of 2,033 compounds were obtained and were subsequently subjected to molecular docking with 3CLpro, PLpro, and RdRp. The docking analyses revealed that a total of 14 compounds displayed better docking scores than the reference compounds and have significant molecular interactions with the active site residues of SARS-CoV-2 virus targeted proteins. Furthermore, the stability of docking-derived complexes was analyzed using molecular dynamics simulations and binding free energy calculations. The analyses revealed two hit compounds against each targeted protein displaying stable behavior, binding affinity, and molecular interactions. Our investigation identified two hit compounds against each targeted proteins displaying stable behavior, higher binding affinity and key residual molecular interactions, with good in silico pharmacokinetic properties, therefore can be considered for further in vitro studies.

Potential Candidates against COVID-19 Targeting RNA-Dependent RNA Polymerase: A Comprehensive Review

Due to the extremely contagious nature of SARS-COV-2, it presents a significant threat to humans worldwide. A plethora of studies are going on all over the world to discover the drug to fight SARS-COV-2. One of the most promising targets is RNA-dependent RNA polymerase (RdRp), responsible for viral RNA replication in host cells. Since RdRp is a viral enzyme with no host cell homologs, it allows the development of selective SARS-COV-2 RdRp inhibitors. A variety of studies used in silico approaches for virtual screening, molecular docking, and repurposing of already existing drugs and phytochemicals against SARS-COV-2 RdRp. This review focuses on collating compounds possessing the potential to inhibit SARS-COV-2 RdRp based on in silico studies to give medicinal chemists food for thought so that the existing drugs can be repurposed for the control and treatment of ongoing COVID-19 pandemic after performing in vitro and in vivo experiments.

Pharmacoinformatics approach based identification of potential Nsp15 endoribonuclease modulators for SARS-CoV-2 inhibition

In the current study, a structure-based virtual screening paradigm was used to screen a small molecular database against the Non-structural protein 15 (Nsp15) endoribonuclease of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). The SARS-CoV-2 is the causative agent of the recent outbreak of coronavirus disease 2019 (COVID-19) which left the entire world locked down inside the home. A multi-step molecular docking study was performed against antiviral specific compounds (~8722) collected from the Asinex antiviral database. The less or non-interacting molecules were wiped out sequentially in the molecular docking. Further, MM-GBSA based binding free energy was estimated for 26 compounds which shows a high affinity towards the Nsp15. The drug-likeness and pharmacokinetic parameters of all 26 compounds were explored, and five molecules were found to have an acceptable pharmacokinetic profile. Overall, the Glide-XP docking score and Prime-MM-GBSA binding free energy of the selected molecules were explained strong interaction potentiality towards the Nsp15 endoribonuclease. The dynamic behavior of each molecule with Nsp15 was assessed using conventional molecular dynamics (MD) simulation. The MD simulation information was strongly favors the Nsp15 and each identified ligand stability in dynamic condition. Finally, from the MD simulation trajectories, the binding free energy was estimated using the MM-PBSA method. Hence, the proposed final five molecules might be considered as potential Nsp15 modulators for SARS-CoV-2 inhibition.