Revealing the Inhibition Mechanism of RNA-Dependent RNA Polymerase (RdRp) of SARS-CoV-2 by Remdesivir and Nucleotide Analogues: A Molecular Dynamics Simulation Study

How fingers affect folding of a thumb: Inter-subdomain cooperation in the folding of SARS-CoV-2 RdRp protein

The RNA-dependent RNA polymerase (RdRp) of SARS-CoV-2 is a critical enzyme essential for the virus's replication and transcription, making it a key therapeutic target. The RdRp protein exhibits a characteristic cupped right-hand shaped structure with two vital subdomains: the fingers and the thumb. Despite being distinct, biophysical experiments suggest that these subdomains cooperate to facilitate RNA accommodation, ensuring RdRp functionality. To investigate the structure-based mechanisms underlying the fingers-thumb interaction in both apo and RNA-bound RdRp, we constructed a coarse-grained structure-based model based on recent cryo-electron microscopy data. The simulations reveal frequent open-to-closed conformational transitions in apo RdRp, akin to a breathing-like motion. These conformational changes are regulated by the fingers-thumb association and the folding dynamics of the thumb subdomain. The thumb adopts a stable fold only when tethered by the fingers-thumb interface; when these subdomains are disconnected, the thumb transitions into an open state. A significant number of open-to-closed transition events were analyzed to generate a transition contact probability map, which highlights a few specific residues at the thumb-fingers interface, distant from the RNA accommodation sites, as essential for inducing the thumb's folding process. Given that thumb subdomain folding is critical for RNA binding and viral replication, the study proposes that these interfacial residues may function as remote regulatory switches and could be targeted for the development of allosteric drugs against SARS-CoV-2 and similar RNA viruses.

Characterization of the Conformational Hotspots of the RNA-Dependent RNA Polymerase Complex Identifies a Unique Structural Malleability of nsp8

Several antiviral therapeutic approaches have been targeted toward the RNA-dependent RNA polymerase (RdRp) complex that is involved in viral genome replication. In SARS-CoV-2, although the RdRp is a multiprotein complex, the focus has been on the ligand binding catalytic core (nonstructural protein nsp12), and not the multiprotein functional dynamics. In this study, we focus on the conformational ensembles of the RdRp complex and their modulation by the presence of RNA, performing comprehensive microsecond-scale atomistic simulations of the apo- and RNA-bound complex. We delineate the differential impact of RNA on the constituent proteins, such as conformational polymorphisms, dominant segment-specific fluctuations, and the switch in dynamical crosstalk within the complex. We distinguish dynamical signatures of nsp7, nsp8, and nsp12 in the apo-state that are reduced in the presence of the RNA and appear to "prime" the complex for activity. Importantly, we identify a unique structural malleability of the nsp8 protein with high conformational heterogeneity in the apo state, especially at three sites (Y71 for nsp8A, and D52 and A66 for nsp8B). Our work highlights the functional implications of the polymorphism of nsp8 structures and reveals possibilities for the development of allosteric inhibitors.

Computational investigation of novel synthetic analogs of C-1'β substituted remdesivir against RNA-dependent RNA-polymerase of SARS-CoV-2

Remdesivir, a C-nucleotide prodrug binds to the viral RNA-dependent-RNA polymerase (RdRp) and inhibits the viral replication by terminating RNA transcription prematurely. It is reported in literature that interaction between the C-1'β-CN moiety of Remdesivir (RDV) and the Ser861 residue in RdRp enzyme, causes a delayed chain termination during the RNA replication process and is one of the important aspect of its mechanism of action. In the pursuance of increasing the biological activity of RDV and enhancing the SAR studies, against RNA viruses, we have designed its fourteen C1'β substituted analogs, 10 -23 bearing 4/5-membered heterocyclic rings. The docking and 100 ns molecular dynamics (MD) simulations of 10-23 to the RdRp protein (PDB ID: 7L1F) revealed important interactions between 2',3'-diol, oxo group of phosphoramidate, nitrogen residues of heterocyclic rings of synthetic molecules with Arg555, Arg553, Ser759, Cys622, Asn691, Asp623 amino acid residues of protein. The docking score of 2-ethylbutyl ((S)-(((2R,3S,4R,5R)-5-(4-aminopyrrolo[2,1-f][1,2,4]triazin-7-yl)-3,4-dihydroxy-5-(1H-1,2,3-triazol-4-yl)tetrahydrofuran-2-yl)methoxy)(phenoxy)phosphoryl)-L-alaninate, 11 was found to be the higher than RDV among 14 new compounds i.e. -5.20 kcal/mol. Out of 3 compounds, 10, 12 and 13 submitted for MD simulations and Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) analysis, trifluoro-oxadiazole derivative, 13 showed higher binding energy as compared to Remdesivir. The predicted ADMET properties of 14 compounds showed their potential for being drug candidates. The present study suggests that substitution at the C1'β position by 4/5-membered rings plays an important role in the interactions between nucleoside/tide and target protein.

Unveiling the intricacies of allosteric regulation in aspartate kinase from the Wolbachia endosymbiont of Brugia Malayi: Mechanistic and therapeutic insights

Aspartate kinase (AK), an enzyme from the Wolbachia endosymbiont of Brugia malayi (WBm), plays a pivotal role in the bacterial cell wall and amino acid biosynthesis, rendering it an attractive candidate for therapeutic intervention. Allosteric inhibition of aspartate kinase is a prevalent mode of regulation across microorganisms and plants, often modulated by end products such as lysine, threonine, methionine, or meso-diaminopimelate. The intricate and diverse nature of microbial allosteric regulation underscores the need for rigorous investigation. This study employs a combined experimental and computational approach to decipher the allosteric regulation of WBmAK. Molecular Dynamics (MD) simulations elucidate that ATP (cofactor) and ASP (substrate) binding induce a closed conformation, promoting enzymatic activity. In contrast, the binding of lysine (allosteric inhibitor) leads to enzyme inactivation and an open conformation. The enzymatic assay demonstrates the optimal activity of WBmAK at 28 °C and a pH of 8.0. Notably, the allosteric inhibition study highlights lysine as a more potent inhibitor compared to threonine. Importantly, this investigation sheds light on the allosteric mechanism governing WBmAK and imparts novel insights into structure-based drug discovery, paving the way for the development of effective inhibitors against filarial pathogens.

