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.

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.

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.

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.

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.

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.

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.

Thermal properties of bulk polyimides: insights from computer modeling versus experiment.

Due to the great importance for many industrial applications it is crucial from the point of view of theoretical description to reproduce thermal properties of thermoplastic polyimides as accurate as possible in order to establish "chemical structure-physical properties" relationships of new materials. In this paper we employ differential scanning calorimetry, dilatometry, and atomistic molecular dynamics (MD) simulations to explore whether the state-of-the-art computer modeling can serve as a precise tool for probing thermal properties of polyimides with highly polar groups. For this purpose the polyimide R-BAPS based on dianhydride 1,3-bis(3',4-dicarboxyphenoxy)benzene (dianhydride R) and diamine 4,4'-bis(4''-aminophenoxy)biphenyl sulphone) (diamine BAPS) was synthesized and extensively studied. Overall, our findings show that the widely used glass-transition temperature Tg evaluated from MD simulations should be employed with great caution for verification of the polyimide computational models against experimental data: in addition to the well-known impact of the cooling rate on the glass-transition temperature, correct definition of Tg requires cooling that starts from very high temperatures (no less than 800 K for considered polyimides) and accurate evaluation of the appropriate cooling rate, otherwise the errors in the measured values of Tg become undefined. In contrast to the glass-transition temperature, the volumetric coefficient of thermal expansion (CTE) does not depend on the cooling rate in the low-temperature domain (T < Tg) so that comparison of computational and experimental values of CTE provides a much safer way for proper validation of the theoretical model when electrostatic interactions are taken into account explicitly. Remarkably, this conclusion is most likely of generic nature: we show that it also holds for the commercial polyimide EXTEM, another polyimide with a similar chemical structure.

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.].

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.

Intrinsic potential of cell membranes: opposite effects of lipid transmembrane asymmetry and asymmetric salt ion distribution.

Using atomic-scale molecular dynamics simulations, we consider the intrinsic cell membrane potential that is found to originate from a subtle interplay between lipid transmembrane asymmetry and the asymmetric distribution of monovalent salt ions on the two sides of the cell membrane. It turns out that both the asymmetric distribution of phosphatidylcholine (PC) and phosphatidylethanolamine (PE) lipids across a membrane and the asymmetric distribution of NaCl and KCl induce nonzero drops in the transmembrane potential. However, these potential drops are opposite in sign. As the PC leaflet faces a NaCl saline solution and the PE leaflet is exposed to KCl, the outcome is that the effects of asymmetric lipid and salt ion distributions essentially cancel one another almost completely. Overall, our study highlights the complex nature of the intrinsic potential of cell membranes under physiological conditions.

Ion dynamics in cationic lipid bilayer systems in saline solutions.

Positively charged lipid bilayer systems are a promising class of nonviral vectors for safe and efficient gene and drug delivery. Detailed understanding of these systems is therefore not only of fundamental but also of practical biomedical interest. Here, we study bilayers comprising a binary mixture of cationic dimyristoyltrimethylammoniumpropane (DMTAP) and zwitterionic (neutral) dimyristoylphosphatidylcholine (DMPC) lipids. Using atomistic molecular dynamics simulations, we address the effects of bilayer composition (cationic to zwitterionic lipid fraction) and of NaCl electrolyte concentration on the dynamical properties of these cationic lipid bilayer systems. We find that, despite the fact that DMPCs form complexes via Na(+) ions that bind to the lipid carbonyl oxygens, NaCl concentration has a rather minute effect on lipid diffusion. We also find the dynamics of Cl(-) and Na(+) ions at the water-membrane interface to differ qualitatively. Cl(-) ions have well-defined characteristic residence times of nanosecond scale. In contrast, the binding of Na(+) ions to the carbonyl region appears to lack a characteristic time scale, as the residence time distributions displayed power-law features. As to lateral dynamics, the diffusion of Na(+) ions within the water-membrane interface consists of two qualitatively different modes of motion: very slow diffusion when ions are bound to DMPC, punctuated by fast rapid jumps when detached from the lipids. Overall, the prolonged dynamics of the Na(+) ions are concluded to be interesting for the physics of the whole membrane, especially considering its interaction dynamics with charged macromolecular surfaces.

Calculation of the electrostatic potential of lipid bilayers from molecular dynamics simulations: methodological issues.

The electrostatic properties of lipid membranes are of profound importance as they are directly associated with membrane potential and, consequently, with numerous membrane-mediated biological phenomena. Here we address a number of methodological issues related to the computation of the electrostatic potential from atomic-scale molecular dynamics simulations of lipid bilayers. We discuss two slightly different forms of Poisson equation that are normally used to calculate the membrane potential: (i) a classical form when the potential and the electric field are chosen to be zero on one of the sides of a simulation box and (ii) an alternative form, when the potential is set to be the same on the opposite sides of a simulation box. Both forms differ by a position-dependent correction term, which has been shown to be proportional to the overall dipole moment of a bilayer system (for neutral systems). For symmetric bilayers we demonstrate that both approaches give essentially the same potential profiles, provided that simulations are long enough (a production run of at least 100 ns is required) and that fluctuations of the center of mass of a bilayer are properly accounted for. In contrast, for asymmetric lipid bilayers, the second approach is no longer appropriate due to a nonzero net dipole moment across a simulation box with a single asymmetric bilayer. We demonstrate that in this case the electrostatic potential can adequately be described by the classical form of Poisson equation, provided that it is employed in conjunction with tin-foil boundary conditions, which exactly balance a nonzero surface charge of a periodically replicated multibilayer system. Furthermore, we show that vacuum boundary conditions give qualitatively similar potential profiles for asymmetric lipid bilayers as compared to the conventional periodic boundaries, but accurate determination of the transmembrane potential difference is then hindered due to detachment of some water dipoles from bulk aqueous solution to vacuum.

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.

Role of phosphatidylglycerols in the stability of bacterial membranes.

An extensive 100-ns molecular dynamics simulation of lipid bilayer composed of mixture of phosphatidylethanolamine (PE) and phosphatidylglycerol (PG) was performed to elucidate the role of PGs to the stability of bacterial membranes. In addition, a control simulation of pure PE over 150 ns was performed. We observed that PGs decrease both the PE headgroup protrusions into the water phase, and the PE headgroup motion along bilayer normal. The above effects are caused by stronger inter-lipid interactions in the mixed bilayer: the number of hydrogen bonds created by PEs is 34% higher in the mixed than in the pure bilayer. Another contribution is due to the numerous ion-mediated inter-lipid links, which strongly enhance interface stability. That provides a plausible mechanism for preventing lipid desorption from the membrane, for example, under the influence of an organic solvent. A more compact and less dynamic interface structure also decreases membrane permeability. That provides a possible mechanism for stabilizing, e.g., bacterial membranes.