Size matters: asphaltenes with enlarged aromatic cores promote heat transfer in organic phase-change materials.

Recent experiments and atomistic computer simulations have shown that asphaltene byproducts of oil refineries can serve as thermal conductivity enhancers for organic phase-change materials such as paraffin and therefore have the potential to improve the performance of paraffin-based heat storage devices. In this work, we explore how the size of the polycyclic aromatic cores of asphaltenes affects the properties of paraffin-asphaltene systems by means of atomistic molecular dynamics simulations. We show that increasing the size of the asphaltene core from 7-8 aromatic rings to ∼20 rings drastically changes the aggregation behavior of asphaltenes. Instead of relatively small, compact aggregates formed by small-core asphaltene molecules, enlarged cores promote the formation of extended single-column structures stabilized in paraffin by asphaltene's aliphatic periphery. Chemical modification of the asphaltenes by removing the periphery leads to the formation of bundles of columns. In contrast to small-core molecules, asphaltenes with enlarged cores do not suppress paraffin crystallization even at high filler concentrations. Remarkably, asphaltenes with enlarged aromatic cores are able to increase the thermal conductivity of liquid paraffin to a greater extent compared to their small-core counterparts. This effect becomes even more pronounced for modified asphaltenes without the aliphatic side groups. Overall, our computational findings suggest that asphaltenes with enlarged aromatic cores can significantly improve the performance of heat storage devices based on organic phase change materials.

Mesoscale computer modeling of asphaltene aggregation in liquid paraffin.

Asphaltenes represent a novel class of carbon nanofillers that are of potential interest for many applications, including polymer nanocomposites, solar cells, and domestic heat storage devices. In this work, we developed a realistic coarse-grained Martini model that was refined against the thermodynamic data extracted from atomistic simulations. This allowed us to explore the aggregation behavior of thousands of asphaltene molecules in liquid paraffin on a microsecond time scale. Our computational findings show that native asphaltenes with aliphatic side groups form small clusters that are uniformly distributed in paraffin. The chemical modification of asphaltenes via cutting off their aliphatic periphery changes their aggregation behavior: modified asphaltenes form extended stacks whose size increases with asphaltene concentration. At a certain large concentration (44 mol. %), the stacks of modified asphaltenes partly overlap, leading to the formation of large, disordered super-aggregates. Importantly, the size of such super-aggregates increases with the simulation box due to phase separation in the paraffin-asphaltene system. The mobility of native asphaltenes is systematically lower than that of their modified counterparts since the aliphatic side groups mix with paraffin chains, slowing down the diffusion of native asphaltenes. We also show that diffusion coefficients of asphaltenes are not very sensitive to the system size: enlarging the simulation box results in some increase in diffusion coefficients, with the effect being less pronounced at high asphaltene concentrations. Overall, our findings provide valuable insight into the aggregation behavior of asphaltenes on spatial and time scales that are normally beyond the scales accessible for atomistic simulations.

Elucidating lipid conformations in the ripple phase: Machine learning reveals four lipid populations.

A new mixed radial-angular, three-particle correlation function method in combination with unsupervised machine learning was applied to examine the emergence of the ripple phase in dipalmitoylphosphatidylcholine (DPPC) lipid bilayers using data from atomistic molecular dynamics simulations of system sizes ranging from 128 to 4096 lipids. Based on the acyl tail conformations, the analysis revealed the presence of four distinct conformational populations of lipids in the ripple phases of the DPPC lipid bilayers. The expected gel-like (ordered; Lo) and fluid-like (disordered; Ld) lipids are found along with their splayed tail equivalents (Lo,s and Ld,s). These lipids differ, based on their gauche distribution and tail packing. The disordered (Ld) and disordered-splayed (Ld,s) lipids spatially cluster in the ripple in the groove side, that is, in an asymmetric manner across the bilayer leaflets. The ripple phase does not contain large numbers of Ld lipids; instead they only exist on the interface of the groove side of the undulation. The bulk of the groove side is a complex coexistence of Lo,Lo,s, and Ld,s lipids. The convex side of the undulation contains predominantly Lo lipids. Thus, the structure of the ripple phase is neither a simple coexistence of ordered and disordered lipids nor a coexistence of ordered interdigitating gel-like (Lo) and ordered-splayed (Lo,s) lipids, but instead a coexistence of an ordered phase and a complex mixed phase. Principal component analysis further confirmed the existence of the four lipid groups.