Ab initio Modeling of Hydrogen Bonding of Remdesivir and Adenosine with Uridine

Remdesivir (RDV) emerged as an effective drug against the SARS-CoV-2 virus pandemic. One of the crucial steps in the mechanism of action of RDV is its incorporation into the growing RNA strand. RDV, an adenosine analogue, forms Watson-Crick (WC) type hydrogen bonds with uridine in the complementary strand and the strength of this interaction will control efficacy of RDV. While there is a plethora of structural and energetic information available about WC H-bonds in natural base pairs, the interaction of RDV with uridine has not been studied yet at the atomic level. In this article, we aim to bridge this gap, to understand RDV and its hydrogen bonding interactions, by employing density functional theory (DFT) at the M06-2X/cc-pVDZ level. The interaction energy, QTAIM analysis, NBO and SAPT2 are performed for RDV, adenosine, and their complex with uridine to gain insights into the nature of hydrogen bonding. The computations show that RDV has similar geometry, energetic, molecular orbitals, and aromaticity as adenosine, suggesting that RDV is an effective adenosine analogue. The important geometrical parameters, such as bond distances and red-shift in the stretching vibrational modes of adenosine, RDV and uridine identify two WC-type H-bonds. The relative strength of these two H-bonds is computed using QTAIM parameters and the computed hydrogen bond energy. Finally, the SAPT2 study is performed at the minima and at non-equilibrium base pair distances to understand the dominant intermolecular physical force. This study, based on a thorough analysis of a variety of computations, suggests that both adenosine and RDV have similar structure, energetic, and hydrogen bonding behaviour.

Molecular dynamics simulation study on the binding mechanism between carbon nanotubes and RNA-dependent RNA polymerase

Carbon nanotubes (CNTs) have potential prospects in disease treatment, so it is of great significance to study CNTs as the possible inhibitors of RNA-dependent RNA polymerase (RdRp). Through the way of using the RdRp of SARS-COV-2 as a model, five armchair single-walled carbon nanotubes (SWCNTs) (namely Dn, which stands for CNTs (n, m = n), n = 3-7) and RdRp have been selected to study the interactions by means of molecular docking and molecular dynamics simulation. After five SWCNT-RdRp complex systems have been subjected to the molecular dynamics simulations of 100 ns, and Molecular Mechanics Poisson - Boltzmann Surface Area (MMPBSA) has been used to calculate the binding free energy, it is found that the binding free energy of the D6 system (-189.541 kJ/mol) is significantly higher than that of the other four systems, and most of the amino acids with strong positive effects on binding are usually basic amino acids. What's more, in the further investigation of the specific interaction mechanism between CNT (6,6) and RdRp, it is revealed that the three amino acid residues LYS545, ARG553 and ARG555 located in the nucleoside triphosphate (NTP) entry channel all have strong effects. In addition, it is also observed that when ARG555 has been inserted into SWCNT, a stable structure will be formed, which will break the original NTP entry channel structure and inhibit virus replication. Therefore, it can be concluded that certain specific types of SWCNT, such as CNT (6,6), could be potential small molecule inhibitors in the treatment of coronavirus.Communicated by Ramaswamy H. Sarma.

RNA-dependent RNA polymerase of SARS-CoV-2 regulate host mRNA translation efficiency by hijacking eEF1A factors

The RNA-dependent RNA polymerase (NSP12) of COVID-19 plays a significant role in the viral infection process, which promotes viral RNA replication by cooperating with NSP7 and NSP8, but little is known about its regulation on the function of host cells. We firstly found that overexpression of NSP12 had little effect on host mRNAs transcription. Using iCLIP technology, we found that NSP12 can bind a series of host RNAs with the conserved binding motif G(C/A/G)(U/G/A)UAG, especially ribosomal RNA. We found that NSP12 could directly bind to eEF1A factor via the NIRAN domain of NSP12 and N-terminal domain of eEF1A. NSP12 colocalized with eEF1A to inhibit type I interferon expression upon virus infection. In order to prove that NSP12 regulates the translation level of host cells, we found that NSP12 significantly affected the translation efficiency of many host mRNAs (such as ISG15, NF-κB2, ILK and SERPINI2) via ribosome profiling experiment, and the genes with significant upregulation in translation efficiency were mainly enriched in positive regulation of ubiquitin-dependent proteasomal process and NIK/NF-κB signaling pathway (such as NF-κB2, ILK), and negative regulation of type I interferon production, protein level of these genes were further confirmed in HEK293T and Calu3 cells upon NSP12 overexpression. These results indicate that NSP12 of SARS-CoV-2 can hijack the eEF1A factor to regulate translation efficiency of host mRNAs, which provides a new idea for us to evaluate the impact of SARS-CoV2 virus on the host and study the potential drug targets.

De Novo Potent Peptide Nucleic Acid Antisense Oligomer Inhibitors Targeting SARS-CoV-2 RNA-Dependent RNA Polymerase via Structure-Guided Drug Design

