The influence of 1-alkanols and external pressure on the lateral pressure profiles of lipid bilayers

Understanding the Molecular Mechanism of Anesthesia: Effect of General Anesthetics and Structurally Similar Non-Anesthetics on the Properties of Lipid Membranes

General anesthesia can be caused by various, chemically very different molecules, while several other molecules, many of which are structurally rather similar to them, do not exhibit anesthetic effects at all. To understand the origin of this difference and shed some light on the molecular mechanism of general anesthesia, we report here molecular dynamics simulations of the neat dipalmitoylphosphatidylcholine (DPPC) membrane as well as DPPC membranes containing the anesthetics diethyl ether and chloroform and the structurally similar non-anesthetics n-pentane and carbon tetrachloride, respectively. To also account for the pressure reversal of anesthesia, these simulations are performed both at 1 bar and at 600 bar. Our results indicate that all solutes considered prefer to stay both in the middle of the membrane and close to the boundary of the hydrocarbon domain, at the vicinity of the crowded region of the polar headgroups. However, this latter preference is considerably stronger for the (weakly polar) anesthetics than for the (apolar) non-anesthetics. Anesthetics staying in this outer preferred position increase the lateral separation between the lipid molecules, giving rise to a decrease of the lateral density. The lower lateral density leads to an increased mobility of the DPPC molecules, a decreased order of their tails, an increase of the free volume around this outer preferred position, and a decrease of the lateral pressure at the hydrocarbon side of the apolar/polar interface, a change that might well be in a causal relation with the occurrence of the anesthetic effect. All these changes are clearly reverted by the increase of pressure. Furthermore, non-anesthetics occur in this outer preferred position in a considerably smaller concentration and hence either induce such changes in a much weaker form or do not induce them at all.

A Membrane Biophysics Perspective on the Mechanism of Alcohol Toxicity

Motivations for understanding the underlying mechanisms of alcohol toxicity range from economical to toxicological and clinical. On the one hand, acute alcohol toxicity limits biofuel yields, and on the other hand, acute alcohol toxicity provides a vital defense mechanism to prevent the spread of disease. Herein the role that stored curvature elastic energy (SCE) in biological membranes might play in alcohol toxicity is discussed, for both short and long-chain alcohols. Structure-toxicity relationships for alcohols ranging from methanol to hexadecanol are collated, and estimates of alcohol toxicity per alcohol molecule in the cell membrane are made. The latter reveal a minimum toxicity value per molecule around butanol before alcohol toxicity per molecule increases to a maximum around decanol and subsequently decreases again. The impact of alcohol molecules on the lamellar to inverse hexagonal phase transition temperature (T_H) is then presented and used as a metric to assess the impact of alcohol molecules on SCE. This approach suggests the nonmonotonic relationship between alcohol toxicity and chain length is consistent with SCE being a target of alcohol toxicity. Finally, in vivo evidence for SCE-driven adaptations to alcohol toxicity in the literature are discussed.

Serotonergic drugs modulate the phase behavior of complex lipid bilayers

Serotonin is an endogenous neurotransmitter involved in both physiological and pathophysiological processes. Traditionally, serotonin acts as a ligand for G protein-coupled receptors (GPCRs) leading to subsequent cell signaling. However, serotonin can also bind to lipid membranes with high affinity and modulate the phase behavior in 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC)/N-palmitoyl-D-erythro-sphingosylphosphorylcholine (PSM)/cholesterol model membranes mimicking the outer leaflet of the plasma membrane. Here, we investigated if serotonergic drugs containing the pharmacophore from serotonin would also modulate phase behavior in lipid membranes in a similar fashion. We used ²H NMR spectroscopy to explore the phase behavior of POPC/PSM/cholesterol (4/4/2 molar ratio) mixtures in the presence of the serotonergic drugs aripiprazole, BRL-54443, BW-723C86, and CP-135807 at a lipid to drug molar ratio of 10:1. POPC and PSM were perdeuterated in the palmitoyl chain, respectively, and prepared in individual samples. Numerical lineshape simulations of the ²H NMR spectra were used to calculate the order parameter profiles and projected lengths of the saturated acyl chains. All serotonergic drugs induce two components in ²H NMR spectra, indicating that they increased the hydrophobic mismatch between the thickness of the coexisting lipid phases leading to larger domain sizes, relatively similarly to serotonin. AFM force indentation and Raman spectral studies, which interrogate membrane mechanical properties, also indicate changes in membrane order in the presence of these drugs. These findings highlight how serotonergic drugs alter membrane phase behavior and could modulate both target and other membrane proteins, possibly explaining the side effects observed for serotonergic and other clinically relevant drugs.

