Lateral pressure profiles in cholesterol-DPPC bilayers

Mechanical properties of anionic asymmetric bilayers from atomistic simulations

Mechanotransduction, the biological response to mechanical stress, is often initiated by activation of mechanosensitive (MS) proteins upon mechanically induced deformations of the cell membrane. A current challenge in fully understanding this process is in predicting how lipid bilayers deform upon the application of mechanical stress. In this context, it is now well established that anionic lipids influence the function of many proteins. Here, we test the hypothesis that anionic lipids could indirectly modulate MS proteins by alteration of the lipid bilayer mechanical properties. Using all-atom molecular dynamics simulations, we computed the bilayer bending rigidity (KC), the area compressibility (KA), and the surface shear viscosity (ηm) of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (PC) lipid bilayers with and without phosphatidylserine (PS) or phosphatidylinositol bisphosphate (PIP2) at physiological concentrations in the lower leaflet. Tensionless leaflets were first checked for each asymmetric bilayer model, and a formula for embedding an asymmetric channel in an asymmetric bilayer is proposed. Results from two different sized bilayers show consistently that the addition of 20% surface charge in the lower leaflet of the PC bilayer with PIP2 has minimal impact on its mechanical properties, while PS reduced the bilayer bending rigidity by 22%. As a comparison, supplementing the PIP2-enriched PC membrane with 30% cholesterol, a known rigidifying steroid lipid, produces a significant increase in all three mechanical constants. Analysis of pairwise splay moduli suggests that the effect of anionic lipids on bilayer bending rigidity largely depends on the number of anionic lipid pairs formed during simulations. The potential implication of bilayer bending rigidity is discussed in the framework of MS piezo channels.

Structure and Dynamics of GPCRs in Lipid Membranes: Physical Principles and Experimental Approaches

Over the past decade, the vast amount of information generated through structural and biophysical studies of GPCRs has provided unprecedented mechanistic insight into the complex signalling behaviour of these receptors. With this recent information surge, it has also become increasingly apparent that in order to reproduce the various effects that lipids and membranes exert on the biological function for these allosteric receptors, in vitro studies of GPCRs need to be conducted under conditions that adequately approximate the native lipid bilayer environment. In the first part of this review, we assess some of the more general effects that a membrane environment exerts on lipid bilayer-embedded proteins such as GPCRs. This is then followed by the consideration of more specific effects, including stoichiometric interactions with specific lipid subtypes. In the final section, we survey a range of different membrane mimetics that are currently used for in vitro studies, with a focus on NMR applications.

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.

Physical properties of model biological lipid bilayers: insights from all-atom molecular dynamics simulations

The physical properties of lipid bilayers are sensitive to the specific type and composition of the lipids that make up the many different types of cell membranes. Studying model bilayers of representative heterogeneous compositions can provide key insights into membrane functionality. In this work, we use atomistic molecular dynamics simulations to characterize key properties in a number of bilayer membranes of varying composition. We first examine basic properties, such as lipid area, volume, and bilayer thickness, of simple, homogeneous bilayers comprising several lipid types, which are prevalent in biological membranes. Such lipids are then used in simulations of heterogeneous systems representative of bacterial, mammalian, and cancer membranes. Our analysis is especially focused on depth-dependent, transmembrane profiles; in particular, we calculate lateral pressure and dipole potential profiles, two fundamental properties which play key roles in a large number of biological functions.

Influence of membrane lipid composition on the structure and activity of γ-secretase

γ-Secretase (GS) is a multi-subunit membrane-embedded aspartyl protease that cleaves more than 80 integral membrane proteins, including the amyloid precursor protein (APP) to produce the amyloid-β (Aβ) peptide. Oligomerization and aggregation of the 42-amino acid length Aβ isoform in the brain has been associated with the development and progression of Alzheimer's disease (AD). Based on recent experimental structural studies and using multiscale computational modeling approaches, the conformational states and protein-membrane interactions of the GS complex embedded in six homogeneous and six heterogeneous lipid bilayers were characterized. In order to identify potential lipid and cholesterol binding sites, GS regions with high lipid/cholesterol occupancy values were analyzed using atomistic and coarse-grained simulations. Long lipid residence times were observed to be correlated with a large number of hydrogen bonds between the charged headgroups and key GS amino acids. This observation provides a plausible explanation for the inhibition of GS by charged lipids observed in previous experimental studies. Computed lateral pressure profiles suggest that higher transmembrane pressures favor active state conformations of the catalytic subunit. A probable mechanism for the regulation of the local stress response in cholesterol-rich multicomponent lipid bilayers is identified. Finally, it is demonstrated that interactions between the nicastrin extracellular domain and lipid headgroups leads to a compact structural conformation of the GS complex. Overall, this study provides valuable insight into the effect of bilayer lipid composition on the GS structural ensemble and its function.