Cooling-Rate Computer Simulations for the Description of Crystallization of Organic Phase-Change Materials.

A molecular-level insight into phase transformations is in great demand for many molecular systems. It can be gained through computer simulations in which cooling is applied to a system at a constant rate. However, the impact of the cooling rate on the crystallization process is largely unknown. To this end, here we performed atomic-scale molecular dynamics simulations of organic phase-change materials (paraffins), in which the cooling rate was varied over four orders of magnitude. Our computational results clearly show that a certain threshold (1.2 × 1011 K/min) in the values of cooling rates exists. When cooling is slower than the threshold, the simulations qualitatively reproduce an experimentally observed abrupt change in the temperature dependence of the density, enthalpy, and thermal conductivity of paraffins upon crystallization. Beyond this threshold, when cooling is too fast, the paraffin's properties in simulations start to deviate considerably from experimental data: the faster the cooling, the larger part of the system is trapped in the supercooled liquid state. Thus, a proper choice of a cooling rate is of tremendous importance in computer simulations of organic phase-change materials, which are of great promise for use in domestic heat storage devices.

How to control interactions of cellulose-based biomaterials with skin: the role of acidity in the contact area.

Being able to control the interactions of biomaterials with living tissues and skin is highly desirable for many biomedical applications. This is particularly the case for cellulose-based materials which provide one of the most versatile platforms for tissue engineering due to their strength, biocompatibility and abundance. Achieving such control, however, requires detailed molecular-level knowledge of the dominant interaction mechanisms. Here, we employed both biased and unbiased atomic-scale molecular dynamics simulations to explore how cellulose crystals interact with model stratum corneum bilayers, ternary mixtures of ceramides, cholesterol, and free fatty acids. Our findings show that acidity in the contact area directly affects binding between cellulose and the stratum corneum bilayer: Protonation of free fatty acids in the bilayer promotes attractive cellulose-bilayer interactions. We identified two major factors that control the cellulose-skin interactions: (i) the electrostatic repulsion between a cellulose crystal and the charged (anionic due to deprotonated fatty acids) surface of a stratum corneum bilayer and (ii) the cellulose-stratum corneum hydrogen bonding. When less than half of the fatty acids in the bilayer are protonated, the first factor dominates and there is no binding to skin. At a larger degree of fatty acid protonation the cellulose-stratum corneum hydrogen bonding prevails yielding a tight binding. Remarkably, we found that ceramide molecules are the key component in hydrogen bonding with cellulose. Overall, our findings highlight the critical role of fatty acid protonation in biomaterial-stratum corneum interactions and can be used for optimizing the surface properties of cellulose-based materials aimed at biomedical applications such as wound dressings.

Correction: Toward realistic computer modeling of paraffin-based composite materials: critical assessment of atomic-scale models of paraffins.

[This corrects the article DOI: 10.1039/C9RA07325F.].

Controlled On-Off Switching of Tight-Binding Hydrogen Bonds between Model Cell Membranes and Acetylated Cellulose Surfaces.

