Inclusion-induced bilayer deformations: effects of monolayer equilibrium curvature

Intrinsic Lipid Curvature and Bilayer Elasticity as Regulators of Channel Function: A Comparative Single-Molecule Study

Perturbations in bilayer material properties (thickness, lipid intrinsic curvature and elastic moduli) modulate the free energy difference between different membrane protein conformations, thereby leading to changes in the conformational preferences of bilayer-spanning proteins. To further explore the relative importance of curvature and elasticity in determining the changes in bilayer properties that underlie the modulation of channel function, we investigated how the micelle-forming amphiphiles Triton X-100, reduced Triton X-100 and the HII lipid phase promoter capsaicin modulate the function of alamethicin and gramicidin channels. Whether the amphiphile-induced changes in intrinsic curvature were negative or positive, amphiphile addition increased gramicidin channel appearance rates and lifetimes and stabilized the higher conductance states in alamethicin channels. When the intrinsic curvature was modulated by altering phospholipid head group interactions, however, maneuvers that promote a negative-going curvature stabilized the higher conductance states in alamethicin channels but destabilized gramicidin channels. Using gramicidin channels of different lengths to probe for changes in bilayer elasticity, we found that amphiphile adsorption increases bilayer elasticity, whereas altering head group interactions does not. We draw the following conclusions: first, confirming previous studies, both alamethicin and gramicidin channels are modulated by changes in lipid bilayer material properties, the changes occurring in parallel yet differing dependent on the property that is being changed; second, isolated, negative-going changes in curvature stabilize the higher current levels in alamethicin channels and destabilize gramicidin channels; third, increases in bilayer elasticity stabilize the higher current levels in alamethicin channels and stabilize gramicidin channels; and fourth, the energetic consequences of changes in elasticity tend to dominate over changes in curvature.

Dependence of protein-induced lipid bilayer deformations on protein shape

Membrane proteins typically deform the surrounding lipid bilayer membrane, which can play an important role in the function, regulation, and organization of membrane proteins. Membrane elasticity theory provides a beautiful description of protein-induced lipid bilayer deformations, in which all physical parameters can be directly determined from experiments. While analytic solutions of protein-induced elastic bilayer deformations are most easily developed for proteins with approximately circular cross sections, structural biology has shown that membrane proteins come in a variety of distinct shapes, with often considerable deviations from a circular cross section. We develop here a boundary value method (BVM) that permits the construction of analytic solutions of protein-induced elastic bilayer deformations for protein shapes with arbitrarily large deviations from a circular cross section, for constant as well as variable boundary conditions along the bilayer-protein interface. We apply this BVM to protein-induced lipid bilayer thickness deformations. Our BVM reproduces available analytic solutions for proteins with circular cross section and yields, for proteins with noncircular cross section, excellent agreement with numerical, finite element solutions. On this basis, we formulate a simple analytic approximation of the bilayer thickness deformation energy associated with general protein shapes and show that, for modest deviations from rotational symmetry, this analytic approximation is in good agreement with BVM solutions. Using the BVM, we survey the dependence of protein-induced elastic bilayer thickness deformations on protein shape, and thus explore how the coupling of protein shape and bilayer thickness deformations affects protein oligomerization and transitions in protein conformational state.

Regulation of Gramicidin Channel Function Solely by Changes in Lipid Intrinsic Curvature

Membrane protein function is regulated by the lipid bilayer composition. In many cases the changes in function correlate with changes in the lipid intrinsic curvature (c 0), and c 0 is considered a determinant of protein function. Yet, water-soluble amphiphiles that cause either negative or positive changes in curvature have similar effects on membrane protein function, showing that changes in lipid bilayer properties other than c 0 are important-and may be dominant. To further investigate the mechanisms underlying the bilayer regulation of protein function, we examined how maneuvers that alter phospholipid head groups effective "size"-and thereby c 0-alter gramicidin (gA) channel function. Using dioleoylphospholipids and planar bilayers, we varied the head groups' physical volume and the electrostatic repulsion among head groups (and thus their effective size). When 1,2-dioleyol-sn-glycero-3-phosphocholine (DOPC), was replaced by 1,2-dioleyol-sn-glycero-3-phosphoethanolamine (DOPE) with a smaller head group (causing a more negative c 0), the channel lifetime (τ) is decreased. When the pH of the solution bathing a 1,2-dioleyol-sn-glycero-3-phosphoserine (DOPS) bilayer is decreased from 7 to 3 (causing decreased head group repulsion and a more negative c 0), τ is decreased. When some DOPS head groups are replaced by zwitterionic head groups, τ is similarly decreased. These effects do not depend on the sign of the change in surface charge. In DOPE:DOPC (3:1) bilayers, pH changes from 5→9 to 5→0 (both increasing head group electrostatic repulsion, thereby causing a less negative c 0) both increase τ. Nor do the effects depend on the use of planar, hydrocarbon-containing bilayers, as similar changes were observed in hydrocarbon-free lipid vesicles. Altering the interactions among phospholipid head groups may alter also other bilayer properties such as thickness or elastic moduli. Such changes could be excluded using capacitance measurements and single channel measurements on gA channels of different lengths. We conclude that changes gA channel function caused by changes in head group effective size can be predicted from the expected changes in c 0.

Endothelial Cell Plasma Membrane Biomechanics Mediates Effects of Pro-Inflammatory Factors on Endothelial Mechanosensors: Vicious Circle Formation in Atherogenic Inflammation

Chronic low-grade vascular inflammation and endothelial dysfunction significantly contribute to the pathogenesis of cardiovascular diseases. In endothelial cells (ECs), anti-inflammatory or pro-inflammatory signaling can be induced by different patterns of the fluid shear stress (SS) exerted by blood flow on ECs. Laminar blood flow with high magnitude is anti-inflammatory, while disturbed flow and laminar flow with low magnitude is pro-inflammatory. Endothelial mechanosensors are the key upstream signaling proteins in SS-induced pro- and anti-inflammatory responses. Being transmembrane proteins, mechanosensors, not only experience fluid SS but also become regulated by the biomechanical properties of the lipid bilayer and the cytoskeleton. We review the apparent effects of pro-inflammatory factors (hypoxia, oxidative stress, hypercholesterolemia, and cytokines) on the biomechanics of the lipid bilayer and the cytoskeleton. An analysis of the available data suggests that the formation of a vicious circle may occur, in which pro-inflammatory cytokines enhance and attenuate SS-induced pro-inflammatory and anti-inflammatory signaling, respectively.

Regulation of Antimicrobial Peptide Activity via Tuning Deformation Fields by Membrane-Deforming Inclusions