Effects of lipid composition on membrane permeation

Passive permeation through lipid membranes is an essential process in biology. In vivo membranes typically consist of mixtures of lamellar and nonlamellar lipids. Lamellar lipids are characterized by their tendency to form lamellar sheet-like structures, which are predominant in nature. Nonlamellar lipids, when isolated, instead form more geometrically complex nonlamellar phases. While mixed lamellar/nonlamellar lipid membranes tend to adopt the ubiquitous lamellar bilayer structure, the presence of nonlamellar lipids is known to have profound effects on key membrane properties, such as internal distributions of stress and elastic properties, which in turn may alter related biological processes. This work focuses on one such process, i.e., permeation, by utilising atomistic molecular dynamics simulations in order to obtain transfer free energy profiles, diffusion profiles and permeation coefficients for a series of thirteen small molecules and drugs. Each permeant is tested on two bilayer membranes of different lipid composition, i.e., purely lamellar and mixed lamellar/nonlamellar. Our results indicate that the presence of nonlamellar lipids reduces permeation for smaller molecules (molecular weight < 100) but facilitates it for the largest ones (molecular weight > 100). This work represents an advancement towards the development of more realistic in silico permeability assays, which may have a substantial future impact in the area of rational drug design.

Influence of lipid composition of model membranes on methacrylate antimicrobial polymer-membrane interactions

Using atomistic molecular dynamics simulations, the role of lipid composition in the interactions of multiple methacrylate antimicrobial polymer agents with model membranes, and the consequent response of the membranes is studied. In our earlier study, methacrylate polymers were observed to induce phase demixing and associated thickness mismatch in a POPE-POPG model microbial membrane. In this work, we probe (1) the role of varying the degree of saturation in lipid acyl chains in the membrane interactions of methacrylate polymers, and (2) whether electrostatics (addition of anionic lipids) can influence the interactions of the polymers with model mammalian membranes. Lipid composition is observed to significantly modify membrane-polymer interactions, leading to differences in both the mode of partitioning and the conformations adopted by the polymers, in addition to impacting membrane properties differently. The results strongly suggest that the oft-cited electrostatic interactions between the antimicrobial agents and the microbial membranes do not fully account for the recognition and subsequent partitioning of the antimicrobial agents. The ability of the methacrylate polymers to sense interfacial lipid packing defects, determined by the PE/PC head groups of lipids, is also found to be influential in their membrane partitioning. Deliberate inclusion of charged anionic lipids into a model mammalian membrane, leading to additional favorable electrostatics, does not reproduce a similar polymer partitioning mechanism to that in its microbial counterpart. The differences observed in the interactions of methacrylate polymers with the various model membranes can be instrumental in extending our understanding of underlying modes of membrane disruption by general antimicrobial agents as well.

Gramicidin A Channel Formation Induces Local Lipid Redistribution II: A 3D Continuum Elastic Model

To change conformation, a protein must deform the surrounding bilayer. In this work, a three-dimensional continuum elastic model for gramicidin A in a lipid bilayer is shown to describe the sensitivity to thickness, curvature stress, and the mechanical properties of the lipid bilayer. A method is demonstrated to extract the gramicidin-lipid boundary condition from all-atom simulations that can be used in the three-dimensional continuum model. The boundary condition affects the deformation dramatically, potentially much more than typical variations in the material stiffness do as lipid composition is changed. Moreover, it directly controls the sensitivity to curvature stress. The curvature stress and hydrophobic surfaces of the all-atom and continuum models are found to be in excellent agreement. The continuum model is applied to estimate the enrichment of hydrophobically matched lipids near the channel in a mixture, and the results agree with single-channel experiments and extended molecular dynamics simulations from the companion article by Beaven et al. in this issue of Biophysical Journal.