Controlling interactions between cellulose-based materials and membranes of living cells is critical in medicine and biotechnology in, for example, wound dressing, tissue engineering, hemodialysis membranes, and drug transport. Cellulose acetylation is a widely used approach to tuning those interactions. Surprisingly, however, detailed interactions of acetylated cellulose and membranes have thus far not been characterized. Using atomistic molecular dynamics (MD) simulations, we show that the key to such control is hydrogen bonds: by tuning the number of hydrogen bonds between tissue (cell membranes) and cellulose, binding can be controlled in a precise manner. We demonstrate that the acetylation of each hydroxymethyl group reduces the free energy of cellulose-membrane binding by an order of magnitude as compared to that of pristine cellulose. Remarkably, this acetylation-induced weakening does not occur gradually and is characterized by a sharp threshold in the degree of substitution, beyond which the microscopic character of lipid-cellulose interactions changes drastically. When the degree of substitution does not exceed 0.125, the cellulose-lipid interactions are mainly driven by hydrogen bonding between cellulose's hydroxyl groups and phosphate groups of lipid molecules. This results in the tight binding of a cellulose crystal and a lipid bilayer. Larger degrees of substitution (here, 0.25 and 0.5) prevent hydrogen bonding, leading to rather weak and unstable cellulose-bilayer binding. In this case, the lipid-cellulose binding is controlled by the interactions of lipid choline groups with hydroxyl(hydroxymethyl) groups and carbonyl groups of acetyl moieties of acetylated cellulose.

Molecular-Level Insight into the Interactions of DNA/Polycation Complexes with Model Cell Membranes.

The interactions of DNA/polycation complexes (polyplexes) with cell membranes are crucial for understanding the molecular mechanisms behind polycation-mediated delivery of nucleic acid therapeutics into the target cells. In this study, we employed both biased and unbiased atomic-scale computer simulations to get an insight into such interactions. To this end, we considered complexes of DNA with linear polyethylenimine (PEI) with various polycation contents, ranging from an almost fully neutralized DNA to a highly overcharged polyplex. Our findings clearly show that the free energy gradually increases when a polyplex approaches the surface of a zwitterionic (neutral) phospholipid membrane from bulk water, implying the lack of attractive polyplex/membrane interactions. Remarkably, overcharging of DNA molecules by polycations enhances the repulsion between the polyplex and the zwitterionic lipid membrane. The observed repulsion is most likely driven by the dehydration of a polyplex upon its partitioning into the zwitterionic lipid membrane as well as by the loss of conformational entropy of PEI chains. We also demonstrate that cationic polymer chains are able to protect DNA from the dehydration as well as from contacts with lipid molecules. Interestingly, the absence of local minima in the free energy profiles does not exclude transient weak adsorption of a polyplex on the zwitterionic membrane surface. We show that such spontaneous adsorption can indeed be initiated by the interactions of loose polycation chains of the polyplex with polar head groups of lipids. Overall, our computational findings contribute considerably to the understanding of the initial stages in polycation-mediated DNA transfection. In particular, we demonstrate that a zwitterionic lipid bilayer represents an energetic barrier for polyplexes, so that a proper model of the cell membrane should account for the anionic surface charge of the membrane (e.g., due to the presence of proteoglycans).

The Devil Is in the Details: What Do We Really Track in Single-Particle Tracking Experiments of Diffusion in Biological Membranes?

Single-particle tracking (SPT) is an experimental technique that allows one to follow the dynamics of individual molecules in biological membranes with unprecedented precision. Given the importance of lipid and membrane protein diffusion in the formation of nanoscale functional complexes, it is critical to understand what exactly is measured in SPT experiments. To clarify this issue, we employed nanoscale computer simulations designed to match SPT experiments that exploit streptavidin-functionalized Au nanoparticles (AuNPs). The results show that lipid labeling interferes critically with the diffusion process; thus, the diffusion measured in SPT is a far more complex process than what has been assumed. It turns out that the influence of AuNP-based labels on the dynamics of probe lipids includes not only the AuNP-induced viscous drag that is the more significant the larger the NP but, more importantly, also the effects related to the interactions of the streptavidin linker with membrane lipids. Due to these effects, the probe lipid moves in a concerted manner as a complex with the linker protein and numerous unlabeled lipids, which can slow down the motion of the probe by almost an order of magnitude. Furthermore, our simulations show that nonlinker streptavidin tetramers on the AuNP surface are able to interact with the membrane lipids, which could potentially lead to multivalent labeling of the NPs by the probe lipids. Our results further demonstrate that in the submicrosecond time domain the motion of the probe lipid is uncorrelated with the motion of the AuNP, showing that there is a 1 µs limit for the temporal resolution of the SPT technique. However, this limit for the temporal resolution depends on the nanoparticle size and increases rapidly with growing AuNPs. Overall, the results provide a molecular-scale framework to accurately interpret SPT data and to design protocols that minimize label-induced artifacts.