Antimicrobial peptides (AMPs) are considered prospective antibiotics. Some AMPs fight bacteria via cooperative formation of pores in their plasma membranes. Most AMPs at their working concentrations can induce lysis of eukaryotic cells as well. Gramicidin A (gA) is a peptide, the transmembrane dimers of which form cation-selective channels in membranes. It is highly toxic for mammalians as being majorly hydrophobic gA incorporates and induces leakage of both bacterial and eukaryotic cell membranes. Both pore-forming AMPs and gA deform the membrane. Here we suggest a possible way to reduce the working concentrations of AMPs at the expense of application of highly-selective amplifiers of AMP activity in target membranes. The amplifiers should alter the deformation fields in the membrane in a way favoring the membrane-permeabilizing states. We developed the statistical model that allows describing the effect of membrane-deforming inclusions on the equilibrium between AMP monomers and cooperative membrane-permeabilizing structures. On the example of gA monomer-dimer equilibrium, the model predicts that amphipathic peptides and short transmembrane peptides playing the role of the membrane-deforming inclusions, even in low concentration can substantially increase the lifetime and average number of gA channels.

Membrane-Mediated Interactions Between Protein Inclusions

Integral or peripheral membrane proteins, or protein oligomers often get close to each other on cell membranes and carry out biological tasks in a collective manner. In addition to electrostatic and van der Waals interactions, those proteins also experience membrane-mediated interactions, which may be necessary for their functionality. The membrane-mediated interactions originate from perturbation of lipid membranes by the presence of protein inclusions, and have been the subject of intensive research in membrane biophysics. Here we review both theoretical and numerical studies of such interactions for membrane proteins and for nanoparticles bound to lipid membranes.

Vectorial insertion of a β-helical peptide into membrane: a theoretical study on polytheonamide B

Spontaneous unidirectional, or vectorial, insertion of transmembrane peptides is a fundamental biophysical process for toxin and viral actions. Polytheonamide B (pTB) is a potent cytotoxic peptide with a β^6.3-helical structure. Previous experimental studies revealed that the pTB inserts into the membrane in a vectorial fashion and forms a channel with its single molecular length long enough to span the membrane. Also, molecular dynamics simulation studies demonstrated that the pTB is prefolded in aqueous solution. These are unique features of pTB because most of the peptide toxins form channels through oligomerization of transmembrane helices. Here, we performed all-atom molecular dynamics simulations to examine the dynamic mechanism of the vectorial insertion of pTB, providing underlying elementary processes of the membrane insertion of a prefolded single transmembrane peptide. We find that the insertion of pTB proceeds with only the local lateral compression of the membrane in three successive phases: "landing," "penetration," and "equilibration" phases. The free energy calculations using the replica-exchange umbrella sampling simulations present an energy cost of 4.3 kcal/mol at the membrane surface for the membrane insertion of pTB from bulk water. The trajectories of membrane insertion revealed that the insertion process can occur in two possible pathways, namely "trapped" and "untrapped" insertions; in some cases, pTB is trapped in the upper leaflet during the penetration phase. Our simulations demonstrated the importance of membrane anchoring by the hydrophobic N-terminal blocking group in the landing phase, leading to subsequent vectorial insertion.

Peptide-induced membrane elastic deformations decelerate gramicidin dimer-monomer equilibration

Gramicidin A (gA) is a hydrophobic pentadecapeptide readily incorporating into a planar bilayer lipid membrane (BLM), thereby inducing a large macroscopic current across the BLM. This current results from ion-channel formation due to head-to-head transbilayer dimerization of gA monomers with rapidly established monomer-dimer equilibrium. Any disturbance of the equilibrium, e.g., by sensitized photoinactivation of a portion of gA monomers, causes relaxation toward a new equilibrium state. According to previous studies, the characteristic relaxation time of the gA-mediated electric current decreases as the current increases upon elevating the gA concentration in the membrane. Here, we report data on the current relaxation kinetics for gA analogs with N-terminal valine replaced by glycine or tyrosine. Surprisingly, the relaxation time increased rather than decreased upon elevation of the total membrane conductance induced by these gA analogs, thus contradicting the classical kinetic scheme. We developed a general theoretical model that accounts for lateral interaction of monomers and dimers mediated by membrane elastic deformations. The modified kinetic scheme of the gramicidin dimerization predicts the reverse dependence of the relaxation time on membrane conductance for gA analogs, with a decreased dimerization constant that is in a good agreement with our experimental data. The equilibration process may be also modulated by incorporation of other peptides ("impurities") into the lipid membrane.

Quantitative Analysis of the Membrane Affinity of Local Anesthetics Using a Model Cell Membrane

Local anesthesia is a drug that penetrates the nerve cell membrane and binds to the voltage gate sodium channel, inhibiting the membrane potential and neurotransmission. It is mainly used in clinical uses to address the pain of surgical procedures in the local area. Local anesthetics (LAs), however, can be incorporated into the membrane, reducing the thermal stability of the membrane as well as altering membrane properties such as fluidity, permeability, and lipid packing order. The effects of LAs on the membrane are not yet fully understood, despite a number of previous studies. In particular, it is necessary to analyze which is the more dominant factor, the membrane affinity or the structural perturbation of the membrane. To analyze the effects of LAs on the cell membrane and compare the results with those from model membranes, morphological analysis and 50% inhibitory concentration (IC50) measurement of CCD-1064sk (fibroblast, human skin) membranes were carried out for lidocaine (LDC) and tetracaine (TTC), the most popular LAs in clinical use. Furthermore, the membrane affinity of the LAs was quantitatively analyzed using a colorimetric polydiacetylene assay, where the color shift represents their distribution in the membrane. Further, to confirm the membrane affinity and structural effects of the membranes, we performed an electrophysiological study using a model protein (gramicidin A, gA) and measured the channel lifetime of the model protein on the free-standing lipid bilayer according to the concentration of each LA. Our results show that when LAs interact with cell membranes, membrane affinity is a more dominant factor than steric or conformational effects of the membrane.

Cyclic Activity of an Osmotically Stressed Liposome in a Finite Hypotonic Environment

A lipid vesicle, or simply called a liposome, represents a synthetic compartment for the examination of transmembrane transport and signaling phenomena. Yet, a liposome is always subjected to size and shape fluctuations due to local and global imbalance of internal and external osmotic pressures. Here, we show that an osmotically stressed liposome placed within a hypotonic spherical bath undergoes cyclic dynamics described by a periodic sequence of swelling and relaxation phases. These two phases are interfaced by the appearance of a transient transmembrane pore through which chemical delivery occurs. An analytical model was formulated for the recurrent differential equations that convey the time-dependent swelling phase of a pulsatory liposome during individual cycles. We demonstrate that the time-dependent swelling phases of the last several cycles of a pulsatory liposome are strongly dependent on the size of the external bath. Furthermore, decreasing the size of the hypotonic medium reduces the number of cycles of a pulsatory liposome. Comparisons and contrasts of an infinite hypotonic bath with finite external baths of varying radii are discussed.

Elastic deformations mediate interaction of the raft boundary with membrane inclusions leading to their effective lateral sorting