The Influence of Lipid Bilayer Physicochemical Properties on Gramicidin A Conformer Preferences

The conformational preferences adopted by gramicidin A (GA) dimers inserted into phospholipid bilayers are reported as a function of the bilayer cholesterol content, temperature, and incubation time. Through use of vesicle capture-freeze drying methodology, GA dimers were captured in lipid bilayers and the conformational preferences of the complex were analyzed using ion mobility-mass spectrometry. Perturbations that affect the physicochemical interactions in the lipid bilayer such as cholesterol incorporation, temperature, and incubation time directly alter the conformer preferences of the complex. Regardless of bilayer cholesterol concentration, the antiparallel double helix (ADH) conformation was observed to be most abundant for GA dimers in bilayers composed of lipids with 12 to 22 carbon acyl chains. Incorporation of cholesterol into lipid bilayers yields increased bilayer thickness and rigidity, and an increased abundance of parallel double helix (PDH) and single-stranded head-to-head (SSHH) dimers were observed. Bilayers prepared using 1,2-dilauroyl-sn-glycero-3-phosphocholine, a lipid with 12 carbon acyl chains, yielded a nascent conformer that decreased in abundance as a function of bilayer cholesterol content. High resolution ion mobility-mass spectrometry data revealed two peaks in the ADH region suggesting that ADH populations are composed of two distinct conformers. The conformer preferences of GA dimers from 1,2-distearoyl-sn-glycero-3-phosphocholine bilayers were significantly different for samples incubated at 4°C vs. 60°C; increased cholesterol content yielded more PDH and SSHH at 60°C. The addition of cholesterol as well as incubating samples of 1,2-distearoyl-sn-glycero-3-phosphocholine at 60°C for 24-72 h yielded an increase in PDH and SSHH abundance.

The Role of Cholesterol in Driving IAPP-Membrane Interactions

Our knowledge of the molecular events underlying type 2 diabetes mellitus-a protein conformational disease characterized by the aggregation of islet amyloid polypeptide (IAPP) in pancreatic β cells-is limited. However, amyloid-mediated membrane damage is known to play a key role in IAPP cytotoxicity, and therefore the effects of lipid composition on modulating IAPP-membrane interactions have been the focus of intense research. In particular, membrane cholesterol content varies with aging and consequently with adverse environmental factors such as diet and lifestyle, but its role in the development of the disease is controversial. In this study, we employ a combination of experimental techniques and in silico molecular simulations to shed light on the role of cholesterol in IAPP aggregation and the related membrane disruption. We show that if anionic POPC/POPS vesicles are used as model membranes, cholesterol has a negligible effect on the kinetics of IAPP fibril growth on the surface of the bilayer. In addition, cholesterol inhibits membrane damage by amyloid-induced poration on membranes, but enhances leakage through fiber growth on the membrane surface. Conversely, if 1:2 DOPC/DPPC raft-like model membranes are used, cholesterol accelerates fiber growth. Next, it enhances pore formation and suppresses fiber growth on the membrane surface, leading to leakage. Our results highlight a twofold effect of cholesterol on the amyloidogenicity of IAPP and help explain its debated role in type 2 diabetes mellitus.

Membrane stress profiles from self-consistent field theory

Using self-consistent field theory (SCFT), we develop an accurate, local expression for the stress profiles in membranes and soft matter interfaces, in general. The bond stresses are expressed in terms of pre-computed chain propagators, which are used to describe the statistical weight of the molecules and therefore require minimal additional calculations. In addition, we overcome the resolution limit of the molecular bond length by including the Irving and Kirkwood bond assignment and recover a constant normal stress profile across an interface. Using this theory, we find that the membrane lateral stress profile contains repulsive (positive) stresses in the regions of the head and tail groups, and attractive (negative) stresses near the hydrophobic/hydrophilic interface. We also verify that the zeroth and first moments of the stress profile correspond to the thermodynamic tension and product of the bending modulus and the spontaneous curvature, respectively.

Exploring the structural relationship between encapsulated antimicrobial peptides and the bilayer membrane mimetic lipidic cubic phase: studies with gramicidin A′