Toward realistic computer modeling of paraffin-based composite materials: critical assessment of atomic-scale models of paraffins.

Paraffin-based composites represent a promising class of materials with numerous practical applications such as e.g. heat storage. Computer modeling of these complex multicomponent systems requires a proper theoretical description of both the n-alkane matrix and the non-alkane filler molecules. The latter can be modeled with the use of a state-of-the-art general-purpose force field such as GAFF, CHARMM, OPLS-AA and GROMOS, while the paraffin matrix is traditionally described in the frame of relatively old, alkane-specific force fields (TraPPE, NERD, and PYS). In this paper we link these two types of models and evaluate the performance of several general-purpose force fields in computer modeling of paraffin by their systematic comparison with earlier alkane-specific models as well as with experimental data. To this end, we have performed molecular dynamics simulations of n-eicosane bulk samples with the use of 10 different force fields: TraPPE, NERD, PYS, OPLS-UA, GROMOS, GAFF, GAFF2, OPLS-AA, L-OPLS-AA, and CHARMM36. For each force field we calculated several thermal, structural and dynamic characteristics of n-eicosane over a wide temperature range. Overall, our findings show that the general-purpose force fields such as CHARMM36, L-OPLS-AA and GAFF/GAFF2 are able to provide a realistic description of n-eicosane samples. While alkane-specific models outperform most general-purpose force fields as far as the temperature dependence of mass density, the coefficient of volumetric thermal expansion in the liquid state, and the crystallization temperature are concerned, L-OPLS-AA, CHARMM36 and GAFF2 force fields provide a better match with experiment for the shear viscosity and the diffusion coefficient in melt. Furthermore, we show that most general-purpose force fields are able to reproduce qualitatively the experimental triclinic crystal structure of n-eicosane at low temperatures.

Phospholipid-Cellulose Interactions: Insight from Atomistic Computer Simulations for Understanding the Impact of Cellulose-Based Materials on Plasma Membranes.

Cellulose is an important biocompatible and nontoxic polymer widely used in numerous biomedical applications. The impact of cellulose-based materials on cells and, more specifically, on plasma membranes that surround cells, however, remains poorly understood. To this end, here, we performed atomic-scale molecular dynamics simulations of phosphatidylcholine (PC) and phosphatidylethanolamine (PE) bilayers interacting with the surface of a cellulose crystal. Both biased umbrella sampling and unbiased simulations clearly show the existence of strong attractive interactions between phospholipids and cellulose: the free energy of the cellulose-bilayer binding was found to be -1.89 and -1.96 kJ/mol per cellulose dimer for PC and PE bilayers, respectively. Although the values are similar, there is a pronounced difference between PC and PE bilayers. The driving force in both cases is the formation of hydrogen bonds. There are two distinct types of hydrogen bonds: (1) between the lipid head groups and the hydroxyl (hydroxymethyl) groups of cellulose, and (2) lipid-water and cellulose-water bonds. The former is the dominant component for PE systems whereas the latter dominates in PC systems. This suggests that achieving controlled binding via new cellulose modifications must pay close attention to the lipid head groups involved. The observed attractive phospholipid-cellulose interactions have a significant effect on bilayer properties: a cellulose crystal induces noticeable structural perturbations on the bilayer leaflet next to the crystal. Given that such perturbations can be undesirable when it comes to the interactions of cellulose-based materials with cell membranes, our computational studies suggest that the impact of cellulose could be reduced through chemical modification of the cellulose surface which prevents cellulose-phospholipid hydrogen bonding.

Toward Understanding Liposome-Based siRNA Delivery Vectors: Atomic-Scale Insight into siRNA-Lipid Interactions.

