The influence of membrane lateral pressures on simple geometric models of protein conformational equilibria

Protein-lipid acyl chain interactions: Depth-dependent changes of segmental mobility of phospholipid in contact with bacteriorhodopsin

Boundary lipids surrounding membrane proteins play an essential role in protein function and structure. These protein-lipid interactions are mainly divided into electrostatic interactions between the polar amino acids of proteins and polar heads of phospholipids, and hydrophobic interactions between protein transmembrane sites and phospholipid acyl chains. Our previous report (Kawatake et al., Biochim. Biophys. Acta 1858 [2016] 2106-2115) covered a method for selectively analyzing boundary lipid interactions and showed differences in membrane protein-peripheral lipid interactions due to differences in their head group. Interactions in the hydrophobic acyl chains of phospholipids are relatively consistent among proteins, but the details of these interactions have not been elucidated. In this study, we reconstituted bacteriorhodopsin as a model protein into phospholipid membranes labeled with 2H and 13C for solid-state NMR measurement to investigate the depth-dependent effect of the head group structure on the lipid bilayer. The results showed that the position of the phospholipid near the carbonyl carbon was affected by the head group in terms of selectivity for protein surfaces, whereas in the deep interior of the bilayer near the leaflet interface, there was little difference between the head groups, indicating that the dependence of their interactions on the head group was much reduced.

Cardiolipin binding enhances KcsA channel gating via both its specific and dianion-monoanion interchangeable sites

KcsA is a potassium channel with a plethora of structural and functional information, but its activity in the KcsA-producing actinomycete membranes remains elusive. To determine lipid species involved in channel-modulation, a surface plasmon resonance (SPR)-based methodology, characterized by immobilization of membrane proteins under a membrane environment, was applied. Dianionic cardiolipin (CL) showed extremely higher affinity for KcsA than monoanionic lipids. The SPR experiments further demonstrated that CL bound not only to the N-terminal M0 helix, a lipid-sensor domain, but to the M0 helix-deleted mutant. In contrast, monoanionic lipids interacted primarily with the M0 helix. This indicates the presence of an alternative CL-binding site, plausibly in the transmembrane domain. Single-channel recordings demonstrated that CL enhanced channel opening in an M0-independent manner. Taken together, the action of monoanionic lipids is exclusively mediated by the M0 helix, while CL binds both the M0 helix and its specific site, further enhancing the channel activity.

Nanoarchitectonic approaches for measuring the catalytic behavior of a membrane anchored enzyme. From Langmuir-Blodgett to a novel Langmuir-Schaefer based nanofilm building device

Self-organized lipid monolayers at the air-water interface (Langmuir films, LF) are commonly used for measuring the catalytic properties of membrane-bound enzymes. This methodology allows to provide a consistent flat topography molecular density, packing defects and thickness. The aim of the present work was to show the methodological advantages of using the horizontal transfer method (Langmuir-Schaefer) with respect to the vertical transfer method (Langmuir-Blodgett) when mounting a device to measure catalytic activity of membrane enzymes. Based on the results obtained we can conclude that it is possible to prepare stable Langmuir-Blodgett (LB) and Langmuir-Schaefer (LS) films from Bovine Erythrocyte Membranes (BEM) preserving the catalytic activity of its native Acetylcholinesterase (BEA). In comparison, the LS films showed Vmax values more similar to the enzyme present in the vesicles of natural membranes. In addition, it was much easier to produce large amounts of transferred areas with the horizontal transfer methodology. It was possible to decrease the time required to mount an assay with numerous activity points, such as building activity curves as a function of substrate concentration. The present results show that LSBEM provides a proof of concept for the development of biosensors based on transferred purified membrane for the screening of new products acting on an enzyme embedded on its natural milieu. In the case of BEA, the application of these enzymatic sensors could have medical interest, providing drug screening tools for the treatment of Alzheimer's disease.

Screening for bilayer-active and likely cytotoxic molecules reveals bilayer-mediated regulation of cell function

A perennial problem encountered when using small molecules (drugs) to manipulate cell or protein function is to assess whether observed changes in function result from specific interactions with a desired target or from less specific off-target mechanisms. This is important in laboratory research as well as in drug development, where the goal is to identify molecules that are unlikely to be successful therapeutics early in the process, thereby avoiding costly mistakes. We pursued this challenge from the perspective that many bioactive molecules (drugs) are amphiphiles that alter lipid bilayer elastic properties, which may cause indiscriminate changes in membrane protein (and cell) function and, in turn, cytotoxicity. Such drug-induced changes in bilayer properties can be quantified as changes in the monomer↔dimer equilibrium for bilayer-spanning gramicidin channels. Using this approach, we tested whether molecules in the Pathogen Box (a library of 400 drugs and drug-like molecules with confirmed activity against tropical diseases released by Medicines for Malaria Venture to encourage the development of therapies for neglected tropical diseases) are bilayer modifiers. 32% of the molecules in the Pathogen Box were bilayer modifiers, defined as molecules that at 10 µM shifted the monomer↔dimer equilibrium toward the conducting dimers by at least 50%. Correlation analysis of the molecules' reported HepG2 cell cytotoxicity to bilayer-modifying potency, quantified as the shift in the gramicidin monomer↔dimer equilibrium, revealed that molecules producing <25% change in the equilibrium had significantly lower probability of being cytotoxic than molecules producing >50% change. Neither cytotoxicity nor bilayer-modifying potency (quantified as the shift in the gramicidin monomer↔dimer equilibrium) was well predicted by conventional physico-chemical descriptors (hydrophobicity, polar surface area, etc.). We conclude that drug-induced changes in lipid bilayer properties are robust predictors of the likelihood of membrane-mediated off-target effects, including cytotoxicity.

Interaction of a Homologous Series of Amphiphiles with P-glycoprotein in a Membrane Environment-Contributions of Polar and Non-Polar Interactions

The transport of drugs by efflux transporters in biomembranes limits their bioavailability and is a major determinant of drug resistance development by cancer cells and pathogens. A large number of chemically dissimilar drugs are transported, and despite extensive studies, the molecular determinants of substrate specificity are still not well understood. In this work, we explore the role of polar and non-polar interactions on the interaction of a homologous series of fluorescent amphiphiles with the efflux transporter P-glycoprotein. The interaction of the amphiphiles with P-glycoprotein is evaluated through effects on ATPase activity, efficiency in inhibition of [¹²⁵I]-IAAP binding, and partition to the whole native membranes containing the transporter. The results were complemented with partition to model membranes with a representative lipid composition, and details on the interactions established were obtained from MD simulations. We show that when the total concentration of amphiphile is considered, the binding parameters obtained are apparent and do not reflect the affinity for P-gp. A new formalism is proposed that includes sequestration of the amphiphiles in the lipid bilayer and the possible binding of several molecules in P-gp's substrate-binding pocket. The intrinsic binding affinity thus obtained is essentially independent of amphiphile hydrophobicity, highlighting the importance of polar interactions. An increase in the lipophilicity and amphiphilicity led to a more efficient association with the lipid bilayer, which maintains the non-polar groups of the amphiphiles in the bilayer, while the polar groups interact with P-gp's binding pocket. The presence of several amphiphiles in this orientation is proposed as a mechanism for inhibition of P-pg function.