Global reports of novel SARS-CoV-2 variants and recurrence cases continue despite substantial vaccination campaigns, raising severe concerns about COVID-19. While repurposed drugs offer some treatment options for COVID-19, notably, nucleoside inhibitors like Remdesivir stand out as curative therapies for COVID-19 that are approved by the US Food and Drug Administration (FDA). The emergence of highly contagious SARS-CoV-2 variants underscores the imperative for antiviral drugs adaptable to evolving viral mutations. RNA-dependent RNA polymerase (RdRp) plays a key role in viral genome replication. Currently, inhibiting viral RdRp function remains a pivotal strategy to tackle the notorious virus. Peptide nucleic acid (PNA) therapy shows promise by effectively targeting specific genome regions, reducing viral replication, and inhibiting infection. In our study, we designed PNA antisense oligomers conjugated with cell-penetrating peptides (CPP) aiming to evaluate their antiviral effects against RdRp target using structure-guided drug design, which involves molecular docking simulations, drug likeliness and pharmacokinetic evaluations, molecular dynamics simulations, and computing binding free energy. The in silico analysis predicts that chemically modified PNAs might act as antisense molecules in order to disrupt ribosome assembly at RdRp's translation start site, and their chemically stable and neutral backbone might enhance sequence-specific RNA binding interaction. Notably, our findings demonstrate that PNA-peptide conjugates might be the most promising inhibitors of SARS-CoV-2 RdRp, with superior binding free energy compared to Remdesivir in the current COVID-19 medication. Specifically, PNA-CPP-1 could bind simultaneously to the active site residues of RdRp protein and sequence-specific RdRp-RNA target in order to control viral replication.

Main and papain-like proteases as prospective targets for pharmacological treatment of coronavirus SARS-CoV-2

The pandemic caused by the coronavirus SARS-CoV-2 led to a global crisis in the world healthcare system. Despite some progress in the creation of antiviral vaccines and mass vaccination of the population, the number of patients continues to grow because of the spread of new SARS-CoV-2 mutations. There is an urgent need for direct-acting drugs capable of suppressing or stopping the main mechanisms of reproduction of the coronavirus SARS-CoV-2. Several studies have shown that the successful replication of the virus in the cell requires proteolytic cleavage of the protein structures of the virus. Two proteases are crucial in replicating SARS-CoV-2 and other coronaviruses: the main protease (Mpro) and the papain-like protease (PLpro). In this review, we summarize the essential viral proteins of SARS-CoV-2 required for its viral life cycle as targets for chemotherapy of coronavirus infection and provide a critical summary of the development of drugs against COVID-19 from the drug repurposing strategy up to the molecular design of novel covalent and non-covalent agents capable of inhibiting virus replication. We overview the main antiviral strategy and the choice of SARS-CoV-2 Mpro and PLpro proteases as promising targets for pharmacological impact on the coronavirus life cycle.

SARS-CoV-2 RNA-Dependent RNA Polymerase Follows Asynchronous Translocation Pathway for Viral Transcription and Replication

Translocation is one essential step for the SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) to exert viral replication and transcription. Although cryo-EM structures of SARS-CoV-2 RdRp are available, the molecular mechanisms of dynamic translocation remain elusive. Herein, we constructed a Markov state model based on extensive molecular dynamics simulations to elucidate the translocation dynamics of the SARS-CoV-2 RdRp. We identified two intermediates that pinpoint the rate-limiting step of translocation and characterize the asynchronous movement of the template-primer duplex. The 3'-terminal nucleotide in the primer strand lags behind due to the uneven distribution of protein-RNA interactions, while the translocation of the template strand is delayed by the hurdle residue K500. Even so, the two strands share the same "ratchet" to stabilize the polymerase in the post-translocation state, suggesting a Brownian-ratchet model. Overall, our study provides intriguing insights into SARS-CoV-2 replication and transcription, which would open a new avenue for drug discoveries.

Computational assessment of herbal medicine-derived compounds as potential inhibitors of SARS-CoV-2 main protease

Since the main protease (Mpro) is crucial for the COVID-19 virus replication and transcription, searching for Mpro inhibitors is one possible treatment option. In our study, 258 small molecules were collected from lung-related herbal medicines, and their structures were optimized with the B3LYP-D3/6-31G* method. After the molecular docking with Mpro, we selected the top 20 compounds for the further geometry optimization with the larger basis sets. After the further molecular docking, the top eight compounds were screened out. Then we performed molecular dynamics simulations and binding free energy calculations to determine stability of the complexes. Our results show that mulberrofuran G, Xambioona, and kuwanon D can bind Mpro well. In quantum chemistry studies, such as ESP and CDFT analyses, the compounds properties are predicted. Additionally, the drug-likeness analyses and ADME studies on these three candidate compounds verified that all of them conform to Libinski's rule and may be drug-like compounds.

Polypharmacology guided drug repositioning approach for SARS-CoV2

Drug repurposing has emerged as an important strategy and it has a great potential in identifying therapeutic applications for COVID-19. An extensive virtual screening of 4193 FDA approved drugs has been carried out against 24 proteins of SARS-CoV2 (NSP1-10 and NSP12-16, envelope, membrane, nucleoprotein, spike, ORF3a, ORF6, ORF7a, ORF8, and ORF9b). The drugs were classified into top 10 and bottom 10 drugs based on the docking scores followed by the distribution of their therapeutic indications. As a result, the top 10 drugs were found to have therapeutic indications for cancer, pain, neurological disorders, and viral and bacterial diseases. As drug resistance is one of the major challenges in antiviral drug discovery, polypharmacology and network pharmacology approaches were employed in the study to identify drugs interacting with multiple targets and drugs such as dihydroergotamine, ergotamine, bisdequalinium chloride, midostaurin, temoporfin, tirilazad, and venetoclax were identified among the multi-targeting drugs. Further, a pathway analysis of the genes related to the multi-targeting drugs was carried which provides insight into the mechanism of drugs and identifying targetable genes and biological pathways involved in SARS-CoV2.

