Predicting the packing parameter for lipids in monolayers with the use of molecular dynamics

Lipid molecules form the backbone of biological membranes. Due to their amphiphilic structure, they can self-organize in a plethora of different structures when in contact with water. The type of self-assembled structure and its curvature depend on so-called shape factor or critical packing parameter, CPP, that can be derived knowing the molecular volume of a lipid (V), optimal surface area (a₀) and critical chain length (l_c) (see Intermolecular and Surface Forces by Jacob N. Israelachvili, Third Edition, 2011). The value of CPP allows not only to predict the type of self-assembled structure but also is a key factor for molecular interactions, which play a great role both in physiological and pathological conditions. The greatest difficulties arise when calculating the a₀ parameter, and although for some typical membrane lipids these values have been determined, there are a number of derivatives for which this parameter, and thus CPP, are unknown. The value of CPP allows not only to predict the type of self-assembled structure but also is a key factor for molecular interactions, which play a great role both in physiological and pathological conditions. So far, the determination of the packing parameter required the use of theoretical models with assumptions deviating from the physical conditions. Here we report a method based on molecular dynamics, which was applied to simulate lipid membranes consisting of cholesterol, oxysterols, sphingolipids, phosphatidylcholines, and phosphatidylethanolamines. For lipid molecules for which CPPs have already been determined, high compliance has been demonstrated. This proves that the method presented herein can be successfully used to determine packing parameters for other membrane lipids and amphiphilic molecules.