Heterosynaptic plasticity in biomembrane memristors controlled by pH

Abstract

In biology, heterosynaptic plasticity maintains homeostasis in synaptic inputs during associative learning and memory, and initiates long-term changes in synaptic strengths that nonspecifically modulate different synapse types. In bioinspired neuromorphic circuits, heterosynaptic plasticity may be used to extend the functionality of two-terminal, biomimetic memristors. In this article, we explore how changes in the pH of droplet interface bilayer aqueous solutions modulate the memristive responses of a lipid bilayer membrane in the pH range 4.97-7.40. Surprisingly, we did not find conclusive evidence for pH-dependent shifts in the voltage thresholds (V*) needed for alamethicin ion channel formation in the membrane. However, we did observe a clear modulation in the dynamics of pore formation with pH in time-dependent, pulsed voltage experiments. Moreover, at the same voltage, lowering the pH resulted in higher steady-state currents because of increased numbers of conductive peptide ion channels in the membrane. This was due to increased partitioning of alamethicin monomers into the membrane at pH 4.97, which is below the pKa (~5.3-5.7) of carboxylate groups on the glutamate residues of the peptide, making the monomers more hydrophobic. Neutralization of the negative charges on these residues, under acidic conditions, increased the concentration of peptide monomers in the membrane, shifting the equilibrium concentrations of peptide aggregate assemblies in the membrane to favor greater numbers of larger, increasingly more conductive pores. It also increased the relaxation time constants for pore formation and decay, and enhanced short-term facilitation and depression of the switching characteristics of the device. Modulating these thresholds globally and independently of alamethicin concentration and applied voltage will enable the assembly of neuromorphic computational circuitry with enhanced functionality.

Impact statement

We describe how to use pH as a modulatory "interneuron" that changes the voltage-dependent memristance of alamethicin ion channels in lipid bilayers by changing the structure and dynamical properties of the bilayer. Having the ability to independently control the threshold levels for pore conduction from voltage or ion channel concentration enables additional levels of programmability in a neuromorphic system. In this article, we note that barriers to conduction from membrane-bound ion channels can be lowered by reducing solution pH, resulting in higher currents, and enhanced short-term learning behavior in the form of paired-pulse facilitation. Tuning threshold values with environmental variables, such as pH, provide additional training and learning algorithms that can be used to elicit complex functionality within spiking neural networks.

Graphical abstract

Supplementary information

The online version contains supplementary material available at 10.1557/s43577-022-00344-z.

Mechanism Study of Proteins under Membrane Environment

Membrane proteins play crucial roles in various physiological processes, including molecule transport across membranes, cell communication, and signal transduction. Approximately 60% of known drug targets are membrane proteins. There is a significant need to deeply understand the working mechanism of membrane proteins in detail, which is a challenging work due to the lack of available membrane structures and their large spatial scale. Membrane proteins carry out vital physiological functions through conformational changes. In the current study, we utilized a coarse-grained (CG) model to investigate three representative membrane protein systems: the TMEM16A channel, the family C GPCRs mGlu2 receptor, and the P4-ATPase phospholipid transporter. We constructed the reaction pathway of conformational changes between the two-end structures. Energy profiles and energy barriers were calculated. These data could provide reasonable explanations for TMEM16A activation, the mGlu2 receptor activation process, and P4-ATPase phospholipid transport. Although they all belong to the members of membrane proteins, they behave differently in terms of energy. Our work investigated the working mechanism of membrane proteins and could give novel insights into other membrane protein systems of interest.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Membrane-Mediated Activity of Local Anesthetics

The activity of local anesthetics (LAs) has been attributed to the inhibition of ion channels, causing anesthesia. However, there is a growing body of research showing that LAs act on a wide range of receptors and channel proteins far beyond simple analgesia. The current concept of ligand recognition may no longer explain the multitude of protein targets influenced by LAs. We hypothesize that LAs can cause anesthesia without directly binding to the receptor proteins just by changing the physical properties of the lipid bilayer surrounding these proteins and ion channels based on LAs' amphiphilicity. It is possible that LAs act in one of the following ways: They 1) dissolve raft-like membrane microdomains, 2) impede nerve impulse propagation by lowering the lipid phase transition temperature, or 3) modulate the lateral pressure profile of the lipid bilayer. This could also explain the numerous additional effects of LAs besides anesthesia. Furthermore, the concepts of membrane-mediated activity and binding to ion channels do not have to exclude each other. If we were to consider LA as the middle part of a continuum between unspecific membrane-mediated activity on one end and highly specific ligand binding on the other end, we could describe LA as the link between the unspecific action of general anesthetics and toxins with their highly specific receptor binding. This comprehensive membrane-mediated model offers a fresh perspective to clinical and pharmaceutical research and therapeutic applications of local anesthetics. SIGNIFICANCE STATEMENT: Local anesthetics, according to the World Health Organization, belong to the most important drugs available to mankind. Their rediscovery as therapeutics and not only anesthetics marks a milestone in global pain therapy. The membrane-mediated mechanism of action proposed in this review can explain their puzzling variety of target proteins and their thus far inexplicable therapeutic effects. The new concept presented here places LAs on a continuum of structures and molecular mechanisms in between small general anesthetics and the more complex molecular toxins.

Partitioning of Aβ Peptide Fragments into Blood-Brain Barrier Mimetic Bilayer

We used all-atom replica-exchange umbrella sampling molecular dynamics simulations to investigate the partitioning of the charged tetrapeptide KLVF and its neutral apolar counterpart VVIA into the blood-brain barrier (BBB)-mimetic bilayer. Our findings allowed us to reconstruct the partitioning mechanism for these two Aβ peptide fragments. Despite dissimilar sequences, their permeation shares significant common features. Computations of free energies and permeabilities show that partitioning of both peptides is highly unfavorable, ruling out passive transport. The peptides experience multiple rotational transitions within the bilayer and typically cause considerable lipid disorder and bilayer thinning. Near the bilayer midplane, they lose almost entirely their solvation shell and the interactions with the lipid headgroups. The peptides cause complex reorganization within the proximal bilayer region. Upon insertion, they induce striking cholesterol influx reversed by its depletion and the influx of DMPC when the peptides reach the midplane. The differences in partitioning mechanisms are due to the much higher polarity of KLVF peptide, the permeation of which is more unfavorable and which exclusively assumes vertical orientations within the bilayer. In contrast, VVIA positions itself flat between the leaflets, causing minor disorder and even thickening of the BBB-mimetic bilayer. Due to the high density of the cholesterol-rich BBB bilayer, the unfavorable work associated with the peptide insertion provides a significant, but not dominant, contribution to the partition free energy, which is still governed by dehydration and loss of peptide-headgroup interactions. Comparison with experiments indicates that KLVF and VVIA permeation is similar to that of proline tetrapeptide, mannitol, or cimetidine, all of which exhibit no passive transport.