Liquid-ordered lipid domains represent a lateral inhomogeneity in cellular membranes. These domains have elastic and physicochemical properties different from those of the surrounding membrane. In particular, their thickness exceeds that of the disordered membrane. Thus, elastic deformations arise at the domain boundary in order to compensate for the thickness mismatch. In equilibrium, the deformations lead to an incomplete register of monolayer ordered domains: the elastic energy is minimal if domains in opposing monolayers lie on the top of each other, and their boundaries are laterally shifted by about 3 nm. This configuration introduces a region, composed of one ordered and one disordered monolayers, with an intermediate bilayer thickness. Besides, a jump in a local monolayer curvature takes place in this intermediate region, concentrating here most of the elastic stress. This region can participate in a lateral sorting of membrane inclusions by offering them an optimal bilayer thickness and local curvature conditions. In the present study, we consider the sorting of deformable lipid inclusions, undeformable peripheral and deeply incorporated peptide inclusions, and undeformable transmembrane inclusions of different molecular geometry. With rare exceptions, all types of inclusions have an affinity to the ordered domain boundary as compared to the bulk phases. The optimal lateral distribution of inclusions allows relaxing the elastic stress at the boundary of domains.

Modeling membrane nanotube morphology: the role of heterogeneity in composition and material properties

Membrane nanotubes are dynamic structures that may connect cells over long distances. Nanotubes are typically thin cylindrical tubes, but they may occasionally have a beaded architecture along the tube. In this paper, we study the role of membrane mechanics in governing the architecture of these tubes and show that the formation of bead-like structures along the nanotubes can result from local heterogeneities in the membrane either due to protein aggregation or due to membrane composition. We present numerical results that predict how membrane properties, protein density, and local tension compete to create a phase space that governs the morphology of a nanotube. We also find that there exists a discontinuity in the energy that impedes two beads from fusing. These results suggest that the membrane-protein interaction, membrane composition, and membrane tension closely govern the tube radius, number of beads, and the bead morphology.

Monolayerwise application of linear elasticity theory well describes strongly deformed lipid membranes and the effect of solvent

The theory of elasticity of lipid membranes is used widely to describe processes of cell membrane remodeling. Classically, the functional of a membrane's elastic energy is derived under assumption of small deformations; the membrane is considered as an infinitely thin film. This functional is quadratic on membrane surface curvature, with half of the splay modulus as its proportionality coefficient; it is generally applicable for small deformations only. Any validity of this functional for the regime of strong deformations should be verified experimentally. Recently, research using molecular dynamics simulations challenged the validity of this classic, linear model, i.e. the constancy of the splay modulus for strongly bent membranes. Here we demonstrate that the quadratic energy functional still can be applied for calculation of the elastic energy of strongly deformed membranes without introducing higher order terms with additional elastic moduli, but only if applied separately for each lipid monolayer. For cylindrical membranes, both classic and monolayerwise models yield equally accurate results. For cylindrical deformations we experimentally show that the elastic energy of lipid monolayers is additive: a low molecular weight solvent leads to an approximately twofold decrease in the membrane bending stiffness. Accumulation of solvent molecules in the inner monolayer of a membrane cylinder can explain these results, as the solvent partially prevents lipid molecules from splaying there. Thus, the linear theory of elasticity can be expanded through the range from weak to strong deformations-its simplicity and physical transparency describe various membrane phenomena.

A new approach for investigating the response of lipid membranes to electrocompression by coupling droplet mechanics and membrane biophysics

A new method for quantifying lipid-lipid interactions within biomimetic membranes undergoing electrocompression is demonstrated by coupling droplet mechanics and membrane biophysics. The membrane properties are varied by altering the lipid packing through the introduction of cholesterol. Pendant drop tensiometry is used to measure the lipid monolayer tension at an oil-water interface. Next, two lipid-coated aqueous droplets are manipulated into contact to form a bilayer membrane at their adhered interface. The droplet geometries are captured from two angles to provide accurate measurements of both the membrane area and the contact angle between the adhered droplets. Combining the monolayer tension and contact angle measurements enables estimations of the membrane tension with respect to lipid composition. Then, the membrane is electromechanically compressed using a transmembrane voltage. Electrostatic pressure, membrane tension and the work necessary for bilayer thinning are tracked, and a model is proposed to capture the mechanics of membrane compression. The results highlight that a previously unaccounted for energetic term is produced during compression, potentially reflecting changes in the lateral membrane structure. This residual energy is eliminated in cases with cholesterol mole fractions of 0.2 and higher, suggesting that cholesterol diminishes these adjustments.

Mechanical properties of lipid bilayers: a note on the Poisson ratio

We investigate the Poisson ratio ν of fluid lipid bilayers, i.e., the question how area strains compare to the changes in membrane thickness (or, equivalently, volume) that accompany them. We first examine existing experimental results on the area- and volume compressibility of lipid membranes. Analyzing them within the framework of linear elasticity theory for homogeneous thin fluid sheets leads us to conclude that lipid membrane deformations are to a very good approximation volume-preserving, with a Poisson ratio that is likely about 3% smaller than the common soft matter limit . These results are fully consistent with atomistic simulations of a DOPC membrane at varying amount of applied lateral stress, for which we instead deduce ν by directly comparing area- and volume strains. To assess the problematic assumption of transverse homogeneity, we also define a depth-resolved Poisson ratio ν(z) and determine it through a refined analysis of the same set of simulations. We find that throughout the membrane's thickness, ν(z) is close to the value derived assuming homogeneity, with only minor variations of borderline statistical significance.

Molecular Mechanism for Gramicidin Dimerization and Dissociation in Bilayers of Different Thickness

Membrane protein functions can be altered by subtle changes in the host lipid bilayer physical properties. Gramicidin channels have emerged as a powerful system for elucidating the underlying mechanisms of membrane protein function regulation through changes in bilayer properties, which are reflected in the thermodynamic equilibrium distribution between nonconducting gramicidin monomers and conducting bilayer-spanning dimers. To improve our understanding of how subtle changes in bilayer thickness alter the gramicidin monomer and dimer distributions, we performed extensive atomistic molecular dynamics simulations and fluorescence-quenching experiments on gramicidin A (gA). The free-energy calculations predicted a nonlinear coupling between the bilayer thickness and channel formation. The energetic barrier inhibiting gA channel formation was sharply increased in the thickest bilayer (1,2-dierucoyl-sn-glycero-3-phosphocholine). This prediction was corroborated by experimental results on gramicidin channel activity in bilayers of different thickness. To further explore the mechanism of channel formation, we performed extensive unbiased molecular dynamics simulations, which allowed us to observe spontaneous gA dimer formation in lipid bilayers. The simulations revealed structural rearrangements in the gA subunits and changes in lipid packing, as well as water reorganization, that occur during the dimerization process. Together, the simulations and experiments provide new, to our knowledge, insights into the process and mechanism of gramicidin channel formation, as a prototypical example of the bilayer regulation of membrane protein function.

Quantitative Characterization of Protein-Lipid Interactions by Free Energy Simulation between Binary Bilayers

