Computational study on the drug resistance mechanism of HCV NS5B RNA-dependent RNA polymerase mutants V494I, V494A, M426A, and M423T to Filibuvir

Small molecule NS5B RdRp non-nucleoside inhibitors for the treatment of HCV infection: A medicinal chemistry perspective

Hepatitis C virus (HCV) infection has become a global health problem with enormous risks. Nonstructural protein 5B (NS5B) RNA-dependent RNA polymerase (RdRp) is a component of HCV, which can promote the formation of the viral RNA replication complex and is also an essential part of the replication complex itself. It plays a vital role in the synthesis of the positive and negative strands of HCV RNA. Therefore, the development of small-molecule inhibitors targeting NS5B RdRp is of great value for treating HCV infection-related diseases. Compared with NS5B RdRp nucleoside inhibitors, non-nucleoside inhibitors have more flexible structures, simpler mechanisms of action, and more predictable efficacy and safety of drugs in humans. Technological advances over the past decade have led to remarkable achievements in developing NS5B RdRp inhibitors. This review will summarize the non-nucleoside inhibitors targeting NS5B RdRp developed in the past decade and describe their structure optimization process and structure-activity relationship.

Insight into the drug resistance mechanisms of GS-9669 caused by mutations of HCV NS5B polymerase via molecular simulation

GS-9669 is a non-nucleos(t)ide inhibitor (NNI) binding to the thumb site II of the Hepatitis C virus (HCV) NS5B polymerase and has advanced into phase II trials. To clarify the drug resistance mechanisms of GS-9669 caused by M423T/I/V, L419M, R422K, and I482L mutations of NS5B polymerase (GT1b) and the receptor-ligand interactions during the binding process, a series of molecular simulation methods including molecular dynamics (MD) simulations and adaptive steered molecular dynamics (ASMD) simulations were performed for the wild-type (WT) and six mutant NS5B/GS-9669 complexes. The calculated results indicate that the binding free energies of the mutant systems are less negative than that of the WT system, indicating that these mutations will indeed cause NS5B to produce different degrees of resistance to GS-9669. The mutation-induced drug resistances are mainly caused by the loss of binding affinities of Leu419 and Trp528 with GS-9669 or the formation of multiple solvent bridges. Moreover, the ASMD calculations show that GS-9669 binds to the thumb II sites of the seven NS5B polymerases in distinct pathways without any obvious energy barriers. Although the recognition methods and binding pathways are distinct, the binding processes of GS-9669 with the WT and mutant NS5B polymerases are basically controlled thermodynamically. This study clearly reveals the resistance mechanisms of GS-9669 caused by M423T/I/V, L419M, R422K, and I482L mutations of HCV NS5B polymerase and provides some valuable clues for further optimization and design of novel NS5B inhibitors.

Computational insights into the molecular mechanisms of differentiated allosteric modulation at the mu opioid receptor by structurally similar bitopic modulators

Targeting the mu opioid receptor (MOR) by applying orthosteric ligands is the most frequently employed method to treat opioid use disorder (OUD). Unfortunately, most of MOR orthosteric ligands produce severe side effects, mainly due to their low selectivity over other opioid receptors. In contrast, some G protein-coupled receptor allosteric modulators have been reported to exhibit high subtype selectivity and can effectively modulate the potency and/or efficacy of orthosteric ligands. Recently, NAQ and its analog NCQ were identified as novel MOR bitopic modulators. Interestingly, NAQ and NCQ were similar in structure but exhibited different efficacy profiles to the MOR. NAQ exhibited an antagonism activity to the MOR while NCQ showed a partial agonism activity to the MOR. In the present study, molecular modeling methods were applied to explore the putative molecular mechanisms of their different functional profiles to the MOR. When NAQ binding with the inactive MOR, the 'address' portion of NAQ interacted with the MOR allosteric site but showed no significant allosteric modulation of the efficacy of the 'message' portion of NAQ. However, when NCQ binding with the inactive and active MOR, the 'address' portion of NCQ seemed to be able to positively modulate the efficacy of the 'message' portion of NCQ at varying levels. Evidentially, the substituents at the 1'- and 4'-positions of the isoquinoline ring of NCQ seemed to play a critical role in the modulatory function of the 'address' portion of NCQ. These findings will be invaluable to develop our next generation of MOR bitopic modulators with high affinity and subtype selectivity to potentially treat OUD.

Molecular mechanism of resveratrol inhibition of Zika virus NS3 helicase: behind the scenes

Aim: Zika virus (ZIKV) still poses a health risk to women and their babies without US FDA-approved vaccines or treatments. Experimentation has proved resveratrol inhibition of ZIKV NS3 helicase without specifying the molecular events during inhibition. Materials & methods: Herein, we leaped forward to study the molecular dynamics of the bound and unbound enzyme, identifying precise binding residues and interactions, and the enzyme's adaptation to support binding, since loop dynamics affect viral RNA replication. Results: Resveratrol stabilizes the P-loop and causes the RNA-binding loop to block the RNA-binding pocket for 200 ns, which is concurrent with experimental evidence that resveratrol binding significantly reduces ATP hydrolysis activity. Conclusion: This study illuminates the structural dynamics of ZIKV helicase and druglikeness of resveratrol, which will advance anti-ZIKV drug development.

Probing the "fingers" domain binding pocket of Hepatitis C virus NS5B RdRp and D559G resistance mutation via molecular docking, molecular dynamics simulation and binding free energy calculations

