Molecular dynamics simulations of HPr under hydrostatic pressure

The histidine-containing protein (HPr) plays an important role in the phosphotransferase system (PTS). The deformations induced on the protein structure at high hydrostatic pressure values (4, 50, 100, 150, and 200 MPa) were previously (H. Kalbitzer, A. Gorler, H. Li, P. Dubovskii, A. Hengstenberg, C. Kowolik, H. Yamada, and K. Akasaka, Protein Science 2000, Vol. 9, pp. 693-703) analyzed by NMR experiments: the nonlinear variations of the amide chemical shifts at high pressure values were supposed to arise from induced shifts in the protein conformational equilibrium. Molecular dynamics (MD) simulations are here performed, to analyze the protein internal mobility at 0.1 MPa, and to relate the nonlinear variations of chemical shifts observed at high pressure, to variations in conformational equilibrium. The global features of the protein structure are only slightly modified along the pressure. Nevertheless, the values of the Voronoi residues volumes show that the residues of (alpha)-helices are more compressed that those belonging to the (beta)-sheet. The (alpha)-helices are also displaying the largest internal mobility and deformation in the simulations. The nonlinearity of the 1H chemical shifts, computed from the MD simulation snapshots, is in qualitative agreement with the nonlinearity of the experimentally observed chemical shifts.