Using a recently developed binary bilayer system (BBS) consisting of two patches of laterally contacting bilayers, umbrella sampling molecular dynamics (MD) simulations were performed for quantitative characterization of protein-lipid interactions. The BBS is composed of 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) with an embedded model membrane protein, a gramicidin A (gA) channel. The calculated free energy difference for the transfer of a gA channel from DLPC (hydrophobic thickness ≈ 21.5 Å) to DMPC (hydrophobic thickness ≈ 25.5 Å) bilayers, ΔG(DLPC → DMPC), is -2.2 ± 0.7 kcal/mol. This value appears at odds with the traditional view that the hydrophobic length of the gA channel is ∼22 Å. To understand this discrepancy, we first note that recent MD simulations by different groups have shown that lipid bilayer thickness profiles in the vicinity of a gA channel differ qualitatively from the deformation profile predicted from continuum elastic bilayer models. Our MD simulations at low and high gA:lipid molar ratios and different membrane compositions indicate that the gA channel's effective hydrophobic length is ∼26 Å. Using this effective hydrophobic length, ΔG(DLPC → DMPC) determined here is in excellent agreement with predictions based on continuum elastic models (-3.0 to -2.2 kcal/mol) where the bilayer deformation energy is approximated as a harmonic function of the mismatch between the channel's effective hydrophobic length and the hydrophobic thickness of the bilayer. The free energy profile for gA in the BBS includes a barrier at the interface between the two bilayers which can be attributed to the line tension at the interface between two bilayers with different hydrophobic thicknesses. This observation implies that translation of a peptide between two different regions of a cell membrane (such as between the liquid ordered and disordered phases) may include effects of a barrier at the interface in addition to the relative free energies of the species far from the interface. The BBS allows for direct transfer free energy calculations between bilayers without a need of a reference medium, such as bulk water, and thus provides an efficient simulation protocol for the quantitative characterization of protein-lipid interactions at all-atom resolution.

Antidepressants are modifiers of lipid bilayer properties

The two major classes of antidepressants, tricyclic antidepressants (TCAs) and selective serotonin reuptake inhibitors (SSRIs), inhibit neurotransmitter reuptake at synapses. They also have off-target effects on proteins other than neurotransmitter transporters, which may contribute to both desired changes in brain function and the development of side effects. Many proteins modulated by antidepressants are bilayer spanning and coupled to the bilayer through hydrophobic interactions such that the conformational changes underlying their function will perturb the surrounding lipid bilayer, with an energetic cost (ΔG _def) that varies with changes in bilayer properties. Here, we test whether changes in ΔG _def caused by amphiphilic antidepressants partitioning into the bilayer are sufficient to alter membrane protein function. Using gramicidin A (gA) channels to probe whether TCAs and SSRIs alter the bilayer contribution to the free energy difference for the gramicidin monomer⇔dimer equilibrium (representing a well-defined conformational transition), we find that antidepressants alter gA channel activity with varying potency and no stereospecificity but with different effects on bilayer elasticity and intrinsic curvature. Measuring the antidepressant partition coefficients using isothermal titration calorimetry (ITC) or cLogP shows that the bilayer-modifying potency is predicted quite well by the ITC-determined partition coefficients, and channel activity is doubled at an antidepressant/lipid mole ratio of 0.02-0.07. These results suggest a mechanism by which antidepressants could alter the function of diverse membrane proteins by partitioning into cell membranes and thereby altering the bilayer contribution to the energetics of membrane protein conformational changes.

Membrane Elastic Deformations Modulate Gramicidin A Transbilayer Dimerization and Lateral Clustering

Gramicidin A (gA) is a short β-helical peptide known to form conducting channels in lipid membranes because of transbilayer dimerization. The gA conducting dimer, being shorter than the lipid bilayer thickness, deforms the membrane in its vicinity, and the bilayer elastic energy contributes to the gA dimer formation energy. Likewise, membrane incorporation of a gA monomer, which is shorter than the lipid monolayer thickness, creates a void, thereby forcing surrounding lipid molecules to tilt to fill it. The energy of membrane deformation was calculated in the framework of the continuum elasticity theory, taking into account splay, tilt, lateral stretching/compression, Gaussian splay deformations, and external membrane tension. We obtained the interaction energy profiles for two gA monomers located either in the same or in the opposite monolayers. The profiles demonstrated the long-range attraction and short-range repulsion behavior of the monomers resulting from the membrane deformation. Analysis of the profile features revealed conditions under which clusters of gA monomers would not dissipate because of diffusion. The calculated dependence of the dimer formation and decay energy barriers on the membrane elastic properties was in good agreement with the available experimental data and suggested an explanation for a hitherto contentious phenomenon.

Clustering and separation of hydrophobic nanoparticles in lipid bilayer explained by membrane mechanics

Small hydrophobic gold nanoparticles with diameter lower than the membrane thickness can form clusters or uniformly distribute within the hydrophobic core of the bilayer. The coexistence of two stable phases (clustered and dispersed) indicates the energy barrier between nanoparticles. We calculated the distance dependence of the membrane-mediated interaction between two adjacent nanoparticles. In our model we consider two deformation modes: the monolayer bending and the hydroxycarbon chain stretching. Existence of an energy barrier between the clustered and the separated state of nanoparticles was predicted. Variation analysis of the membrane mechanical parameters revealed that the energy barrier between two membrane embedded nanoparticles is mainly the consequence of the bending deformation and not change of the thickness of the bilayer in the vicinity of nanoparticles. It is shown, that the forces between the nanoparticles embedded in the biological membrane could be either attractive or repulsive, depending on the mutual distance between them.

Tuning biomimetic membrane barrier properties by hydrocarbon, cholesterol and polymeric additives

The barrier properties of cellular membranes are increasingly attracting attention as a source of inspiration for designing biomimetic membranes. The broad range of potential technological applications makes the use of lipid and lately also polymeric materials a popular choice for constructing biomimetic membranes, where the barrier properties can be controlled by the composition of the membrane constituent elements. Here we investigate the membrane properties reported by the light-induced proton pumping activity of bacteriorhodopsin (bR) reconstituted in three vesicle systems of different membrane composition. Specifically we quantify how the resulting proton influx and efflux rates are influenced by the membrane composition using a variety of membrane modulators. We demonstrate that by adding hydrocarbons to vesicles with reconstituted bR formed from asolectin lipids the resulting transmembrane proton fluxes changes proportional to the carbon chain length when compared against control. We observe a similar proportionality in single-component 1,2-Dioleoyl-sn-glycero-3-phosphocholine model membranes when using cholesterol. Lastly we investigate the effects of adding the amphiphilic di-block co-polymer polybutadiene-polyethyleneoxide (PB₁₂-PEO₁₀) to phospholipid membranes formed from 1,2-Dioleoyl-sn-glycero-3-phosphocholine, 1,2-Dioleoyl-sn-glycero-3-phosphatidylethanolamine, and 1,2-Dioleoyl-sn-glycero-3-phosphatidylserine. The proton pumping activity of bR (measured as a change in extra-vesicular pH) in mixed lipid/PB₁₂-PEO₁₀ lipid systems is up to six-fold higher compared to that observed for bR containing vesicles made from PB₁₂-PEO₁₀ alone. Interestingly, bR inserts with apparent opposite orientation in pure PB₁₂-PEO₁₀ vesicles as compared to pure lipid vesicles. Addition of equimolar amounts of lipids to PB₁₂-PEO₁₀ results in bR orientation similar to that observed for pure lipids. In conclusion our results show how the barrier properties of the membranes can be controlled by the composition of the membrane. In particular the use of mixed lipid-polymer systems may pave the way for constructing biomimetic membranes tailored for optimal properties in various applications including drug delivery systems, biosensors and energy conservation technology.

