Cholesterol-induced modifications in lipid bilayers: a simulation study

We present analysis of new configurational bias Monte Carlo and molecular dynamics simulation data for bilayers of dipalmitoyl phosphatidyl choline and cholesterol for dipalmitoyl phosphatidyl choline:cholesterol ratios of 24:1, 47:3, 11.5:1, 8:1, 7:1, 4:1, 3:1, 2:1, and 1:1, using long molecular dynamics runs and interspersed configurational bias Monte Carlo to ensure equilibration and enhance sampling. In all cases with cholesterol concentrations above 12.5% the area per molecule of the heterogeneous membrane varied linearly with cholesterol fraction. By extrapolation to pure cholesterol, we find the cross-sectional area of cholesterol in these mixtures is approximately 22.3 A(2). From the slope of the area/molecule relationship, we also find that the phospholipid in these mixtures is in a liquid ordered state with an average cross-sectional area per lipid of 50.7 A(2), slightly above the molecular area of a pure phospholipid gel. For lower concentrations of cholesterol, the molecular area rises above the straight line, indicating the "melting" of at least some of the phospholipid into a fluid state. Analysis of the lateral distribution of cholesterol molecules in the leaflets reveals peaks in radial distributions of cholesterols at multiples of approximately 5 A. These peaks grow in size as the simulation progresses, suggesting a tendency for small subunits of one lipid plus one cholesterol, hydrogen bonded together, to act as one composite particle, and perhaps to aggregate with other composites. Our results are consistent with experimentally observed effects of cholesterol, including the condensation effect of cholesterol in phospholipid monolayers and the tendency of cholesterol-rich domains to form in cholesterol-lipid bilayers. We are continuing to analyze this tendency on longer timescales and for larger bilayer patches.