Liposome carriers for delivering small interfering RNA (siRNA) into target cells are of tremendous importance because the siRNA-based therapy offers a completely new approach for treating a wide range of diseases, including cancer and viral infections. In this paper, we employ the state-of-the-art computer simulations to get an atomic-scale insight into the interactions of siRNA with zwitterionic (neutral) lipids. Our computational findings clearly demonstrate that siRNA does adsorb on the surface of a neutral lipid bilayer. The siRNA adsorption, being rather weak and unstable, is driven by attractive interactions of overhanging unpaired nucleotides with choline moieties of lipid molecules. It is the presence of the unpaired terminal nucleotides that underlies a drastic difference between siRNA and DNA; the latter is not able to bind to the zwitterionic lipid bilayer. We also show that adding divalent Ca ions leads to the formation of stable siRNA-lipid system complexes; these complexes are stabilized by Ca-mediated aggregates of siRNA and lipid molecules rather than by the overhanging siRNA nucleotides. Furthermore, the molecular mechanism of interactions between siRNA and the lipid bilayer in the presence of divalent cations seems to involve exchange of Ca ions between the outer mouth of the major groove of siRNA and the lipid/water interface. Overall, our findings contribute significantly to a deeper understanding of the structure and function of liposome carriers used for siRNA delivery and can be used as a theoretical basis for further development of siRNA-based therapeutics.

Molecular-Level Insight into the Interaction of Phospholipid Bilayers with Cellulose.

Molecular-level insight into the interactions of phospholipid molecules with cellulose is crucial for the development of novel cellulose-based materials for wound dressing. Here we employ the state-of-the-art computer simulations to unlock for the first time the molecular mechanisms behind such interactions. To this end, we performed a series of atomic-scale molecular dynamics simulations of phospholipid bilayers on a crystalline cellulose support at various hydration levels of the bilayer leaflets next to the cellulose surface. Our findings clearly demonstrate the existence of strong interactions between polar lipid head groups and the hydrophilic surface of a cellulose crystal. We identified two major types of interactions between phospholipid molecules and cellulose chains: (i) direct attractive interactions between lipid choline groups and oxygens of hydroxyl (hydroxymethyl) groups of cellulose and (ii) hydrogen bonding between phosphate groups of lipids and cellulose's hydroxymethyl/hydroxyl groups. When the hydration level of the interfacial bilayer/support region is low, these interactions lead to a pronounced asymmetry in the properties of the opposite bilayer leaflets. In particular, the mass density profiles of the proximal leaflets are split into two peaks and lipid head groups become more horizontally oriented with respect to the bilayer surface. Furthermore, the lateral mobility of lipids in the leaflets next to the cellulose surface is found to slow down considerably. Most of these cellulose-induced effects are likely due to hydrogen bonding between lipid phosphate groups and hydroxymethyl/hydroxyl groups of cellulose: the lipid phosphate groups are pulled toward the water/lipid interface due to the formation of hydrogen bonds. Overall, our findings shed light on the molecular details of the interactions between phospholipid bilayers and cellulose nanocrystals and can be used for identifying possible strategies for improving the properties of cellulose-based dressing materials via, e.g., chemical modification of their surface.

Glycosylation and Lipids Working in Concert Direct CD2 Ectodomain Orientation and Presentation.

Proteins embedded in the plasma membrane mediate interactions with the cell environment and play decisive roles in many signaling events. For cell-cell recognition molecules, it is highly likely that their structures and behavior have been optimized in ways that overcome the limitations of membrane tethering. In particular, the ligand binding regions of these proteins likely need to be maximally exposed. Here we show by means of atomistic simulations of membrane-bound CD2, a small cell adhesion receptor expressed by human T-cells and natural killer cells, that the presentation of its ectodomain is highly dependent on membrane lipids and receptor glycosylation acting in apparent unison. Detailed analysis shows that the underlying mechanism is based on electrostatic interactions complemented by steric interactions between glycans in the protein and the membrane surface. The findings are significant for understanding the factors that render membrane receptors accessible for binding and signaling.

What Can We Learn about Cholesterol's Transmembrane Distribution Based on Cholesterol-Induced Changes in Membrane Dipole Potential?

