Lipid chains and cholesterol in model membranes: a Monte Carlo Study

The Monte Carlo method has been employed to study the equilibrium properties of a planar array of hydrocarbon chains interacting with a cholesterol molecule. The chains are arranged to model one monolayer of a lipid bilayer and within this monolayer are allowed to move laterally and change conformations by gauche rotations. In the simulation cell there are 90 lipid chains and a single cholesterol molecule. Periodic boundary conditions are imposed upon the cell. The primary results of the calculations are order parameter profiles for the C-C bonds. These are calculated for (i) all chains, (ii) the 6 chains which are nearest neighbors to the cholesterol, and (iii) the 12 chains which are next-nearest neighbors to the cholesterol. Calculations are carried out for C-14, C-16, and C-18 chains. The results show that cholesterol strongly affects the upper portions of the chains, leaving them less able to change conformations. For C-16 and C-18 chains, the chain termini of the cholesterol neighbors are more disordered than the bulk chain termini. The magnitude of the effect depends strongly on the chain length. The results suggest that the changes in the lipid phase transition caused by cholesterol are a consequence of each cholesterol hindering the rotameric freedom of five to seven lipid chains.