Asymmetry of lipid bilayers induced by monovalent salt: Atomistic molecular-dynamics study

Sulfated Nutrition Modifies Nutrient Content and Photosynthetic Pigment Concentration in Cabbage under Salt Stress

Negative effects of salt stress may be counteracted by adequate management of sulfated nutrition. Herein, we applied 3.50, 4.25, and 5.00 mM SO42- in a nutrient solution to counteract salt stress induced by 75 and 150 mM NaCl in cabbage cv. Royal. The increase in NaCl concentration from 75 to 150 mM reduced the contents of macronutrients and micronutrients in the shoot. When increasing from 3.50 to 4.25 mM SO42-, the contents of nitrogen (N), phosphorous (P), potassium (K), calcium (Ca), magnesium (Mg), and sulfur (S) in shoots were enhanced, at both concentrations of NaCl. Increasing from 3.50 to 4.25 mM SO42- enhanced iron (Fe), zinc (Zn), manganese (Mn), and sodium (Na) concentrations with 75 mM NaCl. With 150 mM NaCl, the increase from 3.50 to 4.25 mM SO42- enhanced the contents of Cu and Mn, but also those of Na. Chlorophylls a, b, and total decreased as the concentration of SO42- increased in plants treated with 150 mM NaCl. With 75 mM NaCl, carotenoid concentration had a positive relationship with SO42-. Hence, the 4.25 mM SO42- concentration increased the contents of macronutrients and micronutrients in the presence of 75 mM NaCl, while, with 150 mM NaCl, it improved the contents of macronutrients except K. The chlorophyll a/chlorophyll b ratio remained close to 3 when the plants were treated with 5.00 mM SO42-, regardless of NaCl. Similarly, this level of SO42- increased the concentration of carotenoids, which translated into reductions in the total chlorophyll/carotenoid ratios, indicating a protective effect of the photosynthetic apparatus. It is concluded that higher doses of sulfur favor the accumulation of nutrients and increase the concentration of carotenoids under salt stress.

Efficient Simulations of Solvent Asymmetry Across Lipid Membranes Using Flat-Bottom Restraints

The routinely employed periodic boundary conditions complicate molecular simulations of physiologically relevant asymmetric lipid membranes together with their distinct solvent environments. Therefore, separating the extracellular fluid from its cytosolic counterpart has often been performed using a costly double-bilayer setup. Here, we demonstrate that the lipid membrane and solvent asymmetry can be efficiently modeled with a single lipid bilayer by applying an inverted flat-bottom potential to ions and other solute molecules, thereby restraining them to only interact with the relevant leaflet. We carefully optimized the parameters of the suggested method so that the results obtained using the flat-bottom and double-bilayer approaches become mutually indistinguishable. Then, we apply the flat-bottom approach to lipid bilayers with various compositions and solvent environments, covering ions and cationic peptides to validate the approach in a realistic use case. We also discuss the possible limitations of the method as well as its computational efficiency and provide a step-by-step guide on how to set up such simulations in a straightforward manner.

Tilted State Population of Antimicrobial Peptide PGLa Is Coupled to the Transmembrane Potential

Cationic antimicrobial peptide PGLa gets into close contact with the anionic bacterial cell membrane, facilitating cross-membrane transport phenomena and membrane disruption depending on the concentration. The mechanisms of action are closely associated with the tilted insertion geometry of PGLa. Therefore, we aimed to understand the interaction between the transmembrane potential (TMP) and the orientation of the membrane-bound PGLa helix. Molecular dynamics simulations were performed with TMP, and we found that the PGLa tilt angle relative to the membrane is coupled with the TMP. Elevated TMP increases the population of the tilted state. We observed positive feedback between the tilt angle and the TMP, which occurs due to the electrostatic interaction between the peptidic helix and the Na⁺ cations at the membrane-water interface. These TMP coupled phenomena can contribute to understanding the direct antimicrobial and adjuvant effects of PGLa in combination with regular antibiotics.

The Voltage-Dependent Deactivation of the KvAP Channel Involves the Breakage of Its S4 Helix