Increasing complexity in small-angle X-ray and neutron scattering experiments: from biological membrane mimics to live cells

Small-angle X-ray and neutron scattering are well-established, non-invasive experimental techniques to interrogate global structural properties of biological membrane mimicking systems under physiologically relevant conditions. Recent developments, both in bottom-up sample preparation techniques for increasingly complex model systems, and in data analysis techniques have opened the path toward addressing long standing issues of biological membrane remodelling processes. These efforts also include emerging quantitative scattering studies on live cells, thus enabling a bridging of molecular to cellular length scales. Here, we review recent progress in devising compositional models for joint small-angle X-ray and neutron scattering studies on diverse membrane mimics - with a specific focus on membrane structural coupling to amphiphatic peptides and integral proteins - and live Escherichia coli. In particular, we outline the present state-of-the-art in small-angle scattering methods applied to complex membrane systems, highlighting how increasing system complexity must be followed by an advance in compositional modelling and data-analysis tools.

Discoveries in structure and physiology of mechanically activated ion channels

The ability to sense physical forces is conserved across all organisms. Cells convert mechanical stimuli into electrical or chemical signals via mechanically activated ion channels. In recent years, the identification of new families of mechanosensitive ion channels-such as PIEZO and OSCA/TMEM63 channels-along with surprising insights into well-studied mechanosensitive channels have driven further developments in the mechanotransduction field. Several well-characterized mechanosensory roles such as touch, blood-pressure sensing and hearing are now linked with primary mechanotransducers. Unanticipated roles of mechanical force sensing continue to be uncovered. Furthermore, high-resolution structures representative of nearly every family of mechanically activated channel described so far have underscored their diversity while advancing our understanding of the biophysical mechanisms of pressure sensing. Here we summarize recent discoveries in the physiology and structures of known mechanically activated ion channel families and discuss their implications for understanding the mechanisms of mechanical force sensing.

Membrane Curvature Revisited-the Archetype of Rhodopsin Studied by Time-Resolved Electronic Spectroscopy

G-protein-coupled receptors (GPCRs) comprise the largest and most pharmacologically targeted membrane protein family. Here, we used the visual receptor rhodopsin as an archetype for understanding membrane lipid influences on conformational changes involved in GPCR activation. Visual rhodopsin was recombined with lipids varying in their degree of acyl chain unsaturation and polar headgroup size using 1-palmitoyl-2-oleoyl-sn-glycero- and 1,2-dioleoyl-sn-glycerophospholipids with phosphocholine (PC) or phosphoethanolamine (PE) substituents. The receptor activation profile after light excitation was measured using time-resolved ultraviolet-visible spectroscopy. We discovered that more saturated POPC lipids back shifted the equilibrium to the inactive state, whereas the small-headgroup, highly unsaturated DOPE lipids favored the active state. Increasing unsaturation and decreasing headgroup size have similar effects that combine to yield control of rhodopsin activation, and necessitate factors beyond proteolipid solvation energy and bilayer surface electrostatics. Hence, we consider a balance of curvature free energy with hydrophobic matching and demonstrate how our data support a flexible surface model (FSM) for the coupling between proteins and lipids. The FSM is based on the Helfrich formulation of membrane bending energy as we previously first applied to lipid-protein interactions. Membrane elasticity and curvature strain are induced by lateral pressure imbalances between the constituent lipids and drive key physiological processes at the membrane level. Spontaneous negative monolayer curvature toward water is mediated by unsaturated, small-headgroup lipids and couples directly to GPCR activation upon light absorption by rhodopsin. For the first time to our knowledge, we demonstrate this modulation in both the equilibrium and pre-equilibrium evolving states using a time-resolved approach.

In silico studies of the interactions between propofol and fentanyl using Gaussian accelerated molecular dynamics

Fentanyl is a potent opioid analgesic, which for decades has been used routinely in surgical and therapeutic applications. In addition to its analgesic properties, fentanyl also possesses anesthetic properties, which are not well understood. Fentanyl is used in the general anesthesia process to induce and maintain anesthesia in combination with the general anesthetic propofol, which fentanyl is known to potentiate. As the atomic-level mechanism behind the potentiation of propofol is unclear, we have used classical molecular dynamics simulations to study the interactions of these drugs with the Gloeobacter violaceus ion channel (GLIC). This ion channel has been identified as a target for many anesthetic drugs. We identified multiple binding sites using flooding style and Gaussian accelerated molecular dynamics (GaMD) simulations, showing fentanyl acting as a stabiliser that holds propofol within binding sites. Our extensive GaMD simulations were also able to show the pathway by which propofol blocks the channel pore, which has previously been suggested as a mechanism for ion channel modulation. General anesthesia is a multi-drug process and this study provides the first insight into the interactions between two different drugs in the anesthesia process in a relevant biological environment.Communicated by Ramaswamy H. Sarma.

Degeneracy in molecular scale organization of biological membranes

The scale-rich spatiotemporal organization in biological membranes has its origin in the differential inter- and intra-molecular interactions among their constituents. In this work, we explore the molecular-origin behind that variety and possible degeneracy in lateral organization in membranes. For our study, we post-process microsecond long all-atom molecular dynamics trajectories for three systems that exhibit fluid phase coexistence: (i) PSM/POPC/Chol (0.47/0.32/0.21), (ii) PSM/DOPC/Chol (0.43/0.38/0.19) and (iii) DPPC/DOPC/Chol (0.37/0.36/0.27). To distinguish the liquid ordered and disordered regions at molecular scales, we calculate the degree of non-affineness of individual lipids in their neighbourhood and track their topological rearrangements. Disconnectivity graph analysis with respect to membrane organization shows that the DPPC/DOPC/Chol and PSM/DOPC/Chol systems exhibit funnel-like energy landscapes as opposed to a highly frustrated energy landscape for the more biomimetic PSM/POPC/Chol system. We use these measurements to develop a continuous lattice Hamiltonian and evolve that using Monte Carlo simulated annealing to explore the possibility of structural degeneracy in membrane organization. Our data show that model membranes with lipid constituents that are biomimetic (PSM/POPC/Chol) have the ability to access a large range of membrane sub-structure space (higher degeneracy) as compared to the other two systems, which form only one kind of substructure even with changing composition. Since the spatiotemporal organization in biological membranes dictates the "molecular encounters" and in turn larger scale biological processes such as molecular transport, trafficking and cellular signalling, we posit that this structural degeneracy could enable access to a larger repository to functionally important molecular organization in systems with physiologically relevant compositions.