Lipid based bicontinuous cubic mesophases provide a low-cost, robust membrane mimetic nanomaterial which allows for the incorporation of membrane peptides and proteins.

Effects of Lipid Composition on Bilayer Membranes Quantified by All-Atom Molecular Dynamics

Biological bilayer membranes typically contain varying amounts of lamellar and nonlamellar lipids. Lamellar lipids, such as dioleoylphosphatidylcholine (DOPC), are defined by their tendency to form the lamellar phase, ubiquitous in biology. Nonlamellar lipids, such as dioleoylphosphatidylethanolamine (DOPE), prefer instead to form nonlamellar phases, which are mostly nonbiological. However, nonlamellar lipids mix with lamellar lipids in biomembrane structures that remain overall lamellar. Importantly, changes in the lamellar vs nonlamellar lipid composition are believed to affect membrane function and modulate membrane proteins. In this work, we employ atomistic molecular dynamics simulations to quantify how a range of bilayer properties are altered by variations in the lamellar vs nonlamellar lipid composition. Specifically, we simulate five DOPC/DOPE bilayers at mixing ratios of 1/0, 3/1, 1/1, 1/3, and 0/1. We examine properties including lipid area and bilayer thickness, as well as the transmembrane profiles of electron density, lateral pressure, electric field, and dipole potential. While the bilayer structure is only marginally altered by lipid composition changes, dramatic effects are observed for the lateral pressure, electric field, and dipole potential profiles. Possible implications for membrane function are discussed.

Interaction of multiple biomimetic antimicrobial polymers with model bacterial membranes

Using atomistic molecular dynamics simulations, interaction of multiple synthetic random copolymers based on methacrylates on prototypical bacterial membranes is investigated. The simulations show that the cationic polymers form a micellar aggregate in water phase and the aggregate, when interacting with the bacterial membrane, induces clustering of oppositely charged anionic lipid molecules to form clusters and enhances ordering of lipid chains. The model bacterial membrane, consequently, develops lateral inhomogeneity in membrane thickness profile compared to polymer-free system. The individual polymers in the aggregate are released into the bacterial membrane in a phased manner and the simulations suggest that the most probable location of the partitioned polymers is near the 1-palmitoyl-2-oleoyl-phosphatidylglycerol (POPG) clusters. The partitioned polymers preferentially adopt facially amphiphilic conformations at lipid-water interface, despite lacking intrinsic secondary structures such as α-helix or β-sheet found in naturally occurring antimicrobial peptides.

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.

Atomistic study of lipid membranes containing chloroform: looking for a lipid-mediated mechanism of anesthesia

The molecular mechanism of general anesthesia is still a controversial issue. Direct effect by linking of anesthetics to proteins and indirect action on the lipid membrane properties are the two hypotheses in conflict. Atomistic simulations of different lipid membranes subjected to the effect of small volatile organohalogen compounds are used to explore plausible lipid-mediated mechanisms. Simulations of homogeneous membranes reveal that electrostatic potential and lateral pressure transversal profiles are affected differently by chloroform (anesthetic) and carbon tetrachloride (non-anesthetic). Simulations of structured membranes that combine ordered and disordered regions show that chloroform molecules accumulate preferentially in highly disordered lipid domains, suggesting that the combination of both lateral and transversal partitioning of chloroform in the cell membrane could be responsible of its anesthetic action.

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.

GxxxG motifs, phenylalanine, and cholesterol guide the self-association of transmembrane domains of ErbB2 receptors

GxxxG motifs are common in transmembrane domains of membrane proteins and are often introduced to artificial peptides to inhibit or promote association to stable structures. The transmembrane domain of ErbB2 presents two separate such motifs that are proposed to be connected to stability and activity of the dimer. Using molecular simulations, we show that these sequences play a critical role during the recognition stage, forming transient complexes that lead to stable dimers. In pure phospholipid bilayers association occurs by contacts formed at the C-terminus promoted by the presence of phenylalanine residues. Helices subsequently rotate to eventually pack at short separations favored by lipid entropic contributions. In contrast, at intermediate cholesterol concentrations, a different pathway is followed that involves dimers with a weaker interface toward the N-terminus. However, at high cholesterol content, a switch toward the C-terminus is observed with an overall nonmonotonic change of the dimerization affinity. This conformational switch modulated by cholesterol has important implications on the thermodynamic, structural, and kinetic characteristics of helix-helix association in lipid membranes.