Structure-Based Drug Design of RdRp Inhibitors against SARS-CoV-2

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused a worldwide pandemic since 2019, spreading rapidly and posing a significant threat to human health and life. With over 6 billion confirmed cases of the virus, the need for effective therapeutic drugs has become more urgent than ever before. RNA-dependent RNA polymerase (RdRp) is crucial in viral replication and transcription, catalysing viral RNA synthesis and serving as a promising therapeutic target for developing antiviral drugs. In this article, we explore the inhibition of RdRp as a potential treatment for viral diseases, analysing the structural information of RdRp in virus proliferation and summarizing the reported inhibitors' pharmacophore features and structure-activity relationship profiles. We hope that the information provided by this review will aid in structure-based drug design and aid in the global fight against SARS-CoV-2 infection.

Biochemical interaction of human hemoglobin with ionic liquids of noscapinoids: Spectroscopic and computational approach

In this work, we have developed noscapine based ionic liquids i.e., Noscapine (MeNOS) and 9-Bromonoscapine (MeBrNOS) as cation supported with bis(trifluoromethylsulfonyl)amide (NTf₂^-) as anion. We have reported the mechanism of binding interaction between noscapine based ILs and human hemoglobin (Hb) using various spectroscopic and computational techniques. The corresponding thermodynamics studies showed that the binding is exothermic in nature and major forces responsible for binding are Van der waals and hydrogen bonding interaction. The fluorescence spectra showed that the intensity of Hb decreases in the presence of [MeNOS]NTf₂ and [MeBrNOS]NTf₂ both shows static quenching. The secondary structural changes in Hb were observed and calculated by using CD spectroscopy. Molecular docking studies revealed that both the ILs show strong binding in β1 fragment of the tetrameric structure of Hb, but the binding of [MeNOS]NTf₂ is relatively stronger than [MeBrNOS]NTf₂ and the results are supported by MD simulations.

On the Sensitivity and Affinity of Gold, Silver, and Platinum Surfaces against the SARS-CoV-2 Virus: A Comparative Computational Study

The novel coronavirus disease and its complications have motivated the design of new sensors with the highest sensitivity, and affinity for the detection of the SARS-CoV-2 virus is considered in many research studies. In this research article, we employ full atomistic molecular dynamics (MD) models to study the interactions between the receptor binding domain (RBD) and spike protein of the coronavirus and different metals such as gold (Au), platinum (Pt), and silver (Ag) to analyze their sensitivity against this virus. The comparison between the RBD interactions with ACE2 (angiotensin-converting enzyme 2) and different metals indicates that metals have remarkable effects on the structural features and dynamical properties of the RBD. The binding site of the RBD has more affinity to the surfaces of gold, platinum, and silver than to the other parts of the protein. Moreover, the initial configuration of the RBD relative to the metal surface plays an important role in the stability of metal complexes with the RBD. The binding face of the protein to the metal surface has been changed in the presence of different metals. In other words, the residues of the RBD that participate in RBD interactions with the metals are different irrespective of the initial configurations in which the [Asn, Thr, Tyr], [Ser, Thr, Tyr], and [Asn, Asp, Tyr] residues of the protein have a greater affinity to Ag, Au, and Pt, respectively. The corresponding metals have a considerable affinity to the RBD, which due to strong interactions with the protein can change the secondary structure and structural features. Based on the obtained results during the complexation process between the protein and metals, the helical structure of the protein changes to the bend and antiparallel β-sheets. The calculated binding energies for the RBD complexes with silver, gold, and platinum are -95.03, -138.03, and -133.96 kcal·mol-1, respectively. The adsorption process of the spike protein on the surfaces of different metals represents similar results and indicates that the entire spike protein of the coronavirus forms a more stable complex with the gold surface compared with other metals. Moreover, the RBD of the spike protein has more interactions with the surfaces than with the other parts of the protein. Therefore, it is possible to predict the properties of the coronavirus on the metal surface based on the dynamical behavior of the RBD. Overall, our computational results confirm that the gold surface can be considered as an outstanding substrate for developing new sensors with the highest sensitivity against SARS-CoV-2.

In silico identification of novel benzophenone-coumarin derivatives as SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) inhibitors

In this study, we propose our novel benzophenone-coumarin derivatives (BCDs) as potent inhibitors of the RNA-dependent RNA polymerase (RdRp) of SARS-CoV-2 virus, one of the key targets that are involved in the viral genome replication. We aim to evaluate the in silico antiviral potential of BCDs against this protein target, which involves molecular docking simulations, druglikeliness and pharmacokinetic evaluations, PASS analysis, molecular dynamics simulations, and computing binding free energy. Out of all the BCDs screened through these parameters, BCD-8 was found to be the most efficient and potent inhibitor of SARS-CoV-2 RdRp. During molecular docking simulation, BCD-8 showed an extensive molecular interaction in comparison with that of the standard control used, remdesivir. The druglikeliness and pharmacokinetic analyses also proved the efficiency of BCD-8 as an effective drug without adverse effects. Further, pharmacological potential analysis through PASS depicted the antiviral property of BCD-8. With these findings, we performed molecular dynamics simulations, where BCD-8 edged out remdesivir with its exemplary stable interaction with SARS-CoV-2 RdRp. Furthermore, binding free energy of both BCD-8 and remdesivir was calculated, where BCD-8 showed a lower binding energy and standard deviations in comparison with that of remdesivir. Moreover, being a non-nucleoside analogue, BCD-8 can be used effectively against SARS-CoV-2, whereas nucleoside analogues like remdesivir may become non-functional or less functional due to exonuclease activity of nsp14 of the virus. Therefore, we propose BCD-8 as a SARS-CoV-2 RdRp inhibitor, showing higher predicted efficiency than remdesivir in all the in silico experiments conducted.Communicated by Ramaswamy H. Sarma.

Arjunetin as a promising drug candidate against SARS-CoV-2: molecular dynamics simulation studies

