Atomic-scale Molecular Dynamics Simulations Research Articles

Molecular-Scale Events in Hyperthermal Deposition of Organic Semiconductors Implicated from Experiment and Molecular Simulation

We have examined the dynamics of the adsorption of pentacene molecules on SiO2 and on a pentacene surface at room temperature using a combination of supersonic molecular beam techniques, atomic force microscopy, and atomic-scale molecular dynamics simulations. For incident kinetic energies, Ei = 1.5−6.7 eV, and angles of incidence, θi = 0–75°, we find a significant change in the dynamics of growth of pentacene on SiO2 compared to that on pentacene itself: at incident energies at and above 2 eV, the rate of growth accelerates, moving from the sub-monolayer to the multilayer regime. Molecular-scale insight from simulation attributes a significant contribution to this effect to molecular insertion events into the topmost layer of the pentacene thin film. Insertion can occur on terrace sites at incident kinetic energies as low as 1 eV and is more likely at step edges than on terraces. We discuss possible implications of this mechanism on the properties of deposited thin films.

Role of Cardiolipins in the Inner Mitochondrial Membrane: Insight Gained through Atom-Scale Simulations

Mitochondrial membranes are unique in many ways. Unlike other cellular membranes, they are comprised of two membranes instead of just one, and cardiolipins, one of the abundant lipid species in mitochondrial membranes, are not found in significant amounts elsewhere in the cell. Among other aspects, the exceptional nature of cardiolipins is characterized by their small charged head group connected to typically four hydrocarbon chains. In this work, we present atomic-scale molecular dynamics simulations of the inner mitochondrial membrane modeled as a mixture of cardiolipins (CLs), phosphatidylcholines (PCs), and phosphatidylethanolamines (PEs). For comparison, we also consider pure one-component bilayers and mixed PC-PE, PC-CL, and PE-CL membranes. We find that the influence of CLs on membrane properties depends strongly on membrane composition. This is highlighted by studies of the stability of CL-containing membranes, which indicate that the interactions of CL in ternary lipid bilayers cannot be deduced from the corresponding ones in binary membranes. Moreover, while the membrane properties in the hydrocarbon region are only weakly affected by CLs, the changes at the membrane-water interface turn out to be prominent. The effects at the interface are most evident in membrane properties related to hydrogen bonding and the binding phenomena associated with electrostatic interactions.

Interaction of Ethanol with Biological Membranes: The Formation of Non-bilayer Structures within the Membrane Interior and their Significance

To gain a better understanding of how ethanol affects biological membranes, we have performed a series of atomic-scale molecular dynamics simulations of phospholipid membranes in aqueous solution with ethanol, whose concentration was varied from 2.5 to 30 mol % (lipid-free basis). At concentrations below the threshold value of approximately 12 mol % (30.5 v/v %) ethanol induces expansion of the membrane, accompanied by a drop in the membrane thickness as well as disordering and enhanced interdigitation of lipid acyl chains. These changes become more pronounced with increase in ethanol concentration, but the bilayer structure of the membrane is maintained. Above the threshold concentration the appearance of multiple transient defects in the lipid/water interface eventually gives rise to desorption and assembly of some of the lipids into non-bilayer structures within the membrane interior. These structures, being small and irregular, resemble inverted micelles and have a long-lived character. Furthermore, formation of the non-bilayer structures is accompanied by mixing of lipids that belong to the opposite membrane leaflets, thereby leading to irreversible changes in the membrane structure. Remarkably, this observation of the formation of non-bilayer structures within the membrane interior, being in good agreement with experimental data, is found to be robust with respect to both the simulation conditions (the system size and the presence of salt) and the type of lipids (phosphatidylcholine and phosphatidylethanolamine). We discuss the significance of these non-bilayer structures in relation with the well-known ability of ethanol to promote membrane hemifusion as well as with the possible role of the micelle-like structures as a delivery system for polar solutes and ions. The ethanol-induced "damage" to the bilayer structure also suggests that strong alcoholic beverages ( approximately 40 v/v %) might be potentially hazardous to the epithelial tissues of the human body (such as lips, mouth, throat, gullet, and stomach) that come in direct contact with high-concentrations of ethanol.

Predicted Changes to Alkali Concentration Adjacent to Glass–Crystal Interfaces

