Water dynamics around T0 vs R4 of hemoglobin from local hydrophobicity analysis

Abstract

The local hydration around tetrameric hemoglobin (Hb) in its T0 and R4 conformational substates is analyzed based on molecular dynamics simulations. Analysis of the local hydrophobicity (LH) for all residues at the α1 β2 and α2 β1 interfaces, responsible for the quaternary T → R transition, which is encoded in the Monod–Wyman–Changeux model, as well as comparison with earlier computations of the solvent accessible surface area, makes clear that the two quantities measure different aspects of hydration. Local hydrophobicity quantifies the presence and structure of water molecules at the interface, whereas “buried surface” reports on the available space for solvent. For simulations with Hb frozen in its T0 and R4 states, the correlation coefficient between LH and buried surface is 0.36 and 0.44, respectively, but it increases considerably if the 95% confidence interval is used. The LH with Hb frozen and flexible changes little for most residues at the interfaces but is significantly altered for a few select ones: Thr41 α, Tyr42 α, Tyr140 α, Trp37 β, Glu101 β (for T0) and Thr38 α, Tyr42 α, Tyr140 α (for R4). The number of water molecules at the interface is found to increase by [Formula: see text]% for T0 → R4, which is consistent with earlier measurements. Since hydration is found to be essential to protein function, it is clear that hydration also plays an essential role in allostery.