Analysis of water molecules around GTP in Hras-GTP complex and GDP in Hras-GDP complex by molecular dynamics simulations

We investigate the structures of the Hras-GTP and the Hras-GDP complexes in water solvents in order to understand the mechanism of GTP hydrolysis in the Hras-GTP complex. We performed MD simulations of these complexes in order to study the positions and the orientations of water molecules around the guanosine nucleotides. Using trajectories we calculated the angular distribution of water molecules around the most distant phosphorus from guanosine in our previous work. It was shown that water molecules are distributed evenly in GTP, although unevenly in GDP. This suggests that the trigger of GTP hydrolysis is possibly the attack of water molecule to γ?phosphate from the appropriate direction. In this paper, in order to investigate the role of water molecules in GTP hydrolysis in detail, we calculate the orientation of water molecules. The distribution of the orientation is different between GTP and GDP. In order to investigate the cause of this difference, we examine the hydrogen bonds between water molecules and oxygen atom of the most distant phosphate from guanosine. We find that these hydrogen bonds are formed. We also find that the oxygen atom of hydrogen bond is determined by the position of the water molecule of hydrogen bond.