Without any exaggeration, cholesterol is one of the most important lipid species in eukaryotic cells. Its effects on cellular membranes and functions range from purely mechanistic to complex metabolic ones, besides which it is also a precursor of the sex hormones (steroids) and several vitamins. In this review, we discuss the biophysical effects of cholesterol on the lipid bilayer, in particular the ordering and condensing effects, concentrating on the molecular level or inter-atomic interactions perspective, starting from two-component systems and proceeding to many-component ones e.g., modeling lipid rafts. Particular attention is paid to the roles of the methyl groups in the cholesterol ring system, and their possible biological function. Although our main research methodology is computer modeling, in this review we make extensive comparisons between experiments and different modeling approaches.

Phosphatidylethanolamine (PE) and phosphatidylglycerol (PG) are the main lipid components of the inner bacterial membrane. A computer model for such a membrane was built of palmitoyloleoyl PE (POPE) and palmitoyloleoyl PG (POPG) in the proportion 3:1, and sodium ions (Na+) to neutralize the net negative charge on each POPG (POPE-POPG bilayer). The bilayer was simulated for 25 ns. A final 10-ns trajectory fragment was used for analyses. In the bilayer interfacial region, POPEs and POPGs interact readily with one another via intermolecular hydrogen (H) bonds and water bridges. POPE is the main H-bond donor in either PEPE or PEPG H-bonds; PGPG H-bonds are rarely formed. Almost all POPEs are H-bonded and/or water bridged to either POPE or POPG but PE-PG links are favored. In effect, the atom packing in the near-the-interface regions of the bilayer core is tight. Na+ does not bind readily to lipids, and interlipid links via Na+ are not numerous. Although POPG and POPE comprise one bilayer, their bilayer properties differ. The average surface area per POPG is larger and the average vertical location of the POPG phosphate group is lower than those of POPE. Also, the alkyl chains of POPG are more ordered and less densely packed than the POPE chains. The main conclusion of this study is that in the PE-PG bilayer PE interacts more strongly with PG than with PE. This is a likely molecular-level event behind a regulating mechanism developed by the bacteria to control its membrane permeability and stability consisting in changes of the relative PG/PE concentration in the membrane.

Hydrogen (H-) bonding between water and phosphatidylcholine was studied using a molecular dynamics simulation of a hydrated phosphatidylcholine bilayer membrane in the liquid crystalline phase. A membrane in the liquid-crystalline phase composed of 72 l-α-dimyristoylphosphatidylcholine (DMPC) and 1622 water molecules was generated, starting from the crystal structure of DMPC. At the beginning of the equilibration process, the temperature of the system was raised to 550 K for 20 ps, which was effective in breaking the initial crystalline structure. The thermodynamic and structural parameters became stable after the equilibration period of 1100 ps, and the trajectory of the system obtained during the following 500 ps agreed well with most of the published experimental data. Each DMPC molecule forms 5.3 H-bonds with water, while only 4.5 water molecules are H-bonded to DMPC. The primary targets of water for the formation of H-bonds are the non-ester phosphate oxygens (4.0 H-bonds) and the carbonyl oxygens (∼1.0 H-bonds). Of DMPC's H-bonds, 1.7 are formed with water molecules that are simultaneously H-bonded to two different DMPC oxygens (bridging water). In effect, approximately 70% of the DMPC molecules are linked by water molecules and form clusters of two to seven DMPC molecules. Approximately 70% of the intermolecular water bridges are formed between non-ester phosphate oxygens. The rest are formed between non-ester phosphate and carbonyl oxygens. About half of the intermolecular water bridges are involved in formation of multiple bridges, where two DMPC molecules are linked by more than one parallel bridge. These results suggest a possibility that water bridges are involved in reducing head group mobility and in stabilizing the membrane structure. Non-ester phosphate oxygen of DMPC makes one, two, or three H-bonds with water, but two H-bonds are formed most often (≈60%). In the case where two H-bonds are formed on non-ester phosphate or carbonyl oxygens, the average geometry of H-bonding is planar trigonal (in the case of water oxygen with two H-bonds, geometry is steric tetragonal). When oxygen atoms form three H-bonds, the geometry of H-bonding is steric tetragonal both for non-ester phosphate and water oxygens. On average, H-bonds make nearly right angles with each other when two or three water molecules are bound to the same DMPC oxygen, but the distribution of the angle is broad.