Voltage-gated potassium channels (Kv) allow ion permeation upon changes of the membrane electrostatic potential (Vm). Each subunit of these tetrameric channels is composed of six transmembrane helices, of which the anti-parallel helix bundle S1-S4 constitutes the voltage-sensor domain (VSD) and S5-S6 forms the pore domain. Here, using 82 molecular dynamics (MD) simulations involving 266 replicated VSDs, we report novel responses of the archaebacterial potassium channel KvAP to membrane polarization. We show that the S4 α-helix, which is straight in the experimental crystal structure solved under depolarized conditions (Vm ∼ 0), breaks into two segments when the cell membrane is hyperpolarized (Vm << 0), and reversibly forms a single straight helix following depolarization (Vm = 0). The outermost segment of S4 translates along the normal to the membrane, bringing new perspective to previously paradoxical accessibility experiments that were initially thought to imply the displacement of the whole VSD across the membrane. The novel model is applied through steered and unbiased MD simulations to the recently solved whole structure of KvAP. The simulations show that the resting state involves a re-orientation of the S5 α-helix by ∼ 5-6 degrees in respect to the normal of the bilayer, which could result in the constriction and closure of the selectivity filter. Our findings support the idea that the breakage of S4 under (hyper)polarization is a general feature of Kv channels with a non-swapped topology.

Influence of External NaCl Salt on Membrane Rigidity of Neutral DOPC Vesicles

Sodium chloride (NaCl) is a very common molecule in biotic and abiotic aqueous environments. In both cases, variation of ionic strength is inevitable. In addition to the osmotic variation posed by such perturbations, the question of whether the interactions of monovalent ions Na⁺ and Cl^-, especially with the neutral head groups of phospholipid membranes are impactful enough to change the membrane rigidity, is still not entirely understood. We investigated the dynamics of 1,2-di-(octadecenoyl)-sn-glycero-3-phosphocholine (DOPC) vesicles with zwitterionic neutral head groups in the fluid phase with increasing external salt concentration. At higher salt concentrations, we observe an increase in bending rigidity from neutron spin echo (NSE) spectroscopy and an increase in bilayer thickness from small-angle X-ray scattering (SAXS). We compared different models to distinguish membrane undulations, lipid tail motions, and the translational diffusion of the vesicles. All of the models indicate an increase in bending rigidity by a factor of 1.3-3.6. We demonstrate that even down to t > 10 ns and for Q > 0.07 Å^-1, the observed NSE relaxation spectra are influenced by translational diffusion of the vesicles. For t < 5 ns, the lipid tail motion dominates the intermediate dynamic structure factor. As the salt concentration increases, this contribution diminishes. We introduced a time-dependent analysis for the bending rigidity that highlights only a limited Zilman-Granek time window in which the rigidity is physically meaningful.

Molecular dynamics study of lipid bilayers modeling outer and inner leaflets of plasma membranes of mouse hepatocytes. I. Differences in physicochemical properties between the two leaflets

Outer and inner leaflets of plasma cell membranes have different lipid compositions, and the membrane properties of each leaflet can differ from each other significantly due to these composition differences. However, because of the experimental difficulty in measuring the membrane properties for each leaflet separately, the differences are not well understood at a molecular level. In this study, we constructed two lipid bilayer systems, modeling outer and inner leaflets of plasma membranes of mouse hepatocytes based on experimental composition data. The ion concentration in the interlamellar water phase was also set to match the concentration in extra- and intracellular fluids. The differences in physical properties between the outer and inner leaflets of mouse hepatocyte cell membrane models were investigated by performing 1.2 μs-long all-atomistic molecular dynamics calculations under physiological temperature and pressure conditions (310.15 K and 1 atm). The calculated electron density profiles along the bilayer normal for each model bilayer system captured well the asymmetric feature of the experimental electron density profile across actual cell plasma membranes, indicating that our procedure of modeling the outer and inner leaflets of the cell plasma membranes was satisfactory. We found that compared to the outer leaflet model, the inner leaflet model had a very bulky and soft structure in the lateral direction. To confirm the differences, membrane fluidity was measured from the lateral diffusivity and relaxation times. The fluidity was significantly higher in the inner leaflet model than in the outer leaflet model. We also discuss two topics that are of wide interest in biology, i.e., the interdigitation of acyl tails of lipid molecules between two monolayers and the lateral concentration fluctuation of lipid molecules in the bilayers.

Interaction Profiles and Stability of Rigid and Polymer-Tethered Lipid Bilayer Models at Highly Charged and Highly Adhesive Contacts

