Rotation of Lipids in Membranes: Molecular Dynamics Simulation, 31P Spin-Lattice Relaxation, and Rigid-Body Dynamics

Abstract

Molecular dynamics simulations and 31P-NMR spin-lattice (R1) relaxation rates from 0.022 to 21.1 T of fluid phase dipalmitoylphosphatidylcholine bilayers are compared. Agreement between experiment and direct prediction from simulation indicates that the dominant slow relaxation (correlation) times of the dipolar and chemical shift anisotropy spin-lattice relaxation are ∼10ns and 3ns, respectively. Overall reorientation of the lipid body, consisting of the phosphorus, glycerol, and acyl chains, is well described within a rigid-body model. Wobble, with D⊥=1–2×108s−1, is the primary component of the 10ns relaxation; this timescale is consistent with the tumbling of a lipid-sized cylinder in a medium with the viscosity of liquid hexadecane. The value for D|| the diffusion constant for rotation about the long axis of the lipid body, is difficult to determine precisely because of averaging by fast motions and wobble; it is tentatively estimated to be 1×107s−1. The resulting D||/D⊥≈0.1 implies that axial rotation is strongly modulated by interactions at the lipid/water interface. Rigid-body modeling and potential of mean force evaluations show that the choline group is relatively uncoupled from the rest of the lipid. This is consistent with the ratio of chemical shift anisotropy and dipolar correlation times reported here and the previous observations that 31P-NMR lineshapes are axially symmetric even in the gel phase of dipalmitoylphosphatidylcholine.