Volume 17 Supplement 8
Modeling the functional state of the reverse transcriptase of hepatitis B virus and its application to probing drug-protein interaction
© The Author(s). 2016
Published: 31 August 2016
Herein, the predicted atomic structures of five representative sequence variants of the reverse transcriptase protein (RT) of hepatitis B virus (HBV), sampled from patients with rapid or slow response to tenofovir disoproxil fumarate (TDF) treatment, have been examined to identify structural variations between them in order to assess structural and functional properties of HBV-RT variants associated with the differential responses to TDF treatment.
We utilized a hybrid computational approach to model the atomistic structures of HBV-RT/DNA-RNA/dATP and HBV-RT/DNA-RNA/TFV-DP (tenofovir diphosphate) complexes with the native hybrid DNA-RNA substrate in place. Multi-nanosecond molecular dynamics (MD) simulations of HBV-RT/DNA-RNA/dATP complexes revealed strong coupling of the natural nucleotide substrate, dATP, to the active site of the RT, and the differential involvement of the two putative magnesium cations (Mg2+) at the active site, whereby one Mg2+ directly bridges the interaction between dATP and HBV-RT and the other serves as a coordinator to maintain an optimal configuration of the active site. Solvated interaction energy (SIE) calculated in MD simulations of HBV-RT/DNA-RNA/TFV-DP complexes indicate no differential binding affinity between TFV-DP and HBV-RT variants identified in patients with slow or rapid response to TDF treatment.
The predicted atomic structures accurately represent functional states of HBV-RT. The equivalent interaction between TFV-DP and each examined HBV-RT variants suggests that binding affinity of TFV-DP to HBV-RT is not associated with delayed viral clearance.
KeywordsHybrid structure modeling Molecular dynamics (MD) Solvated interaction energy (SIE) Hepatitis B Tenofovir Drug resistance Reverse transcriptase
During the life cycle of HBV, RT utilizes a single-stranded viral genomic RNA as a template to synthesize a hybrid RNA-DNA duplex, and then converts it to double-stranded DNA (ds-DNA). As this step is critical for the viral genome replication, it makes RT an attractive target for antiviral treatment. Thus, both nucleoside and nucleotide RT inhibitors (NRTIs) and non-nucleoside RT inhibitors (NNRTI) are extensively used as antiretroviral agents against HBV and HIV infection [1, 2].
TDF is a non-selective RT inhibitor with demonstrated in vitro activity against both wild-type and mutant strains of HIV [4–6] and is also effective for treating chronic hepatitis B (CHB) [7–10]. While resistance mutations within the HBV-RT have been identified for NRTI monotherapies, such as adefovir, entecavir, lamivudine, telbivudine and clevudine , there is no evidence of TDF resistance . However, differential therapy responses can be observed among patients and have been previously defined: a delayed viral clearance response in slow responders (SR), or a rapid decline in viral load immediately following treatment initiation in rapid responders (RR) [12, 13].
TFV resistance in the setting of HIV has been extensively studied at the molecular level [14–16]. These studies have shown that the molecular mechanism of TFV resistance is complex and involves reverse transcriptase, RNase and exonuclease activities of HIV-RT. In order to remain functional, resistant mutants of HIV-RT need to maintain their capability of incorporating nucleotide at a sufficient rate, as well as to discriminate against NRTIs or improve the excision rate of incorporated NRTIs .
As the HBV-RT is homologous to HIV-RT , it is reasonable to propose that a similar mechanism could be involved in the differential treatment response to TDF. Whether and/or how the molecular level interaction between TFV-DP and HBV-RT has a part in the aforementioned differential patient responses upon TDF treatment is unclear due to the lack of direct experimental evidence. There is no solved structure of the HBV-RT in any form.
In this study, we set out to address the molecular level interaction between TFV-DP and HBV-RT and the relationship to patient response using a hybrid computational approach. We utilized two X-ray crystal structures of HIV-RT to model the atomic structure of HBV-RT in its functional state with the aim of delineating the detailed drug-protein interactions between TFV-DP and the active sites of genetically distinct variants of HBV-RT.
Results and discussion
Atomic structure models of the HBV-RT in its functional state
While substantial efforts have been devoted to solving the crystal structures of HIV-RT, with or without ligands (e.g. substrates, small inhibitory molecules), no experimentally derived HBV-RT structures have become available. HIV-RT and HBV-RT are recognized as homologous , and possess similar multiple catalytic activities. A computational approach is often applied to structural implications in drug-resistance studies on HBV-RT using HIV-RT as structural template for homology modeling of HBV-RT [17, 18]. Still, the low sequence identity and modeling HBV RT without considering the particular functional state imposes technical difficulties in obtaining structural models of high confidence.
In models of NRTIs bound to HBV-RT, it is reasonable to hypothesize that an effective drug binds to the active site as if it were a native substrate in position to be further covalently incorporated, similar to the HIV-RT/DNA-DNA/TFV-DP complex . Even though we are not to rule out the possibility of a NRTI being able to bind elsewhere than the active site on HBV-RT or HIV-RT, to date there is no evidence for such type of interaction. Moreover, a putative strong binding of the drug to the active site would seem to preclude the likelihood of additional binding sites on the enzyme. To use constraints derived from experimental results to guide the refinement, we first modeled the HBV-RT structure itself using ITASSER [20, 21] with an additional template, the X-ray crystal structure of HIV-RT/DNA-RNA/dATP , and then included the native substrates (the DNA-RNA duplex and dATP) from this template complex in the subsequent refinement (MD simulation).
The C-scores  (data not shown) of our HBV-RT models generated by ITASSER for the variants of interest indicated high confidence; however, we noticed that the exact folding patterns in the fingers domain differ (Additional file 1: Figure S1). Therefore, we chose the HBV-RT4 model (HBV-RT4, Additional file 1: Figure S1C) which has the same folding pattern as the HIV-RT template (Additional file 1: Figure S1F) in order to proceed with the refinement.
The structure of the N-terminal 40 amino acid region was of relatively low accuracy (Additional file 1: Figure S1A-E). Thus, we re-modeled the N-terminal of HBV-RT4 separately using ITASSER, and assembled it with the rest of the protein using ab initio assembly . The N-terminus in the resultant assembled structure of HBV-RT4 is of much more ordered secondary structural elements than it was in the initial model from ITASSER (Additional file 2: Figure S2).
The assembled structure was then combined with the duplex DNA-RNA and dATP from the template HIV-RT/DNA-RNA/dATP crystal strcuture, refined by ~120 ns MD simulation in explicit solvent to obtain the final model of HBV-RT4/DNA-RNA/dATP.
In preparation for incorporation, the dATP situates itself in the active site with good base ring stacking to the penultimate dG, and by making two persistent hydrogen bonds (H-bond) with the uracil from the opposing template RNA strand (Fig. 2a, Additional file 2: Figure S2A); these features persist throughout the un-constrained MD (20 -120 ns). The triphosphate group of the dATP is stably coupled with residue L146 and D205 of HBV-RT by cooperatively chelating one Mg2+ (Fig. 2b). This particular configuration differs from the HIV-RT active site mainly by the absence of a third aspartate residue in the active site; in HIV-RT this third aspartate (D110) participates in chelating the two Mg2+ in the active site of the HIV-RT (Fig. 2c). Without this aspartate, the Mg2+ ions were stabilized over the equilibration stage of the MD by incorporating extra water molecules, and this constellation remained stable throughout in the production run (see next section).
Different roles of Mg2+ ions in the active site
Approximated binding free energies with standard errors of dATP or TNV to HBV-RT/DNA-RNA using Solvated Interaction Energy method (SIE)
dATP, one Mg2+ (Kcal mol-1)
dATP, two Mg2+ (Kcal mol-1)
TFV-DP, two Mg2+ (Kcal mol-1)
−10.43 ± 0.02
−9.93 ± 0.02
−8.91 ± 0.01
−11.51 ± 0.02
−11.21 ± 0.03
−9.22 ± 0.01
−10.59 ± 0.01
−9.69 ± 0.02
−9.56 ± 0.01
−12.51 ± 0.01
−9.50 ± 0.02
−9.05 ± 0.02
−10.51 ± 0.02
−10.59 ± 0.01
−9.51 ± 0.01
Binding free energy of TFV-DP to HBV-RT
Drug resistance often arises from structural alterations consequent to mutations in the target protein, which somehow disrupts or interferes with drug binding. To understand the underlying mechanism of the differential response to TDF treatment, we start by asking whether the TFV-DP binds differently to the RT variants from patients with slow or rapid response to TDF. In building a structural model of HBV-RT4/DNA-RNA/TFV-DP, we started with the HBV-RT4 model generated by ITASSER and AIDA (Additional file 2: Figure S2), then inserted the duplex DNA-RNA taken originally from the HIV-RT/DNA-RNA/dATP template structure , and finally incorporated the TFV-DP from the HIV-RT/DNA-DNA/TFV-DP structure (PDB code: 1 T05) . Models for the other four variants were built analogously from the homology models we obtained using the HBV-RT4 as template, into which the same DNA-RNA duplex and TFV-DP were incorporated. We then pursued MD simulations to refine all the HBV-RT/DNA-RNA/TFV-DP structures. The centroid of the major cluster from the second half of each MD trajectory was then chosen as representative model.
We were then able to compute the binding energies of the dATP and TFV-DP to the HBV-RT/DNA-RNA complex with extensive sampling of the complex structures through MD simulations using the solvated interaction energies (SIE) method [24, 25]. The favorable and consistent binding energies of dATP to the variants observed in both Mg2+ systems (Table 1) indicate strong binding of the natural substrate to the enzyme, hence proper configuration of the active sites for all of our HBV-RT variant models. As for the TFV-DP ligand, the SIE results also show favorable binding to all HBV-RT variants (Table 1), with no substantial differences between them (<1 Kcal mol-1). These findings suggest that RT enzymatic activity and RT’s affinity to bind TFV-DP are not altered by the amino-acid heterogeneity observed among the HBV-RT variants examined herein, which, taken together with the observed stability of ligand-protein interactions and consistency in the structural configuration of the active site in all HBV-RT/DNA-RNA/TFV-DP complexes, suggest that factors other than a decreased binding affinity of HBV-RT to TFV-DP might be influencing therapy response in patients.
Interestingly, evidence of there being no substantial changes in drug-binding affinity in our HBV-RT/DNA-RNA/TFV-DP complexes (Table 1) is in agreement with our previous study . A potential prediction model of response to TDF revealed that genetic diversity outside the RT domain of polymerase contributed robust predictive features strongly associated to SR/RR characteristics of patients, whereas in RT none were found . Moreover, a switch in high-frequency intra-host HBV variants was observed in patient 4, which did not involve changes in RT. Amino-acid substitution S202P (GenBank reference AY AY233278) in the spacer domain of polymerase presented after 4 weeks of TDF treatment, suggesting a possible association between response to TDF and amino-acid heterogeneity in spacer, a suspected intrinsically disordered protein (IDP) . In fact, there is compelling evidence of the involvement of IDPs in evolution of drug resistance in cancer cells [28, 29] and HBV . Nevertheless, an association between amino-acid substitutions in HBV variants RT1–RT5 and differential response to TDF treatment cannot be dismissed solely based on findings presented herein. Drug-binding affinity may not be the only mechanism involved in response rates or resistance to TDF. Amino-acid substitutions observed in HBV-RT could be affecting some other structure-activity-resistance relationships. In addition, intramolecular conformational changes in RT and intermolecular dynamics of the RT complex provide an alternative mechanism for drug resistance to NNRTIs .
While many drug resistant mutants of HBV-RT have been identified, very little is known about the molecular mechanism driving resistance development. Compared to the massive structural and biochemical studies than have been conducted on HIV-RT drug resistant mutants, HBV-RT has not drawn much attention among structural biologists. However, over 350 million people have developed CHB and it remains a serious, global health problem . Therefore, it is important to unravel the drug resistance associated with anti-viral treatment to HBV infections.
To probe the structural basis for HBV-RT drug resistance, researchers rely on structural modeling approaches using crystal structures of HIV-RT as a starting point [17, 18, 31–33]. As protein structures are dynamic and strongly correlated with its functional states, we devised a novel hybrid approach in which we utilized the most relevant X-ray crystal structures (in terms of its functional state) to model the initial structures of HBV-RT and then extensively sampled the conformational space with experimentally derived constraints (crystal contacts between substrate and protein) to select the final representative models (Additional file 3: Figure S3). Without direct experimentally derived structure data on HBV-RT to test the models, this rigorous approach presumably generates structure models of the highest confidence and relevance to the particular biological function of interest. However, as this approach comprises only physics-based methods, it is advisable to apply other methods to evaluate the resulting model from a different perspective. In this study, we have applied a statistical potential based method (ProSA-Web)  to evaluate our model. The resulting Z-score for our HBV-RT4 model confirms good quality in terms of the naturalness defined by sets of X-ray and NMR structures (Additional file 4: Figure S4A), and the local quality plot indicates the majority of the HBV-RT4 model is of high quality measured by a knowledge-based energetic term (Additional file 4: Figure S4B).
Using these models to probe drug-protein interaction, we focused on the TFV-DP binding. The present findings indicate that, despite sequence changes in or around the NRTI binding site, favorable drug-protein interactions persisted. The dynamics of how key molecular groups interact, probed in our studies, reveal new insights into the mechanisms, through which HBV-RT couples to its natural substrate. However, this apparent drug-protein interaction stability leaves us without a good mechanistic rationalization as to why or how HBV-RT amino-acid polymorphisms may contribute to reduced rates of viral clearance responses to TDF treatment, as is often observed in SR patients.
The modeling approach presented in this study can be extended to other drugs and serve as a paradigm to ascertain drug binding and estimate binding free energies. It is important to note that the starting conformation (the binding site and the particular interaction pattern between the drug and protein) would have much effect on the results. Therefore, it is essential to start modeling using all available information for the accurate and reliable prediction of the binding site structure.
Patients and HBV quasispecies (QS)
Five patients from Study GS-US-203-0101, a phase 2 study evaluating TDF and FTC + TDF in treatment naive patients with HBV DNA >1.7 x 107 IU/mL and normal ALT levels for 192 weeks, were selected . 3 patients had a slow response (SR), never achieving HBV DNA <69 IU/mL through 192 weeks of treatment, while 2 patients had a rapid response (RR), achieving HBV DNA <69 IU/mL by Week 96. Patients were matched by HBV DNA, ALT, and HBeAg status at baseline. All patients received TDF monotherapy. Whole-genome sequences of HBV QS were obtained using end-point limiting-dilution real-time PCR coupled with sequencing. Baseline (BL), Week 4 (W4), and Week 40 (W40) time points were evaluated for each patient. Median joining networks (MJN) were performed using Network 4.0 to analyze HBV QS genetic diversity . The high-frequency variants of intra-host HBV QS from each patient were selected to predict the atomic structure of the RT protein domain in HBV polymerase (protein positions 349 – 693, based on GenBank reference sequence AY233278): from SR patients, variants RT1 and RT2 of HBV genotype C (GT C) and RT4 (GT B), and from RR patients, variants RT3 (GT C) and RT5 (GT B).
Structural modeling of HBV-RT variants
Due to the nonexistence of experimental 3D-structures of HBV-RT, a hybrid approach to protein 3D-strucutre prediction was implemented to generate accurate, atomistic structures of the NRTI binding site of HBV-RT. Three major steps were involved in the modeling. First, all the representative HBV-RT sequences were submitted to ITASSER web-server [20, 21], including specification on an homologous HIV-RT structure complex (PDB accession code: 4PQU) as an additional template . The model of the highest C-score  were then chosen to proceed. The N-terminal 1-40 amino-acid of HBV-RT were modeled separately by ITASSER, due to a relatively low local accuracy for the initial model of the entire length. The N-terminal structure was then assembled to the preliminary model by using Ab Initio Domain Assembly Server (AIDA) .
Second, all the HBV-RT variant models were then pooled for inspection to identify the one(s) with the folding pattern most consistent with that seen in the designated HIV-RT structure.
Last, variant HBV-RT4 was chosen for further refinement using multi-nanosecond molecular dynamics (MD) simulation with the native hybrid DNA-RNA duplex and with dATP incorporated into the active site, according to their original orientation in the HIV-RT/DNA-RNA/dATP crystal structure . Pairwise rmsd clustering analysis of the second half of the MD trajectory was then performed to select the centroid of the dominant cluster as the final model of HBV-RT4/DNA-RNA/dATP. The structures of the other variants were then obtained by homology modeling using Modeller 9.13 [37, 38] with the HBV-RT4 as template, followed by same MD simulation procedures, also with the native hybrid DNA-RNA duplexes and the dATPs superimposed to the active sites.
In the case of modeling the HBV-RT4/DNA-RNA/TFV-DP variants, we replaced the dATP with TFV-DP in the initial conformation of the HBV-RT4/DNA-RNA/dATP MD simulation, maintaining the crystal contacts between TFV-DP and corresponding residues that was identified from the HIV-RT4/DNA-DNA/TFV-DP crystal structure (PDB code: 1 T05) .
Molecular dynamic (MD) simulation
To conduct MD simulations, hydrogens were first added by the tleap module of AMBER 14 , and the ionizable side chains of the proteins were assigned to their ionization states at pH 7.0 using the WHATIF webserver . Each system was then solvated with TIP3P water molecules  leaving a minimum distance of 10.0 Å from the protein surface to the edge of the simulation box. Counter-ions were added to neutralize the net charge and reach 100 mM NaCl concentration to mimic physiological conditions. The systems were then minimized and equilibrated -- with extra bonds constraints (2.5 kcal mol−1 Å−2) on the Mg2+ ions and on the ionic bonded oxygen atoms -- for 10–20 ns before free production runs in the isothermal isobaric ensemble (1 atm, 300 K). The short-range non-bonded interactions were evaluated by employing a cut off of 10 Å with a switching function starting at 8.5 Å. The smooth particle mesh Ewald (SPME) algorithm were applied in computing the ong-range electrostatic interactions . All the bonds formed by hydrogen and heavy atoms were fixed to eliminate the most frequent oscillatory motions. The r-RESPA multiple time step method  was applied with a 2 fs time step for bonded, 2 fs for short-range non-bonded interactions and 4 fs for long-range electrostatic interactions. All simulations were performed using the NAMD 2.9 code [44, 45] with the AMBER Parm99 parameter set  containing the force field for nucleic acids and proteins. The parameterization of TFV-DP was conducted using Antechamber  with the general Amber force field (GAFF) [46, 48]. Data were analyzed using the CPPTRAJ utility in AMBER  and custom VMD TCL scripts . Time evolution of the rmsd values for the HBV-RT in each simulation system was monitored for simulation convergence (Additional file 5: Figure S5).
Binding free energy approximation
CHB, chronic hepatitis B; dATP, deoxyadenosine triphosphate; DNA, deoxyribonucleic acid; ds-DNA, double stranded DNA; HBV, hepatitis B virus; MD, molecular dynamics; NNRTI, non-nucleoside reverse transcriptase inhibitors; NRTIs, nucleoside and nucleotide reverse transcriptase inhibitors; rmsd, root mean square deviation; RNA, ribonucleic acid; RR, rapid response; RT, reverse transcriptase; SIE, solvated interaction energy; SR, slow response; TDF, tenofovir disoproxil fumarate; TFV-DP, tenofovir diphosphate (active form of TDF).
We thank Dr. Formijn J van Hemert from University of Amsterdam and Dr. Jianyi Yang from University of Michigan for advises on modeling the HBV-RT structure. We thank Dr. Robert Wohlhueter (Georgia State University) and Dr. Michael Mitchell (Gilead Science, Inc.) for their contribution in editing the manuscript.
This article has been published as part of BMC Bioinformatics Volume 17 Supplement 8, 2016. Selected articles from the 11th International Symposium on Bioinformatics Research and Applications (ISBRA ‘15): bioinformatics. The full contents of the supplement are available online https://bmcbioinformatics.biomedcentral.com/articles/supplements/volume-17-supplement-8.
This study was supported by CDC intramural funding and by APHL postdoctoral fellowship funding (2014-2015) to [X.X]. This information is distributed solely for the purpose of pre-dissemination peer review under applicable information quality guidelines. It has not been formally disseminated by the Centers for Disease Control and Prevention/Agency for Toxic Substances and Disease Registry. It does not represent and should not be construed to represent any agency determination or policy. Publication costs for this work were funded by CDC.
Availability of data and materials
Data and materials are available upon request.
XX, YK and JL designed the experiment; XX implemented the 3D modeling and MD simulations, analyzed the 3D output data and prepared the manuscript; HT and LG-R conducted DNA sequencing of HBV variants; GX conducted clean up and formatting of DNA sequence data; KK, AG and MP provided serum samples and patient clinical/demographic information. All authors discussed the results and implications and commented on the manuscript at all stages.
X.X., H.T., G.X., L.G-R., Y.K. and J.L. declare no competing interests. At the time the work was completed, K.K., A.G., and M.P. were employees and stockholders of Gilead Sciences.
Consent for publication
Ethics approval and consent to participate
All patients provided written consent, and the clinical trial protocol was approved by the Institutional Review Boards of each participating institution.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Menendez-Arias L, Alvarez M, Pacheco B. Nucleoside/nucleotide analog inhibitors of hepatitis B virus polymerase: mechanism of action and resistance. Curr Opin Virol. 2014;8:1–9.View ArticlePubMedGoogle Scholar
- Coppolino G et al. The case of chronic hepatitis B treatment with tenofovir: an update for nephrologists. J Nephrol. 2015;28(4):393–402.View ArticlePubMedGoogle Scholar
- Arimilli MN et al. Synthesis, in vitro biological evaluation and oral bioavailability of 9-[2-(phosphonomethoxy)propyl]adenine (PMPA) prodrugs. Antivir Chem Chemother. 1997;8(6):557–64.Google Scholar
- Robbins BL et al. Anti-human immunodeficiency virus activity and cellular metabolism of a potential prodrug of the acyclic nucleoside phosphonate 9-R-(2-phosphonomethoxypropyl)adenine (PMPA), Bis(isopropyloxymethylcarbonyl)PMPA. Antimicrob Agents Chemother. 1998;42(3):612–7.PubMedPubMed CentralGoogle Scholar
- Miller MD et al. Antiviral activity of tenofovir (PMPA) against nucleoside-resistant clinical HIV samples. Nucleosides Nucleotides Nucleic Acids. 2001;20(4-7):1025–8.View ArticlePubMedGoogle Scholar
- Durand-Gasselin L et al. Nucleotide analogue prodrug tenofovir disoproxil enhances lymphoid cell loading following oral administration in monkeys. Mol Pharm. 2009;6(4):1145–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Perry CM, Simpson D. Tenofovir disoproxil fumarate: in chronic hepatitis B. Drugs. 2009;69(16):2245–56.View ArticlePubMedGoogle Scholar
- Michailidis E et al. Antiviral therapies: focus on hepatitis B reverse transcriptase. Int J Biochem Cell Biol. 2012;44(7):1060–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Jones SA, Hu J. Hepatitis B virus reverse transcriptase: diverse functions as classical and emerging targets for antiviral intervention. Emerg Microbes Infect. 2013;2:e56.View ArticlePubMedPubMed CentralGoogle Scholar
- Jung SK et al. Tenofovir disoproxil fumarate monotherapy for nucleos(t)ide analogue-naive and nucleos(t)ide analogue-experienced chronic hepatitis B patients. Clin Mol Hepatol. 2015;21(1):41–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Kitrinos KM et al. No detectable resistance to tenofovir disoproxil fumarate after 6 years of therapy in patients with chronic hepatitis B. Hepatology. 2014;59(2):434–42.View ArticlePubMedGoogle Scholar
- Marcellin P, Asselah T. Long-term therapy for chronic hepatitis B: hepatitis B virus DNA suppression leading to cirrhosis reversal. J Gastroenterol Hepatol. 2013;28(6):912–23.View ArticlePubMedGoogle Scholar
- Baran B et al. Efficacy of tenofovir in patients with Lamivudine failure is not different from that in nucleoside/nucleotide analogue-naive patients with chronic hepatitis B. Antimicrob Agents Chemother. 2013;57(4):1790–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Das K et al. Structural basis for the role of the K65R mutation in HIV-1 reverse transcriptase polymerization, excision antagonism, and tenofovir resistance. J Biol Chem. 2009;284(50):35092–100.View ArticlePubMedPubMed CentralGoogle Scholar
- Hachiya A et al. K70Q adds high-level tenofovir resistance to “Q151M Complex” HIV reverse transcriptase through the enhanced discrimination mechanism. PLoS One. 2011;6(1):e16242.View ArticlePubMedPubMed CentralGoogle Scholar
- White KL et al. Molecular mechanisms of tenofovir resistance conferred by human immunodeficiency virus type 1 reverse transcriptase containing a diserine insertion after residue 69 and multiple thymidine analog-associated mutations. Antimicrob Agents Chemother. 2004;48(3):992–1003.View ArticlePubMedPubMed CentralGoogle Scholar
- van Hemert FJ, Berkhout B, Zaaijer HL. Differential binding of tenofovir and adefovir to reverse transcriptase of hepatitis B virus. PLoS One. 2014;9(9):e106324.View ArticlePubMedPubMed CentralGoogle Scholar
- Das K et al. Molecular modeling and biochemical characterization reveal the mechanism of hepatitis B virus polymerase resistance to lamivudine (3TC) and emtricitabine (FTC). J Virol. 2001;75(10):4771–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Tuske S et al. Structures of HIV-1 RT-DNA complexes before and after incorporation of the anti-AIDS drug tenofovir. Nat Struct Mol Biol. 2004;11(5):469–74.View ArticlePubMedGoogle Scholar
- Zhang Y. I-TASSER server for protein 3D structure prediction. BMC Bioinformatics. 2008;9:40.View ArticlePubMedPubMed CentralGoogle Scholar
- Roy A, Kucukural A, Zhang Y. I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc. 2010;5(4):725–38.View ArticlePubMedPubMed CentralGoogle Scholar
- Das K et al. Structures of HIV-1 RT-RNA/DNA ternary complexes with dATP and nevirapine reveal conformational flexibility of RNA/DNA: insights into requirements for RNase H cleavage. Nucleic Acids Res. 2014;42(12):8125–37.View ArticlePubMedPubMed CentralGoogle Scholar
- Xu D et al. AIDA: ab initio domain assembly server. Nucleic Acids Res. 2014;42(Web Server issue):W308–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Naim M et al. Solvated interaction energy (SIE) for scoring protein-ligand binding affinities. 1. Exploring the parameter space. J Chem Inf Model. 2007;47(1):122–33.View ArticlePubMedGoogle Scholar
- Sulea T, Cui Q, Purisima EO. Solvated interaction energy (SIE) for scoring protein-ligand binding affinities. 2. Benchmark in the CSAR-2010 scoring exercise. J Chem Inf Model. 2011;51(9):2066–81.View ArticlePubMedGoogle Scholar
- Lara, J., et al., P0633 : Association of hepatitis B virus genetic diversity with tenofovir disoproxil fumarate response rates in immune tolerant chronic hepatitis B patients. Journal of Hepatology. 2015;62: p. S556.Google Scholar
- Campo DS et al. Coordinated evolution of the hepatitis B virus polymerase. In Silico Biol. 2011;11(5-6):175–82.PubMedGoogle Scholar
- Eldar A, Elowitz MB. Functional roles for noise in genetic circuits. Nature. 2010;467(7312):167–73.View ArticlePubMedPubMed CentralGoogle Scholar
- Sharma, S.V., et al., A Chromatin-Mediated Reversible Drug-Tolerant State in Cancer Cell Subpopulations. Cell. 2010;141(1): p. 69-80.Google Scholar
- Schauer GD et al. Mechanism of allosteric inhibition of HIV-1 reverse transcriptase revealed by single-molecule and ensemble fluorescence. Nucleic Acids Res. 2014;42(18):11687–96.View ArticlePubMedPubMed CentralGoogle Scholar
- Daga PR, Duan J, Doerksen RJ. Computational model of hepatitis B virus DNA polymerase: molecular dynamics and docking to understand resistant mutations. Protein Sci. 2010;19(4):796–807.View ArticlePubMedPubMed CentralGoogle Scholar
- Walsh AW et al. Mechanistic characterization and molecular modeling of hepatitis B virus polymerase resistance to entecavir. PLoS One. 2010;5(2):e9195.View ArticlePubMedPubMed CentralGoogle Scholar
- Lin X et al. A single amino acid in the reverse transcriptase domain of hepatitis B virus affects virus replication efficiency. J Virol. 2001;75(23):11827–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Wiederstein M, Sippl MJ. ProSA-web: interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res. 2007;35(Web Server issue):W407–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Chan HL et al. Effects of tenofovir disoproxil fumarate in hepatitis B e antigen-positive patients with normal levels of alanine aminotransferase and high levels of hepatitis B virus DNA. Gastroenterology. 2014;146(5):1240–8.View ArticlePubMedGoogle Scholar
- Bandelt HJ, Forster P, Rohl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16(1):37–48.View ArticlePubMedGoogle Scholar
- Eswar N et al. Comparative protein structure modeling using Modeller. Curr Protoc Bioinformatics. 2006;Chapter 5:Unit 5 6.PubMedGoogle Scholar
- Webb B, Sali A. Comparative protein structure modeling using MODELLER. Curr Protoc Bioinformatics. 2014;47:5 6 1–5 6 32.View ArticleGoogle Scholar
- Case DA et al. The Amber biomolecular simulation programs. J Comput Chem. 2005;26(16):1668–88.View ArticlePubMedPubMed CentralGoogle Scholar
- Vriend G. WHAT IF: a molecular modeling and drug design program. J Mol Graph. 1990;8(1):52–6. 29.View ArticlePubMedGoogle Scholar
- Jorgensen WL et al. Comparison of simple potential functions for simulating liquid water. J Chem Phys. 1983;79(2):10.View ArticleGoogle Scholar
- Essmann U et al. A smooth particle mesh Ewald method. J Chem Phys. 1995;103:8577–93.View ArticleGoogle Scholar
- Tuckerman M, Berne BJ, Martyna GJ. Reversible multiple time scale molecular-dynamics. J Chem Phys. 1992;97(3):1990–2001.View ArticleGoogle Scholar
- Kale L et al. NAMD2: Greater scalability for parallel molecular dynamics. J Comput Phys. 1999;151(1):283–312.View ArticleGoogle Scholar
- Phillips JC et al. Scalable molecular dynamics with NAMD. J Comput Chem. 2005;26:1781–802.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang J et al. Development and testing of a general amber force field. J Comput Chem. 2004;25(9):1157–74.View ArticlePubMedGoogle Scholar
- Wang J et al. Automatic atom type and bond type perception in molecular mechanical calculations. J Mol Graph Model. 2006;25(2):247–60.View ArticlePubMedGoogle Scholar
- Hornak V et al. Comparison of multiple Amber force fields and development of improved protein backbone parameters. Proteins. 2006;65(3):712–25.View ArticlePubMedPubMed CentralGoogle Scholar
- Roe DR, Cheatham TE. PTRAJ and CPPTRAJ: software for processing and analysis of molecular dynamics trajectory data. J Chem Theory Comput. 2013;9(7):3084–95.View ArticlePubMedGoogle Scholar
- Humphrey W, Dalke A, Schulten K. VMD: visual molecular dynamics. J Mol Graph Model. 1996;14(1):33–8.View ArticleGoogle Scholar
- Aqvist J, Medina C, Samuelsson JE. A new method for predicting binding affinity in computer-aided drug design. Protein Eng. 1994;7(3):385–91.View ArticlePubMedGoogle Scholar
- Jones-Hertzog DK, Jorgensen WL. Binding affinities for sulfonamide inhibitors with Human Thrombin Using Monte Carlo Simulations with a linear response method. J Med Chem. 1997;40(10):1539–49.View ArticlePubMedGoogle Scholar