Cholesterol is abundant in the plasma membranes of animal cells and is known to regulate a variety of membrane properties. Despite decades of research, the transmembrane distribution of cholesterol is still a matter of debate. Here we consider this outstanding issue through atomistic simulations of asymmetric lipid membranes, whose composition is largely consistent with eukaryotic plasma membranes. We show that the membrane dipole potential changes in a cholesterol-dependent manner. Remarkably, moving cholesterol from the extracellular to the cytosolic leaflet increases the dipole potential on the cytosolic side, and vice versa. Biologically this implies that by altering the dipole potential, cholesterol can provide a driving force for cholesterol molecules to favor the cytosolic leaflet, in order to compensate for the intramembrane field that arises from the resting potential.

Adsorption of Synthetic Cationic Polymers on Model Phospholipid Membranes: Insight from Atomic-Scale Molecular Dynamics Simulations.

Although synthetic cationic polymers represent a promising class of effective antibacterial agents, the molecular mechanisms behind their antimicrobial activity remain poorly understood. To this end, we employ atomic-scale molecular dynamics simulations to explore adsorption of several linear cationic polymers of different chemical structure and protonation (polyallylamine (PAA), polyethylenimine (PEI), polyvinylamine (PVA), and poly-l-lysine (PLL)) on model bacterial membranes (4:1 mixture of zwitterionic phosphatidylethanolamine (PE) and anionic phosphatidylglycerol (PG) lipids). Overall, our findings show that binding of polycations to the anionic membrane surface effectively neutralizes its charge, leading to the reorientation of water molecules close to the lipid/water interface and to the partial release of counterions to the water phase. In certain cases, one has even an overcharging of the membrane, which was shown to be a cooperative effect of polymer charges and lipid counterions. Protonated amine groups of polycations are found to interact preferably with head groups of anionic lipids, giving rise to formation of hydrogen bonds and to a noticeable lateral immobilization of the lipids. While all the above findings are mostly defined by the overall charge of a polymer, we found that the polymer architecture also matters. In particular, PVA and PEI are able to accumulate anionic PG lipids on the membrane surface, leading to lipid segregation. In turn, PLL whose charge twice exceeds charges of PVA/PEI does not induce such lipid segregation due to its considerably less compact architecture and relatively long side chains. We also show that partitioning of a polycation into the lipid/water interface is an interplay between its protonation level (the overall charge) and hydrophobicity of the backbone. Therefore, a possible strategy in creating highly efficient antimicrobial polymeric agents could be in tuning these polycation's properties through proper combination of protonated and hydrophobic blocks.

Atomic-Scale Molecular Dynamics Simulations of DNA-Polycation Complexes: Two Distinct Binding Patterns.

Synthetic cationic polymers represent a promising class of delivery vectors for gene therapy. Here, we employ atomistic molecular dynamics simulations to gain insight into the structure and properties of complexes of DNA with four linear polycations: polyethylenimine (PEI), poly-l-lysine (PLL), polyvinylamine (PVA), and polyallylamine (PAA). These polycations differ in their polymer geometries, protonation states, and hydrophobicities of their backbone chains. Overall, our results demonstrate for the first time the existence of two distinct patterns of binding of DNA with polycations. For PEI, PLL, and PAA, the complex is stabilized by the electrostatic attraction between protonated amine groups of the polycation and phosphate groups of DNA. In contrast, PVA demonstrates an alternative binding pattern as it gets embedded into the DNA major groove. It is likely that both the polymer topology and affinity of the backbone chain of PVA to the DNA groove are responsible for such behavior. The differences in binding patterns can have important biomedical implications: embedding PVA into a DNA groove makes it less sensitive to changes in the aqueous environment (pH level, ionic strength, etc.) and could therefore hinder the intracellular release of genetic material from a delivery vector, leading to lower transfection activity.

Molecular mechanism of calcium-induced adsorption of DNA on zwitterionic phospholipid membranes.

