Molecular dynamics simulations for water and ions in protein crystals

The spatial and temporal properties of water and ions in bionanoporous materials-protein crystals-have been investigated using molecular dynamics simulations. Three protein crystals are considered systematically with different morphologies and chemical topologies: tetragonal lysozyme, orthorhombic lysozyme, and tetragonal thermolysin. It is found that the thermal fluctuations of C(alpha) atoms in the secondary structures of protein molecules are relatively weak due to hydrogen bonding. The solvent-accessible surface area per residue is nearly identical in the three protein crystals; the hydrophobic and hydrophilic residues in each crystal possess approximately the same solvent-accessible surface area. Water distributes heterogeneously and has different local structures within the biological nanopores of the three protein crystals. The mobility of water and ions in the crystals is enhanced as the porosity increases and also by the fluctuations of protein atoms particularly in the two lysozyme crystals. Anisotropic diffusion is found preferentially along the pore axis, as experimentally observed. The anisotropy of the three crystals increases in the order: tetragonal thermolysin < tetragonal lysozyme < orthorhombic lysozyme.