Abstract
Cholesterol molecules were put into a computer-modeled hydrated bilayer of dimyristoyl phosphatidyl choline molecules, and molecular dynamics simulations were run to characterize the effect of this important molecule on membrane structure and dynamics. The effect was judged by observing differences in order parameters, tilt angles, and the fraction of gauche bonds along the hydrocarbon chains between lipids adjacent to cholesterol molecules and comparing them with those further away. It was observed that cholesterol causes an increase in the fraction of trans dihedrals and motional ordering of chains close to the rigid steroid ring system with a decrease in the kink population. The hydrogen-bonding interactions between cholesterol and lipid molecules were determined from radial distribution calculations and showed the cholesterol hydroxyl groups either solvated by water, or forming hydrogen bond contacts with the oxygens of lipid carbonyl and phosphate groups. The dynamics and conformation of the cholesterol molecules were investigated and it was seen that they had a smaller tilt with respect to the bilayer normal than the lipid chains and furthermore that the hydrocarbon tail of the cholesterol was conformationally flexible.