Short-Chained Alcohols Make Membrane Surfaces Conducive for Melittin Action: Implication for the Physiological Role of Alcohols in Cells

Alcohols are a part of cellular metabolism, but their physiological roles are not well understood. We investigated the effects of short-chain alcohols on Daphnia pulex and model membranes mimicking the lipid composition of eukaryotic inner mitochondrial membranes. We also studied the synergistic effects of alcohols with the bee venom membrane-active peptide, melittin, which is structurally similar to endogenous membrane-active peptides. The alcohols, from ethanol to octanol, gradually decreased the heart rate and the mitochondrial ATP synthesis of daphnia; in contrast, in combination with melittin, which exerted no sizeable effect, they gradually increased both the heart rate and the ATP synthesis. Lipid packing and the order parameter of oriented films, monitored by EPR spectroscopy of the spin-labeled probe 5-doxylstrearic acid, revealed gradual alcohol-assisted bilayer to non-bilayer transitions in the presence of melittin; further, while the alcohols decreased, in combination with melittin they increased the order parameter of the film, which is attributed to the alcohol-facilitated association of melittin with the membrane. A ¹H-NMR spectroscopy of the liposomes confirmed the enhanced induction of a non-bilayer lipid phase that formed around the melittin, without the permeabilization of the liposomal membrane. Our data suggest that short-chain alcohols, in combination with endogenous peptides, regulate protein functions via modulating the lipid polymorphism of membranes.

Spontaneous Membrane Nanodomain Formation in the Absence or Presence of the Neurotransmitter Serotonin

Detailed knowledge on the formation of biomembrane domains, their structure, composition, and physical characteristics is scarce. Despite its frequently discussed importance in signaling, e.g., in obtaining localized non-homogeneous receptor compositions in the plasma membrane, the nanometer size as well as the dynamic and transient nature of domains impede their experimental characterization. In turn, atomistic molecular dynamics (MD) simulations combine both, high spatial and high temporal resolution. Here, using microsecond atomistic MD simulations, we characterize the spontaneous and unbiased formation of nano-domains in a plasma membrane model containing phosphatidylcholine (POPC), palmitoyl-sphingomyelin (PSM), and cholesterol (Chol) in the presence or absence of the neurotransmitter serotonin at different temperatures. In the ternary mixture, highly ordered and highly disordered domains of similar composition coexist at 303 K. The distinction of domains by lipid acyl chain order gets lost at lower temperatures of 298 and 294 K, suggesting a phase transition at ambient temperature. By comparison of domain ordering and composition, we demonstrate how the domain-specific binding of the neurotransmitter serotonin results in a modified domain lipid composition and a substantial downward shift of the phase transition temperature. Our simulations thus suggest a novel mode of action of neurotransmitters possibly of importance in neuronal signal transmission.

Serotonin Alters the Phase Equilibrium of a Ternary Mixture of Phospholipids and Cholesterol

Unsaturated and saturated phospholipids tend to laterally segregate, especially in the presence of cholesterol. Small molecules such as neurotransmitters, toxins, drugs etc. possibly modulate this lateral segregation. The small aromatic neurotransmitter serotonin (5-HT) has been found to bind to membranes. We studied the lipid structure and packing of a ternary membrane mixture consisting of palmitoyl-oleoyl-phosphatidylcholine, palmitoyl-sphingomyelin, and cholesterol at a molar ratio of 4/4/2 in the absence and in the presence of 5-HT, using a combination of solid-state 2H NMR, atomic force microscopy, and atomistic molecular dynamics (MD) simulations. Both NMR and MD report formation of a liquid ordered (L o ) and a liquid disordered (L d ) phase coexistence with small domains. Lipid exchange between the domains was fast such that single component 2H NMR spectra are detected over a wide temperature range. A drastic restructuring of the domains was induced when 5-HT is added to the membranes at a 9 mol% concentration relative to the lipids. 2H NMR spectra of all components of the mixture showed two prominent contributions indicative of molecules of the same kind residing both in the disordered and the ordered phase. Compared to the data in the absence of 5-HT, the lipid chain order in the disordered phase was further decreased in the presence of 5-HT. Likewise, addition of serotonin increased lipid chain order within the ordered phase. These characteristic lipid chain order changes were confirmed by MD simulations. The 5-HT-induced larger difference in lipid chain order results in more pronounced differences in the hydrophobic thickness of the individual membrane domains. The correspondingly enlarged hydrophobic mismatch between ordered and disordered phases is assumed to increase the line tension at the domain boundary, which drives the system into formation of larger size domains. These results not only demonstrate that small membrane binding molecules such as neurotransmitters have a profound impact on essential membrane properties. It also suggests a mechanism by which the interaction of small molecules with membranes can influence the function of membrane proteins and non-cognate receptors. Altered membrane properties may modify lateral sorting of membrane protein, membrane protein conformation, and thus influence their function as suspected for neurotransmitters, local anesthetics, and other small drug molecules.