Atomic scale molecular dynamics simulation is used to investigate how alkali concentrations in glass adjacent to interfaces with crystal surfaces, might differ from those within the interior of a bulk glass. A sodium, lithium alkali-borosilicate glass and its interfaces with the (1 0 0) and (1 1 0) surfaces of MgO, CaO, and SrO are considered. These six interfaces were formed by melt-quenching the glass in the presence of the crystal surface, providing glass alkali ions with the opportunity to either migrate towards or away from the crystal surface. The resulting alkali contents adjacent to the interface are found to differ depending on crystal composition and surface orientation.

Membrane Potential and Electrostatics of Phospholipid Bilayers with Asymmetric Transmembrane Distribution of Anionic Lipids

It is well-established that native plasma membranes are characterized by an asymmetric distribution of charged (anionic) lipids across the membrane. To clarify how the asymmetry can affect membrane electrostatics, we have performed extensive atomic-scale molecular dynamics simulations of asymmetric lipid membranes composed of zwitterionic (phosphatidylcholine (PC) or phosphatidylethanolamine (PE)) and anionic (phosphatidylserine (PS)) leaflets. It turns out that the asymmetry in transmembrane distribution of anionic lipids gives rise to a nonzero potential difference between the two sides of the membrane. This potential arises from the difference in surface charges of the two leaflets. The magnitude of the intrinsic membrane potential was found to be 238 mV and 198 mV for PS/PC and PS/PE membranes, respectively. Remarkably, this potential is of the same sign as the membrane potential in cells. Our findings, being in reasonable agreement with available experimental data, lend support to the idea that the transmembrane lipid asymmetry typical of most living cells contributes to the membrane potential.

Influence of Ethanol on Lipid Membranes: From Lateral Pressure Profiles to Dynamics and Partitioning

We have combined experiments with atomic-scale molecular dynamics simulations to consider the influence of ethanol on a variety of lipid membrane properties. We first employed isothermal titration calorimetry together with the solvent-null method to study the partitioning of ethanol molecules into saturated and unsaturated membrane systems. The results show that ethanol partitioning is considerably more favorable in unsaturated bilayers, which are characterized by their more disordered nature compared to their saturated counterparts. Simulation studies at varying ethanol concentrations propose that the partitioning of ethanol depends on its concentration, implying that the partitioning is a nonideal process. To gain further insight into the permeation of alcohols and their influence on lipid dynamics, we also employed molecular dynamics simulations to quantify kinetic events associated with the permeation of alcohols across a membrane, and to characterize the rotational and lateral diffusion of lipids and alcohols in these systems. The simulation results are in agreement with available experimental data and further show that alcohols have a small but non-vanishing effect on the dynamics of lipids in a membrane. The influence of ethanol on the lateral pressure profile of a lipid bilayer is found to be prominent: ethanol reduces the tension at the membrane-water interface and reduces the peaks in the lateral pressure profile close to the membrane-water interface. The changes in the lateral pressure profile are several hundred atmospheres. This supports the hypothesis that anesthetics may act by changing the lateral pressure profile exerted on proteins embedded in membranes.

Significance of Cholesterol Methyl Groups

Cholesterol is an indispensable molecule in mammalian cell membranes. To truly understand its role in the functions of membranes, it is essential to unravel cholesterol's structure-function relationship determined by underlying molecular interactions. For this purpose, we elaborate on this issue by considering the previously proposed idea that cholesterol's effects on a number of physical properties of membranes have been optimized during the evolution by removal of its excess methyl groups from the alpha-face of cholesterol, thus "smoothening" the structure. Consequently, the methyl groups still attached to cholesterol are one of the most intriguing structural features of the molecule. An obvious question arises: Why do these methyl groups still exist, and could cholesterol properties be further optimized by their removal? Because of the nature of the biosynthetic pathways of cholesterol, and the evidence of decreased interactions between sterols and lipid acyl chains when methyl groups are present, it seems plausible that removal of the methyl groups might indeed lead to stronger ordering and condensing effects of the cholesterol molecule. Atomic-scale molecular dynamics simulations of numerous modified sterols embedded in saturated lipid bilayers demonstrate, however, that the issue is more subtle. The analysis reveals a complex interplay between the lipid acyl chains and the structural details of cholesterol. Changes in cholesterol structure typically do not improve its performance in terms of promoting membrane order. This view is substantiated by a detailed analysis of the simulation data. In particular, it highlights the importance of the methyl group C18 for cholesterol properties. The C18 group resides between the third and fourth ring of cholesterol on its "rough" beta-side, and the results provide compelling evidence that C18 is crucial for the proper orientation of the sterol. More generally, the data provide insight into the role of the methyl groups of cholesterol.

