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.

Anionic palmitoyloleoylphosphatidylglycerol (POPG) is one of the most abundant lipids in nature, yet its atomic-scale properties have not received significant attention. Here we report extensive 150-ns molecular dynamics simulations of a pure POPG lipid membrane with sodium counterions. It turns out that the average area per lipid of the POPG bilayer under physiological conditions is approximately 19% smaller than that of a bilayer built from its zwitterionic phosphatidylcholine analog, palmitoyloleoylphosphatidylcholine. This suggests that there are strong attractive interactions between anionic POPG lipids, which overcome the electrostatic repulsion between negative charges of PG headgroups. We demonstrate that interlipid counterion bridges and strong intra- and intermolecular hydrogen bonding play a key role in this seemingly counterintuitive behavior. In particular, the substantial strength and stability of ion-mediated binding between anionic lipid headgroups leads to complexation of PG molecules and ions and formation of large PG-ion clusters that act in a concerted manner. The ion-mediated binding seems to provide a possible molecular-level explanation for the low permeability of PG-containing bacterial membranes to organic solvents: highly polar interactions at the water/membrane interface are able to create a high free energy barrier for hydrophobic molecules such as benzene.

A molecular dynamics (MD) simulation of a fully hydrated, liquid-crystalline dimyristoylphosphatidylcholine (DMPC)-Chol bilayer membrane containing approximately 22 mol% Chol was carried out for 4.3 ns. The bilayer reached thermal equilibrium after 2.3 ns of MD simulation. A 2.0-ns trajectory generated during 2.3-4.3 ns of MD simulation was used for analyses to determine the effects of Chol on the membrane/water interfacial region. In this region, 70% of Chol molecules are linked to DMPC molecules via short-distance interactions, where the Chol hydroxyl group (OH-Chol) is 1) charge paired to methyl groups of the DMPC choline moiety ( approximately 34%), via the hydroxyl oxygen atom (Och); 2) water bridged to carbonyl ( approximately 19%) and nonester phosphate ( approximately 14%) oxygen atoms, via both Och and the hydroxyl hydrogen atom (Hch); and 3) directly hydrogen (H) bonded to carbonyl ( approximately 11%) and nonester phosphate ( approximately 5%) oxygen atoms, via Hch ( approximately 17% of DMPC-Chol links are multiple). DMPC's gamma-chain carbonyl oxygen atom is involved in 44% of water bridges and 51% of direct H bonds formed between DMPC and Chol. On average, a Chol molecule forms 0.9 links with DMPC molecules, while a DMPC molecule forms 2.2 and 0.3 links with DMPC and Chol molecules, respectively. OH-Chol makes hydrogen bonds with 1.1 water molecules, preferentially via Hch. The average number of water molecules H bonded to the DMPC headgroup is increased by 7% in the presence of Chol. These results indicate that inclusion of Chol decreases interlipid links and increases hydration in the polar region of the membrane.

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

There is no comprehensive model for the dynamics of cellular membranes. Even mechanisms of basic dynamic processes, such as lateral diffusion of lipids, are poorly understood. Our atomic-scale molecular dynamics simulations support a novel, concerted mechanism for lipid diffusion. We find that a lipid and its nearest neighbors move in unison, forming loosely defined clusters. What is more, the motions of lipids are correlated over tens of nanometers: the lateral displacements of lipids in a given monolayer produce striking two-dimensional flow patterns. These flow patterns should have wide implications, affecting, for example, the formation of membrane domains, protein functionality, and action of lipases and drugs on membranes.

The structural and dynamical properties of lipid membranes rich in phospholipids and cholesterol are known to be strongly affected by the unsaturation of lipid acyl chains. We show that not only unsaturation but also the position of a double bond has a pronounced effect on membrane properties. We consider how cholesterol interacts with phosphatidylcholines comprising two 18-carbon long monounsaturated acyl chains, where the position of the double bond is varied systematically along the acyl chains. Atomistic molecular dynamics simulations indicate that when the double bond is not in contact with the cholesterol ring, and especially with the C18 group on its rough beta-side, the membrane properties are closest to those of the saturated bilayer. However, any interaction between the double bond and the ring promotes membrane disorder and fluidity. Maximal disorder is found when the double bond is located in the middle of a lipid acyl chain, the case most commonly found in monounsaturated acyl chains of phospholipids. The results suggest a cholesterol-mediated lipid selection mechanism in eukaryotic cell membranes. With saturated lipids, cholesterol promotes the formation of highly ordered raft-like membrane domains, whereas domains rich in unsaturated lipids with a double bond in the middle remain highly fluid despite the presence of cholesterol.

No abstract