The cut-off effect of n-alcohols in lipid rafts: A lipid-dependent phenomenon☆

n-Aliphatic alcohols act as anesthetics only up to a certain chain length, beyond which its biological activity disappears. This is known as the 'cut-off' phenomenon. Although the most accepted explanation is based on action sites in membrane proteins, it is not well understood why alcohols alter their functions. The structural dependence of these protein receptors to lipid domains known as 'lipid rafts', suggests a new approach to tackle the puzzling phenomenon. In this work, by performing molecular dynamic simulations (MDS) to explore the lipid role, we provide relevant molecular details about the membrane-alcohol interaction at the cut-off point regime. Since the high variability of the cut-off points found on protein receptors in neurons may be a consequence of differences in the lipid composition surrounding such proteins, our results could have a clear-cut importance.

Microscopic Picture of Calcium-Assisted Lipid Demixing and Membrane Remodeling Using Multiscale Simulations

The specificity of anionic phospholipids-calcium ion interaction and lipid demixing has been established as a key regulatory mechanism in several cellular signaling processes. The mechanism and implications of this calcium-assisted demixing have not been elucidated from a microscopic point of view. Here, we present an overview of atomic interactions between calcium and phospholipids that can drive nonideal mixing of lipid molecules in a model lipid bilayer composed of zwitterionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)) and anionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine (POPS)) lipids with computer simulations at multiple resolutions. Lipid nanodomain formation and growth were driven by calcium-enabled lipid bridging of the charged phosphatidylserine (PS) headgroups, which were favored against inter-POPS dipole interactions. Consistent with several experimental studies of calcium-associated membrane sculpting, our analyses also suggest modifications in local membrane curvature and cross-leaflet couplings as a response to such induced lateral heterogeneity. In addition, reverse mapping to a complementary atomistic description revealed structural insights in the presence of anionic nanodomains, at timescales not accessed by previous computational studies. This work bridges information across multiple scales to reveal a mechanistic picture of calcium ion's impact on membrane biophysics.

Membrane Curvature Sensing by Amphipathic Helices: Insights from Implicit Membrane Modeling

Sensing and generation of lipid membrane curvature, mediated by the binding of specific proteins onto the membrane surface, play crucial roles in cell biology. A number of mechanisms have been proposed, but the molecular understanding of these processes is incomplete. All-atom molecular dynamics simulations have offered valuable insights but are extremely demanding computationally. Implicit membrane simulations could provide a viable alternative, but current models apply only to planar membranes. In this work, the implicit membrane model 1 is extended to spherical and tubular membranes. The geometric change from planar to curved shapes is straightforward but insufficient for capturing the full curvature effect, which includes changes in lipid packing. Here, these packing effects are taken into account via the lateral pressure profile. The extended implicit membrane model 1 is tested on the wild-types and mutants of the antimicrobial peptide magainin, the ALPS motif of arfgap1, α-synuclein, and an ENTH domain. In these systems, the model is in qualitative agreement with experiments. We confirm that favorable electrostatic interactions tend to weaken curvature sensitivity in the presence of strong hydrophobic interactions but may actually have a positive effect when those are weak. We also find that binding to vesicles is more favorable than binding to tubes of the same diameter and that the long helix of α-synuclein tends to orient along the axis of tubes, whereas shorter helices tend to orient perpendicular to it. Adoption of a specific orientation could provide a mechanism for coupling protein oligomerization to tubule formation.

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.

Effect of general anesthetics on the properties of lipid membranes of various compositions

Computer simulations of four lipid membranes of different compositions, namely neat DPPC and PSM, and equimolar DPPC-cholesterol and PSM-cholesterol mixtures, are performed in the presence and absence of the general anesthetics diethylether and sevoflurane both at 1 and 600 bar. The results are analyzed in order to identify membrane properties that are potentially related to the molecular mechanism of anesthesia, namely that change in the same way in any membrane with any anesthetics, and change oppositely with increasing pressure. We find that the lateral lipid density satisfies both criteria: it is decreased by anesthetics and increased by pressure. This anesthetic-induced swelling is attributed to only those anesthetic molecules that are located close to the boundary of the apolar phase. This lateral expansion is found to lead to increased lateral mobility of the lipids, an effect often thought to be related to general anesthesia; to an increased fraction of the free volume around the outer preferred position of anesthetics; and to the decrease of the lateral pressure in the nearby range of the ester and amide groups, a region into which anesthetic molecules already cannot penetrate. All these changes are reverted by the increase of pressure. Another important finding of this study is that cholesterol has an opposite effect on the membrane properties than anesthetics, and, correspondingly, these changes are less marked in the presence of cholesterol. Therefore, changes in the membrane that can lead to general anesthesia are expected to occur in the membrane domains of low cholesterol content.