Understanding interaction force versus distance profiles of supported lipid bilayers (SLBs) is relevant to a number of areas, which rely on these model systems, including, e.g., characterization of ligand/receptor interactions or bacterial adhesion. Here, the stability of 4 different SLB architectures was compared using the surface forces apparatus (SFA) and atomic force microscopy (AFM). Specifically, the outer envelope of the bilayer systems remained constant as 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC). The inner layer was varied between DPPC and 1,2-dipalmitoyl-3-trimethylammonium-propane (DPTAP) both on mica, and self-assembled monolayers (SAMs) of hexadecanethiol and the polymer-tethered diphytanylglycerol-tetraethylene glycol-lipoid acid (DPhyTL) on smooth gold surfaces. In that same order these gave an increasing strength of interaction between the inner layer and the supporting substrate and hence improved stability under highly adhesive conditions. Detachment profiles from highly charged and highly adhesive contacts were characterized, and approach characteristics were fitted to DLVO models. We find increasing stability under highly adhesive loads, approaching the hydrophobic limit of the adhesive energy between the inner and outer layers for the SAM-based systems. For all four SLBs we further compare AFM surface topographies, which strongly depend on preparation conditions, and the DLVO fitting of the SFA approach curves finds a strong charge regulation behavior during interaction, dependent on the particular model system. In addition, we find undulation characteristics during approach and separation. The increased stability of the complex architectures on a gold support makes these model systems an ideal starting point for studying more complex strongly adhesive/interacting systems, including, for example, ligand/receptor interactions, biosensing interactions, or cell/surface interactions.

Coupling of Membrane Nanodomain Formation and Enhanced Electroporation near Phase Transition

Biological cells are enveloped by a heterogeneous lipid bilayer that prevents the uncontrolled exchange of substances between the cell interior and its environment. In particular, membranes act as a continuous barrier for salt and macromolecules to ensure proper physiological functions within the cell. However, it has been shown that membrane permeability strongly depends on temperature and, for phospholipid bilayers, displays a maximum at the transition between the gel and fluid phase. Here, extensive molecular dynamics simulations of dipalmitoylphosphatidylcholine bilayers were employed to characterize the membrane structure and dynamics close to phase transition, as well as its stability with respect to an external electric field. Atomistic simulations revealed the dynamic appearance and disappearance of spatially related nanometer-sized thick ordered and thin interdigitating domains in a fluid-like bilayer close to the phase transition temperature (Tm). These structures likely represent metastable precursors of the ripple phase that vanished at increased temperatures. Similarly, a two-phase bilayer with coexisting gel and fluid domains featured a thickness minimum at the interface because of splaying and interdigitating lipids. For all systems, application of an external electric field revealed a reduced bilayer stability with respect to pore formation for temperatures close to Tm. Pore formation occurred exclusively in thin interdigitating membrane nanodomains. These findings provide a link between the increased membrane permeability and the structural heterogeneity close to phase transition.

Molecular Dynamics Simulations of Membrane Proteins: An Overview

Simulations of membrane proteins have been rising in popularity in the past decade. Advancements in technology and force fields made it possible to simulate behavior of membrane proteins. Membrane protein simulations can now be used as supporting evidence for experimental findings, for elucidating protein mechanisms, and validating protein crystal structures. Unrelated to experimental data, these simulations can also serve to investigate larger scale processes like protein sorting, protein-membrane interactions, and more. In this review, the history as well as the state-of-the-art methodologies in membrane protein simulations will be summarized. An emphasis will be put on how to set up the system and on the current models for the different components of the simulation system. An overview of the available tools for membrane protein simulation will be given, and current limitations and prospects will also be discussed.

Direct measurement of surface charge distribution in phase separating supported lipid bilayers

The local surface charge density of the cell membrane influences regulation and localization of membrane proteins. The local surface charge density could, until recently, not be measured directly under physiological conditions, and it was largely a hypothetical yet very important parameter. Here we use unsaturated lipids of a distinct charge (DOTAP, DOPC, and DOPG) and a neutral fully saturated lipid (DPPC) to create model membranes with phase separating domains of a defined charge. We then apply quantitative surface charge microscopy (QSCM) to investigate the local surface charge density; this is a technique based on a scanning ion conductance microscope (SICM) capable of measuring surface charge density with nanoscale lateral resolution. We are able to clearly distinguish lipid domains from charge and topography in all three model membranes. The measured surface charge densities furthermore reveal that disordered domains formed by charged lipids are in fact not only impure, but also incorporate uncharged saturated lipids. We estimate that at least 30% of disordered domains in DOPG : DPPC and DOTAP : DPPC will be DPPC. These ratios could present a limit for the formation of charged domains in lipid membranes.

