Can the protonation state of histidine residues be determined from molecular dynamics simulations?
Uranga, Jon ; Mikulskis, Paulius LU ; Genheden, Samuel LU and Ryde, Ulf LU
(2012)
In Computational and Theoretical Chemistry
1000.
p.75-84
- Abstract
- Histidine (His) residues in proteins can attain three different protonation states at normal pH. This constitutes a prominent problem when adding protons to a protein crystal structure, e.g. in order to perform molecular simulations. Typically, the His protonation is deduced from the hydrogen-bond pattern in crystal structures. Here, we study whether it is possible to detect erroneous His protonation state by analysing short molecular dynamics (MD) trajectories. We systematically vary the His protonation state and measure the root-mean-squared deviation (RMSD) of the His residues and nearby residues relative to the starting structure, as well as the distribution of the dihedral angle that determines the rotation of the His side chain. We... (More)
- Histidine (His) residues in proteins can attain three different protonation states at normal pH. This constitutes a prominent problem when adding protons to a protein crystal structure, e.g. in order to perform molecular simulations. Typically, the His protonation is deduced from the hydrogen-bond pattern in crystal structures. Here, we study whether it is possible to detect erroneous His protonation state by analysing short molecular dynamics (MD) trajectories. We systematically vary the His protonation state and measure the root-mean-squared deviation (RMSD) of the His residues and nearby residues relative to the starting structure, as well as the distribution of the dihedral angle that determines the rotation of the His side chain. We study three proteins, hisactophilin with 31 solvent-exposed His residues, galectin-3, for which an experimental assignment is available for two of the His residues, and trypsin, for which the hydrogen-bond analysis is quite conclusive. The results show that improper protonation states have larger RMSD values and larger widths of the dihedral distribution, compared to the correct protonation states. Unfortunately, the variation among different His residues in the same and different proteins is so large that it is hard to define unambiguous thresholds between proper and improper protonation states. Therefore, simulations of all three protonation states are needed for conclusive results. For trypsin, we could obtain a conclusive assignment for all three His residues, which was better than the simple hydrogen-bond analysis. For galectin-3, the MD trajectories confirmed the results of hydrogen-bond analysis and experiments. They also gave additional, more uncertain information for some of the residues. However, for the solvent-exposed His residues in hisactophilin, no unambiguous conclusions regarding the protonation states could be reached. On the other hand, this indicates that protein structures are quite insensitive to the protonation state of the His residues, besides those that involve direct hydrogen bonds to the His side chain. (C) 2012 Elsevier B.V. All rights reserved. (Less)
Please use this url to cite or link to this publication:
https://lup.lub.lu.se/record/3366316
- author
- Uranga, Jon
; Mikulskis, Paulius
LU
; Genheden, Samuel
LU
and Ryde, Ulf
LU
- organization
- publishing date
- 2012
- type
- Contribution to journal
- publication status
- published
- subject
- keywords
- Histidine protonation, Molecular dynamics simulations, Side-chain, rotation, Dihedral distribution
- in
- Computational and Theoretical Chemistry
- volume
- 1000
- pages
- 75 - 84
- publisher
- Elsevier
- external identifiers
-
- wos:000312421100011
- scopus:84868683184
- ISSN
- 2210-2728
- DOI
- 10.1016/j.comptc.2012.09.025
- language
- English
- LU publication?
- yes
- additional info
- The information about affiliations in this record was updated in December 2015. The record was previously connected to the following departments: Theoretical Chemistry (S) (011001039)
- id
- 0af7deda-ae7c-45e9-9fd9-d7dc0a7efda0 (old id 3366316)
- date added to LUP
- 2016-04-01 11:16:14
- date last changed
- 2023-01-10 07:09:59
@article{0af7deda-ae7c-45e9-9fd9-d7dc0a7efda0,
abstract = {{Histidine (His) residues in proteins can attain three different protonation states at normal pH. This constitutes a prominent problem when adding protons to a protein crystal structure, e.g. in order to perform molecular simulations. Typically, the His protonation is deduced from the hydrogen-bond pattern in crystal structures. Here, we study whether it is possible to detect erroneous His protonation state by analysing short molecular dynamics (MD) trajectories. We systematically vary the His protonation state and measure the root-mean-squared deviation (RMSD) of the His residues and nearby residues relative to the starting structure, as well as the distribution of the dihedral angle that determines the rotation of the His side chain. We study three proteins, hisactophilin with 31 solvent-exposed His residues, galectin-3, for which an experimental assignment is available for two of the His residues, and trypsin, for which the hydrogen-bond analysis is quite conclusive. The results show that improper protonation states have larger RMSD values and larger widths of the dihedral distribution, compared to the correct protonation states. Unfortunately, the variation among different His residues in the same and different proteins is so large that it is hard to define unambiguous thresholds between proper and improper protonation states. Therefore, simulations of all three protonation states are needed for conclusive results. For trypsin, we could obtain a conclusive assignment for all three His residues, which was better than the simple hydrogen-bond analysis. For galectin-3, the MD trajectories confirmed the results of hydrogen-bond analysis and experiments. They also gave additional, more uncertain information for some of the residues. However, for the solvent-exposed His residues in hisactophilin, no unambiguous conclusions regarding the protonation states could be reached. On the other hand, this indicates that protein structures are quite insensitive to the protonation state of the His residues, besides those that involve direct hydrogen bonds to the His side chain. (C) 2012 Elsevier B.V. All rights reserved.}},
author = {{Uranga, Jon and Mikulskis, Paulius and Genheden, Samuel and Ryde, Ulf}},
issn = {{2210-2728}},
keywords = {{Histidine protonation; Molecular dynamics simulations; Side-chain; rotation; Dihedral distribution}},
language = {{eng}},
pages = {{75--84}},
publisher = {{Elsevier}},
series = {{Computational and Theoretical Chemistry}},
title = {{Can the protonation state of histidine residues be determined from molecular dynamics simulations?}},
url = {{https://lup.lub.lu.se/search/files/2522738/3412377.pdf}},
doi = {{10.1016/j.comptc.2012.09.025}},
volume = {{1000}},
year = {{2012}},
}