Fluorinated Alcohols' Effects on Lipid Bilayer Properties

Fluorinated alcohols (fluoroalcohols) have physicochemical properties that make them excellent solvents of peptides, proteins, and other compounds. Like other alcohols, fluoroalcohols also alter membrane protein function and lipid bilayer properties and stability. Thus, the questions arise: how potent are fluoroalcohols as lipid-bilayer-perturbing compounds, could small residual amounts that remain after adding compounds dissolved in fluoroalcohols alter lipid bilayer properties sufficiently to affect membranes and membrane protein function, and do they behave like other alcohols? To address these questions, we used a gramicidin-based fluorescence assay to determine the bilayer-modifying potency of selected fluoroalcohols: trifluoroethanol (TFE), HFIP, and perfluoro-tert-butanol (PFTB). These fluoroalcohols alter bilayer properties in the low (PFTB) to high (TFE) mM range. Using the same assay, we determined the bilayer partitioning of the alcohols. When referenced to the aqueous concentrations, the fluoroalcohols are more bilayer perturbing than their nonfluorinated counterparts, with the largest fluoroalcohol, PFTB, being the most potent and the smallest, TFE, the least. When referenced to the mole fractions in the membrane, however, the fluoroalcohols have equal or lesser bilayer-perturbing potency than their nonfluorinated counterparts, with TFE being more bilayer perturbing than PFTB. We compared the fluoroalcohols' molecular level bilayer interactions using atomistic molecular dynamics simulations and showed how, at higher concentrations, they can cause bilayer breakdown using absorbance measurements and ³¹P nuclear magnetic resonance.

Lipid-mediated mode of action of local anesthetics on lipid pores induced by polyenes, peptides and lipopeptides

The effects of local anesthetics (LAs), namely, lidocaine (LDC), prilocaine (PLC), mepivacaine (MPV), bupivacaine (BPV), procaine (PC), and tetracaine (TTC), on the steady-state transmembrane conductance induced by the cis-side addition of the antifungal polyene macrolide antibiotic, nystatin (NYS), in planar lipid bilayers were studied. The addition of TTC to model membranes comprising DOPC and cholesterol (33 mol%) led to a nearly twenty-fold increase in the steady-state NYS-induced membrane conductance. BPV slightly enhanced the channel-forming activity of polyene. LDC, PLC, MPV, and PC did not affect the NYS-induced transmembrane current. We concluded that the effects of LAs on the channel-forming activity of NYS were in agreement with their effects on the elastic properties of model membranes. The ability of aminoamide LAs to promote calcein leakage from large unilamellar DOPC-vesicles was decreased in the following order: BPV >> LDC ≈ PLC ≈ MPV. LDC, PLC, and MPV produced a graded leakage of fluorescent marker from liposomes, up to 10-13%. A initial sharp jump in fluorescence after the introduction of BPV was attributed to the solubilization of liposomes and the formation of mixed DOPC:BPV-micelles. Differential scanning microcalorimetry (DSC) of large unilamellar DPPC-vesicles showed that the main transition temperature (Tm) is continuously decreased upon increasing concentrations of TTC. A sharp drop in the enthalpy of the transition at higher TTC concentrations indicated a formation of anesthetic/lipid mixed micelles. In contrast to TTC, PC slightly decreased Tm, broadened the DSC signal and did not provoke vesicle-to-micelle transition. Both the calcein leakage and DSC data together with the results of measurements of threshold voltages that are required to cause the lipid bilayer breakdown might indicate an alteration in the curvature lipid packing stress, induced by BPV and TTC. The data presented here lend support to a lipid-mediated mode of LAs action on NYS pores via an alteration in curvature stress near the trans-mouth. Similar results were obtained for several lipid pores, formed by polyene amphotericin B, lipopeptide syringomycin E, and the peptides magainin and melittin. This finding further developed the concept of non-specific regulation of lipid pores by LAs. In conclusion, the combination of nystatin with LAs could be a novel treatment for efficient therapy of superficial and mucosal candidiasis.

Alcohol Interactions with Lipid Bilayers

