title: Multiple Gaussian network modes alignment reveals dynamically variable regions: the hemoglobin case.
abstract: Gaussian network model (GNM) modes of motion are calculated to a dataset of hemoglobin (Hb) structures and modes with dynamics similarity to the T state are multiply aligned. The sole criterion for the alignment is the mode shape itself and not sequence or structural similarity. Standard deviation (SD) of the GNM value score along the alignment is calculated, regions with high SD are defined as dynamically variable. The analysis shows that the α1β1/α2β2 interface is a dynamically variable region but not the α1β2/α2β1 and the α1α2/β1β2 interfaces. The results are in accordance with the T→R2 transition of Hb. We suggest that dynamically variable regions are regions that are likely to undergo structural change in the protein upon binding, conformational transition, or any other relevant chemical event. The represented technique of multiple dynamics-based alignment of modes is novel and may offer a new insight in proteins' dynamics to function relation.
authors.affiliation: Department of Computer Sciences and Mathematics, Ariel University, Ariel, 40700, Israel.
authors.forename: Meir
authors.initials: M
authors.lastname: Davis
authors.forename: Dror
authors.initials: D
authors.lastname: Tobi
mesh_terms: D000465:Algorithms; D000595:Amino Acid Sequence; D014644:Genetic Variation / Q000235:genetics*; D006454:Hemoglobins / Q000032:analysis / Q000737:chemistry* / Q000235:genetics*; D056004:Molecular Dynamics Simulation; D016011:Normal Distribution; D017434:Protein Structure, Tertiary; D016415:Sequence Alignment
chemical_list: D006454:Hemoglobins
keywords: Gaussian network model; dynamics alignment; hemoglobin; multiple dynamics alignment; normal mode analysis; protein dynamics