The NS5B RdRp polymerase is a prominent enzyme for the replication of Hepatitis C virus (HCV). During the HCV replication, the template RNA binding takes place in the "fingers" sub-domain of NS5B. The "fingers" domain is a new emerging allosteric site for the HCV drug development. The inhibitors of the "fingers" sub-domain adopt a new antiviral mechanism called RNA intervention. The details of essential amino acid residues, binding mode of the ligand, and the active site intermolecular interactions of RNA intervention reflect that this mechanism is ambiguous in the experimental study. To elucidate these details, we performed molecular docking analysis of the fingers domain inhibitor quercetagetin (QGN) with NS5B polymerase. The detailed analysis of QGN-NS5B intermolecular interactions was carried out and found that QGN interacts with the binding pocket amino acid residues Ala97, Ala140, Ile160, Phe162, Gly283, Gly557, and Asp559; and also forms π⋯π stacking interaction with Phe162 and hydrogen bonding interaction with Gly283. These are found to be the essential interactions for the RNA intervention mechanism. Among the strong hydrogen bonding interactions, the QGN⋯Ala140 is a newly identified important hydrogen bonding interaction by the present work and this interaction was not resolved by the previously reported crystal structure. Since D559G mutation at the fingers domain was reported for reducing the inhibition percentage of QGN to sevenfold, we carried out molecular dynamics (MD) simulation for wild and D559G mutated complexes to study the stability of protein conformation and intermolecular interactions. At the end of 50 ns MD simulation, the π⋯π stacking interaction of Phe162 with QGN found in the wild-type complex is altered into T-shaped π stacking interaction, which reduces the inhibition strength. The origin of the D559G resistance mutation was studied using combined MD simulation, binding free energy calculations and principal component analysis. The results were compared with the wild-type complex. The mutation D559G reduces the binding affinity of the QGN molecule to the fingers domain. The free energy decomposition analysis of each residue of wild-type and mutated complexes revealed that the loss of non-polar energy contribution is the origin of the resistance. Communicated by Ramaswamy H. Sarma.

Exploring the binding mechanisms of diaminopimelic acid analogs to meso-diaminopimelate dehydrogenase by molecular modeling

Meso-Diaminopimelic acid (meso-2,6-diamino-heptanedioic acid, DAP) is an important component of the cell wall of many bacteria. Meso-diaminopimelate dehydrogenase (m-Ddh) is a critical enzyme in the process of converting tetrahydrodipicolinate to DAP. Here, we are proposing that DAP analogs targeting m-Ddh may be considered as potential antibiotics. Four DAP analogs without significant structural change from DAP have been obtained and their inhibitory potencies against m-Ddh from the P. gingivalis strain W83 show significant differences from that of DAP. However, their inhibitory mechanisms as for how simple structural change influences the inhibitory potency remain unknown. Therefore, we employed molecular modeling methods to obtain insight into the inhibitory mechanisms of DAP and analogs with m-Ddh. The predicted binding mode of DAP was highly consistent with the experimental structural data and disclosed the important roles played by the binding pocket residues. According to our predictions, the isoxazoline ring of compounds 1 and 2 and the double bonds in compounds 3 and 4 had distinct influences on these compounds' binding to m-Ddh. This enriched understanding of the inhibitory mechanisms of DAP and these four analogs to m-Ddh has provided new and relevant information for future rational development of potent inhibitors targeting m-Ddh.

Computational study on the drug resistance mechanism of HCV NS3 protease to BMS-605339

NS3 protease plays a vital role in the replication of the hepatitis C virus (HCV). BMS-605339 is a novel linear tetra-peptide α-ketoamide inhibitor of NS3 protease and shows specificity for HCV NS3 protease genotype 1a and genotype 1b. Mutation at the key site 168 of the HCV NS3 protease can induce resistance to BMS-605339, which greatly affects the antiviral therapy efficacy to hepatitis C. In the present study, we employed molecular dynamics simulations, free energy calculations, and free energy decomposition to explore the drug resistance mechanism of BMS-605339 due to the three representative mutations D168C/Y/V. The free energy decomposition analysis indicates that the decrease in the binding affinity is mainly attributed to the decrease in both van der Waals and electrostatic interactions. After detailed analysis of our calculated results, we observed that the break of the salt bridge between residues 155 and 168 caused by the mutations D168C/Y/V is the original reason for the decrease in the binding ability between BMS-605339 and the mutant NS3 proteases. The obtained results will reveal the drug resistance mechanism between BMS-605339 and the mutant NS3 proteases, and provide valuable clue for designing novel and more potent drugs to HCV NS3 protease.

Molecular modeling study on the drug resistance mechanism of NS5B polymerase to TMC647055

NS5B polymerase plays an important role in viral replication machinery. TMC647055 (TMC) is a novel and potent non-nucleoside inhibitor of the HCV NS5B polymerase. However, mutations that result in drug resistance to TMC have been reported. In this study, we used molecular dynamics (MD) simulations, binding free energy calculations, and free energy decomposition to investigate the drug resistance mechanism of HCV to TMC resulting from L392I, P495T, P495S, and P495L mutations in NS5B polymerase. From the calculated results we determined that the decrease in the binding affinity between TMC and NS5B(L392I) polymerase is mainly caused by the extra methyl group at the CB atom of Ile. The polarity of the side-chain of residue 495 has no distinct influence on residue 495 binding with TMC, whereas the smaller size of the side-chain of residue 495 causes a substantial decrease in the van der Walls interaction between TMC and residue 495. Moreover, the longer length of the side-chain of residue 495 has a significant effect on the electrostatic interaction between TMC and Arg-503. Finally, we performed the same calculations and detailed analysis on other 3 mutations (L392V, P495V, and P495I). The results further confirmed our conclusions. The computational results not only reveal the drug resistance mechanism between TMC647055 and NS5B polymerase, but also provide valuable information for the rational design of more potent non-nucleoside inhibitors targeting HCV NS5B polymerase.