We investigate the structural changes to lipid membrane that ensue from the addition of aliphatic alcohols with various alkyl tail lengths. Small angle neutron diffraction from flat lipid bilayers that are hydrated through water vapor has been employed to eliminate possible artefacts of the membrane curvature and the alcohol's membrane-water partitioning. We have observed clear changes to membrane structure in both transversal and lateral directions. Most importantly, our results suggest the alteration of the membrane-water interface. The water encroachment has shifted in the way that alcohol loaded bilayers absorbed more water molecules when compared to the neat lipid bilayers. The experimental results have been corroborated by molecular dynamics simulations to reveal further details. Namely, the order parameter profiles have been fruitful in correlating the mechanical model of structural changes to the effect of anesthesia.

High-Throughput Simulations Reveal Membrane-Mediated Effects of Alcohols on MscL Gating

The mechanosensitive channels of large conductance (MscL) are bacterial membrane proteins that serve as last resort emergency release valves in case of severe osmotic downshock. Sensing bilayer tension, MscL channels are sensitive to changes in the bilayer environment and are, therefore, an ideal test case for exploring membrane protein coupling. Here, we use high-throughput coarse-grained molecular dynamics simulations to characterize MscL gating kinetics in different bilayer environments under the influence of alcohols. We performed over five hundred simulations to obtain sufficient statistics to reveal the subtle effects of changes in the membrane environment on MscL gating. MscL opening times were found to increase with the addition of the straight-chain alcohols ethanol, octanol, and to some extent dodecanol but not with hexadecanol. Increasing concentration of octanol increased the impeding effect, but only up to 10-20 mol %. Our in silico predictions were experimentally confirmed using reconstituted MscL in a liposomal fluorescent efflux assay. Our combined data reveal that the effect of alcohols on MscL gating arises not through specific binding sites but through a combination of the alcohol-induced changes to a number of bilayer properties and their alteration of the MscL-bilayer interface. Our work provides a key example of how extensive molecular simulations can be used to predict the functional modification of membrane proteins by subtle changes in their bilayer environment.

Nanosecond pulsed electric field (nsPEF) enhance cytotoxicity of cisplatin to hepatocellular cells by microdomain disruption on plasma membrane

Previous studies showed nanosecond pulsed electric field (nsPEF) can ablate solid tumors including hepatocellular carcinoma (HCC) but its effect on cell membrane is not fully understood. We hypothesized nsPEF disrupt the microdomains on outer-cellular membrane with direct mechanical force and as a result the plasma membrane permeability increases to facilitate the small molecule intake. Three HCC cells were pulsed one pulse per minute, an interval longer than nanopore resealing time. The cationized ferritin was used to mark up the electronegative microdomains, propidium iodide (PI) for membrane permeabilization, energy dispersive X-ray spectroscopy (EDS) for the negative cell surface charge and cisplatin for inner-cellular cytotoxicity. We demonstrated that the ferritin marked-microdomain and negative cell surface charge were disrupted by nsPEF caused-mechanical force. The cell uptake of propidium and cytotoxicity of DNA-targeted cisplatin increased with a dose effect. Cisplatin gains its maximum inner-cellular cytotoxicity when combining with nsPEF stimulation. We conclude that nsPEF disrupt the microdomains on the outer cellular membrane directly and increase the membrane permeabilization for PI and cisplatin. The microdomain disruption and membrane infiltration changes are caused by the mechanical force from the changes of negative cell surface charge.

Modelling of noble anaesthetic gases and high hydrostatic pressure effects in lipid bilayers

Our objective was to study molecular processes that might be responsible for inert gas narcosis and high-pressure nervous syndrome. The classical molecular dynamics trajectories (200 ns) of dioleoylphosphatidylcholine (DOPC) bilayers simulated by the Berger force field were evaluated for water and the atomic distribution of noble gases around DOPC molecules in the pressure range of 1-1000 bar and at a temperature of 310 K. Xenon and argon have been tested as model gases for general anaesthetics, and neon has been investigated for distortions that are potentially responsible for neurological tremors in hyperbaric conditions. The analysis of stacked radial pair distribution functions of DOPC headgroup atoms revealed the explicit solvation potential of the gas molecules, which correlates with their dimensions. The orientational dynamics of water molecules at the biomolecular interface should be considered as an influential factor, while excessive solvation effects appearing in the lumen of membrane-embedded ion channels could be a possible cause of inert gas narcosis. All the noble gases tested exhibit similar order parameter patterns for both DOPC acyl chains, which are opposite of the patterns found for the order parameter curve at high hydrostatic pressures in intact bilayers. This finding supports the 'critical volume' hypothesis of anaesthesia pressure reversal. The irregular lipid headgroup-water boundary observed in DOPC bilayers saturated with neon in the pressure range of 1-100 bar could be associated with the possible manifestation of neurological tremors at the atomic scale. The non-immobiliser neon also demonstrated the highest momentum impact on the normal component of the DOPC diffusion coefficient representing the monolayer undulation rate, which indicates that enhanced diffusivity rather than atomic size is the key factor.