Effect of head group orientation on phospholipid assembly

The relationship between bilayer stability and lipid head group orientation is reported. In this work, molecular-dynamics simulations are performed to analyze the structure-property relationship of lipid biomembranes, taking into account coarse-grained model lipid interactions. The work explains the molecular scale mechanism of the phase behavior of lipid systems due to ion-lipid or anesthetic-lipid interactions, where reorientations of dipoles play a key role in modifying lipid phases and thereby alter biomembrane function. Our study demonstrates that simple dipolar reorientation is indeed sufficient in tuning a bilayer to a randomly flipped nonbilayer lamellar phase. This study may be used to assess the impact of changes in lipid phase characteristics on biomembrane structure due to the presence of anesthetics and ions.

How does ytterbium chloride interact with DMPC bilayers? A computational and experimental study

Lanthanide salts have been studied for many years, primarily in Nuclear Magnetic Resonance (NMR) experiments of mixed lipid–protein systems and more recently to study lipid flip-flop in model membrane systems.

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.

Computer modelling studies of the bilayer/water interface

This review summarises high resolution studies on the interface of lamellar lipid bilayers composed of the most typical lipid molecules which constitute the lipid matrix of biomembranes. The presented results were obtained predominantly by computer modelling methods. Whenever possible, the results were compared with experimental results obtained for similar systems. The first and main section of the review is concerned with the bilayer-water interface and is divided into four subsections. The first describes the simplest case, where the interface consists only of lipid head groups and water molecules and focuses on interactions between the lipid heads and water molecules; the second describes the interface containing also mono- and divalent ions and concentrates on lipid-ion interactions; the third describes direct inter-lipid interactions. These three subsections are followed by a discussion on the network of direct and indirect inter-lipid interactions at the bilayer interface. The second section summarises recent computer simulation studies on the interactions of antibacterial membrane active compounds with various models of the bacterial outer membrane. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Free energy landscapes of sodium ions bound to DMPC–cholesterol membrane surfaces at infinite dilution

Typical configuration of two DMPC lipids and one cholesterol molecule solvating one sodium ion, together with two additional water molecules.

Molecular-level insight into the interactions of DNA with phospholipid bilayers: barriers and triggers

A zwitterionic phospholipid bilayer represents a repulsive barrier for DNA binding; this barrier can be overcome through adsorption of divalent cations to the bilayer surface.

Membrane pore formation in atomistic and coarse-grained simulations

Biological cells and their organelles are protected by ultra thin membranes. These membranes accomplish a broad variety of important tasks like separating the cell content from the outer environment, they are the site for cell-cell interactions and many enzymatic reactions, and control the in- and efflux of metabolites. For certain physiological functions e.g. in the fusion of membranes and also in a number of biotechnological applications like gene transfection the membrane integrity needs to be compromised to allow for instance for the exchange of polar molecules across the membrane barrier. Mechanisms enabling the transport of molecules across the membrane involve membrane proteins that form specific pores or act as transporters, but also so-called lipid pores induced by external fields, stress, or peptides. Recent progress in the simulation field enabled to closely mimic pore formation as supposed to occur in vivo or in vitro. Here, we review different simulation-based approaches in the study of membrane pores with a focus on lipid pore properties such as their size and energetics, poration mechanisms based on the application of external fields, charge imbalances, or surface tension, and on pores that are induced by small molecules, peptides, and lipids. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Effect of Force Field Parameters on Sodium and Potassium Ion Binding to Dipalmitoyl Phosphatidylcholine Bilayers

The behavior of electrolytes in molecular dynamics simulations of zwitterionic phospholipid bilayers is very sensitive to the force field parameters used. Here, several 200 ns molecular dynamics of simulations of dipalmitoyl phosphotidylcholine (PC) bilayers in 0.2 M sodium or potassium chloride using various common force field parameters for the cations are presented. All employed parameter sets give a larger number of Na(+) ions than K(+) ions that bind to the lipid heads, but depending on the parameter choice quite different results are seen. A wide range of coordination numbers for the Na(+) and K(+) ions is also observed. These findings have been analyzed and compared to published experimental data. Some simulations produce aggregates of potassium chloride, indicating (in accordance with published simulations) that these force fields do not reproduce the delicate balance between salt and solvated ions. The differences between the force fields can be characterized by one single parameter, the electrostatic radius of the ion, which is correlated to σMO (M represents Na(+)/K(+)), the Lennard-Jones radius. When this parameter exceeds a certain threshold, binding to the lipid heads is no longer observed. One would, however, need more accurate experimental data to judge or rank the different force fields precisely. Still, reasons for the poor performance of some of the parameter sets are clearly demonstrated, and a quality control procedure is provided.