Model answers to lipid membrane questions

Ever since it was discovered that biological membranes have a core of a bimolecular sheet of lipid molecules, lipid bilayers have been a model laboratory for investigating physicochemical and functional properties of biological membranes. Experimental and theoretical models help the experimental scientist to plan experiments and interpret data. Theoretical models are the theoretical scientist's preferred toys to make contact between membrane theory and experiments. Most importantly, models serve to shape our intuition about which membrane questions are the more fundamental and relevant ones to pursue. Here we review some membrane models for lipid self-assembly, monolayers, bilayers, liposomes, and lipid-protein interactions and illustrate how such models can help answering questions in modern lipid cell biology.

Atomistic simulation of lipid and DiI dynamics in membrane bilayers under tension

Membrane tension modulates cellular processes by initiating changes in the dynamics of its molecular constituents. To quantify the precise relationship between tension, structural properties of the membrane, and the dynamics of lipids and a lipophilic reporter dye, we performed atomistic molecular dynamics (MD) simulations of DiI-labeled dipalmitoylphosphatidylcholine (DPPC) lipid bilayers under physiological lateral tensions ranging from -2.6 mN m(-1) to 15.9 mN m(-1). Simulations showed that the bilayer thickness decreased linearly with tension consistent with volume-incompressibility, and this thinning was facilitated by a significant increase in acyl chain interdigitation at the bilayer midplane and spreading of the acyl chains. Tension caused a significant drop in the bilayer's peak electrostatic potential, which correlated with the strong reordering of water and lipid dipoles. For the low tension regime, the DPPC lateral diffusion coefficient increased with increasing tension in accordance with free-area theory. For larger tensions, free area theory broke down due to tension-induced changes in molecular shape and friction. Simulated DiI rotational and lateral diffusion coefficients were lower than those of DPPC but increased with tension in a manner similar to DPPC. Direct correlation of membrane order and viscosity near the DiI chromophore, which was just under the DPPC headgroup, indicated that measured DiI fluorescence lifetime, which is reported to decrease with decreasing lipid order, is likely to be a good reporter of tension-induced decreases in lipid headgroup viscosity. Together, these results offer new molecular-level insights into membrane tension-related mechanotransduction and into the utility of DiI in characterizing tension-induced changes in lipid packing.

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.

On the Validation of Molecular Dynamics Simulations of Saturated and cis-Monounsaturated Phosphatidylcholine Lipid Bilayers: A Comparison with Experiment

Molecular dynamics simulations of fully hydrated pure bilayers of four widely studied phospholipids, 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), and 2-oleoyl-1-palmitoyl-sn-glycero-3-phosphocholine (POPC) using a recent revision of the GROMOS96 force field are reported. It is shown that the force field reproduces the structure and the hydration of bilayers formed by each of the four lipids with high accuracy. Specifically, the solvation and the orientation of the dipole of the phosphocholine headgroup and of the ester carbonyls show that the structure of the primary hydration shell in the simulations closely matches experimental findings. This work highlights the need to reproduce a broad range of properties beyond the area per lipid, which is poorly defined experimentally, and to consider the effect of system size and sampling times well beyond those commonly used.

Screening for the drug-phospholipid interaction: correlation to phospholipidosis

Phospholipid bilayers represent a complex, anisotropic environment fundamentally different from bulk oil or octanol, for instance. Even "simple" drug association to phospholipid bilayers can only be fully understood if the slab-of-hydrocarbon approach is abandoned and the complex, anisotropic properties of lipid bilayers reflecting the chemical structures and organization of the constituent phospholipids are considered. The interactions of drugs with phospholipids are important in various processes, such as drug absorption, tissue distribution, and subcellular distribution. In addition, drug-lipid interactions may lead to changes in lipid-dependent protein activities, and further, to functional and morphological changes in cells, a prominent example being the phospholipidosis (PLD) induced by cationic amphiphilic drugs. Herein we briefly review drug-lipid interactions in general and the significance of these interactions in PLD in particular. We also focus on a potential causal connection between drug-induced PLD and steatohepatitis, which is induced by some cationic amphiphilic drugs.