Sensing molecular organizational changes through the catalytic activity of acetylcholinesterase from erythrocyte membranes in Langmuir-Blodgett films

Langmuir films prepared from bovine erythrocyte membranes (LF_BEM) were studied and transferred to alkylated glasses (Langmuir-Blodgett films, LB_BEM) in order to assess the effects of membrane molecular packing on Bovine Erythrocyte Acetylcholinesterase (BEA) catalytic activity. Surface pressure (π) vs Area isotherms showed three 2D-transitions at ~7, ~18 and ~44 mN/m and a collapse pressure at π_c = 49 mN/m. The 0-12-0 mN/m compression-decompression cycles resulted reversible while those 0-40-0 mN/m exhibited a significant hysteresis. Taken together, EFM, BAM and AFM images and the stability of the film after 3C-D cycles, we can suggest that over the air-water interface as well as over the silanized glass substrate the surface is mostly covered by a monolayer with a few particles dispersed. Acetylthiocholine hydrolysis was assayed with BEA in bovine erythrocyte membrane suspensions (S_BEM) and in LB_BEM packed at 10 (LB_BEM,10) and 35 mN/m (LB_BEM,35), which gave the following kinetic parameters: V_max = 3.41 ± 0.15, 0.021 ± 0.002 and 0.030 ± 0.003 nmol.min^-1·μg prot^-1 and K_M = 0.11 ± 0.02, 0.047 ± 0.017 and 0.026 ± 0.017 mM, respectively. Although from S_BEM to LB_BEM we lost active enzyme, the catalytic efficiency (V_max/K_M) increased ~750 times. Eugenol and 1,8-cineol inhibited BEA catalytic activity in LB_BEM,35. Our results demonstrate the transmission of information between the membrane and the environment within the subphase immediately below the membrane, where anchored proteins are hosted. This was reflected by the membrane packing-induced modulation of BEA catalytic activity. Furthermore, LB_BEM provides a proof of concept for the development of biosensors to screen new green pesticides acting through BEA interaction.

Archaeal Glycolipid S-TGA-1 Is Crucial for Trimer Formation and Photocycle Activity of Bacteriorhodopsin

Although it has been demonstrated that membrane proteins (MPs) require lipids to ensure their structural and functional integrity, details on how lipid-MP interactions regulate MPs are still unclear. Recently, we developed a concise method for quantitatively evaluating lipid-MP interactions and applied it to bacteriorhodopsin (bR), a halobacterial MP that forms trimers and acts as a light-driven proton pump. Consequently, we found that the halobacterial glycolipid, S-TGA-1, has the highest affinity for bR, among other lipids. In this study, we examined the effects of S-TGA-1 on bR via visible circular dichroism spectroscopy, flash photolysis, and proton influx measurement. The results showed that S-TGA-1 efficiently promotes trimer formation, photocycle, and proton pumping in bR. Our data also suggested that the bR photocycle is restored as a consequence of the trimerization induced by the lipid. This study demonstrates clearly that lipids specifically interacting with MPs can have significant impacts on MP structure and/or function. The methodology adopted in our studies can be applied to other MPs and will help elucidate the physiological functions of lipids in terms of lipid-MP interactions, thus accelerating "lipid chemical biology" studies.

Molecular basis of force-from-lipids gating in the mechanosensitive channel MscS

Prokaryotic mechanosensitive (MS) channels open by sensing the physical state of the membrane. As such, lipid-protein interactions represent the defining molecular process underlying mechanotransduction. Here, we describe cryo-electron microscopy (cryo-EM) structures of the E. coli small-conductance mechanosensitive channel (MscS) in nanodiscs (ND). They reveal a novel membrane-anchoring fold that plays a significant role in channel activation and establish a new location for the lipid bilayer, shifted ~14 Å from previous consensus placements. Two types of lipid densities are explicitly observed. A phospholipid that ‘hooks’ the top of each TM2-TM3 hairpin and likely plays a role in force sensing, and a bundle of acyl chains occluding the permeation path above the L105 cuff. These observations reshape our understanding of force-from-lipids gating in MscS and highlight the key role of allosteric interactions between TM segments and phospholipids bound to key dynamic components of the channel.

Interaction of Antimicrobial Lipopeptides with Bacterial Lipid Bilayers

The resistance of pathogens to traditional antibiotics is currently a global issue of enormous concern. As the discovery and development of new antibiotics become increasingly challenging, synthetic antimicrobial lipopeptides (AMLPs) are now receiving renewed attention as a new class of antimicrobial agents. In contrast to traditional antibiotics, AMLPs act by physically disrupting the cell membrane (rather than targeting specific proteins), thus reducing the risk of inducing bacterial resistance. In this study, we use microsecond-timescale atomistic molecular dynamics simulations to quantify the interaction of a short AMLP (C16-KKK) with model bacterial lipid bilayers. In particular, we investigate how fundamental transmembrane properties change in relation to a range of lipopeptide concentrations. A number of structural, mechanical, and dynamical features are found to be significantly altered in a non-linear fashion. At 10 mol% concentration, lipopeptides have a condensing effect on bacterial bilayers, characterized by a decrease in the area per lipid and an increase in the bilayer order. Higher AMLP concentrations of 25 and 40 mol% destabilize the membrane by disrupting the bilayer core structure, inducing membrane thinning and water leakage. Important transmembrane properties such as the lateral pressure and dipole potential profiles are also affected. Potential implications on membrane function and associated proteins are discussed.

Electrophysiological-mechanical coupling in the neuronal membrane and its role in ultrasound neuromodulation and general anaesthesia

The current understanding of the role of the cell membrane is in a state of flux. Recent experiments show that conventional models, considering only electrophysiological properties of a passive membrane, are incomplete. The neuronal membrane is an active structure with mechanical properties that modulate electrophysiology. Protein transport, lipid bilayer phase, membrane pressure and stiffness can all influence membrane capacitance and action potential propagation. A mounting body of evidence indicates that neuronal mechanics and electrophysiology are coupled, and together shape the membrane potential in tight coordination with other physical properties. In this review, we summarise recent updates concerning electrophysiological-mechanical coupling in neuronal function. In particular, we aim at making the link with two relevant yet often disconnected fields with strong clinical potential: the use of mechanical vibrations-ultrasound-to alter the electrophysiogical state of neurons, e.g., in neuromodulation, and the theories attempting to explain the action of general anaesthetics. STATEMENT OF SIGNIFICANCE: General anaesthetics revolutionised medical practice; now an apparently unrelated technique, ultrasound neuromodulation-aimed at controlling neuronal activity by means of ultrasound-is poised to achieve a similar level of impact. While both technologies are known to alter the electrophysiology of neurons, the way they achieve it is still largely unknown. In this review, we argue that in order to explain their mechanisms/effects, the neuronal membrane must be considered as a coupled mechano-electrophysiological system that consists of multiple physical processes occurring concurrently and collaboratively, as opposed to sequentially and independently. In this framework the behaviour of the cell membrane is not the result of stereotypical mechanisms in isolation but instead emerges from the integrative behaviour of a complexly coupled multiphysics system.

