title: Determination of the protonation state of the Asp dyad: conventional molecular dynamics versus thermodynamic integration.
abstract: The protonation state of the Asp dyad is important as it can reveal enzymatic mechanisms, and the information this provides can be used in the development of drugs for proteins such as memapsin 2 (BACE-1), HIV-1 protease, and rennin. Conventional molecular dynamics (MD) simulations have been successfully used to determine the preferred protonation state of the Asp dyad. In the present work, we demonstrate that the results obtained from conventional MD simulations can be greatly influenced by the particular force field applied or the values used for control parameters. In principle, free-energy changes between possible protonation states can be used to determine the protonation state. We show that protonation state prediction by the thermodynamic integration (TI) method is insensitive to force field version or to the cutoff for calculating nonbonded interactions (a control parameter). In the present study, the protonation state of the Asp dyad predicted by TI calculations was the same regardless of the force field and cutoff value applied. Contrary to the intuition that conventional MD is more efficient, our results clearly show that the TI method is actually more efficient and more reliable for determining the protonation state of the Asp dyad.
authors.affiliation: School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, 13 Hangkong Road, Wuhan, Hubei, 430030, People's Republic of China.
authors.forename: Jinfeng
authors.initials: J
authors.lastname: Huang
authors.affiliation: School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, 13 Hangkong Road, Wuhan, Hubei, 430030, People's Republic of China.
authors.forename: Yali
authors.initials: Y
authors.lastname: Zhu
authors.affiliation: School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, 13 Hangkong Road, Wuhan, Hubei, 430030, People's Republic of China.
authors.forename: Bin
authors.initials: B
authors.lastname: Sun
authors.affiliation: Institute of Theoretical and Simulational Chemistry, The Academy of Fundamental and Interdisciplinary Science, Harbin Institute of Technology, 2 Yikuang Street, Harbin, 150080, People's Republic of China.
authors.forename: Yuan
authors.initials: Y
authors.lastname: Yao
authors.affiliation: School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, 13 Hangkong Road, Wuhan, Hubei, 430030, People's Republic of China. junjun.liu@hust.edu.cn.
authors.forename: Junjun
authors.initials: J
authors.lastname: Liu
mesh_terms: D000465:Algorithms; D053829:Amyloid Precursor Protein Secretases / Q000737:chemistry; D016282:Aspartic Acid Endopeptidases / Q000737:chemistry; D006860:Hydrogen Bonding; D008958:Models, Molecular; D056004:Molecular Dynamics Simulation; D011522:Protons; D013816:Thermodynamics
chemical_list: D011522:Protons; D053829:Amyloid Precursor Protein Secretases; D016282:Aspartic Acid Endopeptidases; C509280:BACE1 protein, human
keywords: Asp dyad; BACE-1; Molecular dynamics; Protonation state; Thermodynamic integration