Skip to content
Licensed Unlicensed Requires Authentication Published by De Gruyter April 17, 2020

Free energy profile of permeation of Entecavir through Hepatitis B virus capsid studied by molecular dynamics calculation

  • Kazushi Fujimoto , Motohiro Fukai , Ryo Urano , Wataru Shinoda , Tetsuya Ishikawa , Katsumi Omagari , Yasuhito Tanaka , Atsushi Nakagawa and Susumu Okazaki EMAIL logo

Abstract

Entecavir, triphosphorylated in liver cells, is an antiviral reagent against Hepatitis B virus (HBV). The reagent inhibits reverse transcription of RNA inside the virus capsid. In the present study, free energy profile of an Entecavir triphosphate (ETVTP) molecule has been calculated when it passes through pores of the capsid along two- and three-fold rotational symmetry axes in order to investigate permeation pathway of the reagent to the inside of the capsid. The calculations have been done based on thermodynamic integration (TI) method combined with all-atomistic molecular dynamic (MD) calculations. A free energy minimum of −19 kJ/mol was found at the entrance of the pore from the outside along the three-fold symmetry axis. This stabilization is from the interaction of negatively charged ETVTP with positively charged capsid methionine residues. This excess free energy concentrates of the reagent at the entrance of the pore by a factor of about 2000. A free energy barrier of approximately 13 kJ/mol was also found near the exit of the pore to the inside of the capsid due to narrow space of the pore surrounded by hydrophobic wall made by proline residues and negatively charged wall by aspartic acid residues. There, ETVTP is partially dehydrated in order to pass through the narrow space, which causes the great free energy loss. Further, the negatively charged residues produce repulsive forces on the ETVTP molecule. In contrast, in the case of the pore along the two-fold symmetry axis, the calculated free energy profile showed shallower free energy minimum, −4 kJ/mol at the entrance in spite of the similarly high barrier, 7 kJ/mol, near the exit of the pore.


Article note

A collection of invited papers based on presentations at the 36th International Conference of Solution Chemistry (ICSC-36), held in Xining, China, 4–8 August 2019.


Acknowledgment

We would like to thank Mr. S. Tsurumaki for preparing the figures of this paper. This work was supported by MEXT as “Priority Issue on Post-K computer” (Building Innovative Drug Discovery Infrastructure Through Functional Control of Biomolecular Systems) and “General Issue on Post-K computer” (Absorption Mechanism of Inhibitor into Capsid by using All-Atom Molecular Dynamics of Hepatitis B Virus in Solution) using computational resources of the K computer provided by the Funder Id: http://dx.doi.org/10.13039/501100006321, RIKEN, Funder Id: http://dx.doi.org/10.13039/501100006264, Center for Computational Science through the HPCI System Research project (Project ID: hp170255 and hp170205, respectively). Calculations were also performed at the Research Center for Computational Science, Okazaki, Japan, at the Information Technology Center of Nagoya University, and at the Institute for Solid State Physics, the University of Tokyo.

References

[1] WHO. Global Hepatitis Report, 2017 (2017).Search in Google Scholar

[2] D. Knipe, P. Howley. Fields Virology: Sixth ed, 05, Lippincott Williams & Wilkins, Philadelphia, PA, USA (2013).Search in Google Scholar

[3] P. Karayiannis. J. Antimicrob. Chemother. 51, 761 (2003).10.1093/jac/dkg163Search in Google Scholar

[4] D. J. Tenney, R. E. Rose, C. J. Baldick, S. M. Levine, K. A. Pokornowski, A. W. Walsh, J. Fang, C.-F. Yu, S. Zhang, C. E. Mazzucco, B. Eggers, M. Hsu, M. J. Plym, P. Poundstone, J. Yang, R. J. Colonno. Antimicrob. Agents Chemother. 51, 902 (2007).10.1128/AAC.00833-06Search in Google Scholar

[5] D. J. Tenney, R. E. Rose, C. J. Baldick, K. A. Pokornowski, B. J. Eggers, J. Fang, M. J. Wichroski, D. Xu, J. Yang, R. B. Wilber, R. J. Colonno. Hepatology 49, 1503 (2009).10.1002/hep.22841Search in Google Scholar