Spontaneous adsorption of coiled-coil model peptides K and E to a mixed lipid bilayer

A molecular description of the lipid-protein interactions underlying the adsorption of proteins to membranes is crucial for understanding, for example, the specificity of adsorption or the binding strength of a protein to a bilayer, or for characterizing protein-induced changes of membrane properties. In this paper, we extend an automated in silico assay (DAFT) for binding studies and apply it to characterize the adsorption of the model fusion peptides E and K to a mixed phospholipid/cholesterol membrane using coarse-grained molecular dynamics simulations. In addition, we couple the coarse-grained protocol to reverse transformation to atomistic resolution, thereby allowing to study molecular interactions with high detail. The experimentally observed differential binding of the peptides E and K to membranes, as well as the increased binding affinity of helical over unstructered peptides, could be well reproduced using the polarizable Martini coarse-grained (CG) force field. Binding to neutral membranes is shown to be dominated by initial binding of the positively charged N-terminus to the phospholipid headgroup region, followed by membrane surface-aligned insertion of the peptide at the interface between the hydrophobic core of the membrane and its polar headgroup region. Both coarse-grained and atomistic simulations confirm a before hypothesized snorkeling of lysine side chains for the membrane-bound state of the peptide K. Cholesterol was found to be enriched in peptide vicinity, which is probably of importance for the mechanism of membrane fusion. The applied sequential multiscale method, using coarse-grained simulations for the slow adsorption process of peptides to membranes followed by backward transformation to atomistic detail and subsequent atomistic simulations of the preformed peptide-lipid complexes, is shown to be a versatile approach to study the interactions of peptides or proteins with biomembranes.

The effect of anaesthetics on the properties of a lipid membrane in the biologically relevant phase: a computer simulation study

Phospholipid membranes containing four different general anaesthetic molecules are simulated in the biologically relevant L_α phase at atmospheric and high pressures.

Diffusion and spectroscopy of water and lipids in fully hydrated dimyristoylphosphatidylcholine bilayer membranes

Microscopic structure and dynamics of water and lipids in a fully hydrated dimyristoylphosphatidylcholine phospholipid lipid bilayer membrane in the liquid-crystalline phase have been analyzed with all-atom molecular dynamics simulations based on the recently parameterized CHARMM36 force field. The diffusive dynamics of the membrane lipids and of its hydration water, their reorientational motions as well as their corresponding spectral densities, related to the absorption of radiation, have been considered for the first time using the present force field. In addition, structural properties such as density and pressure profiles, a deuterium-order parameter, surface tension, and the extent of water penetration in the membrane have been analyzed. Molecular self-diffusion, reorientational motions, and spectral densities of atomic species reveal a variety of time scales playing a role in membrane dynamics. The mechanisms of lipid motion strongly depend on the time scale considered, from fast ballistic translation at the scale of picoseconds (effective diffusion coefficients of the order of 10(-5) cm(2)/s) to diffusive flow of a few lipids forming nanodomains at the scale of hundreds of nanoseconds (diffusion coefficients of the order of 10(-8) cm(2)/s). In the intermediate regime of sub-diffusion, collisions with nearest neighbors prevent the lipids to achieve full diffusion. Lipid reorientations along selected directions agree well with reported nuclear magnetic resonance data and indicate two different time scales, one about 1 ns and a second one in the range of 2-8 ns. We associated the two time scales of reorientational motions with angular distributions of selected vectors. Calculated spectral densities corresponding to lipid and water reveal an overall good qualitative agreement with Fourier transform infrared spectroscopy experiments. Our simulations indicate a blue-shift of the low frequency spectral bands of hydration water as a result of its interaction with lipids. We have thoroughly analyzed the physical meaning of all spectral features from lipid atomic sites and correlated them with experimental data. Our findings include a "wagging of the tails" frequency around 30 cm(-1), which essentially corresponds to motions of the tail-group along the instantaneous plane formed by the two lipid tails, i.e., in-plane oscillations are clearly of bigger importance than those along the normal-to-the plane direction.

High hydrostatic pressure effects investigated by neutron scattering on lipid multilamellar vesicles