Effect of NaCl and KCl on Phosphatidylcholine and Phosphatidylethanolamine Lipid Membranes: Insight from Atomic-Scale Simulations for Understanding Salt-Induced Effects in the Plasma Membrane

To gain a better understanding of how monovalent salt under physiological conditions affects plasma membranes, we have performed 200 ns atomic-scale molecular dynamics simulations of phosphatidylcholine (PC) and phosphatidylethanolamine (PE) lipid bilayers. These two systems provide representative models for the outer and inner leaflets of the plasma membrane, respectively. The implications of cation-lipid interactions in these lipid systems have been considered in two different aqueous salt solutions, namely NaCl and KCl, and the sensitivity of the results on the details of interactions used for ions is determined by repeating the simulations with two distinctly different force fields. We demonstrate that the main effect of monovalent salt on a phospholipid membrane is determined by cations binding to the carbonyl region of a membrane, while chloride anions mostly stay in the water phase. It turns out that the strength and character of the cation-lipid interactions are quite different for different types of lipids and cations. PC membranes and Na+ ions demonstrate strongest interactions, leading to notable membrane compression. This finding was confirmed by both force fields (Gromacs and Charmm) employed for the ions. The binding of potassium ions to PC membranes (and the overall effect of KCl), in turn, was found to be much weaker mainly due to the larger size of a K+ ion compared to Na+. Furthermore, the effect of KCl on PC membranes was found to be force-field sensitive: The binding of a potassium ion was not observed at all in simulations performed with the Gromacs force-field, which seems to exaggerate the size of a K+ ion. As far as PE lipid bilayers are concerned, they are found to be influenced by monovalent salt to a significantly lesser extent compared to PC bilayers, which is a direct consequence of the ability of PE lipids to form both intra- and intermolecular hydrogen bonds and hence to adopt a more densely packed bilayer structure. Whereas for NaCl we observed weak binding of Na+ cations to the PE lipid-water interface, in the case of KCl we witnessed almost complete lack of cation binding. Overall, our findings indicate that monovalent salt ions affect lipids in the inner and outer leaflets of plasma cell membranes in substantially different ways.

Lateral Diffusion in Lipid Membranes through Collective Flows

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.

Control of Calcium Oxalate Crystal Growth by Face-Specific Adsorption of an Osteopontin Phosphopeptide

Mineral-associated proteins have been proposed to regulate many aspects of biomineralization, including the location, type, orientation, shape, and texture of crystals. To understand how proteins achieve this exquisite level of control, we are studying the interaction between the phosphoprotein osteopontin (OPN) and the biomineral calcium oxalate monohydrate (COM). In the present study, we have synthesized peptides corresponding to amino acids 220-235 of rat bone OPN (pSHEpSTEQSDAIDpSAEK), one of several highly phosphorylated, aspartic-, and glutamic acid-rich sequences found in the protein. To investigate the role of phosphorylation in interaction with crystals, peptides containing no (P0), one (P1), or all three (P3) phosphates were prepared. Using a novel combination of confocal microscopy and scanning electron microscopy, we show that these peptides adsorb preferentially to {100} faces of COM and inhibit growth of these faces in a phosphorylation-dependent manner. To characterize the mechanism of adsorption of OPN peptides to COM, we have performed the first atomic-scale molecular-dynamics simulation of a protein-crystal interaction. P3 adsorbs to the {100} face much more rapidly than P1, which in turn adsorbs more rapidly than P0. In all cases, aspartic and glutamic acid, not phosphoserine, are the amino acids in closest contact with the crystal surface. These studies have identified a COM face-specific adsorption motif in OPN and delineated separate roles for carboxylate and phosphate groups in inhibition of crystal growth by mineral-associated phosphoproteins. We propose that the formation of close-range, stable, and face-specific interactions is a key factor in the ability of phosphoproteins to regulate biomineralization processes.