Stem and bark of the tree Terminalia arjuna Wight & Arn. (Combretaceae) has been documented to exhibit therapeutic properties like cardiotonic, anticancer, antiviral, antibacterial, antifungal, hypercholesterolemia, hypolipidemic, and anti-coagulant. Our previous studies have shown that, ethanolic extract of T. arjuna bark exhibits radical scavenging anti-oxidant activity and also effectively inhibited catalase activity. In this study, oleanane triterpenoids type compounds viz., oleanolic acid, arjunolic acid, arjunolitin, arjunetin were isolated from ethanolic bark extract as bio-active compound and their structures were elucidated using ¹H, ¹³C NMR, HR-ESIMS, IR. Of the various compounds, Arjunetin showed significant inhibition of catalase activity as compared to the other compounds. Based on the structural similarity between arjunetin and current antiviral drugs, we propose that arjunetin might exhibit antiviral activity. Molecular docking and molecular dynamics studies showed that arjunetin binds to the binds to key targets of SARS-CoV-2 namely, 3CLpro, PLpro, and RdRp) with a higher binding energy values (3CLpro, -8.4 kcal/mol; PLpro, -7.6 kcal/mol and RdRp, -8.1 kcal/mol) as compared with FDA approved protease inhibitor drugs to Lopinavir (3CLpro, -7.2 kcal/mole and PLpro -7.7 kcal/mole) and Remdesivir (RdRp -7.6 kcal/mole). To further investigate this, we performed 200-500 ns molecular dynamics simulation studies. The results transpired that the binding affinity of Arjunetin is higher than Remdesivir in the RNA binding cavity of RdRp. Based on structural similarity between arjunetin and Saikosaponin (a known antiviral agents) and based on our molecular docking and molecular dynamic simulation studies, we propose that arjunetin can be a promising drug candidate against Covid-19.Communicated by Ramaswamy H. Sarma.

Cheminformatics-Based Discovery of Potential Chemical Probe Inhibitors of Omicron Spike Protein

During the past two decades, the world has witnessed the emergence of various SARS-CoV-2 variants with distinct mutational profiles influencing the global health, economy, and clinical aspects of the COVID-19 pandemic. These variants or mutants have raised major concerns regarding the protection provided by neutralizing monoclonal antibodies and vaccination, rates of virus transmission, and/or the risk of reinfection. The newly emerged Omicron, a genetically distinct lineage of SARS-CoV-2, continues its spread in the face of rising vaccine-induced immunity while maintaining its replication fitness. Efforts have been made to improve the therapeutic interventions and the FDA has issued Emergency Use Authorization for a few monoclonal antibodies and drug treatments for COVID-19. However, the current situation of rapidly spreading Omicron and its lineages demands the need for effective therapeutic interventions to reduce the COVID-19 pandemic. Several experimental studies have indicated that the FDA-approved monoclonal antibodies are less effective than antiviral drugs against the Omicron variant. Thus, in this study, we aim to identify antiviral compounds against the Spike protein of Omicron, which binds to the human angiotensin-converting enzyme 2 (ACE2) receptor and facilitates virus invasion. Initially, docking-based virtual screening of the in-house database was performed to extract the potential hit compounds against the Spike protein. The obtained hits were optimized by DFT calculations to determine the electronic properties and molecular reactivity of the compounds. Further, MD simulation studies were carried out to evaluate the dynamics of protein–ligand interactions at an atomistic level in a time-dependent manner. Collectively, five compounds (AKS-01, AKS-02, AKS-03, AKS-04, and AKS-05) with diverse scaffolds were identified as potential hits against the Spike protein of Omicron. Our study paves the way for further in vitro and in vivo studies.

Practical and Highly Efficient Synthesis of Remdesivir from GS-441524

A three-step sequence for preparing remdesivir, an important anti-SARS-CoV-2 drug, is described. Employing N,N-dimethylformamide dimethyl acetal (DMF-DMA) as a protecting agent, this synthesis started from (2R,3R,4S,5R)-2-(4-aminopyrrolo[2,1-f][1,2,4]triazin-7-yl)-3,4-dihydroxy-5-(hydroxymethyl)tetrahydro-furan-2-carbonitrile (GS-441524) and consisted of three reactions, including protection, phosphoramidation, and deprotection. The advantages of this approach are as follows: (1) the protecting group could be removed under a mild deprotection condition, which avoided the generation of the degraded impurity; (2) high stereoselectivity was achieved in the phosphorylated reaction; (3) this synthesis could be performed successively without purification of intermediates. Moreover, the overall yield of this approach on a gram scale could be up to 85% with an excellent purity of 99.4% analyzed by high-performance liquid chromatography (HPLC).

Mechanism of reaction of RNA-dependent RNA polymerase from SARS-CoV-2

We combine molecular dynamics, statistical mechanics, and hybrid quantum mechanics/molecular mechanics simulations to describe mechanistically the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA-dependent RNA polymerase (RdRp). Our study analyzes the binding mode of both natural triphosphate substrates as well as remdesivir triphosphate (the active form of drug), which is bound preferentially over ATP by RdRp while being poorly recognized by human RNA polymerase II (RNA Pol II). A comparison of incorporation rates between natural and antiviral nucleotides shows that remdesivir is incorporated more slowly into the nascent RNA compared with ATP, leading to an RNA duplex that is structurally very similar to an unmodified one, arguing against the hypothesis that remdesivir is a competitive inhibitor of ATP. We characterize the entire mechanism of reaction, finding that viral RdRp is highly processive and displays a higher catalytic rate of incorporation than human RNA Pol II. Overall, our study provides the first detailed explanation of the replication mechanism of RdRp.

Whole genome analysis and homology modeling of SARS-CoV-2 Indian isolate reveals potent FDA approved drug choice for treating COVID-19