New Continuum Approaches for Determining Protein-Induced Membrane Deformations

The influence of the membrane on transmembrane proteins is central to a number of biological phenomena, notably the gating of stretch activated ion channels. Conversely, membrane proteins can influence the bilayer, leading to the stabilization of particular membrane shapes, topological changes that occur during vesicle fission and fusion, and shape-dependent protein aggregation. Continuum elastic models of the membrane have been widely used to study protein-membrane interactions. These mathematical approaches produce physically interpretable membrane shapes, energy estimates for the cost of deformation, and a snapshot of the equilibrium configuration. Moreover, elastic models are much less computationally demanding than fully atomistic and coarse-grained simulation methodologies; however, it has been argued that continuum models cannot reproduce the distortions observed in fully atomistic molecular dynamics simulations. We suggest that this failure can be overcome by using chemically and geometrically accurate representations of the protein. Here, we present a fast and reliable hybrid continuum-atomistic model that couples the protein to the membrane. We show that the model is in excellent agreement with fully atomistic simulations of the ion channel gramicidin embedded in a POPC membrane. Our continuum calculations not only reproduce the membrane distortions produced by the channel but also accurately determine the channel's orientation. Finally, we use our method to investigate the role of membrane bending around the charged voltage sensors of the transient receptor potential cation channel TRPV1. We find that membrane deformation significantly stabilizes the energy of insertion of TRPV1 by exposing charged residues on the S4 segment to solution.

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.

Membrane-mediated action of the endocannabinoid anandamide on membrane proteins: implications for understanding the receptor-independent mechanism

Endocannabinoids are amphiphilic molecules that play crucial neurophysiological functions acting as lipid messengers. Antagonists and knockdown of the classical CB1 and CB2 cannabinoid receptors do not completely abolish many endocannabinoid activities, supporting the idea of a mechanism independent of receptors whose mode of action remains unclear. Here we combine gramicidin A (gA) single channel recordings and membrane capacitance measurements to investigate the lipid bilayer-modifying activity of endocannabinoids. Single channel recordings show that the incorporation of endocannabinoids into lipid bilayers reduces the free energy necessary for gramicidin channels to transit from the monomeric to the dimeric conformation. Membrane capacitance demonstrates that the endocannabinoid anandamide has limited effects on the overall structure of the lipid bilayers. Our results associated with the theory of membrane elastic deformation reveal that the action of endocannabinoids on membrane proteins can involve local adjustments of the lipid/protein hydrophobic interface. The current findings shed new light on the receptor-independent mode of action of endocannabinoids on membrane proteins, with important implications towards their neurobiological function.

Revisiting the curvature-mediated interactions between proteins in biological membranes

Proteins embedded in soft biological membranes experience a long-range force mediated by elastic curvature deformations. The classical linearized Helfrich-Canham Hamiltonian based derivations reveal the nature of the force between a pair of proteins to be repulsive in the zero-temperature limit and the interaction potential is inversely proportional to the fourth power of the distance separating the inclusions. Such a result is the starting point to understand many-body interactions between proteins in biological membranes and the study of their clustering or, more broadly, self-organization. A key observation regarding this widely quoted result is that any two (mechanically rigid) proteins will experience an identical force. In other words, there is no specificity in the currently employed continuum models that purport to explain protein interactions. In this work we argue that each protein has a unique mechanical signature based on its interaction with the surrounding lipid bilayer membrane and cannot be treated as a non-specific rigid object. We modify the classical Helfrich-Canham theory of curvature elasticity to incorporate protein-membrane specificity, discuss the estimation of the new model parameters via atomistic simulations and re-evaluate the curvature-mediated force between proteins. We find that the incorporation of protein-specificity can reduce the interaction force by several orders of magnitude. Our result may provide at least one plausible reason behind why in some computational and experimental studies, a net attractive force between proteins is in evidence.

Continuum descriptions of membranes and their interaction with proteins: Towards chemically accurate models

Biological membranes deform in response to resident proteins leading to a coupling between membrane shape and protein localization. Additionally, the membrane influences the function of membrane proteins. Here we review contributions to this field from continuum elastic membrane models focusing on the class of models that couple the protein to the membrane. While it has been argued that continuum models cannot reproduce the distortions observed in fully-atomistic molecular dynamics simulations, we suggest that this failure can be overcome by using chemically accurate representations of the protein. We outline our recent advances along these lines with our hybrid continuum-atomistic model, and we show the model is in excellent agreement with fully-atomistic simulations of the nhTMEM16 lipid scramblase. We believe that the speed and accuracy of continuum-atomistic methodologies will make it possible to simulate large scale, slow biological processes, such as membrane morphological changes, that are currently beyond the scope of other computational approaches. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Bilayer-thickness-mediated interactions between integral membrane proteins

Hydrophobic thickness mismatch between integral membrane proteins and the surrounding lipid bilayer can produce lipid bilayer thickness deformations. Experiment and theory have shown that protein-induced lipid bilayer thickness deformations can yield energetically favorable bilayer-mediated interactions between integral membrane proteins, and large-scale organization of integral membrane proteins into protein clusters in cell membranes. Within the continuum elasticity theory of membranes, the energy cost of protein-induced bilayer thickness deformations can be captured by considering compression and expansion of the bilayer hydrophobic core, membrane tension, and bilayer bending, resulting in biharmonic equilibrium equations describing the shape of lipid bilayers for a given set of bilayer-protein boundary conditions. Here we develop a combined analytic and numerical methodology for the solution of the equilibrium elastic equations associated with protein-induced lipid bilayer deformations. Our methodology allows accurate prediction of thickness-mediated protein interactions for arbitrary protein symmetries at arbitrary protein separations and relative orientations. We provide exact analytic solutions for cylindrical integral membrane proteins with constant and varying hydrophobic thickness, and develop perturbative analytic solutions for noncylindrical protein shapes. We complement these analytic solutions, and assess their accuracy, by developing both finite element and finite difference numerical solution schemes. We provide error estimates of our numerical solution schemes and systematically assess their convergence properties. Taken together, the work presented here puts into place an analytic and numerical framework which allows calculation of bilayer-mediated elastic interactions between integral membrane proteins for the complicated protein shapes suggested by structural biology and at the small protein separations most relevant for the crowded membrane environments provided by living cells.

A general mechanism for drug promiscuity: Studies with amiodarone and other antiarrhythmics

