Acyl chain order parameter profiles in phospholipid bilayers: computation from molecular dynamics simulations and comparison with 2H NMR experiments

Order parameters from deuterium NMR are often used to validate or calibrate molecular dynamics simulations. This paper gives a short overview of the literature in which experimental order parameters from 2H NMR are compared to those calculated from MD simulations. The different ways in which order parameters from experiment are used to calibrate and validate simulations are reviewed. In the second part of this review, a case study of cholesterol in a DMPC bilayer is presented. It is concluded that the agreement between experimental data and simulation is favorable in the hydrophobic region of the membrane, for both the phospholipids and cholesterol. In the interfacial region the agreement is less satisfactory, probably because of the high polarity of this region which makes the correct computation of the electrostatics more complex.