The effects of high hydrostatic pressure on the structure and dynamics of model membrane systems were investigated using neutron scattering. Diffraction experiments show shifts of the pre- and main-phase transitions of multilamellar vesicles of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) to higher temperatures with increased pressure which are close to results observed previously by other techniques, namely (10.4 ± 1.0) K kbar(-1) and (20.0 ± 0.5) K kbar(-1) for the two transitions. Backscattering spectroscopy reveals that the mean square displacements in the liquid phase are about 10% smaller at 300 bar and about 20% smaller at 600 bar compared to atmospheric pressure, whereas in the gel phase below the main phase transition the mean square displacements show a smaller difference in the dynamics of the three pressure values within the studied pressure range.

Inclusion of lateral pressure/curvature stress effects in implicit membrane models

Implicit membrane models usually treat the membrane as a hydrophobic slab and neglect lateral pressure/curvature stress effects. As a result, they cannot distinguish, for example, PE from PC lipids. Here, the implicit membrane model IMM1 is extended to include these effects using a combination of classical thermodynamics and membrane elasticity theory. The proposed model is tested by molecular dynamics simulation of the peptides alamethicin, melittin, cyclotide kalata B1, 18A, and KKpL15. The lateral pressure term stabilizes interfacial binding due to the negative pressure at the hydrocarbon-water interface. In agreement with experiment, increase in the peptide/lipid molar ratio shifts the equilibrium from the interfacial to the transmembrane orientation. Simulations of mixed DOPC/DOPE bilayers show that increase of the DOPE mole fraction in general stabilizes interfacial orientations and destabilizes transmembrane orientations. The extent of the stabilization or destabilization varies depending on the exact position of the peptides. The computational results are in good agreement with experiments.

Molecular mechanisms of drug action: an emerging view

Volatile anesthetics serve as useful probes of a conserved biological process that is essential to the proper functioning of the central nervous system. A kinetic and thermodynamic analysis of their unusual pharmacological and physiological characteristics has led to a general, predictive theory in which small molecules that adsorb to membranes modulate ion channel function by altering physical properties of membrane bilayers. A kinetic model that is both parsimonious and falsifiable has been developed to test this mechanism. This theory leads to predictions about the structure, function, origin, and evolution of synapses, the etiology of several diseases and disease symptoms affecting the brain, and the mechanism of action of several drugs that are used therapeutically. Neuronal membranes may offer an appealing drug target, given the large number of compounds that adsorb to interfaces and hence membranes.

Biophysical changes induced by xenon on phospholipid bilayers

Structural and dynamic changes in cell membrane properties induced by xenon, a volatile anesthetic molecule, may affect the function of membrane-mediated proteins, providing a hypothesis for the mechanism of general anesthetic action. Here, we use molecular dynamics simulation and differential scanning calorimetry to examine the biophysical and thermodynamic effects of xenon on model lipid membranes. Our results indicate that xenon atoms preferentially localize in the hydrophobic core of the lipid bilayer, inducing substantial increases in the area per lipid and bilayer thickness. Xenon depresses the membrane gel-liquid crystalline phase transition temperature, increasing membrane fluidity and lipid head group spacing, while inducing net local ordering effects in a small region of the lipid carbon tails and modulating the bilayer lateral pressure profile. Our results are consistent with a role for nonspecific, lipid bilayer-mediated mechanisms in producing xenon's general anesthetic action.

Diffusive nature of xenon anesthetic changes properties of a lipid bilayer: molecular dynamics simulations

Effects of general anesthesia can be controllable by the ambient pressure. We perform molecular dynamics simulations for a 1-palmitoyl-2-oleoyl phosphatidylethanolamine lipid bilayer with or without xenon molecules by changing the pressure to elucidate the mechanism of the pressure reversal of general anesthesia. According to the diffusive nature of xenon molecules in the lipid bilayer, a decrease in the orientational order of the lipid tails, an increase in the area and volume per lipid molecule, and an increase in the diffusivity of lipid molecules are observed. We show that the properties of the lipid bilayer with xenon molecules at high pressure come close to those without xenon molecules at 0.1 MPa. Furthermore, we find that xenon molecules are concentrated in the middle of the lipid bilayer at high pressures by the pushing effect and that the diffusivity of xenon molecules is suppressed. These results suggest that the pressure reversal originates from a jamming and suppression of the diffusivity of xenon molecules in lipid bilayers.

Fluorescence anisotropy of diphenylhexatriene and its cationic Trimethylamino derivative in liquid dipalmitoylphosphatidylcholine liposomes: opposing responses to isoflurane

Background

The mechanism of action of volatile general anesthetics has not yet been resolved. In order to identify the effects of isoflurane on the membrane, we measured the steady-state anisotropy of two fluorescent probes that reside at different depths. Incorporation of anesthetic was confirmed by shifting of the main phase transition temperature.

Results