Amiodarone is a widely prescribed antiarrhythmic drug used to treat the most prevalent type of arrhythmia, atrial fibrillation (AF). At therapeutic concentrations, amiodarone alters the function of many diverse membrane proteins, which results in complex therapeutic and toxicity profiles. Other antiarrhythmics, such as dronedarone, similarly alter the function of multiple membrane proteins, suggesting that a multipronged mechanism may be beneficial for treating AF, but raising questions about how these antiarrhythmics regulate a diverse range of membrane proteins at similar concentrations. One possible mechanism is that these molecules regulate membrane protein function by altering the common environment provided by the host lipid bilayer. We took advantage of the gramicidin (gA) channels' sensitivity to changes in bilayer properties to determine whether commonly used antiarrhythmics--amiodarone, dronedarone, propranolol, and pindolol, whose pharmacological modes of action range from multi-target to specific--perturb lipid bilayer properties at therapeutic concentrations. Using a gA-based fluorescence assay, we found that amiodarone and dronedarone are potent bilayer modifiers at therapeutic concentrations; propranolol alters bilayer properties only at supratherapeutic concentration, and pindolol has little effect. Using single-channel electrophysiology, we found that amiodarone and dronedarone, but not propranolol or pindolol, increase bilayer elasticity. The overlap between therapeutic and bilayer-altering concentrations, which is observed also using plasma membrane-like lipid mixtures, underscores the need to explore the role of the bilayer in therapeutic as well as toxic effects of antiarrhythmic agents.

Bilayer Effects of Antimalarial Compounds

Because of the perpetual development of resistance to current therapies for malaria, the Medicines for Malaria Venture developed the Malaria Box to facilitate the drug development process. We tested the 80 most potent compounds from the box for bilayer-mediated effects on membrane protein conformational changes (a measure of likely toxicity) in a gramicidin-based stopped flow fluorescence assay. Among the Malaria Box compounds tested, four compounds altered membrane properties (p< 0.05); MMV007384 stood out as a potent bilayer-perturbing compound that is toxic in many cell-based assays, suggesting that testing for membrane perturbation could help identify toxic compounds. In any case, MMV007384 should be approached with caution, if at all.

Cell interaction with graphene microsheets: near-orthogonal cutting versus parallel attachment

Recent experiments indicate that graphene microsheets can either undergo a near-orthogonal cutting or a parallel attachment mode of interaction with cell membranes. Here we perform a theoretical analysis to characterize the deformed membrane microstructure and investigate how these two interaction modes are influenced by the splay, tilt, compression, tension, bending and adhesion energies of the membrane. Our analysis indicates that, driven by the membrane splay and tension energies, a two-dimensional microsheet such as graphene would adopt a near-perpendicular configuration with respect to the membrane in the transmembrane penetration mode, whereas the membrane bending and tension energies would lead to parallel attachment in the absence of cross membrane penetration. These interaction modes may have broad implications in applications involving drug delivery, cell encapsulation and protection, and the measurement of the dynamic cell response.

Volatile anesthetics inhibit sodium channels without altering bulk lipid bilayer properties

Although general anesthetics are clinically important and widely used, their molecular mechanisms of action remain poorly understood. Volatile anesthetics such as isoflurane (ISO) are thought to alter neuronal function by depressing excitatory and facilitating inhibitory neurotransmission through direct interactions with specific protein targets, including voltage-gated sodium channels (Na(v)). Many anesthetics alter lipid bilayer properties, suggesting that ion channel function might also be altered indirectly through effects on the lipid bilayer. We compared the effects of ISO and of a series of fluorobenzene (FB) model volatile anesthetics on Na(v) function and lipid bilayer properties. We examined the effects of these agents on Na(v) in neuronal cells using whole-cell electrophysiology, and on lipid bilayer properties using a gramicidin-based fluorescence assay, which is a functional assay for detecting changes in lipid bilayer properties sensed by a bilayer-spanning ion channel. At clinically relevant concentrations (defined by the minimum alveolar concentration), both the FBs and ISO produced prepulse-dependent inhibition of Na(v) and shifted the voltage dependence of inactivation toward more hyperpolarized potentials without affecting lipid bilayer properties, as sensed by gramicidin channels. Only at supra-anesthetic (toxic) concentrations did ISO alter lipid bilayer properties. These results suggest that clinically relevant concentrations of volatile anesthetics alter Na(v) function through direct interactions with the channel protein with little, if any, contribution from changes in bulk lipid bilayer properties. Our findings further suggest that changes in lipid bilayer properties are not involved in clinical anesthesia.

Regulation of ion channel function by the host lipid bilayer examined by a stopped-flow spectrofluorometric assay

To examine the function of ligand-gated ion channels in a defined membrane environment, we developed a robust sequential-mixing fluorescence-based stopped-flow assay. Channel activity is determined using a channel-permeable quencher (e.g., thallium, Tl(+)) of a water-soluble fluorophore (8-aminonaphthalene-1,3,6-trisulfonic acid) encapsulated in large unilamellar vesicles in which the channel of interest has been reconstituted, which allows for rapid solution changes. To validate the method, we explored the activation of wild-type KcsA channel, as well as it's noninactivating (E71A) KcsA mutant, by extravesicular protons (H(+)). For both channel types, the day-to-day variability in the reconstitution yield (as judged from the time course of fluorescence quenching) is <10%. The activation curve for E71A KcsA is similar to that obtained previously using single-channel electrophysiology, and the activation curves for wild-type and E71A KcsA are indistinguishable, indicating that channel activation and inactivation are separate processes. We then investigated the regulation of KcsA activation by changes in lipid bilayer composition. Increasing the acyl chain length (from C18:1 to C22:1 in diacylphosphatidylcholine), but not the mole fraction of POPG (>0.25) in the bilayer-forming phospholipid mixture, alters KcsA H(+) gating. The bilayer-thickness-dependent shift in the activation curve is suggestive of a decrease in an apparent H(+) affinity and cooperativity. The control over bilayer environment and time resolution makes this method a powerful assay for exploring ligand activation and inactivation of ion channels, and how channel gating varies with changes in the channels' lipid bilayer environment or other regulatory processes.

Exploring the local elastic properties of bilayer membranes using molecular dynamics simulations

Membrane mechanical elastic properties regulate a variety of cellular processes involving local membrane deformation, such as ion channel function and vesicle fusion. In this work, we used molecular dynamics simulations to estimate the local elastic properties of a membrane. For this, we calculated the energy needed to extract a DOPE lipid molecule, modified with a linker chain, from a POPC bilayer membrane using the umbrella sampling technique. Although the extraction energy entails several contributions related not only to elastic deformation but also to solvation, careful analysis of the potential of mean force (PMF) allowed us to dissect the elastic contribution. With this information, we calculated an effective linear spring constant of 44 ± 4 kJ·nm(-2)·mol(-1) for the DOPC membrane, in agreement with experimental estimates. The membrane deformation profile was determined independently during the stretching process in molecular detail, allowing us to fit this profile to a previously proposed continuum elastic model. Through this approach, we calculated an effective membrane spring constant of 42 kJ·nm(-2)·mol(-1), which is in good agreement with the PMF calculation. Furthermore, the solvation energy we derived from the data is shown to match the solvation energy estimated from critical micelle formation constants. This methodology can be used to determine how changes in lipid composition or the presence of membrane modifiers can affect the elastic properties of a membrane at a local level.