Coronaviruses have caused enough devastation in the last two decades. These viruses have some rare features while sharing some common features. Novel coronavirus disease (nCoV-19) caused an outbreak with a fatality rate of 5%. It emerged from China and spread into many countries. The present research focused on genome analysis of Indian nCoV-19 Isolate and its translational product subjected to homology modeling and its subsequent molecular simulations to find out potent FDA approved drug for treating COVID-19. Phylogenetic analysis of SARS-CoV-2 Indian isolate shows close resemblance with 17 countries SARS-CoV-2 isolates. Homology modeling of four non-structural proteins translational product of Indian SARS-CoV-2 genome shows high similarity and allowed regions with the existing PDB deposited SARS-CoV-2 target proteins. Finally, these four generated proteins show more affinity with cobicistat, remdesivir and indinavir out of 14 screened FDA approved drugs in molecular docking which is further proven by molecular dynamics simulation and MMGBSA analysis of target ligand complex with best simulation trajectories. Overall our present research findings is that three proposed drugs namely cobicistat, remdesivir and indinavir showed higher interaction with the model SARS-CoV-2 viral target proteins from the Indian nCoV-19 isolate. These compounds could be used as a starting point for the creation of active antiviral drugs to combat the deadly COVID-19 virus during global pandemic and its subsequent viral infection waves across the globe.Communicated by Ramaswamy H. Sarma.

Synthesis of Nucleoside-like Molecules from a Pyrolysis Product of Cellulose and Their Computational Prediction as Potential SARS-CoV-2 RNA-Dependent RNA Polymerase Inhibitors

(1R,5S)-1-Hydroxy-3,6-dioxa-bicyclo[3.2.1]octan-2-one, available by an efficient catalytic pyrolysis of cellulose, has been applied as a chiral building block in the synthesis of seven new nucleoside analogues, with structural modifications on the nucleobase moiety and on the carboxyl- derived unit. The inverted configuration by Mitsunobu reaction used in their synthesis was verified by 2D-NOESY correlations, supported by the optimized structure employing the DFT methods. An in silico screening of these compounds as inhibitors of SARS-CoV-2 RNA-dependent RNA polymerase has been carried out in comparison with both remdesivir, a mono-phosphoramidate prodrug recently approved for COVID-19 treatment, and its ribonucleoside metabolite GS-441524. Drug-likeness prediction and data by docking calculation indicated compound 6 [=(3S,5S)-methyl 5-(hydroxymethyl)-3-(6-(4-methylpiperazin-1-yl)-9H-purin-9-yl)tetrahydrofuran-3-carboxylate] as the best candidate. Furthermore, molecular dynamics simulation showed a stable interaction of structure 6 in RNA-dependent RNA polymerase (RdRp) complex and a lower average atomic fluctuation than GS-441524, suggesting a well accommodation in the RdRp binding pocket.

Applications of density functional theory in COVID-19 drug modeling

The rapidly evolving Coronavirus 2019 (COVID-19) pandemic has led to millions of deaths around the world, highlighting the pressing need to develop effective antiviral pharmaceuticals. Recent efforts with computer-aided rational drug discovery have allowed detailed examination of drug–macromolecule interactions primarily by molecular mechanics (MM) techniques. Less widely applied in COVID-19 drug modeling is density functional theory (DFT), a quantum mechanics (QM) method that enables electronic structure calculations and elucidations of reaction mechanisms. Here, we review recent advances in applying DFT in molecular modeling studies of COVID-19 pharmaceuticals. We start by providing an overview of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) drugs and targets, followed by a brief introduction to DFT. We then provide a discussion of different approaches by which DFT has been applied. Finally, we discuss essential factors to consider when incorporating DFT in future drug modeling research.

WaterMap and Molecular Dynamic Simulation-Guided Discovery of Potential PAK1 Inhibitors Using Repurposing Approaches

p21-Activated kinase 1 (PAK1) is positioned at the nexus of several oncogenic signaling pathways. Currently, there are no approved inhibitors for disabling the transfer of phosphate in the active site directly, as they are limited by lower affinity, and poor kinase selectivity. In this work, a repurposing study utilizing FDA-approved drugs from the DrugBank database was pursued with an initial selection of 27 molecules out of ∼2162 drug molecules, based on their docking energies and molecular interaction patterns. From the molecules that were considered for WaterMap analysis, seven molecules, namely, Mitoxantrone, Labetalol, Acalabrutinib, Sacubitril, Flubendazole, Trazodone, and Niraparib, ascertained the ability to overlap with high-energy hydration sites. Considering many other displaced unfavorable water molecules, only Acalabrutinib, Flubendazole, and Trazodone molecules highlighted their prominence in terms of binding affinity gains through ΔΔG that ranges between 6.44 and 2.59 kcal/mol. Even if Mitoxantrone exhibited the highest docking score and greater interaction strength, it did not comply with the WaterMap and molecular dynamics simulation results. Moreover, detailed MD simulation trajectory analyses suggested that the drug molecules Flubendazole, Niraparib, and Acalabrutinib were highly stable, observed from their RMSD values and consistent interaction pattern with Glu315, Glu345, Leu347, and Asp407 including the hydrophobic interactions maintained in the three replicates. However, the drug molecule Trazodone displayed a loss of crucial interaction with Leu347, which was essential to inhibit the kinase activity of PAK1. The molecular orbital and electrostatic potential analyses elucidated the reactivity and strong complementarity potentials of the drug molecules in the binding pocket of PAK1. Therefore, the CADD-based reposition efforts, reported in this work, helped in the successful identification of new PAK1 inhibitors that requires further investigation by in vitro analysis.

Incorporation efficiency and inhibition mechanism of 2'-substituted nucleotide analogs against SARS-CoV-2 RNA-dependent RNA polymerase

