Molecular Dynamics Study on Protein and it's Water Structure at High Pressure

Abstract

We have performed NPT molecular dynamics simulations (Langevin Piston Method) on two types of solvated proteins-‘denaturation-unfavorable’ protein (insulin) and ‘denaturation-favorable protein’ (ribonuclease A) at high pressure (from 1 bar up to 20 kbar). The method is based on the extended system formalism introduced by Andersen, where the deterministic equations of motion for the piston degree of freedom are replaced by Langevin equation. We report the structural changes of proteins (ribonuclease A and insulin) and water molecules through radius of gyration, solvent accessible surface area, hydrogen bond pattern, and the topology of water clusters connected by the hydrogen bonded circular network. The solvent accessibility of ribonuclease A is mainly decreased by hydrophilic residues rather than hydrophobic residues under high pressure. From the results of hydrogen bond analysis, we have found that α-helix is more stable than β-sheet under high pressure. In addition, from the analysis of the water cluster, we have observed that for ribonuclease A, 5-membered ring structure is more favorable than 6-membered ring at higher pressure. However, for insulin, the ratio of 5 to 6-ring is constant over the pressure ranges for which we have performed MD simulation. This indicates that the water structure around insulin does not change under high pressure.