Phytochemicals perturb membranes and promiscuously alter protein function

A wide variety of phytochemicals are consumed for their perceived health benefits. Many of these phytochemicals have been found to alter numerous cell functions, but the mechanisms underlying their biological activity tend to be poorly understood. Phenolic phytochemicals are particularly promiscuous modifiers of membrane protein function, suggesting that some of their actions may be due to a common, membrane bilayer-mediated mechanism. To test whether bilayer perturbation may underlie this diversity of actions, we examined five bioactive phenols reported to have medicinal value: capsaicin from chili peppers, curcumin from turmeric, EGCG from green tea, genistein from soybeans, and resveratrol from grapes. We find that each of these widely consumed phytochemicals alters lipid bilayer properties and the function of diverse membrane proteins. Molecular dynamics simulations show that these phytochemicals modify bilayer properties by localizing to the bilayer/solution interface. Bilayer-modifying propensity was verified using a gramicidin-based assay, and indiscriminate modulation of membrane protein function was demonstrated using four proteins: membrane-anchored metalloproteases, mechanosensitive ion channels, and voltage-dependent potassium and sodium channels. Each protein exhibited similar responses to multiple phytochemicals, consistent with a common, bilayer-mediated mechanism. Our results suggest that many effects of amphiphilic phytochemicals are due to cell membrane perturbations, rather than specific protein binding.

Membrane tension influences the spike propagation between voltage-gated ion channel clusters of excitable membranes

Ion channels of excitable membranes are known to be sensitive to various kinds of stimuli, but the case of simultaneous occurrence of different stimuli is poorly understood. Here, we theoretically analyze the influence of membrane tension on the dynamics of voltage-gated ion channels of excitable membranes. To do so, we develop a modification of the well-known Hodgkin-Huxley model to study numerically the spike generation and propagation in a single and two coupled excitable cells. We find that these cells can use membrane tension to trigger sub-threshold spike propagation, to suppress spike propagation and to alter the intensity of the signal transmission. These effects indicate that cells could use membrane tension to regulate cell-to-cell communication.

Simulating G protein-coupled receptors in native-like membranes: from monomers to oligomers

G protein-coupled receptors (GPCRs) are one of the most relevant superfamilies of transmembrane proteins as they participate in an important variety of biological events. Recently, the scientific community is witnessing an advent of a GPCR crystallization age along with impressive improvements achieved in the field of computer simulations during the last two decades. Computer simulation techniques such as molecular dynamics (MD) simulations are now frequent tools to study the dynamic behavior of GPCRs and, more importantly, to model the complex membrane environment where these proteins spend their lifetime. Thanks to these tools, GPCRs can be simulated not only longer but also in a more "physiological" fashion. In this scenario, scientists are taking advantage of such advances to approach certain phenomena such as GPCR oligomerization occurring only at timescales not reachable until now. Thus, despite current MD simulations having important limitations today, they have become an essential tool to study key biophysical properties of GPCRs and GPCR oligomers.

Regulation of channel function due to physical energetic coupling with a lipid bilayer

Regulation of membrane protein functions due to hydrophobic coupling with a lipid bilayer has been investigated. An energy formula describing interactions between lipid bilayer and integral ion channels with different structures, which is based on the screened Coulomb interaction approximation, has been developed. Here the interaction energy is represented as being due to charge-based interactions between channel and lipid bilayer. The hydrophobic bilayer thickness channel length mismatch is found to induce channel destabilization exponentially while negative lipid curvature linearly. Experimental parameters related to channel dynamics are consistent with theoretical predictions. To measure comparable energy parameters directly in the system and to elucidate the mechanism at an atomistic level we performed molecular dynamics (MD) simulations of the ion channel forming peptide-lipid complexes. MD simulations indicate that peptides and lipids experience electrostatic and van der Waals interactions for short period of time when found within each other's proximity. The energies from these two interactions are found to be similar to the energies derived theoretically using the screened Coulomb and the van der Waals interactions between peptides (in ion channel) and lipids (in lipid bilayer) due to mainly their charge properties. The results of in silico MD studies taken together with experimental observable parameters and theoretical energetic predictions suggest that the peptides induce ion channels inside lipid membranes due to peptide-lipid physical interactions. This study provides a new insight helping better understand of the underlying mechanisms of membrane protein functions in cell membrane leading to important biological implications.

Small-molecule photostabilizing agents are modifiers of lipid bilayer properties

Small-molecule photostabilizing or protective agents (PAs) provide essential support for the stability demands on fluorescent dyes in single-molecule spectroscopy and fluorescence microscopy. These agents are employed also in studies of cell membranes and model systems mimicking lipid bilayer environments, but there is little information about their possible effects on membrane structure and physical properties. Given the impact of amphipathic small molecules on bilayer properties such as elasticity and intrinsic curvature, we investigated the effects of six commonly used PAs--cyclooctatetraene (COT), para-nitrobenzyl alcohol (NBA), Trolox (TX), 1,4-diazabicyclo[2.2.2]octane (DABCO), para-nitrobenzoic acid (pNBA), and n-propyl gallate (nPG)--on bilayer properties using a gramicidin A (gA)-based fluorescence quench assay to probe for PA-induced changes in the gramicidin monomer↔dimer equilibrium. The experiments were done using fluorophore-loaded large unilamellar vesicles that had been doped with gA, and changes in the gA monomer↔dimer equilibrium were assayed using a gA channel-permeable fluorescence quencher (Tl⁺). Changes in bilayer properties caused by, e.g., PA adsorption at the bilayer/solution interface that alter the equilibrium constant for gA channel formation, and thus the number of conducting gA channels in the large unilamellar vesicle membrane, will be detectable as changes in the rate of Tl⁺ influx-the fluorescence quench rate. Over the experimentally relevant millimolar concentration range, TX, NBA, and pNBA, caused comparable increases in gA channel activity. COT, also in the millimolar range, caused a slight decrease in gA channel activity. nPG increased channel activity at submillimolar concentrations. DABCO did not alter gA activity. Five of the six tested PAs thus alter lipid bilayer properties at experimentally relevant concentrations, which becomes important for the design and analysis of fluorescence studies in cells and model membrane systems. We therefore tested combinations of COT, NBA, and TX; the combinations altered the fluorescence quench rate less than would be predicted assuming their effects on bilayer properties were additive. The combination of equimolar concentrations of COT and NBA caused minimal changes in the fluorescence quench rate.

Phosphoinositides alter lipid bilayer properties