Systematic Parameterization of Monovalent Ions Employing the Nonbonded Model

Monovalent ions play fundamental roles in many biological processes in organisms. Modeling these ions in molecular simulations continues to be a challenging problem. The 12-6 Lennard-Jones (LJ) nonbonded model is widely used to model monovalent ions in classical molecular dynamics simulations. A lot of parameterization efforts have been reported for these ions with a number of experimental end points. However, some reported parameter sets do not have a good balance between the two Lennard-Jones parameters (the van der Waals (VDW) radius and potential well depth), which affects their transferability. In the present work, via the use of a noble gas curve we fitted in former work (J. Chem. Theory Comput. 2013, 9, 2733), we reoptimized the 12-6 LJ parameters for 15 monovalent ions (11 positive and 4 negative ions) for three extensively used water models (TIP3P, SPC/E, and TIP4P(EW)). Since the 12-6 LJ nonbonded model performs poorly in some instances for these ions, we have also parameterized the 12-6-4 LJ-type nonbonded model (J. Chem. Theory Comput. 2014, 10, 289) using the same three water models. The three derived parameter sets focused on reproducing the hydration free energies (the HFE set) and the ion-oxygen distance (the IOD set) using the 12-6 LJ nonbonded model and the 12-6-4 LJ-type nonbonded model (the 12-6-4 set) overall give improved results. In particular, the final parameter sets showed better agreement with quantum mechanically calculated VDW radii and improved transferability to ion-pair solutions when compared to previous parameter sets.

Ionic strength and composition govern the elasticity of biological membranes. A study of model DMPC bilayers by force- and transmission IR spectroscopy

Infrared (IR) spectroscopy was used to quantify the ion mixture effect of seawater (SW), particularly the contribution of Mg(2+) and Ca(2+) as dominant divalent cations, on the thermotropic phase behaviour of 1,2-dimyristoyl-sn-glycero-3-posphocholine (DMPC) bilayers. The changed character of the main transition at 24 °C from sharp to gradual in films and the 1 °C shift of the main transition temperature in dispersions reflect the interactions of lipid headgroups with the ions in SW. Force spectroscopy was used to quantify the nanomechanical hardness of a DMPC supported lipid bilayer (SLB). Considering the electrostatic and ion binding equilibrium contributions while systematically probing the SLB in various salt solutions, we showed that ionic strength had a decisive influence on its nanomechanics. The mechanical hardness of DMPC SLBs in the liquid crystalline phase linearly increases with the increasing fraction of all ion-bound lipids in a series of monovalent salt solutions. It also linearly increases in the gel phase but almost three times faster (the corresponding slopes are 4.9 nN/100 mM and 13.32 nN/100 mM, respectively). We also showed that in the presence of divalent ions (Ca(2+) and Mg(2+)) the bilayer mechanical hardness was unproportionally increased, and that was accompanied with the decrease of Na(+) ion and increase of Cl(-) ion bound lipids. The underlying process is a cooperative and competitive ion binding in both the gel and the liquid crystalline phase. Bilayer hardness thus turned out to be very sensitive to ionic strength as well as to ionic composition of the surrounding medium. In particular, the indicated correlation helped us to emphasize the colligative properties of SW as a naturally occurring complex ion mixture.

A methodological combined framework for roadmapping biosensor research: a fault tree analysis approach within a strategic technology evaluation frame

