Diffusion of water and sodium counter-ions in nanopores of a β-lactoglobulin crystal: a molecular dynamics study.

The dynamics of water and sodium counter-ions (Na+) in a C2221 orthorhombic β-lactoglobulin crystal is investigated by means of 5 ns molecular dynamics simulations. The effect of the fluctuation of the protein atoms on the motion of water and sodium ions is studied by comparing simulations in a rigid and in a flexible lattice. The electrostatic interactions of sodium ions with the positively charged LYS residues inside the crystal channels significantly influence the ionic motion. According to our results, water molecules close to the protein surface undergo an anomalous diffusive motion. On the other hand, the motion of water molecules further away from the protein surface is normal diffusive. Protein fluctuations affect the diffusion constant of water, which increases from 0.646 ± 0.108 to 0.887 ± 0.41 nm2 ns−1, when protein fluctuations are taken into account. The pore size (0.63–1.05 nm) and the water diffusivities are in good agreement with previous experimental results. The dynamics of sodium ions is disordered. LYS residues inside the pore are the main obstacles to the motion of sodium ions. However, the simulation time is still too short for providing a precise description of anomalous diffusion of sodium ions. The results are not only of interest for studying ion and water transport through biological nanopores, but may also elucidate water–protein and ion–protein interactions in protein crystals.