Hydration of Proteins: A Comparison of Experimental Residence Times of Water Molecules Solvating the Bovine Pancreatic Trypsin Inhibitor with Theoretical Model Calculations

Abstract

A 1 ns trajectory from a molecular dynamics study of 1·4 ns total length was used for a detailed analysis of the residence times of water molecules located near 227 selected bovine pancreatic trypsin inhibitor (BPTI) atoms. The simulation was performed using the GROMOS force field, with apolar hydrogen atoms treated as united atoms, and the SPC/E water model. The system consisted of 568 BPTI atoms and 2371 water molecules. The theoretical results are in good agreement with experimental data available from nuclear magnetic resonance spectroscopy. The residence times of individual water molecules coming near a give BPTI atoms, as obtained from the simulation, vary greatly and range betweem 10 and 500 ps. The effective residence time, calculated using a correlation function technique from the presence of all individual water molecules visiting the hydration shell of a given BPTI atom, never exceeds 200 ps. The average residence time near backbone and side-chain atoms is approximately 39 ps and 24 ps, respectively. The shortest residence times, on average, are found near charged atoms (19 ps), whereas near non-polar and polar side-chain atoms the residence times are 25 ps and 36 ps, respectively. There is no apparent correlation between the residence times of the hydration water molecules of solvent-accessible residues and their location in different regular or non-regular secondary structures.