Hydration dynamics of a lipid membrane: Hydrogen bond networks and lipid-lipid associations

Dynamics of hydration layers of a dimyristoylphosphatidylcholine (DMPC) bilayer are investigated using an all atom molecular dynamics simulation. Based upon the geometric criteria, continuously residing interface water molecules which form hydrogen bonds solely among themselves and then concertedly hydrogen bonded to carbonyl, phosphate, and glycerol head groups of DMPC are identified. The interface water hydrogen bonded to lipids shows slower relaxation rates for translational and rotational dynamics compared to that of the bulk water and is found to follow sub-diffusive and non-diffusive behaviors, respectively. The mean square displacements and the reorientational auto-correlation functions are slowest for the interfacial waters hydrogen bonded to the carbonyl oxygen since these are buried deep in the hydrophobic core among all interfacial water studied. The intermittent hydrogen bond auto-correlation functions are calculated, which allows breaking and reformations of the hydrogen bonds. The auto-correlation functions for interfacial hydrogen bonded networks develop humps during a transition from cage-like motion to eventual power law behavior of t-3/2. The asymptotic t-3/2 behavior indicates translational diffusion dictated dynamics during hydrogen bond breaking and formation irrespective of the nature of the chemical confinement. Employing reactive flux correlation analysis, the forward rate constant of hydrogen bond breaking and formation is calculated which is used to obtain Gibbs energy of activation of the hydrogen bond breaking. The relaxation rates of the networks buried in the hydrophobic core are slower than the networks near the lipid-water interface which is again slower than bulk due to the higher Gibbs energy of activation. Since hydrogen bond breakage follows a translational diffusion dictated mechanism, chemically confined hydrogen bond networks need an activation energy to diffuse through water depleted hydrophobic environments. Our calculations reveal that the slow relaxation rates of interfacial waters in the vicinity of lipids are originated from the chemical confinement of concerted hydrogen bond networks. The analysis suggests that the networks in the hydration layer of membranes dynamically facilitate the water mediated lipid-lipid associations which can provide insights on the thermodynamic stability of soft interfaces relevant to biological systems in the future.