Direct-acting antiviral agents (DAAs) are known to interfere with various intracellular stages of the hepatitis C virus (HCV) life cycle and have demonstrated efficacy in treating HCV infection. However, DAA monotherapy can lead to drug resistance due to mutations. This paper explores the impact of DAA therapy on HCV dynamics using a multiscale age-structured partial differential equation (PDE) model that incorporates intracellular viral RNA replication within infected cells and two strains of viruses representing a drug-sensitive strain and a drug-resistant mutant variant, respectively. We derived an equivalent ordinary differential equation (ODE) model from the PDE model to simplify mathematical analysis and numerical simulations. We studied the dynamics of the two virus strains before treatment and investigated the impact of mutations on the evolution kinetics of drug-sensitive and drug-resistant viruses, as well as the competition between the two strains during treatment. We also explored the role of DAAs in blocking HCV RNA replication and releasing new virus particles from cells. During treatment, mutations do not significantly influence the dynamics of various virus strains; however, they can generate low-level HCV that may be completely inhibited due to their poor fitness. The fitness of the mutant strain compared to the drug-sensitive strain determines which strain dominates the virus population. We also investigated the prevalence and drug resistance evolution of HCV variants during DAA treatment.
more »
« less
Binding of GS-461203 and Its Halogen Derivatives to HCV Genotype 2a RNA Polymerase Drug Resistance Mutants
Hepatitis C Virus (HCV) is reported to develop GS-461203 resistance because of multiple mutations within the RNA-dependent RNA Polymerase (RdRp) of HCV. The lack of a high-resolution structure of these RdRp mutants in complex with GS-461203 hinders efforts to understand the drug resistance. Here we decipher the binding differences of GS-461203 in the wild type and mutated systems T179A or M289L of HCV RdRp Genotype 2a using homology modeling, molecular docking, and molecular dynamics simulation. Key residues responsible for GS-461203 binding were identified to be Arg48, Arg158, Asp318, Asp319, and Asp220, and that mutations T179A or M289L have caused conformational changes of GS-461203 in the RdRp active site. The affinities of GS-461203 were reduced in T179A system, but it became slightly stronger in the M289L system. Furthermore, we designed two new analogues of GS-461203 which encouragingly induced more stable interactions than GS-461203, and thus resulted in much better binding energies. This present study reveals how a single mutation, T179A or M289L, will modulate GS-461203 binding in HCV RdRp Genotype 2a, while introducing two novel analogues to overcome the drug resistance which may be good candidate for further experimental verification.
more »
« less
- Award ID(s):
- 1904797
- PAR ID:
- 10324150
- Date Published:
- Journal Name:
- Scientia Pharmaceutica
- Volume:
- 90
- Issue:
- 2
- ISSN:
- 2218-0532
- Page Range / eLocation ID:
- 26
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract 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.more » « less
-
null (Ed.)The novel coronavirus disease 19 (Covid-19) which is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been a pandemic across the world, which necessitate the need for the antiviral drug discovery. One of the potential protein targets for coronavirus treatment is RNA-dependent RNA polymerase. It is the key enzyme in the viral replication machinery, and it does not exist in human beings, therefore its targeting has been considered as a strategic approach. Here we describe the identification of potential hits from Indonesian Herbal and ZINC databases. The pharmacophore modeling was employed followed by molecular docking and dynamics simulation for 40 ns. 151 and 14480 hit molecules were retrieved from Indonesian herbal and ZINC databases, respectively. Three hits that were selected based on the structural analysis were stable during 40 ns, while binding energy prediction further implied that ZINC1529045114, ZINC169730811, and 9-Ribosyl-trans-zeatin had tighter binding affinities compared to Remdesivir. The ZINC169730811 had the strongest affinity toward RdRp compared to the other two hits including Remdesivir and its binding was corroborated by electrostatic, van der Waals, and nonpolar contribution for solvation energies. The present study offers three hits showing tighter binding to RdRp based on MM-PBSA binding energy prediction for further experimental verification.more » « less
-
Abstract Understanding how mutations render a drug ineffective is a problem of immense relevance. Often the mechanism through which mutations cause drug resistance can be explained purely through thermodynamics. However, the more perplexing situation is when two proteins have the same drug binding affinities but different residence times. In this work, we demonstrate how all‐atom molecular dynamics simulations using recent developments grounded in statistical mechanics can provide a detailed mechanistic rationale for such variances. We discover dissociation mechanisms for the anti‐cancer drug Imatinib (Gleevec) against wild‐type and the N368S mutant of Abl kinase. We show how this point mutation triggers far‐reaching changes in the protein's flexibility and leads to a different, much faster, drug dissociation pathway. We believe that this work marks an efficient and scalable approach to obtain mechanistic insight into resistance mutations in biomolecular receptors that are hard to explain using a structural perspective.more » « less
-
null (Ed.)Background: Estrogen Receptors (ER) are members of the nuclear intracellular receptorsfamily. ER once activated by estrogen, it binds to DNA via translocating into the nucleus and regulatesthe activity of various genes. Withaferin A (WA) - an active compound of a medicinal plant Withaniasomnifera was reported to be a very effective anti-cancer agent and some of the recent studies hasdemonstrated that WA is capable of arresting the development of breast cancer via targeting estrogenreceptor. Objective: The present study is aimed at understanding the molecular level interactions of ER and Tamoxifenin comparison to Withaferin A using In-silico approaches with emphasis on Withaferin Abinding capability with ER in presence of point mutations which are causing de novo drug resistance toexisting drugs like Tamoxifen. Methods: Molecular modeling and docking studies were performed for the Tamoxifen and WithaferinA with the Estrogen receptor. Molecular docking simulations of estrogen receptor in complex withTamoxifen and Withaferin A were also performed. Results: Amino acid residues, Glu353, Arg394 and Leu387 was observed as crucial for binding andstabilizing the protein-ligand complex in case of Tamoxifen and Withaferin-A. The potential ofWithaferin A to overcome the drug resistance caused by the mutations in estrogen receptor to the existingdrugs such as Tamoxifen was demonstrated. Conclusion: In-silico analysis has elucidated the binding mode and molecular level interactions whichare expected to be of great help in further optimizing Withaferin A or design / discovery of futurebreast cancer inhibitors targeting estrogen receptor.more » « less
An official website of the United States government