Biosensor technology began in the 1960s to revolutionize instrumentation and measurement. Despite the glucose sensor market success that revolutionized medical diagnostics, and artificial pancreas promise currently the approval stage, the industry is reluctant to capitalize on other relevant university-produced knowledge and innovation. On the other hand, the scientific literature is extensive and persisting, while the number of university-hosted biosensor groups is growing. Considering the limited marketability of biosensors compared to the available research output, the biosensor field has been used by the present authors as a suitable paradigm for developing a methodological combined framework for "roadmapping" university research output in this discipline. This framework adopts the basic principles of the Analytic Hierarchy Process (AHP), replacing the lower level of technology alternatives with internal barriers (drawbacks, limitations, disadvantages), modeled through fault tree analysis (FTA) relying on fuzzy reasoning to count for uncertainty. The proposed methodology is validated retrospectively using ion selective field effect transistor (ISFET) - based biosensors as a case example, and then implemented prospectively membrane biosensors, putting an emphasis on the manufacturability issues. The analysis performed the trajectory of membrane platforms differently than the available market roadmaps that, considering the vast industrial experience in tailoring and handling crystallic forms, suggest the technology path of biomimetic and synthetic materials. The results presented herein indicate that future trajectories lie along with nanotechnology, and especially nanofabrication and nano-bioinformatics, and focused, more on the science-path, that is, on controlling the natural process of self-assembly and the thermodynamics of bioelement-lipid interaction. This retained the nature-derived sensitivity of the biosensor platform, pointing out the differences between the scope of academic research and the market viewpoint.

Macro- versus microscopic view on the electrokinetics of a water-membrane interface

Electrophoresis is an experimental method widely used to study electrostatic properties of interfaces. Here, we question the validity of the macroscopic theory for the planar geometry by Helmholtz and Smoluchowski by considering a POPC bilayer in an aqueous solution with 500 mM NaCl, using molecular dynamics simulations. We find that POPC shows positive electrophoretic mobility due to adsorption of sodium ions at the lipid headgroups. The theory assumes that the region in which the water density undergoes a transition from the bulk value to zero (interfacial width) is small compared to the Debye screening length. This separation of length scale is not fullfilled in the present case. Hence, contrasting the theory, we observe that the surface is not sharply defined, continuum hydrodynamics is not applicable, the effective viscosity in the double layer is increased compared to the bulk, and the zeta potential is dominated by the dipole potential. Our results might have widespread implications for interpretation of electrokinetic studies in general.

A theoretical model for calculating voltage sensitivity of ion channels and the application on Kv1.2 potassium channel

Voltage sensing confers conversion of a change in membrane potential to signaling activities underlying the physiological processes. For an ion channel, voltage sensitivity is usually experimentally measured by fitting electrophysiological data to Boltzmann distributions. In our study, a two-state model of the ion channel and equilibrium statistical mechanics principle were used to test the hypothesis of empirically calculating the overall voltage sensitivity of an ion channel on the basis of its closed and open conformations, and determine the contribution of individual residues to the voltage sensing. We examined the theoretical paradigm by performing experimental measurements with Kv1.2 channel and a series of mutants. The correlation between the calculated values and the experimental values is at respective level, R(2) = 0.73. Our report therefore provides in silico prediction of key conformations and has identified additional residues critical for voltage sensing.

Force spectroscopy reveals the effect of different ions in the nanomechanical behavior of phospholipid model membranes: the case of potassium cation

How do metal cations affect the stability and structure of phospholipid bilayers? What role does ion binding play in the insertion of proteins and the overall mechanical stability of biological membranes? Investigators have used different theoretical and microscopic approaches to study the mechanical properties of lipid bilayers. Although they are crucial for such studies, molecular-dynamics simulations cannot yet span the complexity of biological membranes. In addition, there are still some experimental difficulties when it comes to testing the ion binding to lipid bilayers in an accurate way. Hence, there is a need to establish a new approach from the perspective of the nanometric scale, where most of the specific molecular phenomena take place. Atomic force microscopy has become an essential tool for examining the structure and behavior of lipid bilayers. In this work, we used force spectroscopy to quantitatively characterize nanomechanical resistance as a function of the electrolyte composition by means of a reliable molecular fingerprint that reveals itself as a repetitive jump in the approaching force curve. By systematically probing a set of bilayers of different composition immersed in electrolytes composed of a variety of monovalent and divalent metal cations, we were able to obtain a wealth of information showing that each ion makes an independent and important contribution to the gross mechanical resistance and its plastic properties. This work addresses the need to assess the effects of different ions on the structure of phospholipid membranes, and opens new avenues for characterizing the (nano)mechanical stability of membranes.

Concentration dependence of NaCl ion distributions around DPPC lipid bilayers

