Polyunsaturated lipids are an essential component of biological membranes, influencing order
and dynamics of lipids, protein−lipid interaction, and membrane transport properties. To gain an atomic
level picture of the impact of polyunsaturation on membrane properties, quantum mechanical (QM) and
empirical force field based calculations have been undertaken. The QM calculations of the torsional energy
surface for rotation about vinyl−methylene bonds reveal low barriers to rotation, indicating an intrinsic
propensity toward flexibility. Based on QM and experimental data, empirical force field parameters were
developed for polyunsaturated lipids and applied in a 16 ns molecular dynamics (MD) simulation of a
1-stearoyl-2-docosahexaenoyl-sn-glyerco-3-phosphocholine (SDPC) lipid bilayer. The simulation results
are in good agreement with experimental data, suggesting an unusually high degree of conformational
flexibility of polyunsaturated hydrocarbon chains in membranes. The detailed analysis of chain conformation
and dynamics by simulations is aiding the interpretation of experimental data and is useful for understanding
the unique role of polyunsaturated lipids in biological membranes. The complete force field is included as
Supporting Information and is available from http://www.pharmacy.umaryland.edu/faculty/amackere/research.html