Direct and Fast Assessment of Antimicrobial Surface Activity Using Molecular Dynamics Simulation and Time-Lapse Imaging.

With the alarming rise of antimicrobial resistance, studies on bacteria-surface interactions are both relevant and timely. Scanning electron microscopy and colony forming unit counting are commonly used techniques but require sophisticated sample preparation and long incubation time. Here, we present a direct method based on molecular dynamics simulation of nanostructured surfaces providing in silico predictions, complemented with time-lapse fluorescence imaging to study live interactions of bacteria at the membrane-substrate level. We evaluate its effectiveness in predicting and statistically analyzing the temporal evolution and spatial distribution of prototypical bacteria with costained nucleoids and membranes (E. coli) on surfaces with nanopillars. We observed cell reorientation, clustering, membrane damage, growth inhibition, and in the extreme case of hydrocarbon-coated nanopillars, this was followed by cell disappearance, validating the obtained simulation results. Contrary to commonly used experimental methods, microscopy data are fast processed, in less than 1 h. In particular, the bactericidal effects can be straightforwardly detected and correlated with surface morphology and/or wettability.

Interaction of Graphene Nanoparticles and Lipid Membranes Displaying Different Liquid Orderings: A Molecular Dynamics Study.

Understanding the effects of graphene-based nanomaterials on lipid membranes is fundamental to determine their environmental impact and the efficiency of their biomedical use. By means of molecular dynamics simulations of simple model lipid bilayers, we analyze in detail the different interaction modes. We have studied bilayers consisting of lipid species (including cholesterol) which display different internal liquid orderings. Nanometric graphene layers can be transiently adsorbed onto the lipid membrane and/or inserted in its hydrophobic region. Once inserted, graphene nanometric flakes display a diffusive dynamics in the membrane plane, they adopt diverse orientations depending on their size and oxidation degree, and they show a particular aversion to be placed close to cholesterol molecules in the membrane. Addition of graphene to phase-segregated ternary membranes is also investigated in the context of the lipid raft model for the lipid organization of biological membranes. Our simulation results show that graphene layers can be inserted indistinctly in the ordered and disordered regions. Once inserted, nanometric flakes migrate to disordered and cholesterol-poor lipid phases.

Effects of fullerene on lipid bilayers displaying different liquid ordering: a coarse-grained molecular dynamics study.

BACKGROUND: The toxic effects and environmental impact of nanomaterials, and in particular of Fullerene particles, are matters of serious concern. It has been reported that fullerene molecules enter the cell membrane and occupy its hydrophobic region. Understanding the effects of carbon-based nanoparticles on biological membranes is therefore of critical importance to determine their exposure risks. METHODS: We report on a systematic coarse-grained molecular dynamics study of the interaction of fullerene molecules with simple model cell membranes. We have analyzed bilayers consisting of lipid species with different degrees of unsaturation and a variety of cholesterol fractions. Addition of fullerene particles to phase-segregated ternary membranes is also investigated in the context of the lipid raft model for the organization of the cell membrane. RESULTS: Fullerene addition to lipid membranes modifies their structural properties like thickness, area and internal ordering of the lipid species, as well as dynamical aspects such as molecular diffusion and cholesterol flip-flop. Interestingly, we show that phase-segregating ternary lipid membranes accumulate fullerene molecules preferentially in the liquid-disordered domains promoting phase-segregation and domain alignment across the membrane. CONCLUSIONS: Lipid membrane internal ordering determines the behavior and distribution of fullerene particle, and this, in turn, determines the influence of fullerene on the membrane. Lipid membranes are good solvents of fullerene molecules, and in particular those with low internal ordering. GENERAL SIGNIFICANCE: Preference of fullerene molecules to be dissolved in the more disordered hydrophobic regions of a lipid bilayer and the consequent alteration of its phase behavior may have important consequences on the activity of biological cell membranes and on the bioconcentration of fullerene in living organisms.

Lipid Vesicle Interaction with Hydrophobic Surfaces: A Coarse-Grained Molecular Dynamics Study.