Membrane trauma and Na+ leak from Nav1.6 channels

During brain trauma, white matter experiences shear and stretch forces that, without severing axons, nevertheless trigger their secondary degeneration. In central nervous system (CNS) trauma models, voltage-gated sodium channel (Nav) blockers are neuroprotective. This, plus the rapid tetrodotoxin-sensitive Ca2+ overload of stretch-traumatized axons, points to "leaky" Nav channels as a pivotal early lesion in brain trauma. Direct effects of mechanical trauma on neuronal Nav channels have not, however, been tested. Here, we monitor immediate responses of recombinant neuronal Nav channels to stretch, using patch-clamp and Na+-dye approaches. Trauma constituted either bleb-inducing aspiration of cell-attached oocyte patches or abrupt uniaxial stretch of cells on an extensible substrate. Nav1.6 channel transient current displayed irreversible hyperpolarizing shifts of steady-state inactivation [availability(V)] and of activation [g(V)] and, thus, of window current. Left shift increased progressively with trauma intensity. For moderately intense patch trauma, a approximately 20-mV hyperpolarizing shift was registered. Nav1.6 voltage sensors evidently see lower energy barriers posttrauma, probably because of the different bilayer mechanics of blebbed versus intact membrane. Na+ dye-loaded human embryonic kidney (HEK) cells stably transfected with alphaNav1.6 were subjected to traumatic brain injury-like stretch. Cytoplasmic Na+ levels abruptly increased and the trauma-induced influx had a significant tetrodotoxin-sensitive component. Nav1.6 channel responses to cell and membrane trauma are therefore consistent with the hypothesis that mechanically induced Nav channel leak is a primary lesion in traumatic brain injury. Nav1.6 is the CNS node of Ranvier Nav isoform. When, during head trauma, nodes experienced bleb-inducing membrane damage of varying intensities, nodal Nav1.6 channels should immediately "leak" over a broadly left-smeared window current range.

The molecular-scale arrangement and mechanical strength of phospholipid/cholesterol mixed bilayers investigated by frequency modulation atomic force microscopy in liquid

Cholesterols play key roles in controlling molecular fluidity in a biological membrane, yet little is known about their molecular-scale arrangements in real space. In this study, we have directly imaged lipid-cholesterol complexes in a model biological membrane consisting of dipalmitoylphosphatidylcholine (DPPC) and cholesterols by frequency modulation atomic force microscopy (FM-AFM) in phosphate buffer solution. FM-AFM images of a DPPC/cholesterol bilayer in the liquid-ordered phase showed higher energy dissipation values compared to those measured on a nanoscale DPPC domain in the gel phase, reflecting the increased molecular fluidity due to the insertion of cholesterols. Molecular-resolution FM-AFM images of a DPPC/cholesterol bilayer revealed the existence of a rhombic molecular arrangement (lattice constants: a = 0.46 nm, b = 0.71 nm) consisting of alternating rows of DPPC and cholesterols as well as the increased defect density and reduced molecular ordering. The mechanical strength of a DPPC/cholesterol bilayer was quantitatively evaluated by measuring a loading force required to penetrate the membrane with an AFM tip. The result revealed the significant decrease of mechanical strength upon insertion of cholesterols. Based on the molecular-scale arrangement found in this study, we propose a model to explain the reduced mechanical strength in relation to the formation of lipid-ion networks.

Solvent-free lipid bilayer model using multiscale coarse-graining

The multiscale coarse-graining (MS-CG) approach developed in our previous work is extended here to model solvent-free lipid bilayers. The water (solvent) molecules are completely integrated out of the coarse-grained (CG) effective force field. The MS-CG potential, a sum of pairwise central terms, accurately approximates the many-body potential of mean force in the coarse-grained coordinates. It thus incorporates both energetic and entropic contributions. To improve the stability and elastic properties of the MS-CG simulated bilayer, an additional constraint was adopted: the partial virial associated with CG bilayer sites was matched to its corresponding atomistic value for each configuration of the system. The resulting solvent-free MS-CG model reproduces a liquid-state lipid bilayer with accurate structural and elastic properties. Finally, the solvent-free MS-CG model is used to simulate a very large, flat bilayer and two liposome geometries, demonstrating its greatly enhanced computational efficiency.