We study the coordination of excess NaCl to zwitterionic DPPC lipid bilayers using molecular dynamics simulations. We find that Na ions directly coordinate with the DPPC lipid carbonyl groups. As the number of excess ions increases, the number of coordinated ions increases, until it reaches a plateau at a ratio near 1 ion per every four lipids at 310 K, and 1 ion per every six lipids at 323 K. The area per lipid decreases as the number of excess ions is increased. For low number of ions per lipids (1:16 and 1:8), most Na ions are bound to the lipid carbonyls, while the Cl form an ionic cloud around the lipid choline groups. As a result of the Na binding, the lipid has an effective positive charge density. The residence time of Na ions bound to the lipid is longer than 40 ns, while Cl ions exchange faster than the nanoseconds timescale. We find that the bound Na ions replace ordered water around the carbonyls. The net linear charge density near the carbonyl groups stays positive, regardless of the presence of excess salt in the solution.

NaCl interactions with phosphatidylcholine bilayers do not alter membrane structure but induce long-range ordering of ions and water

It is generally accepted that ions interact directly with lipids in biological membranes. Decades of biophysical studies on pure lipid bilayer systems have shown that only certain types of ions, most significantly large anions and multivalent cations, can fundamentally alter the structure and dynamics of lipid bilayers. It has long been accepted that at physiological concentrations NaCl ions do not alter the physical behavior or structure of bilayers composed solely of zwitterionic phosphatidylcholine (PC) lipids. Recent X-ray scattering experiments have reaffirmed this dogma, showing that below 1 M concentration, NaCl does not significantly alter bilayer structure. However, despite this history, there is an ongoing controversy within the molecular dynamics (MD) simulation community regarding NaCl/PC interactions. In particular, the CHARMM and GROMOS force fields show dramatically different behavior, including the effect on bilayer structure, surface potential, and the ability to form stable, coordinated ion-lipid complexes. Here, using long-timescale, constant-pressure simulations under the newest version of the CHARMM force field, we find that Na⁺ and Cl⁻ associate with PC head groups in a POPC bilayer with approximately equal, though weak, affinity, and that the salt has a negligible effect on bilayer structure, consistent with earlier CHARMM results and more recent X-ray data. The results suggest that interpretation of simulations where lipids interact with charged groups of any sort, including charged proteins, must be carefully scrutinized.

Interaction of salicylate and a terpenoid plant extract with model membranes: reconciling experiments and simulations

We investigate the effects of two structurally similar small cyclic molecules: salicylic acid and perillic acid on a zwitterionic model lipid bilayer, and show that both molecules might have biological activity related to membrane thinning. Salicylic acid is a nonsteroidal antiinflammatory drug, some of the pharmacological properties of which arise from its interaction with the lipid bilayer component of the plasma membrane. Prior simulations show that salicylate orders zwitterionic lipid membranes. However, this is in conflict with Raman scattering and vesicle fluctuation analysis data, which suggest the opposite. We show using extensive molecular dynamics simulations, cumulatively >2.5 μs, that salicylic acid indeed disorders membranes with concomitant membrane thinning and that the conflict arose because prior simulations suffered from artifacts related to the sodium-ion induced condensation of zwitterionic lipids modeled by the Berger force field. Perillic acid is a terpenoid plant extract that has antiinfective and anticancer properties, and is extensively used in eastern medicine. We found that perillic acid causes large-scale membrane thinning and could therefore exert its antimicrobial properties via a membrane-lytic mechanism reminiscent of antimicrobial peptides. Being more amphipathic, perillic acid is more potent in disrupting lipid headgroup packing, and significantly modifies headgroup dipole orientation. Like salicylate, the membrane thinning effect of perillic acid is masked by the presence of sodium ions. As an alternative to sodium cations, we advocate the straightforward solution of using larger countercations like potassium or tetra-methyl-ammonium that will maintain electroneutrality but not interact strongly with, and thus not condense, the lipid bilayer.

Ca(2+) bridging of apposed phospholipid bilayers

In an effort to provide insight into the mechanism of Ca(2+)-induced fusion of lipid vesicles, molecular dynamics simulations in the isobaric-isothermal ensemble are used to investigate interactions of Ca(2+) with apposed lipid bilayers in close proximity. Simulations reveal the formation of a Ca(2+)-phospholipid "anhydrous complex" between apposed bilayers, whereas similar calculations performed with Na(+) display only complexation between neighboring lipids within the same bilayer. The binding of Ca(2+) to apposed phospholipids brings large regions of the bilayers into close contact (<4 Å), displacing water from phospholipid head groups in the process and creating regions of local dehydration. Dehydration of the apposed bilayers leads to ordering of the phospholipid tails, which is partially disrupted by the presence of Ca(2+)-phospholipid bridges.