Interaction of DNA with zwitterionic phospholipids is an important long-standing problem in the field of liposome-based gene delivery. Although it is well-established that divalent cations can promote formation of stable DNA-phospholipid complexes, the underlying molecular mechanism remains largely unknown. Here we employ computer simulations to gain atomistically resolved insight into the kinetics of calcium-induced adsorption of DNA on zwitterionic phosphatidylcholine membranes as well as into the structure and stability of the resulting complexes. Overall, our findings show that calcium ions play a dual role in DNA-phospholipid systems. First, binding of divalent cations to the lipid-water interface turns the surface of the zwitterionic membrane positively charged, promoting thereby the initial electrostatic attraction of a polyanionic DNA molecule. Second, we show that calcium ions are crucial for stabilizing the DNA-lipid membrane complex as they bridge together phosphate groups of DNA and lipid molecules. In contrast to previous hypotheses, we demonstrate that direct interactions between choline groups of phospholipids and DNA phosphates play only a rudimentary role as they are relatively short-lived and unstable: typical residence times for such interactions are 2 orders of magnitude smaller than those for Ca-mediated bridges between DNA and lipid phosphate groups. The results of our study can serve as a basis for a deeper understanding of molecular mechanisms behind noncovalent binding of DNA and DNA-based nanodevices to complex surfaces such as cell membranes.

Atomistic simulations of anionic Au144(SR)60 nanoparticles interacting with asymmetric model lipid membranes.

Experimental observations indicate that the interaction between nanoparticles and lipid membranes varies according to the nanoparticle charge and the chemical nature of their protecting side groups. We report atomistic simulations of an anionic Au nanoparticle (AuNP(-)) interacting with membranes whose lipid composition and transmembrane distribution are to a large extent consistent with real plasma membranes of eukaryotic cells. To this end, we use a model system which comprises two cellular compartments, extracellular and cytosolic, divided by two asymmetric lipid bilayers. The simulations clearly show that AuNP(-) attaches to the extracellular membrane surface within a few tens of nanoseconds, while it avoids contact with the membrane on the cytosolic side. This behavior stems from several factors. In essence, when the nanoparticle interacts with lipids in the extracellular compartment, it forms relatively weak contacts with the zwitterionic head groups (in particular choline) of the phosphatidylcholine lipids. Consequently, AuNP(-) does not immerse deeply in the leaflet, enabling, e.g., lateral diffusion of the nanoparticle along the surface. On the cytosolic side, AuNP(-) remains in the water phase due to Coulomb repulsion that arises from negatively charged phosphatidylserine lipids interacting with AuNP(-). A number of structural and dynamical features resulting from these basic phenomena are discussed. We close the article with a brief discussion of potential implications.

Electroporation of asymmetric phospholipid membranes.

As plasma membranes of animal cells are known to be asymmetric, the transmembrane lipid asymmetry, being essential for many membranes' properties and functions, should be properly accounted for in model membrane systems. In this paper, we employ atomic-scale molecular dynamics simulations to explore electroporation phenomena in asymmetric model membranes comprised of phosphatidylcholine (PC) and phosphatidylethanolamine (PE) lipid monolayers that mimic the outer and inner leaflets of plasma membranes, respectively. Our findings clearly demonstrate that the molecular mechanism of electroporation in asymmetric phospholipid membranes differs considerably from the picture observed for their single-component symmetric counterparts: The initial stages of electric-field-induced formation of a water-filled pore turn out to be asymmetric and occur mainly on the PC side of the PC/PE membrane. In particular, water molecules penetrate in the membrane interior mostly from the PC side, and the reorientation of lipid head groups, being crucial for stabilizing the hydrophilic pore, also takes place in the PC leaflet. In contrast, the PE lipid head groups do not enter the central region of the membrane until the water pore becomes rather large and partly stabilized by PC head groups. Such behavior implies that the PE leaflet is considerably more robust against an electric field most likely due to interlipid hydrogen bonding. We also show that an electric field induces asymmetric changes in the lateral pressure profile of PC/PE membranes, decreasing the cohesion between lipid molecules predominantly in the PC membrane leaflet. Overall, our simulations provide compelling evidence that the transmembrane lipid asymmetry can be essential for understanding electroporation phenomena in living cells.