Computer modelling of dynamics of Ser92X deoxymyoglobin mutants

The hydrogen bond between His93 and Ser92, recently discovered in crystal structures of myoglobins (Mbs), may contribute to the oxygen storage capacity of the heme proteins through a stabilization of the proximal ligand. The possible influence of this H-bond on the geometry of the heme proximal side and ligand binding properties of Mb were computationally studied using model proteins with point mutations affecting this bond. The results of the computer modelling of Ser92X (X = Ala, Ile, Thr, Val) mutants of human (H) and sperm whale (SW) Mbs are presented. The OPLS-AMBER-CHARMM forcefield was used in the calculations. Several 10-50 ps molecular dynamics simulations (300 K, in vacuo) were performed. Our results show that the Ser92X mutants are stable molecules. In the wild types and Ser92Thr mutants, the H-bond studied is observed only for a relatively short period of time. It is expected that in both HMb and SW Mb molecules the impact of the proximal histidine interaction with the Ser92(F7) residue on the iron reactivity is rather low. However, the limited torsional flexibility of the proximal histidine imidazole ring was found in hydrogen bonding mutants. This effect may be attributed to the specific long range electrostatic interactions.