Active surfaces are presently tailored to cause specific effects on living cells, which can be useful in many fields. Their development requires the understanding of the molecular mechanisms of interaction between lipid-enveloped entities and solid surfaces. Here, using coarse-grained molecular dynamics simulations, we have analyzed the different interaction modes of coated substrates with lipid vesicles that mimic biological envelopes. For neutral and hydrophobically functionalized substrates, three action modes on contacting vesicles have been obtained including intact, partially broken, and completely destroyed vesicles. The molecular mechanisms for each interaction pathway and the corresponding energy balances have been analyzed in detail. Interestingly, we have shown that any specific action mode can be obtained by appropriately tailoring the wetting characteristics of the surface coating. In particular, we have shown that surfaces that are simultaneously hydrophobic and oleophilic are optimal to fully disrupt the contacting vesicle lipid bilayer.

Electroporation of heterogeneous lipid membranes.

Electroporation is the basis for the transfection of genetic material and for drug delivery to cells, including electrochemotherapy for cancer. By means of molecular dynamics many aspects of membrane electroporation have been unveiled at the molecular detail in simple, homogeneous, lipid bilayers. However, the correspondence of these findings \with the process happening in cell membranes requires, at least, the consideration of laterally structured membranes. Here, I present a systematic molecular dynamics study of bilayers composed of different liquid-ordered and liquid-disordered lipid phases subjected to a transversal electric field. The simulations reveal two significant results. First, the electric field mainly affects the properties of the disordered phases, so that electroporation takes place in these membrane regions. Second, the smaller the disordered domains are, the faster they become electroporated. These findings may have a relevant significance in the experimental application of cell electroporation in vivo since it implies that electro-induced and pore-mediated transport processes occur in particularly small disordered domains of the plasma membrane, thus locally affecting only specific regions of the cell.

Direct mapping of nanoscale compositional connectivity on intact cell membranes.

Lateral segregation of cell membranes is accepted as a primary mechanism for cells to regulate a diversity of cellular functions. In this context, lipid rafts have been conceptualized as organizing principle of biological membranes where underlying cholesterol-mediated selective connectivity must exist even at the resting state. However, such a level of nanoscale compositional connectivity has been challenging to prove. Here we used single-molecule near-field scanning optical microscopy to visualize the nanolandscape of raft ganglioside GM1 after tightening by its ligand cholera toxin (CTxB) on intact cell membranes. We show that CTxB tightening of GM1 is sufficient to initiate a minimal raft coalescence unit, resulting in the formation of cholesterol-dependent GM1 nanodomains < 120 nm in size. This particular arrangement appeared independent of cell type and GM1 expression level on the membrane. Simultaneous dual color high-resolution images revealed that GPI anchored and certain transmembrane proteins were recruited to regions proximal (< 150 nm) to CTxB-GM1 nanodomains without physical intermixing. Together with in silico experiments, our high-resolution data conclusively demonstrate the existence of raft-based interconnectivity at the nanoscale. Such a linked state on resting cell membranes constitutes thus an obligatory step toward the hierarchical evolution of large-scale raft coalescence upon cell activation.

Size-controlled nanopores in lipid membranes with stabilizing electric fields.

Molecular dynamics (MD) has been shown to be a useful tool for unveiling many aspects of pore formation in lipid membranes under the influence of an applied electric field. However, the study of the structure and transport properties of electropores by means of MD has been hampered by difficulties in the maintenance of a stable electropore in the typically small simulated membrane patches. We describe a new simulation scheme in which an initially larger porating field is systematically reduced after pore formation to lower stabilizing values to produce stable, size-controlled electropores, which can then be characterized at the molecular level. A new method allows the three-dimensional modeling of the irregular shape of the pores obtained as well as the quantification of its volume. The size of the pore is a function of the value of the stabilizing field. At lower fields the pore disappears and the membrane recovers its normal shape, although in some cases long-lived, fragmented pores containing unusual lipid orientations in the bilayer are observed.

Effect of integral proteins in the phase stability of a lipid bilayer: application to raft formation in cell membranes.

The existence of lipid rafts is a controversial issue. The affinity of cholesterol for saturated lipids is manifested in macroscopic phase separation in model membranes, and is believed to be the thermodynamic driving force for raft formation. However, there is no clear reason to explain the small (nanometric) size of raft domains in cell membranes. In a recent paper Yethiraj and Weisshaar [Biophys. J. 93, 3113 (2007)] proposed that the effect of neutral integral membrane proteins may prevent from the formation of large lipid domains. In this paper we extend this approach by studying the effect of the protein size, as well as the lipid-protein interaction. Depending on these factors, two different mechanisms for nanodomain stabilization are shown to be possible for static proteins. The application of these results to a biological context is discussed.