In liquid crystalline dipalmitoylphosphatidylcholine liposomes, isoflurane (7-25 mM in the bath) increases trimethylammonium-diphenylhexatriene fluorescence anisotropy by ~0.02 units and decreases diphenylhexatriene anisotropy by the same amount.

Conclusions

The anisotropy data suggest that isoflurane decreases non-axial dye mobility in the headgroup region, while increasing it in the tail region. We propose that these results reflect changes in the lateral pressure profile of the membrane.

Alcohol's effects on lipid bilayer properties

Alcohols are known modulators of lipid bilayer properties. Their biological effects have long been attributed to their bilayer-modifying effects, but alcohols can also alter protein function through direct protein interactions. This raises the question: Do alcohol's biological actions result predominantly from direct protein-alcohol interactions or from general changes in the membrane properties? The efficacy of alcohols of various chain lengths tends to exhibit a so-called cutoff effect (i.e., increasing potency with increased chain length, which that eventually levels off). The cutoff varies depending on the assay, and numerous mechanisms have been proposed such as: limited size of the alcohol-protein interaction site, limited alcohol solubility, and a chain-length-dependent lipid bilayer-alcohol interaction. To address these issues, we determined the bilayer-modifying potency of 27 aliphatic alcohols using a gramicidin-based fluorescence assay. All of the alcohols tested (with chain lengths of 1-16 carbons) alter the bilayer properties, as sensed by a bilayer-spanning channel. The bilayer-modifying potency of the short-chain alcohols scales linearly with their bilayer partitioning; the potency tapers off at higher chain lengths, and eventually changes sign for the longest-chain alcohols, demonstrating an alcohol cutoff effect in a system that has no alcohol-binding pocket.

New insights into the molecular mechanisms of general anaesthetics

This paper provides new insights of how general anaesthetic research should be carried out in the future by an analysis of what we know, what we do not know and what we would like to know. I describe previous hypotheses on the mechanism of action of general anaesthetics (GAs) involving membranes and protein receptors. I provide the reasons why the GABA type A receptor, the NMDA receptor and the glycine receptor are strong candidates for the sites of action of GAs. I follow with a review on attempts to provide a mechanism of action, and how future research should be conducted with the help of physical and chemical methods.

Modulation of KvAP unitary conductance and gating by 1-alkanols and other surface active agents

The actions of alcohols and anesthetics on ion channels are poorly understood. Controversy continues about whether bilayer restructuring is relevant to the modulatory effects of these surface active agents (SAAs). Some voltage-gated K channels (Kv), but not KvAP, have putative low affinity alcohol-binding sites, and because KvAP structures have been determined in bilayers, KvAP could offer insights into the contribution of bilayer mechanics to SAA actions. We monitored KvAP unitary conductance and macroscopic activation and inactivation kinetics in PE:PG/decane bilayers with and without exposure to classic SAAs (short-chain 1-alkanols, cholesterol, and selected anesthetics: halothane, isoflurane, chloroform). At levels that did not measurably alter membrane specific capacitance, alkanols caused functional changes in KvAP behavior including lowered unitary conductance, modified kinetics, and shifted voltage dependence for activation. A simple explanation is that the site of SAA action on KvAP is its entire lateral interface with the PE:PG/decane bilayer, with SAA-induced changes in surface tension and bilayer packing order combining to modulate the shape and stability of various conformations. The KvAP structural adjustment to diverse bilayer pressure profiles has implications for understanding desirable and undesirable actions of SAA-like drugs and, broadly, predicts that channel gating, conductance and pharmacology may differ when membrane packing order differs, as in raft versus nonraft domains.

Coarse-grain modelling of DMPC and DOPC lipid bilayers

Our recently developed coarse-grain model for dimyristoylphosphatidylcholine (DMPC) has been improved and extended to dioleylphosphatidylcholine (DOPC), a more typical constituent of real biological membranes. Single-component DMPC and DOPC bilayers have been simulated using microsecond-long molecular dynamics. We investigated properties that are difficult or impossible to access experimentally, such as the pressure distribution, the spontaneous curvature and the diffusion pattern of individual lipid molecules. Moreover, we studied the dipole potential, a basic physical feature of paramount biological importance that cannot be currently modelled by other coarse-grain approaches. In fact, a complete representation of the system electrostatics and a realistic description of the water component make our method unique amongst the existing coarse-grain membrane models. The spontaneous permeation of water, a phenomenon out of reach of standard atomistic models, was also observed and quantified; this was possible thanks to the efficiency of our model, which is about two orders of magnitude less computationally expensive than atomic-level counterparts. Results are generally in good agreement with the literature data. Further model extensions and future applications are proposed.