Phosphatidylinositol-4,5-bisphosphate (PIP₂), which constitutes ∼1% of the plasma membrane phospholipid, plays a key role in membrane-delimited signaling. PIP₂ regulates structurally and functionally diverse membrane proteins, including voltage- and ligand-gated ion channels, inwardly rectifying ion channels, transporters, and receptors. In some cases, the regulation is known to involve specific lipid–protein interactions, but the mechanisms by which PIP₂ regulates many of its various targets remain to be fully elucidated. Because many PIP₂ targets are membrane-spanning proteins, we explored whether the phosphoinositides might alter bilayer physical properties such as curvature and elasticity, which would alter the equilibrium between membrane protein conformational states—and thereby protein function. Taking advantage of the gramicidin A (gA) channels’ sensitivity to changes in lipid bilayer properties, we used gA-based fluorescence quenching and single-channel assays to examine the effects of long-chain PIP₂s (brain PIP₂, which is predominantly 1-stearyl-2-arachidonyl-PIP₂, and dioleoyl-PIP₂) on bilayer properties. When premixed with dioleoyl-phosphocholine at 2 mol %, both long-chain PIP₂s produced similar changes in gA channel function (bilayer properties); when applied through the aqueous solution, however, brain PIP₂ was a more potent modifier than dioleoyl-PIP₂. Given the widespread use of short-chain dioctanoyl-phosphoinositides, we also examined the effects of diC₈-phosphoinositol (PI), PI(4,5)P₂, PI(3,5)P₂, PI(3,4)P₂, and PI(3,4,5)P₃. The diC₈ phosphoinositides, except for PI(3,5)P₂, altered bilayer properties with potencies that decreased with increasing head group charge. Nonphosphoinositide diC₈ phospholipids generally were more potent bilayer modifiers than the polyphosphoinositides. These results show that physiological increases or decreases in plasma membrane PIP₂ levels, as a result of activation of PI kinases or phosphatases, are likely to alter lipid bilayer properties, in addition to any other effects they may have. The results further show that exogenous PIP₂, as well as structural analogues that differ in acyl chain length or phosphorylation state, alters lipid bilayer properties at the concentrations used in many cell physiological experiments.

Dynamic Heterogeneous Dielectric Generalized Born (DHDGB): An implicit membrane model with a dynamically varying bilayer thickness

An extension to the heterogeneous dielectric generalized Born (HDGB) implicit membrane formalism is presented to allow dynamic membrane deformations in response to membrane-inserted biomolecules during molecular dynamic simulations. The flexible membrane is implemented through additional degrees of freedom that represent the membrane deformation at the contact points of a membrane-inserted solute with the membrane. The extra degrees of freedom determine the dielectric and non-polar solvation free energy profiles that are used to obtain the solvation free energy in the presence of the membrane and are used to calculate membrane deformation free energies according to an elastic membrane model. With the dynamic HDGB (DHDGB) model the membrane is able to deform in response to the insertion of charged molecules thereby avoiding the overestimation of insertion free energies with static membrane models. The DHDGB model also allows the membrane to respond to the insertion of membrane-spanning solutes with hydrophobic mismatch. The model is tested with the membrane insertion of amino acid side chain analogs, arginine-containing helices, the WALP23 peptide, and the gramicidin A channel.

Effects of proteoliposome composition and draw solution types on separation performance of aquaporin-based proteoliposomes: implications for seawater desalination using aquaporin-based biomimetic membranes

Aquaporins are a large family of water transport proteins in cell membranes. Their high water permeability and solute rejection make them potential building blocks for high-performance biomimetic membranes for desalination. In the current study, proteoliposomes were prepared using AquaporinZ from Escherichia coli cells, and their separation properties were characterized by stopped-flow measurements. The current study systematically investigated the effect of proteoliposome composition (lipid type, protein-to-lipid ratio (PLR), and the addition of cholesterol) on water permeability and NaCl retention. Among the various lipids investigated, 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC)-based proteoliposomes were found to have excellent osmotic water permeability and NaCl reflection coefficient values. Increasing the PLR of DOPC proteoliposomes up to 1:200 increased their osmotic water permeability. However, further increase in the PLR reduced the osmotic water permeability probably due to the occurrence of defects in the proteoliposomes, whereas the addition of cholesterol improved their osmotic water permeation likely due to defects sealing. The current study also investigated the effect of major dissolved ions in seawater (e.g., Mg(2+) and SO(4)(2-)) on the stability of proteoliposomes, and design criteria for aquaporin-based biomimetic membranes are proposed in the context of desalination.

UV-visible and infrared methods for investigating lipid-rhodopsin membrane interactions

We describe experimental UV-visible and Fourier transform infrared (FTIR) spectroscopic methods for characterizing lipid-protein interactions for rhodopsin in a membrane bilayer environment. The combination of FTIR and UV-visible difference spectroscopy is used to monitor the structural and functional changes during rhodopsin activation. Investigations of how membrane lipids stabilize various rhodopsin photoproducts are analogous to mutating the protein in terms of gain or loss of function. Interpretation of the results entails a flexible surface model for explaining membrane lipid-protein interactions through material properties relevant to biological activity.

Assessing smectic liquid-crystal continuum models for elastic bilayer deformations

For four decades, since W. Helfrich's pioneering study of smectic A liquid crystals in 1973, continuum elastic models (CEMs) have been employed as tools to understand the energetics of protein-induced lipid bilayer deformations. Among the assumptions underlying this use is that all relevant protein-lipid interactions can be included in the continuum representation of the protein-bilayer interactions through the physical parameters determined for protein-free bilayers and the choice of boundary conditions at the protein/bilayer interface. To better understand this assumption, we review the general structure of CEMs, examine how different choices of boundary conditions and physical moduli profiles alter the predicted bilayer thickness profiles around gramicidin A (gA) and mitochondrial voltage-dependent anion channels (VDAC), respectively, and compare these profiles with those obtained from all-atom molecular dynamics simulations. We find that the profiles differ qualitatively in the first lipid shell around the channels, indicating that the CEMs do not capture accurately the consequences of the protein-induced local changes in lipid bilayer dynamics. Therefore, one needs to be careful when interpreting the results of CEM-based analyses of lipid bilayer-membrane protein interactions.

Rhodopsin-lipid interactions studied by NMR

The biophysical properties of the lipid matrix are known to influence function of integral membrane proteins. We report on a sample preparation method for reconstitution of membrane proteins which uses porous anodic aluminum oxide (AAO) filters with 200-nm-wide pores of high density. The substrate permits formation of tubular, single membranes that line the inner surface of pores. One square centimeter of filter with a thickness of 60μm yields on the order of 500cm(2) of solid-supported single bilayer surface, sufficient for NMR studies. The tubular bilayers are free of detergent, fully hydrated, and accessible for ligands from one side of the membrane. The use of AAO filters greatly improves reproducibility of the reconstitution process such that the influence of protein on lipid order parameters can be studied with high resolution. As an example, results for the G protein-coupled receptor of class A, bovine rhodopsin, are shown. By (2)H NMR order parameter measurements, it is detected that rhodopsin insertion elastically deforms membranes near the protein. Furthermore, by (1)H saturation-transfer NMR under conditions of magic angle spinning, we demonstrate detection of preferences in interactions of rhodopsin with particular lipid species. It is assumed that function of integral membrane proteins depends on both protein-induced elastic deformations of the lipid matrix and preferences for interaction of the protein with particular lipid species in the first layer of lipids surrounding the protein.