General and specific interactions of the phospholipid bilayer with P-type ATPases

Protein structure and function are modulated via interactions with their environment, representing both the surrounding aqueous media and lipid membranes that have an active role in shaping the structural topology of membrane proteins. Compared to a decade ago, there is now an abundance of crystal structural data on membrane proteins, which together with their functional studies have enhanced our understanding of the salient features of lipid-protein interactions. It is now important to recognize that membrane proteins are regulated by both (1) general lipid-protein interactions, where the general physicochemical properties of the lipid environment affect the conformational flexibility of a membrane protein, and (2) by specific lipid-protein interactions, where lipid molecules directly interact via chemical interactions with specific lipid-binding sites located on the protein. However, due to local differences in membrane composition, thickness, and lipid packing, local membrane physical properties and hence the associated lipid-protein interactions also differ due to membrane location, even for the same protein. Such a phenomenon has been shown to be true for one family of integral membrane ion pumps, the P2-type adenosine triphosphatases (ATPases). Despite being highly homologous, individual members of this family have distinct structural and functional activity and are an excellent candidate to highlight how the local membrane physical properties and specific lipid-protein interactions play a vital role in facilitating the structural rearrangements of these proteins necessary for their activity. Hence in this review, we focus on both the general and specific lipid-protein interactions and will mostly discuss the structure-function relationships of the following P2-type ATPases, Na⁺,K⁺-ATPase (NKA), gastric H⁺,K⁺-ATPase (HKA), and sarco(endo)plasmic reticulum Ca²⁺-ATPase (SERCA), in concurrence with their lipid environment.

A concise method for quantitative analysis of interactions between lipids and membrane proteins

Although interactions between lipids and membrane proteins (MPs) have been considered crucially important for understanding the functions of lipids, lack of useful and convincing experimental methods has hampered the analysis of the interactions. Here, we developed a surface plasmon resonance (SPR)-based concise method for quantitative analysis of lipid-MP interactions, coating the sensor chip surface with self-assembled monolayer (SAM) with C₆-chain. To develop this method, we used bacteriorhodopsin (bR) as an MP, and examined its interaction with various types of lipids. The merits of using C₆-SAM-modified sensor chip are as follows: (1) alkyl-chains of SAM confer a better immobilization of MPs because of the efficient preconcentration due to hydrophobic contacts; (2) SAM provides immobilized MPs with a partial membranous environment, which is important for the stabilization of MPs; and (3) a thinner C₆-SAM layer (1 nm) compared with MP size forces the MP to bulge outward from the SAM surface, allowing extraneously injected lipids to be accessible to the hydrophobic transmembrane regions. Actually, the amount of bR immobilized on C₆-SAM is 10 times higher than that on a hydrophilic CM5 sensor chip, and AFM observations confirmed that bR molecules are exposed on the SAM surface. Of the lipids tested, S-TGA-1, a halobacterium-derived glycolipid, had the highest specificity to bR with a nanomolar dissociation constant. This is consistent with the reported co-crystal structure that indicates the formation of several intermolecular hydrogen bonds. Therefore, we not only reproduced the specific lipid-bR recognition, but also succeeded in its quantitative evaluation, demonstrating the validity and utility of this method.

Flexible lipid nanomaterials studied by NMR spectroscopy

Advances in solid-state nuclear magnetic resonance spectroscopy inform the emergence of material properties from atomistic-level interactions in membrane lipid nanostructures.

Dual Role of the C-Terminal Domain in Osmosensing by Bacterial Osmolyte Transporter ProP

ProP is a member of the major facilitator superfamily, a proton-osmolyte symporter, and an osmosensing transporter. ProP proteins share extended cytoplasmic carboxyl terminal domains (CTDs) implicated in osmosensing. The CTDs of the best characterized, group A ProP orthologs, terminate in sequences that form intermolecular, antiparallel α-helical coiled coils (e.g., ProPEc, from Escherichia coli). Group B orthologs lack that feature (e.g., ProPXc, from Xanthomonas campestris). ProPXc was expressed and characterized in E. coli to further elucidate the role of the coiled coil in osmosensing. The activity of ProPXc was a sigmoid function of the osmolality in cells and proteoliposomes. ProPEc and ProPXc attained similar activities at the same expression level in E. coli. ProPEc transports proline and glycine betaine with comparable high affinities at low osmolality. In contrast, proline weakly inhibited high-affinity glycine-betaine uptake via ProPXc. The KM for proline uptake via ProPEc increases dramatically with the osmolality. The KM for glycine-betaine uptake via ProPXc did not. Thus, ProPXc is an osmosensing transporter, and the C-terminal coiled coil is not essential for osmosensing. The role of CTD-membrane interaction in osmosensing was examined further. As for ProPEc, the ProPXc CTD co-sedimented with liposomes comprising E. coli phospholipid. Molecular dynamics simulations illustrated association of the monomeric ProPEc CTD with the membrane surface. Comparison with the available NMR structure for the homodimeric coiled coil formed by the ProPEc-CTD suggested that membrane association and homodimeric coiled-coil formation by that peptide are mutually exclusive. The membrane fluidity in liposomes comprising E. coli phospholipid decreased with increasing osmolality in the range relevant for ProP activation. These data support the proposal that ProP activates as cellular dehydration increases cytoplasmic cation concentration, releasing the CTD from the membrane surface. For group A orthologs, this also favors α-helical coiled-coil formation that stabilizes the transporter in an active form.

Lipid-Protein Interactions Are Unique Fingerprints for Membrane Proteins

Cell membranes contain hundreds of different proteins and lipids in an asymmetric arrangement. Our current understanding of the detailed organization of cell membranes remains rather elusive, because of the challenge to study fluctuating nanoscale assemblies of lipids and proteins with the required spatiotemporal resolution. Here, we use molecular dynamics simulations to characterize the lipid environment of 10 different membrane proteins. To provide a realistic lipid environment, the proteins are embedded in a model plasma membrane, where more than 60 lipid species are represented, asymmetrically distributed between the leaflets. The simulations detail how each protein modulates its local lipid environment in a unique way, through enrichment or depletion of specific lipid components, resulting in thickness and curvature gradients. Our results provide a molecular glimpse of the complexity of lipid-protein interactions, with potentially far-reaching implications for our understanding of the overall organization of real cell membranes.

The role of phospholipid molecular species in determining the physical properties of yeast membranes

In most eukaryotes, including Saccharomyces cerevisiae, glycerophospholipids are the main membrane lipid constituents. Besides serving as general membrane ‘building blocks’, glycerophospholipids play an important role in determining the physical properties of the membrane, which are crucial for proper membrane function. To ensure optimal physical properties, membrane glycerophospholipid composition and synthesis are tightly regulated. This review will summarize our current knowledge of factors and processes determining the membrane glycerophospholipid composition of the reference eukaryote S. cerevisiae at the level of molecular species. Extrapolating from relevant model membrane data, we also discuss how modulation of the molecular species composition can regulate membrane physical properties.