The ongoing pandemic caused by SARS-CoV-2 emphasizes the need for effective therapeutics. Inhibition of SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) by nucleotide analogs provides a promising antiviral strategy. One common group of RdRp inhibitors, 2'-modified nucleotides, are reported to exhibit different behaviors in the SARS-CoV-2 RdRp transcription assay. Three of these analogs, 2'-O-methyl UTP, Sofosbuvir, and 2'-methyl CTP, act as effective inhibitors in previous biochemical experiments, while Gemcitabine and ara-UTP show no inhibitory activity. To understand the impact of the 2'-modification on their inhibitory effects, we conducted extensive molecular dynamics simulations and relative binding free energy calculations using the free energy perturbation method on SARS-CoV-2 replication-transcription complex (RTC) with these five nucleotide analogs. Our results reveal that the five nucleotide analogs display comparable binding affinities to SARS-CoV-2 RdRp and they can all be added to the nascent RNA chain. Moreover, we examine how the incorporation of these nucleotide triphosphate (NTP) analogs will impact the addition of the next nucleotide. Our results indicate that 2'-O-methyl UTP can weaken the binding of the subsequent NTP and consequently lead to partial chain termination. Additionally, Sofosbuvir and 2'-methyl CTP can cause immediate termination due to the strong steric hindrance introduced by their bulky 2'-methyl groups. In contrast, nucleotide analogs with smaller substitutions, such as the fluorine atoms and the ara-hydroxyl group in Gemcitabine and ara-UTP, have a marginal impact on the polymerization process. Our findings are consistent with experimental observations, and more importantly, shed light on the detailed molecular mechanism of SARS-CoV-2 RdRp inhibition by 2'-substituted nucleotide analogs, and may facilitate the rational design of antiviral agents to inhibit SARS-CoV-2 RdRp.

"Bucket brigade" using lysine residues in RNA-dependent RNA polymerase of SARS-CoV-2

The RNA-dependent RNA polymerase (RdRp) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a promising drug target for coronavirus disease 2019 (COVID-19) because it plays the most important role in the replication of the RNA genome. Nucleotide analogs such as remdesivir and favipiravir are thought to interfere with the RNA replication by RdRp. More specifically, they are expected to compete with nucleoside triphosphates, such as ATP. However, the process in which these drug molecules and nucleoside triphosphates are taken up by RdRp remains unknown. In this study, we performed all-atom molecular dynamics simulations to clarify the recognition mechanism of RdRp for these drug molecules and ATP that were at a distance. The ligand recognition ability of RdRp decreased in the order of remdesivir, favipiravir, and ATP. We also identified six recognition paths. Three of them were commonly found in all ligands, and the remaining three paths were ligand-dependent ones. In the common two paths, it was observed that the multiple lysine residues of RdRp carried the ligands to the binding site like a "bucket brigade." In the remaining common path, the ligands directly reached the binding site. Our findings contribute to the understanding of the efficient ligand recognition by RdRp at the atomic level.

Co-crystallization and structure determination: An effective direction for anti-SARS-CoV-2 drug discovery

Safer and more-effective drugs are urgently needed to counter infections with the highly pathogenic SARS-CoV-2, cause of the COVID-19 pandemic. Identification of efficient inhibitors to treat and prevent SARS-CoV-2 infection is a predominant focus. Encouragingly, using X-ray crystal structures of therapeutically relevant drug targets (PLpro, Mpro, RdRp, and S glycoprotein) offers a valuable direction for anti-SARS-CoV-2 drug discovery and lead optimization through direct visualization of interactions. Computational analyses based primarily on MMPBSA calculations have also been proposed for assessing the binding stability of biomolecular structures involving the ligand and receptor. In this study, we focused on state-of-the-art X-ray co-crystal structures of the abovementioned targets complexed with newly identified small-molecule inhibitors (natural products, FDA-approved drugs, candidate drugs, and their analogues) with the assistance of computational analyses to support the precision design and screening of anti-SARS-CoV-2 drugs.

Recent Developments in Free Energy Calculations for Drug Discovery

The grand challenge in structure-based drug design is achieving accurate prediction of binding free energies. Molecular dynamics (MD) simulations enable modeling of conformational changes critical to the binding process, leading to calculation of thermodynamic quantities involved in estimation of binding affinities. With recent advancements in computing capability and predictive accuracy, MD based virtual screening has progressed from the domain of theoretical attempts to real application in drug development. Approaches including the Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA), Linear Interaction Energy (LIE), and alchemical methods have been broadly applied to model molecular recognition for drug discovery and lead optimization. Here we review the varied methodology of these approaches, developments enhancing simulation efficiency and reliability, remaining challenges hindering predictive performance, and applications to problems in the fields of medicine and biochemistry.

A review on the interaction of nucleoside analogues with SARS-CoV-2 RNA dependent RNA polymerase

The outbreaks of severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) in 2019, have highlighted the concerns about the lack of potential vaccines or antivirals approved for inhibition of CoVs infection. SARS-CoV-2 RNA dependent RNA polymerase (RdRp) which is almost preserved across different viral species can be a potential target for development of antiviral drugs, including nucleoside analogues (NA). However, ExoN proofreading activity of CoVs leads to their protection from several NAs. Therefore, potential platforms based on the development of efficient NAs with broad-spectrum efficacy against human CoVs should be explored. This study was then aimed to present an overview on the development of NAs-based drug repurposing for targeting SARS-CoV-2 RdRp by computational analysis. Afterwards, the clinical development of some NAs including Favipiravir, Sofosbuvir, Ribavirin, Tenofovir, and Remdesivir as potential inhibitors of RdRp, were surveyed. Overall, exploring broad-spectrum NAs as promising inhibitors of RdRp may provide useful information about the identification of potential antiviral repurposed drugs against SARS-CoV-2.

Molecular dynamics and in silico mutagenesis on the reversible inhibitor-bound SARS-CoV-2 main protease complexes reveal the role of lateral pocket in enhancing the ligand affinity