Why is the sn-2 chain of monounsaturated glycerophospholipids usually unsaturated whereas the sn-1 chain is saturated? Studies of 1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (SOPC) and 1-oleoyl-2-stearoyl-sn-glycero-3-phosphatidylcholine (OSPC) membranes with and without cholesterol.

Despite the large number of possible glycerol-based phospholipids, biological membranes contain only a small number of them. For example, double bonds in acyl chains are preferably located in the sn-2 chain. The question that emerges is: Why? We have addressed this question through atomistic simulations by considering pure one-component bilayers comprising monounsaturated glycerophospholipids [1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (SOPC) and 1-oleoyl-2-stearoyl-sn-glycero-3-phosphatidylcholine (OSPC)] and membranes of these lipids mixed with cholesterol. By considering the cases in which an individual double bond is in either the sn-1 or the sn-2 chain, we elucidated how membrane properties depend on this intrinsic feature. We found small but systematic differences in all structural and dynamic membrane properties that we considered. It turns out that the differences are driven by two factors: the mismatch in the acyl chain lengths and the interaction of the double bond in the acyl chains with the cholesterol off-plane methyl groups. The results highlight the fact that unsaturated sn-2 chains lead to more disordered membranes than systems with unsaturated sn-1 chains. The differences between the two isomers are enhanced when cholesterol is present as a result of the interaction of the off-plane cholesterol methyl groups with the double-bond carbon segments in the lipid acyl chains.

Nonequilibrium patterns in phase-separating ternary membranes.

We present a nonequilibrium approach for the study of a two-dimensional phase-separating ternary mixture. When the component that promotes phase separation is dynamically exchanged with the medium, the separation process is halted and actively maintained finite-size segregation domains appear in the system. In addition to this effect, already reported in our earlier work [J. Gómez, F. Sagués, and R. Reigada, Phys. Rev. E 77, 021907 (2008)], the use of a generic Ginzburg-Landau formalism and the inclusion of thermal fluctuations provide a more dynamic description of the resulting domain organization. Its size, shape, and stability properties are studied. Larger and more circular and stable domains are formed when decreasing the recycling rate, increasing the mobility of the exchanged component, and the mixture is quenched deeper. We expect this outcome to be of applicability in raft phenomenology in plasmatic cell membranes.

Interaction modes between nanosized graphene flakes and liposomes: Adsorption, insertion and membrane fusion.

BACKGROUND: Understanding the effects of graphene-based nanomaterials on lipid membranes is critical to determine their environmental impact and their efficiency in the biomedical context. Graphene has been reported to favourably interact with biological and model lipid membranes. METHODS: We report on a systematic coarse-grained molecular dynamics study of the interaction modes of graphene nanometric flakes with POPC/cholesterol liposome membranes. We have simulated graphene layers with a variety of sizes and oxidation degrees, and we have analyzed the trajectories, the interaction modes, and the energetics of the observed phenomena. RESULTS: Three interaction modes are reported. Graphene can be transiently adsorbed onto the liposome membrane and/or inserted in its hydrophobic region. Inserted nanosheets prefer a perpendicular orientation, and tilt in order to maximize the contact with phospholipid tails while avoiding the contact with cholesterol molecules. When placed between two liposomes, graphene facilitates their fusion in a single vesicle. CONCLUSIONS: Graphene can be temporary adsorbed on the liposome before insertion. Bilayer curvature has an influence on the orientation of inserted graphene particles. Cholesterol molecules are depleted from the surrounding of graphene particles. Graphene layers may catalyse membrane fusion by bypassing the energy barrier required in stalk formation. GENERAL SIGNIFICANCE: Nanometric graphene layers can be adsorbed/inserted in lipid-based membranes in different manners and affect the cholesterol distribution in the membrane, implying important consequences on the structure and functionality of biological cell membranes, and on the bioaccumulation of graphene in living organisms. The graphene-mediated mechanism opens new possibilities for vesicle fusion in the experimental context.

A spatial toggle switch drives boundary formation in development.