Effects of High Pressure on Phospholipid Bilayers

The response of lipid membranes to changes in external pressure is important for many biological processes, and it can also be exploited for technological applications. In this work, we employ all-atom molecular dynamics simulations to characterize the changes in the physical properties of phospholipid bilayers brought about by high pressure (1000 bar). In particular, we study how the response differs, in relation to different chain unsaturation levels, by comparing monounsaturated 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) and biunsaturated dioleoyl-phosphatidylcholine (DOPC) bilayers. Various structural, mechanical, and dynamical features are found to be altered by the pressure increase in both bilayers. Notably, for most properties, including bilayer area and thickness, lipid order parameters, lateral pressure profile, and curvature frustration energy, we observe significantly more pronounced effects for monounsaturated POPC than biunsaturated DOPC. Possible biological implications of the results obtained are discussed, especially in relation to how different lipids can control the structure and function of membrane proteins.

Cholesterol Changes the Mechanisms of Aβ Peptide Binding to the DMPC Bilayer

Using isobaric-isothermal all-atom replica-exchange molecular dynamics (REMD) simulations, we investigated the equilibrium binding of Aβ10-40 monomers to the zwitterionic dimyristoylphosphatidylcholine (DMPC) bilayer containing cholesterol. Our previous REMD simulations, which studied binding of the same peptide to the cholesterol-free DMPC bilayer, served as a control, against which we measured the impact of cholesterol. Our findings are as follows. First, addition of cholesterol to the DMPC bilayer partially expels the Aβ peptide from the hydrophobic core and promotes its binding to bilayer polar headgroups. Using thermodynamic and energetics analyses, we argued that Aβ partial expulsion is not related to cholesterol-induced changes in lateral pressure within the bilayer but is caused by binding energetics, which favors Aβ binding to the surface of the densely packed cholesterol-rich bilayer. Second, cholesterol has a protective effect on the DMPC bilayer structure against perturbations caused by Aβ binding. More specifically, cholesterol reduces bilayer thinning and overall depletion of bilayer density beneath the Aβ binding footprint. Third, we found that the Aβ peptide contains a single cholesterol binding site, which involves hydrophobic C-terminal amino acids (Ile31-Val36), Phe19, and Phe20 from the central hydrophobic cluster, and cationic Lys28 from the turn region. This binding site accounts for about 76% of all Aβ-cholesterol interactions. Because cholesterol binding site in the Aβ10-40 peptide does not contain the GXXXG motif featured in cholesterol interactions with the transmembrane domain C99 of the β-amyloid precursor protein, we argued that the binding mechanisms for Aβ and C99 are distinct reflecting their different conformations and positions in the lipid bilayer. Fourth, cholesterol sharply reduces the helical propensity in the bound Aβ peptide. As a result, cholesterol largely eliminates the emergence of helical structure observed upon Aβ transition from a water environment to the cholesterol-free DMPC bilayer. We explain this effect by the formation of hydrogen bonds between cholesterol and the Aβ backbone, which prevent helix formation. Taken together, we expect that our simulations will advance understanding of a molecular-level mechanism behind the role of cholesterol in Alzheimer's disease.

How do mechanosensitive channels sense membrane tension?

Mechanosensitive (MS) channels provide protection against hypo-osmotic shock in bacteria whereas eukaryotic MS channels fulfil a multitude of important functions beside osmoregulation. Interactions with the membrane lipids are responsible for the sensing of mechanical force for most known MS channels. It emerged recently that not only prokaryotic, but also eukaryotic, MS channels are able to directly sense the tension in the membrane bilayer without any additional cofactor. If the membrane is solely viewed as a continuous medium with specific anisotropic physical properties, the sensitivity towards tension changes can be explained as result of the hydrophobic coupling between membrane and transmembrane (TM) regions of the channel. The increased cross-sectional area of the MS channel in the active conformation and elastic deformations of the membrane close to the channel have been described as important factors. However, recent studies suggest that molecular interactions of lipids with the channels could play an important role in mechanosensation. Pockets in between TM helices were identified in the MS channel of small conductance (MscS) and YnaI that are filled with lipids. Less lipids are present in the open state of MscS than the closed according to MD simulations. Thus it was suggested that exclusion of lipid fatty acyl chains from these pockets, as a consequence of increased tension, would trigger gating. Similarly, in the eukaryotic MS channel TRAAK it was found that a lipid chain blocks the conducting path in the closed state. The role of these specific lipid interactions in mechanosensation are highlighted in this review.

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.

Effect of head group orientation on phospholipid assembly

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

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.

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.

The Environment Shapes the Inner Vestibule of LeuT

Human neurotransmitter transporters are found in the nervous system terminating synaptic signals by rapid removal of neurotransmitter molecules from the synaptic cleft. The homologous transporter LeuT, found in Aquifex aeolicus, was crystallized in different conformations. Here, we investigated the inward-open state of LeuT. We compared LeuT in membranes and micelles using molecular dynamics simulations and lanthanide-based resonance energy transfer (LRET). Simulations of micelle-solubilized LeuT revealed a stable and widely open inward-facing conformation. However, this conformation was unstable in a membrane environment. The helix dipole and the charged amino acid of the first transmembrane helix (TM1A) partitioned out of the hydrophobic membrane core. Free energy calculations showed that movement of TM1A by 0.30 nm was driven by a free energy difference of ~15 kJ/mol. Distance measurements by LRET showed TM1A movements, consistent with the simulations, confirming a substantially different inward-open conformation in lipid bilayer from that inferred from the crystal structure.

Crystal structures of the bacterial small amino acid transporter LeuT provided structural evidence for the alternating access model. Thereby, these structures shaped our understanding of the mechanisms underlying substrate translocation by neurotransmitter transporters. However, it has been questioned, if the crystallized inward-open conformation of LeuT can exist in the membrane environment. Here we show that, while stable in detergent micelles, the inward-open conformation of LeuT is of high energy and undergoes structural readjustments. We use a multi-faceted approach including molecular dynamics simulations, scintillation proximity assays, free energy calculations and apply for the first time lanthanide resonance energy transfer measurements to verify the in silico predictions. In silico and in vitro approaches using the same conditions allowed us to combine the macroscopic experimental data with microscopic all atom results from simulations to identify the underlying driving forces: partitioning of charged and polar groups from the hydrophobic membrane interior to the hydrophilic environment. We propose that the inward-facing state shows a much smaller movement of TM1A, but large enough to create an access path to the S1 substrate binding site from the vestibule.

CHARMM-GUI HMMM Builder for Membrane Simulations with the Highly Mobile Membrane-Mimetic Model