[6] A. Zlotnick, N. Cheng, J. F. Conway, F. P. Booy, A. C. Steven, S. J. Stahl, P. T. Wing-Field. Biochemistry 35, 7412 (1996).10.1021/bi9604800Search in Google Scholar

[7] S. A. Wynne, R. A. Crowther, A. G. W. Leslie. Mol. Cell 3, 771 (1999).10.1016/S1097-2765(01)80009-5Search in Google Scholar

[8] D. Frenkel, B. Smit. Understanding Molecular Simulation From Algorithms to Applications, 2nd ed. Academic Press, Cambridge, MA, USA (2002).Search in Google Scholar

[9] R. B. Best, X. Zhu, J. Shim, P. E. M. Lopes, J. Mittal, M. Feig, A. D. MacKerell. J. Chem. Theory Comput. 8, 3257 (2012). PMID: 23341755.10.1021/ct300400xSearch in Google Scholar PubMed PubMed Central

[10] D. Beglov, B. Roux. J. Chem. Phys. 100, 9050 (1994).10.1063/1.466711Search in Google Scholar

[11] W. L. Jorgensen, J. Chandrasekhar, J. D. Madura, R. W. Impey, M. L. Klein. J. Chem. Phys. 79, 926 (1983).10.1063/1.445869Search in Google Scholar

[12] E. J. Denning, U. Deva Priyakumar, L. Nilsson, A. D. Mackerell Jr. J. Comput. Chem. 32, 1929 (2011).10.1002/jcc.21777Search in Google Scholar PubMed PubMed Central

[13] C. G. Mayne, J. Saam, K. Schulten, E. Tajkhorshid, J. C. Gumbart. J. Comput. Chem. 34, 2757 (2013).10.1002/jcc.23422Search in Google Scholar PubMed PubMed Central

[14] U. Essmann, L. Perera, M. L. Berkowitz, T. Darden, H. Lee, L. G. Pedersen. J. Chem. Phys. 103, 8577 (1995).10.1063/1.470117Search in Google Scholar

[15] W. G. Hoover. Phys. Rev. A 31, 1695 (1985).10.1103/PhysRevA.31.1695Search in Google Scholar

[16] M. Parrinello, A. Rahman. Phys. Rev. Lett. 45, 1196 (1980).10.1103/PhysRevLett.45.1196Search in Google Scholar

[17] B. Hess, H. Bekker, H. J. C. Berendsen, J. Fraaije. J. Comput. Chem. 18, 1463 (1997).10.1002/(SICI)1096-987X(199709)18:12<1463::AID-JCC4>3.0.CO;2-HSearch in Google Scholar

[18] B. Hess. J. Chem. Theory Comput. 4, 116 (2008).10.1021/ct700200bSearch in Google Scholar

[19] S. Miyamoto, P. A. Kollman. J. Comput. Chem. 13, 952 (1992).10.1002/jcc.540130805Search in Google Scholar

[20] B. Hess, C. Kutzner, D. van der Spoel, E. Lindah. J. Chem. Theory Comput. 4, 435 (2008).10.1021/ct700301qSearch in Google Scholar

[21] B. Hess, G. Groenhof, A. E. Mark, D. van der Spoel, E. Lindahl, H. J. C. Berendsen. J. Comp. Chem. 26, 1701 (2005).10.1002/jcc.20291Search in Google Scholar

[22] E. Lindahl, B. Hess, D. van der Spoel. J. Mol. Model. 7, 306 (2001).10.1007/s008940100045Search in Google Scholar

[23] D. van der Spoel, H. J. C. Berendsen, R. van Drunen. Comp. Phys. Comm. 91, 43 (1995).10.1016/0010-4655(95)00042-ESearch in Google Scholar


Supplementary Material

The online version of this article offers supplementary material (https://doi.org/10.1515/pac-2020-0109).


Published Online: 2020-04-17
Published in Print: 2020-10-25

©2020 IUPAC & De Gruyter. This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License. For more information, please visit: http://creativecommons.org/licenses/by-nc-nd/4.0/

Downloaded on 23.4.2024 from https://www.degruyter.com/document/doi/10.1515/pac-2020-0109/html
Scroll to top button