Direct imaging of salt effects on lipid bilayer ordering at sub-molecular resolution

The interactions of salts with lipid bilayers are known to alter the properties of membranes and therefore influence their structure and dynamics. Sodium and calcium cations penetrate deeply into the headgroup region and bind to the lipids, whereas potassium ions only loosely associate with lipid molecules and mostly remain outside of the headgroup region. We investigated a dipalmitoylphosphatidylcholine (DPPC) bilayer in the gel phase in the presence of all three cations with a concentration of Ca²+ ions an order of magnitude smaller than the Na+ and K+ ions. Our findings indicate that the area per unit cell does not significantly change in these three salt solutions. However the lipid molecules do re-order non-isotropically under the influence of the three different cations. We attribute this reordering to a change in the highly directional intermolecular interactions caused by a variation in the dipole-dipole bonding arising from a tilt of the headgroup out of the membrane plane. Measurements in different NaCl concentrations also show a non-isotropic re-ordering of the lipid molecules.

Characterization of perfluorooctylbromide-based nanoemulsion particles using atomistic molecular dynamics simulations

Perfluorocarbon-based nanoemulsion particles have arisen as promising platforms for the cellular delivery of imaging and therapeutic agents to specific targets. However, current knowledge of the agent delivery mechanism is limited to qualitative and phenomenological models. Lack of detail at the molecular level has hence delayed optimizing or customizing nanoemulsion particles for therapeutic and imaging applications. Here we report the first atomistic structural details of a perfluorooctylbromide-based (PFOB-based) nanoemulsion particle (NEP) with a 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) lipid emulsifier. Newly developed PFOB force-field parameters were used in molecular dynamics simulations to model the PFOB-NEP interface in a planar configuration. These PFOB force field parameters were developed and tested to reproduce the characteristics of bulk PFOB as well as PFOB at interfaces with water and emulsifying phospholipids. The modeled PFOB-NEP interface demonstrated significant intercalation of PFOB into the emulsifying lipid monolayer and consequent changes in the structural, electrostatic, and mechanical properties of the POPC monolayer and PFOB. This intercalation provides an explanation for experimental data demonstrating melittin tryptophan fluorescence quenching upon binding to the nanoemulsion particles through the observation of direct contact between the melittin tryptophan and the PFOB bromine. Additionally, the atomistic details of the PFOB-NEP interface structure provided by our simulations are used to suggest the influence of each component on PFOB-NEP delivery function which will be tested in future coarse-grained simulations.

Examining the origins of the hydration force between lipid bilayers using all-atom simulations

Using 237 all-atom double bilayer simulations, we examined the thermodynamic and structural changes that occur as a phosphatidylcholine lipid bilayer stack is dehydrated. The simulated system represents a micropatch of lipid multilayer systems that are studied experimentally using surface force apparatus, atomic force microscopy and osmotic pressure studies. In these experiments, the hydration level of the system is varied, changing the separation between the bilayers, in order to understand the forces that the bilayers feel as they are brought together. These studies have found a curious, strongly repulsive force when the bilayers are very close to each other, which has been termed the "hydration force," though the origins of this force are not clearly understood. We computationally reproduce this repulsive, relatively free energy change as bilayers come together and make qualitative conclusions as to the enthalpic and entropic origins of the free energy change. This analysis is supported by data showing structural changes in the waters, lipids and salts that have also been seen in experimental work. Increases in solvent ordering as the bilayers are dehydrated are found to be essential in causing the repulsion as the bilayers come together.

Fusogenic activity of cationic lipids and lipid shape distribution

Addition of co-lipids into cationic lipid formulations is considered as promoting cell delivery of DNA by enhancing fusion processes with cell membranes. Here, by combining FRET and confocal microscopy, we demonstrate that some cationic lipids do not require a co-lipid to fuse efficiently with cells. These cationic lipids are able to self-organize into bilayers that are stable enough to form liposomes, while presenting some destabilizing properties reminiscent of the conically shaped fusogenic co-lipid, DOPE. We therefore analyzed the resident lipid structures in cationic bilayers by molecular dynamics simulations, clustering the individual lipid structures into populations of similarly shaped molecules, as opposed to the classical approach of using the static packing parameter to define the lipid shapes. Comparison of fusogenic properties with these lipid populations suggests that the ratio of cylindrical versus conical lipid populations correlates with the ability to fuse with cell membranes.