Slow diffusion of the lipids in conventional all-atom simulations of membrane systems makes it difficult to sample large rearrangements of lipids and protein-lipid interactions. Recently, Tajkhorshid and co-workers developed the highly mobile membrane-mimetic (HMMM) model with accelerated lipid motion by replacing the lipid tails with small organic molecules. The HMMM model provides accelerated lipid diffusion by one to two orders of magnitude, and is particularly useful in studying membrane-protein associations. However, building an HMMM simulation system is not easy, as it requires sophisticated treatment of the lipid tails. In this study, we have developed CHARMM-GUI HMMM Builder (http://www.charmm-gui.org/input/hmmm) to provide users with ready-to-go input files for simulating HMMM membrane systems with/without proteins. Various lipid-only and protein-lipid systems are simulated to validate the qualities of the systems generated by HMMM Builder with focus on the basic properties and advantages of the HMMM model. HMMM Builder supports all lipid types available in CHARMM-GUI and also provides a module to convert back and forth between an HMMM membrane and a full-length membrane. We expect HMMM Builder to be a useful tool in studying membrane systems with enhanced lipid diffusion.

Impact of membrane lipid composition on the structure and stability of the transmembrane domain of amyloid precursor protein

Cleavage of the amyloid precursor protein (APP) by γ-secretase is a crucial first step in the evolution of Alzheimer's disease. To discover the cleavage mechanism, it is urgent to predict the structures of APP monomers and dimers in varying membrane environments. We determined the structures of the C9923-55 monomer and homodimer as a function of membrane lipid composition using a multiscale simulation approach that blends atomistic and coarse-grained models. We demonstrate that the C9923-55 homodimer structures form a heterogeneous ensemble with multiple conformational states, each stabilized by characteristic interpeptide interactions. The relative probabilities of each conformational state are sensitive to the membrane environment, leading to substantial variation in homodimer peptide structure as a function of membrane lipid composition or the presence of an anionic lipid environment. In contrast, the helicity of the transmembrane domain of monomeric C991-55 is relatively insensitive to the membrane lipid composition, in agreement with experimental observations. The dimer structures of human EphA2 receptor depend on the lipid environment, which we show is linked to the location of the structural motifs in the dimer interface, thereby establishing that both sequence and membrane composition modulate the complete energy landscape of membrane-bound proteins. As a by-product of our work, we explain the discrepancy in structures predicted for C99 congener homodimers in membrane and micelle environments. Our study provides insight into the observed dependence of C99 protein cleavage by γ-secretase, critical to the formation of amyloid-β protein, on membrane thickness and lipid composition.

Dipole-Potential-Mediated Effects on Ion Pump Kinetics

The kinetics of conformational changes of P-type ATPases necessary for the occlusion or deocclusion of transported ions are known to be sensitive to the composition of the surrounding membrane, e.g., phospholipid content, mole percentage of cholesterol, and the presence of lipid-bound anions. Research has shown that many membrane components modify the dipole potential of the lipid head-group region. Based on the observation that occlusion/deocclusion reactions of ion pumps perturb the membrane surrounding the protein, a mechanism is suggested whereby dipole potential modifiers induce preferential stabilization or destabilization of occluded or nonoccluded states of the protein, leading to changes in the forward and backward rate constants for the transition. The mechanism relies on the assumption that conformational changes of the protein are associated with changes in its hydrophobic thickness that requires a change in local lipid packing density to allow hydrophobic matching with the membrane. The changes in lipid packing density cause changes in local lipid dipole potential that are responsible for the dependence of conformational kinetics on dipole potential modifiers. The proposed mechanism has the potential to explain effects of lipid composition on the kinetics of any membrane protein undergoing significant changes in its membrane cross-sectional area during its activity.

Lipids modulate the conformational dynamics of a secondary multidrug transporter

Direct interactions with lipids have emerged as key determinants of the folding, structure and function of membrane proteins, but an understanding of how lipids modulate protein dynamics is still lacking. Here, we systematically explored the effects of lipids on the conformational dynamics of the proton-powered multidrug transporter LmrP from Lactococcus lactis, using the pattern of distances between spin-label pairs previously shown to report on alternating access of the protein. We uncovered, at the molecular level, how the lipid headgroups shape the conformational-energy landscape of the transporter. The model emerging from our data suggests a direct interaction between lipid headgroups and a conserved motif of charged residues that control the conformational equilibrium through an interplay of electrostatic interactions within the protein. Together, our data lay the foundation for a comprehensive model of secondary multidrug transport in lipid bilayers.

Lateral pressure-mediated protein partitioning into liquid-ordered/liquid-disordered domains

We have studied the contributions of stored elastic energies in liquid-ordered (Lo) and liquid-disordered (Ld) domains to transmembrane proteins using the lateral pressure concept. In particular we applied previously reported experimental data for the membrane thickness, intrinsic curvature and bending elasticities of coexisting Lo/Ld domains to calculate whether proteins of simple geometric shapes would preferentially diffuse into Lo or Ld domains and form oligomers of a certain size. For the studied lipid mixture we generally found that proteins with convex shapes prefer sorting to Ld phases and the formation of large clusters. Lo domains in turn would be enriched in monomers of concave shaped proteins. We further observed that proteins which are symmetric with respect to the bilayer center prefer symmetric Lo or Ld domains, while asymmetric proteins favor a location in domains with Lo/Ld asymmetry. In the latter case we additionally retrieved a strong dependence on protein directionality, thus providing a mechanism for transmembrane protein orientation.

Exploring the mechanism of general anesthesia: kinetic analysis of GABAA receptor electrophysiology

A kinetic model of the effect of agonist and anesthetics on ligand-gated ion channels, developed in earlier work, is further refined and used to predict traces observed in fast-perfusion electrophysiological studies on recombinant GABAA receptors under a wide range of agonist and/or anesthetic concentrations. The model incorporates only three conformational states (resting, open, and desensitized) but allows for the modulation of the conformational free energy landscape connecting these states resulting from adsorption of agonist and/or anesthetic to the bilayer in which the protein is embedded. The model is shown to reproduce the diverse and complex features of experimental traces remarkably well, including both anesthetic-induced and agonist-induced traces, as well as the modulation of agonist-induced traces by anesthetic, either coapplied or continuously present. The solutions to the kinetic equations, which give the time-dependence of each of the nine protein states (three ligation states for each of the three conformations), describe the flow of probability among these states and thus reveal the kinetic underpinnings of the traces. Many of the parameters in the model, such as the desorption rate constants of anesthetic and agonist, are directly related to model-independent experimental measurements and thus can serve as a definitive test of its validity.

Membranes Do Not Tell Proteins How To Fold

Which properties of the membrane environment are essential for the folding and oligomerization of transmembrane proteins? Because the lipids that surround membrane proteins in situ spontaneously organize into bilayers, it may seem intuitive that interactions with the bilayer provide both hydrophobic and topological constraints that help the protein to achieve a stable and functional three-dimensional structure. However, one may wonder whether folding is actually driven by the membrane environment or whether the folded state just reflects an adaptation of integral proteins to the medium in which they function. Also, apart from the overall transmembrane orientation, might the asymmetry inherent in biosynthesis processes cause proteins to fold to out-of-equilibrium, metastable topologies? Which of the features of a bilayer are essential for membrane protein folding, and which are not? To which extent do translocons dictate transmembrane topologies? Recent data show that many membrane proteins fold and oligomerize very efficiently in media that bear little similarity to a membrane, casting doubt on the essentiality of many bilayer constraints. In the following discussion, we argue that some of the features of bilayers may contribute to protein folding, stability and regulation, but they are not required for the basic three-dimensional structure to be achieved. This idea, if correct, would imply that evolution has steered membrane proteins toward an accommodation to biosynthetic pathways and a good fit into their environment, but that their folding is not driven by the latter or dictated by insertion apparatuses. In other words, the three-dimensional structure of membrane proteins is essentially determined by intramolecular interactions and not by bilayer constraints and insertion pathways. Implications are discussed.

Exploring the Effects on Lipid Bilayer Induced by Noble Gases via Molecular Dynamics Simulations

Noble gases seem to have no significant effect on the anesthetic targets due to their simple, spherical shape. However, xenon has strong narcotic efficacy and can be used clinically, while other noble gases cannot. The mechanism remains unclear. Here, we performed molecular dynamics simulations on phospholipid bilayers with four kinds of noble gases to elucidate the difference of their effects on the membrane. Our results showed that the sequence of effects on membrane exerted by noble gases from weak to strong was Ne, Ar, Kr and Xe, the same order as their relative narcotic potencies as well as their lipid/water partition percentages. Compared with the other three kinds of noble gases, more xenon molecules were distributed between the lipid tails and headgroups, resulting in membrane’s lateral expansion and lipid tail disorder. It may contribute to xenon’s strong anesthetic potency. The results are well consistent with the membrane mediated mechanism of general anesthesia.

Osmotic Stress

Escherichia coli and Salmonella encounter osmotic pressure variations in natural environments that include host tissues, food, soil, and water. Osmotic stress causes water to flow into or out of cells, changing their structure, physics, and chemistry in ways that perturb cell functions. E. coli and Salmonella limit osmotically induced water fluxes by accumulating and releasing electrolytes and small organic solutes, some denoted compatible solutes because they accumulate to high levels without disturbing cell functions. Osmotic upshifts inhibit membrane-based energy transduction and macromolecule synthesis while activating existing osmoregulatory systems and specifically inducing osmoregulatory genes. The osmoregulatory response depends on the availability of osmoprotectants (exogenous organic compounds that can be taken up to become compatible solutes). Without osmoprotectants, K+ accumulates with counterion glutamate, and compatible solute trehalose is synthesized. Available osmoprotectants are taken up via transporters ProP, ProU, BetT, and BetU. The resulting compatible solute accumulation attenuates the K+ glutamate response and more effectively restores cell hydration and growth. Osmotic downshifts abruptly increase turgor pressure and strain the cytoplasmic membrane. Mechanosensitive channels like MscS and MscL open to allow nonspecific solute efflux and forestall cell lysis. Research frontiers include (i) the osmoadaptive remodeling of cell structure, (ii) the mechanisms by which osmotic stress alters gene expression, (iii) the mechanisms by which transporters and channels detect and respond to osmotic pressure changes, (iv) the coordination of osmoregulatory programs and selection of available osmoprotectants, and (v) the roles played by osmoregulatory mechanisms as E. coli and Salmonella survive or thrive in their natural environments.

Lipid insertion domain unfolding regulates protein orientational transition behavior in a lipid bilayer

We have used coarse-grained (CG) and united atom (UA) molecular dynamics simulations to explore the mechanisms of protein orientational transition of a model peptide (Aβ42) in a phosphatidylcholine/cholesterol (PC/CHO) lipid bilayer. We started with an inserted state of Aβ42 containing a folded (I) or unfolded (II) K28-A42 lipid insertion domain (LID), which was stabilized by the K28-snorkeling and A42-anchoring to the PC polar groups in the lipid bilayer. After a UA-to-CG transformation and a 1000ns-CG simulation for enhancing the sampling of protein orientations, we discovered two transitions: I-to-"deep inserted" state with disrupted K28-snorkeling and II-to-"deep surface" state with disrupted A42-anchoring. The new states remained stable after a CG-to-UA transformation and a 200ns-UA simulation relaxation. Significant changes in the cholesterol-binding domain of Aβ42 and protein-induced membrane disruptions were evident after the transitions. We propose that the conformation of the LID regulates protein orientational transitions in the lipid membrane.

Lipid-protein interactions: Lessons learned from stress

Biological membranes are essential for normal function and regulation of cells, forming a physical barrier between extracellular and intracellular space and cellular compartments. These physical barriers are subject to mechanical stresses. As a consequence, nature has developed proteins that are able to transpose mechanical stimuli into meaningful intracellular signals. These proteins, termed Mechanosensitive (MS) proteins provide a variety of roles in response to these stimuli. In prokaryotes these proteins form transmembrane spanning channels that function as osmotically activated nanovalves to prevent cell lysis by hypoosmotic shock. In eukaryotes, the function of MS proteins is more diverse and includes physiological processes such as touch, pain and hearing. The transmembrane portion of these channels is influenced by the physical properties such as charge, shape, thickness and stiffness of the lipid bilayer surrounding it, as well as the bilayer pressure profile. In this review we provide an overview of the progress to date on advances in our understanding of the intimate biophysical and chemical interactions between the lipid bilayer and mechanosensitive membrane channels, focusing on current progress in both eukaryotic and prokaryotic systems. These advances are of importance due to the increasing evidence of the role the MS channels play in disease, such as xerocytosis, muscular dystrophy and cardiac hypertrophy. Moreover, insights gained from lipid-protein interactions of MS channels are likely relevant not only to this class of membrane proteins, but other bilayer embedded proteins as well. This article is part of a Special Issue entitled: Lipid-protein interactions.

Thermosensing via transmembrane protein-lipid interactions

Cell membranes are composed of a lipid bilayer containing proteins that cross and/or interact with lipids on either side of the two leaflets. The basic structure of cell membranes is this bilayer, composed of two opposing lipid monolayers with fascinating properties designed to perform all the functions the cell requires. To coordinate these functions, lipid composition of cellular membranes is tailored to suit their specialized tasks. In this review, we describe the general mechanisms of membrane-protein interactions and relate them to some of the molecular strategies organisms use to adjust the membrane lipid composition in response to a decrease in environmental temperature. While the activities of all biomolecules are altered as a function of temperature, the thermosensors we focus on here are molecules whose temperature sensitivity appears to be linked to changes in the biophysical properties of membrane lipids. This article is part of a Special Issue entitled: Lipid-protein interactions.