Herein we introduce a multicellular network motif that performs as a spatial toggle switch and explains how boundary formation can be faithfully accomplished in developmental processes. Importantly, we show that expression and activity patterns of proteins must be simultaneously characterized for a proper understanding and description of the underlying mechanism. Our in silico experiments, in agreement with in vivo results, evaluate different genetic backgrounds and shed light on the dynamics of boundary formation. In addition, we provide an estimation of relevant biological parameters and a robustness analysis.

Interplay of unsaturated phospholipids and cholesterol in membranes: effect of the double-bond position.

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

Use of an enhanced bulk diffusion-based algorithm for phase separation of a ternary mixture.

Phase separation kinetics of two-dimensional ternary mixtures have been studied by means of a Monte Carlo approach. Standard Kawasaki kinetics are impractical to study the late stages of the segregation process at deep quenches. An extension of the accelerated algorithm for binary mixtures proposed by Marko and Barkema [Phys. Rev. E 52, 2522 (1995)] is presented to overcome this problem in a three-component system discretized with two coupled lattices. We study the domain growth and the scaling behavior over a wide range of quench depths. Computer performances of the Kawasaki and the accelerated schemes are compared.

Influence of cis double-bond parametrization on lipid membrane properties: how seemingly insignificant details in force-field change even qualitative trends.

We have employed atomistic molecular dynamics simulations to investigate the effect of double-bond parametrization on lipid membrane properties. As models, we use one-component membranes composed of glycerol-based phosphatidylcholines (PCs) with monounsaturated acyl chains, and we complement these studies by additional PC/cholesterol simulations. We compare differences between double-bond parametrizations by varying the position of the double bond systematically along the lipid hydrocarbon chains. The results give rise for concern: They indicate that the double-bond description may change not only the quantitative but also the qualitative nature of membrane behavior. In particular, we find that the double-bond description which accounts for skew states in the vicinity of a double bond predicts a maximum in membrane disorder, when the double bond resides at the middle of an acyl chain, in agreement with experiments. The more commonly used description which does not accommodate skew states, however, predicts membrane disorder to decrease monotonically as the double bond is shifted from the glycerol backbone to the end of an acyl chain. The results highlight the importance of properly describing double bonds especially in many-component membranes, where the interplay of different molecule types is difficult to predict on intuitive grounds.

Actively maintained lipid nanodomains in biomembranes.

Lipid rafts, defined as domains rich in cholesterol and sphingolipids, are involved in many important plasma membrane functions. Recent studies suggest that the way cells handle membrane cholesterol is fundamental in the formation of such lateral heterogeneities. We propose to model the plasma membrane as a nonequilibrium phase-separating system where cholesterol is dynamically incorporated and released. The model shows how cellular regulation of membrane cholesterol may determine the nanoscale lipid organization when the lipid mixture is close to a phase separation boundary, providing a plausible mechanism for raft formation in vivo.

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.

Self-organized propagation patterns from dynamic self-assembly in monolayers.

Propagation of localized orientational waves, as imaged by Brewster angle microscopy, is induced by low intensity linearly polarized light inside axisymmetric smectic-C confined domains in a photosensitive molecular thin film at the air/water interface (Langmuir monolayer). Results from numerical simulations of a model that couples photoreorientational effects and long-range elastic forces are presented. Differences are stressed between our scenario and the paradigmatic wave phenomena in excitable chemical media.

Alteration of interleaflet coupling due to compounds displaying rapid translocation in lipid membranes.

The spatial coincidence of lipid domains at both layers of the cell membrane is expected to play an important role in many cellular functions. Competition between the surface interleaflet tension and a line hydrophobic mismatch penalty are conjectured to determine the transversal behavior of laterally heterogeneous lipid membranes. Here, by a combination of molecular dynamics simulations, a continuum field theory and kinetic equations, I demonstrate that the presence of small, rapidly translocating molecules residing in the lipid bilayer may alter its transversal behavior by favoring the spatial coincidence of similar lipid phases.

Equilibrium microphase separation in the two-leaflet model of lipid membranes.

Because of the coupling between local lipid composition and the thickness of the membrane, microphase separation in two-component lipid membranes can take place; such effects may underlie the formation of equilibrium nanoscale rafts. Using a kinetic description, this phenomenon is analytically and numerically investigated. The phase diagram is constructed through the stability analysis for linearized kinetic equations, and conditions for microphase separation are discussed. Simulations of the full kinetic model reveal the development of equilibrium membrane nanostructures with various morphologies from the initial uniform state.