Ion Transport through Chemically Induced Pores in Protein-Free Phospholipid Membranes

We address the possibility of being able to induce the trafficking of salt ions and other solutes across cell membranes without the use of specific protein-based transporters or pumps. On the basis of realistic atomic-scale molecular dynamics simulations, we demonstrate that transmembrane ionic leakage can be initiated by chemical means, in this instance through addition of dimethyl sulfoxide (DMSO), a solvent widely used in cell biology. Our results provide compelling evidence that the small amphiphilic solute DMSO is able to induce transient defects (water pores) in membranes and to promote a subsequent diffusive pore-mediated transport of salt ions. The findings are consistent with available experimental data and offer a molecular-level explanation for the experimentally observed activities of DMSO solvent as an efficient penetration enhancer and a cryoprotectant, as well as an analgesic. Our findings suggest that transient pore formation by chemical means could emerge as an important general principle for therapeutics.

Effect of Double Bond Position on Lipid Bilayer Properties: Insight through Atomistic Simulations

Phosphatidylcholines (PCs) are among the most common phospholipids in plasma membranes. Their structural and dynamic properties are known to be strongly affected by unsaturation of lipid hydrocarbon chains, yet the role of the exact positions of the double bonds is poorly understood. In this work, we shed light on this matter through atomic-scale molecular dynamics simulations of eight different one-component lipid bilayers comprised of PCs with 18 carbons in their acyl chains. By introducing a single double bond in each acyl chain and varying its position in a systematic manner, we elucidate the effects of a double bond on various membrane properties. Studies in the fluid phase show that a number of membrane properties depend on the double bond position. In particular, when the double bond in an acyl chain is located close to the membrane-water interface, the area per lipid is considerably larger than that found for a saturated lipid. Further, when the double bond is shifted from the interfacial region toward membrane center, the area per lipid is observed to increase and have a maximum when the double bond is in the middle of the chain. Beyond this point, the surface area decreases systematically as the double bond approaches membrane center. These changes in area per lipid are accompanied by corresponding changes in membrane thickness and ordering of the chains. Further changes are observed in the tilt angles of the chains, membrane hydration together with changes in the number of gauche conformations, and direct head group interactions. All of these effects can be associated with changes in acyl chain conformations and local effects of the double bond on the packing of the surrounding atoms.

Modulating the Structure and Properties of Cell Membranes: The Molecular Mechanism of Action of Dimethyl Sulfoxide