The 2019 novel coronavirus pandemic caused by SARS-CoV-2 remains a serious health threat to humans and there is an urgent need to develop therapeutics against this deadly virus. Recent scientific evidences have suggested that the main protease (Mpro) enzyme in SARS-CoV-2 can be an ideal drug target due to its crucial role in the viral replication and transcription processes. Therefore, there are ongoing research efforts to identify drug candidates against SARS-CoV-2 Mpro that resulted in hundreds of X-ray crystal structures of ligand-bound Mpro complexes in the Protein Data Bank (PDB) describing the interactions of different fragment chemotypes within different sites of the Mpro. In this work, we performed rigorous molecular dynamics (MD) simulation of 62 reversible ligand-Mpro complexes in the PDB to gain mechanistic insights about their interactions at the atomic level. Using a total of over 3 µs long MD trajectories, we characterized different pockets in the apo Mpro structure, and analyzed the dynamic interactions and binding affinity of ligands within those pockets. Our results identified the key residues that stabilize the ligands in the catalytic sites and other pockets of Mpro. Our analyses unraveled the role of a lateral pocket in the catalytic site in Mpro that is critical for enhancing the ligand binding to the enzyme. We also highlighted the important contribution from HIS163 in the lateral pocket towards ligand binding and affinity against Mpro through computational mutation analyses. Further, we revealed the effects of explicit water molecules and Mpro dimerization in the ligand association with the target. Thus, comprehensive molecular-level insights gained from this work can be useful to identify or design potent small molecule inhibitors against SARS-CoV-2 Mpro.

Modeling the Binding Mechanism of Remdesivir, Favilavir, and Ribavirin to SARS-CoV-2 RNA-Dependent RNA Polymerase

Recent efforts to repurpose drugs to combat COVID-19 have identified Remdesivir as a candidate. It acts on the RNA-dependent, RNA polymerase (RdRp) of the SARS-CoV-2 virus, a protein complex responsible for mediating replication of the virus's genome. However, its exact action mechanism, and that of other nucleotide analogue inhibitors, is not known. In this study, we examine at the molecular level the interaction of this drug and that of similar nucleotide analogue inhibitors, ribavirin and favilavir, by relying on atomistic molecular simulations and advanced sampling. By analyzing the binding free energies of these different drugs, it is found that all of them bind strongly at the active site. Surprisingly, however, ribavirin and favilavir do not bind the nucleotide on the complementary strand as effectively and seem to act by a different mechanism than remdesivir. Remdesivir exhibits similar binding interactions to the natural base adenine. Moreover, by analyzing remdesivir at downstream positions of the RNA, we also find that, consistent with a "delayed" termination mechanism, additional nucleotides can be incorporated after remdesivir is added, and its highly polar 1'-cyano group induces a set of conformational changes that can affect the normal RdRp complex function. By analyzing the fluctuations of residues that are altered by remdesivir binding, and comparing them to those induced by lethal point mutations, we find a possible secondary mechanism in which remdesivir destabilizes the protein complex and its interactions with the RNA strands.

SARS-CoV-2 protein drug targets landscape: a potential pharmacological insight view for the new drug development

Introduction: Protein drug targets play a significant choice in different stages of the drug discovery process. There is an urgent need to understand the drug discovery approaches and protein drug targets (PDT) of SARS-CoV-2, with structural insights for the development of SARS-CoV-2 drugs through targeted therapeutic approach.Areas covered: We have described the protein as a drug target class and also discussed various drug discovery approaches for SARS-CoV-2 involving the protein drug targets such as drug repurposing study, designing of viral entry inhibitors, viral replication inhibitors, and different enzymes of the virus. We have performed comprehensive literature search from the popular databases such as PubMed Google scholar, Web of Science, and Scopus. Finally, we have illustrated the structural landscape of different significant viral proteins (3 CLpro or Mpro, PLpro, RdRp, helicase, S protein) and host proteins as drug targets (cathepsin L, furin, TMPRSS2, ACE2).Expert opinion: The structural landscape of PDT with their binding pockets, and significant residues involved in binding has been discussed further to better understand the PDT and the structure-based drug discovery for SARS-CoV-2. This attempt will increase more therapeutic options, and combination therapies with a multi-target strategy.

RNA dependent RNA polymerase (RdRp) as a drug target for SARS-CoV2

RNA-dependent RNA polymerase (RdRp), also called nsp12, is considered a promising but challenging drug target for inhibiting replication and hence, the growth of various RNA-viruses. In this report, a computational study is performed to offer insights on the binding of Remdesivir and Galidesivir with SARS-CoV2 RdRp with natural substrate, ATP, as the control. It was observed that Remdesivir and Galidesivir exhibited similar binding energies for their best docked poses, −6.6 kcal/mole and −6.2 kcal/mole, respectively. ATP also displayed comparative and strong binding free energy of −6.3 kcal/mole in the catalytic site of RdRp. However, their binding locations within the active site are distinct. Further, the interaction of catalytic site residues (Asp760, Asp761, and Asp618) with Remdesivir and Galidesivir is comprehensively examined. Conformational changes of RdRp and bound molecules are demonstrated using 100 ns explicit solvent simulation of the protein-ligand complex. Simulation suggests that Galidesivir binds at the non-catalytic location and its binding strength is relatively weaker than ATP and Remdesivir. Remdesivir also binds at the catalytic site and showed high potency to inhibit the function of RdRp. Binding of co-factor units nsp7 and nsp8 with RdRp (nsp12) complexed with Remdesivir and Galidesivir was also examined. MMPBSA binding energy for all three complexes has been computed across the 100 ns simulation trajectory. Overall, this study suggests, Remdesivir has anti-RdRp activity via binding at a catalytic site. In contrast, Galidesivir may not have direct anti-RdRp activity but it can induce a conformational change in the RNA polymerase.