Dimethyl sulfoxide (DMSO) is a small amphiphilic molecule which is widely employed in cell biology as an effective penetration enhancer, cell fusogen, and cryoprotectant. Despite the vast number of experimental studies, the molecular basis of its action on lipid membranes is still obscure. A recent simulation study employing coarse-grained models has suggested that DMSO induces pores in the membrane (Notman, R.; Noro, M.; O'Malley, B.; Anwar, J. J. Am. Chem. Soc. 2006, 128, 13982-13983). We report here the molecular mechanism for DMSO's interaction with phospholipid membranes ascertained from atomic-scale molecular dynamics simulations. DMSO is observed to exhibit three distinct modes of action, each over a different concentration range. At low concentrations, DMSO induces membrane thinning and increases fluidity of the membrane's hydrophobic core. At higher concentrations, DMSO induces transient water pores into the membrane. At still higher concentrations, individual lipid molecules are desorbed from the membrane followed by disintegration of the bilayer structure. The study provides further evidence that a key aspect of DMSO's mechanism of action is pore formation, which explains the significant enhancement in permeability of membranes to hydrophilic molecules by DMSO as well as DMSO's cryoprotectant activity. The reduction in the rigidity and the general disruption of the membrane induced by DMSO are considered to be prerequisites for membrane fusion processes. The findings also indicate that the choice of DMSO concentration for a given application is critical, as the concentration defines the specific mode of the solvent's action. Knowledge of the distinct modes of action of DMSO and associated concentration dependency should enable optimization of current application protocols on a rational basis and also promote new applications for DMSO.

What Happens if Cholesterol Is Made Smoother: Importance of Methyl Substituents in Cholesterol Ring Structure on Phosphatidylcholine–Sterol Interaction

Although sterols constitute one of the most important molecular species in cells, the reasons for their structure-function relationships in lipid membranes are not well understood. The main objective of this work is to elucidate the recently suggested possibility that the ordering and condensing effects of sterols on phospholipid membranes are related to the smoothness of a sterol. We focus on cholesterol, which has two methyl groups attached to its β-face, and compare its properties to those of demethylated cholesterol (Dchol), from which the two methyl groups have been removed. Atomic-scale molecular dynamics simulations of lipid membranes comprised of saturated lipids and sterols, either cholesterol or Dchol, provide compelling evidence that despite its smoother structure, the ordering and condensing effects of Dchol are less effective than those of cholesterol. The ordering capability of both cholesterol and Dchol is highly asymmetric with respect to their ring structure, but whereas cholesterol favors the α-face, Dchol favors the β-face. The origin and implications of this difference are analyzed in detail. The picture that emerges from this study supports a view that the two methyl groups at the steroid ring system of cholesterol play an important role in cholesterol-lipid interactions by reducing sterol tilt in the bilayer and hence allowing for an optimal orientation for cholesterol.

Lipid Transmembrane Asymmetry and Intrinsic Membrane Potential: Two Sides of the Same Coin

To better understand the origin of intrinsic membrane potential in cells, we have performed atomic-scale molecular dynamics simulations of an asymmetric lipid membrane comprised of zwitterionic (neutral) phosphatidylcholine (PC) and phosphatidylethanolamine (PE) leaflets in the absence of salt. It turns out that the asymmetry in distribution of zwitterionic lipids across the membrane gives rise to a nonzero intrinsic potential difference of about 100 mV between the two sides of the membrane. This potential arises from the difference in dipole moments of the two leaflets of the asymmetric PC/PE membrane. Our findings suggest that the transmembrane lipid asymmetry typical of most living cells can contribute to the membrane potential.

Free Pyrene Probes in Gel and Fluid Membranes: Perspective through Atomistic Simulations

We consider the properties of free pyrene probes inside gel- and fluidlike phospholipid membranes and unravel their influence on membrane properties. For this purpose, we employ atomic-scale molecular dynamics simulations at several temperatures for varying pyrene concentrations. Molecular dynamics simulations show that free pyrene molecules prefer to be located in the hydrophobic acyl chain region close to the glycerol group of lipid molecules. Their orientation is shown to depend on the phase of the membrane. In the fluid phase, pyrenes favor orientations where they are standing upright in parallel to the membrane normal, while, in the gel phase, the orientation is affected by the tilt of lipid acyl chains. Pyrenes are found to locally perturb membrane structure, while the nature of perturbations in the gel and fluid phases is completely different. In the gel phase, pyrenes break the local packing of lipids and decrease the ordering of lipid acyl chains around them, while, in the fluid phase, pyrenes increase the ordering of nearby acyl chains, thus having an opposite effect. Interestingly, this proposes a similarity to effects induced by cholesterol on structural membrane properties above and below the gel-fluid transition temperature. Further studies express a view that the orientational ordering of pyrene is not a particularly good measure of the acyl chain ordering of lipids. While pyrene ordering provides the correct qualitative behavior of acyl chain ordering in the fluid phase, its capability to predict the correct temperature dependence is limited.

Ion Leakage through Transient Water Pores in Protein-Free Lipid Membranes Driven by Transmembrane Ionic Charge Imbalance

We have employed atomic-scale molecular dynamics simulations to address ion leakage through transient water pores in protein-free phospholipid membranes. Our results for phospholipid membranes in aqueous solution with NaCl and KCl salts show that the formation of transient water pores and the consequent ion leakage can be induced and be driven by a transmembrane ionic charge imbalance, an inherent feature in living cells. These processes take place if the gradient is large enough to develop a sufficiently significant potential difference across the membrane. The transport of cations and anions through the water pores is then seen; it discharges the transmembrane potential, considerably reduces the size of a water pore, and makes the water pore metastable, leading eventually to its sealing. The ion transport is found to be sensitive to the type of ions. It turns out that Na + and Cl − ions leak through a membrane at approximately the same ratio despite the fact that Na + ions are expected to experience a lower potential barrier for the permeation through the pore. This is because of strong interactions of sodium ions with the carbonyl region of a phospholipid membrane as well as with lipid headgroups forming pore “walls,” considerably slowing down the permeation of sodium ions. In contrast, we observed a pronounced selectivity of a phospholipid membrane to the permeation of potassium ions as compared to chloride ions: Potassium ions, being larger than sodium ions, interact only weakly with phospholipid headgroups, so that these interactions are not able to compensate for a large difference in free-energy barriers for permeation of K + and Cl − ions. These findings are found to be robust to a choice of force-field parameters for ions (tested by Gromacs and Charmm force-fields for ions). What is more, a potassium ion is found to be able to permeate a membrane along an alternate, “water-defect-mediated” pathway without actual formation of a pore. The “water-defect-mediated” leakage involves formation of a single water defect only and is found to be at least one order of magnitude faster than the pore-mediated ion leakage.

Coarse-grained model for phospholipid/cholesterol bilayer employing inverse Monte Carlo with thermodynamic constraints

The authors introduce a coarse-grained (CG) model for a lipid membrane comprised of phospholipids and cholesterol at different molar concentrations, which allows them to study systems that are approximately 100 nm in linear size. The systems are studied in the fluid phase above the main transition temperature. The effective interactions for the CG model are extracted from atomic-scale molecular dynamics simulations using the inverse Monte Carlo (IMC) technique, an approach similar to the one the authors used earlier to construct another CG bilayer model [T. Murtola et al., J. Chem. Phys. 121, 9156 (2004)]. Here, the authors improve their original CG model by employing a more accurate description of the molecular structure for the phospholipid molecules. Further, they include a thermodynamic constraint in the IMC procedure to yield area compressibilities in line with experimental data. The more realistic description of the molecular structure of phospholipids and a more accurate representation of the interaction between cholesterols and phospholipid tails are shown to improve the behavior of the model significantly. In particular, the new model predicts the formation of denser transient regions in a pure phospholipid system, a finding that the authors have verified through large scale atomistic simulations. They also find that the model predicts the formation of cholesterol-rich and cholesterol-poor domains at intermediate cholesterol concentrations, in agreement with the original model and the experimental phase diagram. However, the domains observed here are much more distinct compared to the previous model. Finally, the authors also explore the limitations of the model, discussing its advantages and disadvantages.

Conformational analysis of lipid molecules by self-organizing maps

The authors have studied the use of the self-organizing map (SOM) in the analysis of lipid conformations produced by atomic-scale molecular dynamics simulations. First, focusing on the methodological aspects, they have systematically studied how the SOM can be employed in the analysis of lipid conformations in a controlled and reliable fashion. For this purpose, they have used a previously reported 50 ns atomistic molecular dynamics simulation of a 1-palmitoyl-2-linoeayl-sn-glycero-3-phosphatidylcholine (PLPC) lipid bilayer and analyzed separately the conformations of the headgroup and the glycerol regions, as well as the diunsaturated fatty acid chain. They have elucidated the effect of training parameters on the quality of the results, as well as the effect of the size of the SOM. It turns out that the main conformational states of each region in the molecule are easily distinguished together with a variety of other typical structural features. As a second topic, the authors applied the SOM to the PLPC data to demonstrate how it can be used in the analysis that goes beyond the standard methods commonly used to study the structure and dynamics of lipid membranes. Overall, the results suggest that the SOM method provides a relatively simple and robust tool for quickly gaining a qualitative understanding of the most important features of the conformations of the system, without a priori knowledge. It seems plausible that the insight given by the SOM could be applied to a variety of biomolecular systems and the design of coarse-grained models for these systems.

Transient Ordered Domains in Single-Component Phospholipid Bilayers

We report evidence of dense, ordered nanodomains in single-component fluid lipid bilayers. Our atomic-scale molecular dynamics simulations suggest that the area available to a lipid acyl chain exhibits large fluctuations, resulting in denser and sparser domains. The sizes of the dense domains can be up to approximately 10 nm, and their lifetimes are of the order of approximately 10 ns. In addition, our simulations suggest that domains of lipids with highly ordered acyl chains form predominantly within the dense regions, their sizes ranging from a few chains up to a few nanometers, and with lifetimes between approximately 10 ps-10 ns. These observations shed light on the origin of experimentally observed fluctuations, as well as on the mechanisms of phase transitions in lipid membranes.