Energetics of hydrophobic matching in lipid-protein interactions

Probing membrane deformation energy by KcsA potassium channel gating under varied membrane thickness and tension

This study investigated how membrane thickness and tension modify the gating of KcsA potassium channels when simultaneously varied. The KcsA channel undergoes global conformational changes upon gating: expansion of the cross-sectional area and longitudinal shortening upon opening. Thus, membranes impose differential effects on the open and closed conformations, such as hydrophobic mismatches. Here, the single-channel open probability was recorded in the contact bubble bilayer, by which variable thickness membranes under a defined tension were applied. A fully open channel in thin membranes turned to sporadic openings in thick membranes, where the channel responded moderately to tension increase. Quantitative gating analysis prompted the hypothesis that tension augmented the membrane deformation energy when hydrophobic mismatch was enhanced in thick membranes.

Energetics of the Transmembrane Peptide Sorting by Hydrophobic Mismatch

Hydrophobic mismatch between a lipid membrane and embedded transmembrane peptides or proteins plays a role in their lateral localization and function. Earlier studies have resolved numerous mechanisms through which the peptides and membrane proteins adapt to mismatch, yet the energetics of lateral sorting due to hydrophobic mismatch have remained elusive due to the lack of suitable computational or experimental protocols. Here, we pioneer a molecular dynamics simulation approach to study the sorting of peptides along a membrane thickness gradient. Peptides of different lengths tilt and diffuse along the membrane to eliminate mismatch with a rate directly proportional to the magnitude of mismatch. We extract the 2-dimensional free energy profiles as a function of local thickness and peptide orientation, revealing the relative contributions of sorting and tilting, and suggesting their thermally accessible regimes. Our approach can readily be applied to study other membrane systems of biological interest where hydrophobic mismatch, or membrane thickness in general, plays a role.

In Search of a Molecular View of Peptide-Lipid Interactions in Membranes

Lipid bilayer membranes are often represented as a continuous nonpolar slab with a certain thickness bounded by two more polar interfaces. Phenomena such as peptide binding to the membrane surface, folding, insertion, translocation, and diffusion are typically interpreted on the basis of this view. In this Perspective, I argue that this membrane representation as a hydrophobic continuum solvent is not adequate to understand peptide-lipid interactions. Lipids are not small compared to membrane-active peptides: their sizes are similar. Therefore, peptide diffusion needs to be understood in terms of free volume, not classical continuum mechanics; peptide solubility or partitioning in membranes cannot be interpreted in terms of hydrophobic mismatch between membrane thickness and peptide length; peptide folding and translocation, often involving cationic peptides, can only be understood if realizing that lipids adapt to the presence of peptides and the membrane may undergo considerable lipid redistribution in the process. In all of those instances, the detailed molecular interactions between the peptide residues and the lipid components are essential to understand the mechanisms involved.

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.

Length matters: Functional flip of the short TatA transmembrane helix

The twin arginine translocase (Tat) exports folded proteins across bacterial membranes. The putative pore-forming or membrane-weakening component (TatA_d in B. subtilis) is anchored to the lipid bilayer via an unusually short transmembrane α-helix (TMH), with less than 16 residues. Its tilt angle in different membranes was analyzed under hydrophobic mismatch conditions, using synchrotron radiation circular dichroism and solid-state NMR. Positive mismatch (introduced either by reconstitution in short-chain lipids or by extending the hydrophobic TMH length) increased the helix tilt of the TMH as expected. Negative mismatch (introduced either by reconstitution in long-chain lipids or by shortening the TMH), on the other hand, led to protein aggregation. These data suggest that the TMH of TatA is just about long enough for stable membrane insertion. At the same time, its short length is a crucial factor for successful translocation, as demonstrated here in native membrane vesicles using an in vitro translocation assay. Furthermore, when reconstituted in model membranes with negative spontaneous curvature, the TMH was found to be aligned parallel to the membrane surface. This intrinsic ability of TatA to flip out of the membrane core thus seems to play a key role in its membrane-destabilizing effect during Tat-dependent translocation.

Fluorescence microscopy imaging of a neurotransmitter receptor and its cell membrane lipid milieu

Hampered by the diffraction phenomenon, as expressed in 1873 by Abbe, applications of optical microscopy to image biological structures were for a long time limited to resolutions above the ∼200 nm barrier and restricted to the observation of stained specimens. The introduction of fluorescence was a game changer, and since its inception it became the gold standard technique in biological microscopy. The plasma membrane is a tenuous envelope of 4 nm-10 nm in thickness surrounding the cell. Because of its highly versatile spectroscopic properties and availability of suitable instrumentation, fluorescence techniques epitomize the current approach to study this delicate structure and its molecular constituents. The wide spectral range covered by fluorescence, intimately linked to the availability of appropriate intrinsic and extrinsic probes, provides the ability to dissect membrane constituents at the molecular scale in the spatial domain. In addition, the time resolution capabilities of fluorescence methods provide complementary high precision for studying the behavior of membrane molecules in the time domain. This review illustrates the value of various fluorescence techniques to extract information on the topography and motion of plasma membrane receptors. To this end I resort to a paradigmatic membrane-bound neurotransmitter receptor, the nicotinic acetylcholine receptor (nAChR). The structural and dynamic picture emerging from studies of this prototypic pentameric ligand-gated ion channel can be extrapolated not only to other members of this superfamily of ion channels but to other membrane-bound proteins. I also briefly discuss the various emerging techniques in the field of biomembrane labeling with new organic chemistry strategies oriented to applications in fluorescence nanoscopy, the form of fluorescence microscopy that is expanding the depth and scope of interrogation of membrane-associated phenomena.

Membranolytic Mechanism of Amphiphilic Antimicrobial β-Stranded [KL]n Peptides

Amphipathic peptides can act as antibiotics due to membrane permeabilization. KL peptides with the repetitive sequence [Lys-Leu]_n-NH₂ form amphipathic β-strands in the presence of lipid bilayers. As they are known to kill bacteria in a peculiar length-dependent manner, we suggest here several different functional models, all of which seem plausible, including a carpet mechanism, a β-barrel pore, a toroidal wormhole, and a β-helix. To resolve their genuine mechanism, the activity of KL peptides with lengths from 6–26 amino acids (plus some inverted LK analogues) was systematically tested against bacteria and erythrocytes. Vesicle leakage assays served to correlate bilayer thickness and peptide length and to examine the role of membrane curvature and putative pore diameter. KL peptides with 10–12 amino acids showed the best therapeutic potential, i.e., high antimicrobial activity and low hemolytic side effects. Mechanistically, this particular window of an optimum β-strand length around 4 nm (11 amino acids × 3.7 Å) would match the typical thickness of a lipid bilayer, implying the formation of a transmembrane pore. Solid-state ¹⁵N- and ¹⁹F-NMR structure analysis, however, showed that the KL backbone lies flat on the membrane surface under all conditions. We can thus refute any of the pore models and conclude that the KL peptides rather disrupt membranes by a carpet mechanism. The intriguing length-dependent optimum in activity can be fully explained by two counteracting effects, i.e., membrane binding versus amyloid formation. Very short KL peptides are inactive, because they are unable to bind to the lipid bilayer as flexible β-strands, whereas very long peptides are inactive due to vigorous pre-aggregation into β-sheets in solution.

Fluorescence evidence of annexin A6 translocation across membrane in model matrix vesicles during apatite formation

Matrix vesicles (MVs) are 100-300 nm spherical structures released by mineralization competent cells to initiate formation of apatite, the mineral component in bones. Among proteins present in MVs, annexin A6 (AnxA6) is thought to be ubiquitously distributed in the MVs' lumen, on the surface of the internal and external leaflets of the membrane and also inserted in the lipid bilayer. To determine the molecular mechanism(s) that lead to the different locations of AnxA6, we hypothesized the occurrence of a pH drop during the mineralization. Such a change would induce the AnxA6 protonation, which in turn, and because of its isoelectric point of 5.41, would change the protein hydrophobicity facilitating its insertion into the MVs' bilayer. The various distributions of AnxA6 are likely to disturb membrane phospholipid organization. To examine this possibility, we used fluorescein as pH reporter, and established that pH decreased inside MVs during apatite formation. Then, 4-(14-phenyldibenzo[a,c]phenazin-9(14H)-yl)-phenol, a vibration-induced emission fluorescent probe, was used as a reporter of changes in membrane organization occurring with the varying mode of AnxA6 binding. Proteoliposomes containing AnxA6 and 1,2-Dimyristoyl-sn-glycero-3phosphocholine (DMPC) or 1,2-Dimyristoyl-sn-glycero-3phosphocholine: 1,2-Dipalmitoyl-sn-glycero-3-phosphoserine (DMPC:DPPS 9:1), to mimic the external and internal MV membrane leaflet, respectively, served as biomimetic models to investigate the nature of AnxA6 binding. Addition of Anx6 to DMPC at pH 7.4 and 5.4, or DMPC:DPPS (9:1) at pH 7.4 induced a decrease in membrane fluidity, consistent with AnxA6 interactions with the bilayer surface. In contrast, AnxA6 addition to DMPC:DPPS (9:1) at pH 5.4 increased the fluidity of the membrane. This latest result was interpreted as reflecting the insertion of AnxA6 into the bilayer. Taken together, these findings point to a possible mechanism of AnxA6 translocation in MVs from the surface of the internal leaflet into the phospholipid bilayer stimulated upon acidification of the MVs' lumen during formation of apatite.

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.

Amphiphilic Membrane Environments Regulate Enzymatic Behaviors of Salmonella Outer Membrane Protease

The role of an amphiphilic environment in the functional regulation of integral membrane proteins is well appreciated but how specific amphiphilic surrounding influences the conformational plasticity and function of a protein is less obvious. We focus on the Salmonella phosphoglycerate transport system (pgt)-encoded outer membrane protease E (PgtE), which plays an important role in tissue infiltration and survival of Salmonella enterica. Despite our understanding of its physiological functions, elucidation of its enzymatic behavior in response to the immediate amphiphilic surrounding is lacking. We monitor the proteolytic activity of PgtE reconstituted in Zwittergent 3-12 detergent micelles or a 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) bilayer and examine factors that influence its activity. We find, to our surprise, that PgtE, which is thought to elicit a rapid response toward various substrates, showed hysteretic enzymatic behavior, characterized by a prominent lag phase prior to achieving the exponential steady state in its detergent-stabilized form as well as in the outer membrane embedded native state in live bacteria. The lag phase was abolished under three conditions: preformation of an inactive detergent-stabilized PgtE-substrate complex without lipopolysaccharide (LPS), LPS-bound detergent-stabilized PgtE that had reached steady state velocity, or PgtE reconstituted into a POPC bilayer environment. Interestingly, detergent- and bilayer-stabilized PgtE showed comparable steady-state activity. And strikingly, lipopolysaccharide (LPS) becomes nonessential for the activation of PgtE when the protein is reconstituted in the phospholipid bilayer, contrasting a long-standing notion that LPS is required for proteases belonging to the omptin family to be proteolytically active. These findings suggest intriguing biological nuances for the proteolytic function of PgtE that were not well appreciated previously and offer new perspectives that may generally be applicable for omptins.

Thermodynamics-Based Molecular Modeling of α-Helices in Membranes and Micelles

The Folding of Membrane-Associated Peptides (FMAP) method was developed for modeling α-helix formation by linear peptides in micelles and lipid bilayers. FMAP 2.0 identifies locations of α-helices in the amino acid sequence, generates their three-dimensional models in planar bilayers or spherical micelles, and estimates their thermodynamic stabilities and tilt angles, depending on temperature and pH. The method was tested for 723 peptides (926 data points) experimentally studied in different environments and for 170 single-pass transmembrane (TM) proteins with available crystal structures. FMAP 2.0 detected more than 95% of experimentally observed α-helices with an average error in helix end determination of around 2, 3, 4, and 5 residues per helix for peptides in water, micelles, bilayers, and TM proteins, respectively. Helical and nonhelical residue states were predicted with an accuracy from 0.86 to 0.96, and the Matthews correlation coefficient was from 0.64 to 0.88 depending on the environment. Experimental micelle- and membrane-binding energies and tilt angles of peptides were reproduced with a root-mean-square deviation of around 2 kcal/mol and 7°, respectively. The TM and non-TM states of hydrophobic and pH-triggered α-helical peptides in various lipid bilayers were reproduced in more than 95% of cases. The FMAP 2.0 web server (https://membranome.org/fmap) is publicly available to explore the structural polymorphism of antimicrobial, cell-penetrating, fusion, and other membrane-binding peptides, which is important for understanding the mechanisms of their biological activities.

Unusual mode of dimerization of retinitis pigmentosa-associated F220C rhodopsin

Mutations in the G protein-coupled receptor (GPCR) rhodopsin are a common cause of autosomal dominant retinitis pigmentosa, a blinding disease. Rhodopsin self-associates in the membrane, and the purified monomeric apo-protein opsin dimerizes in vitro as it transitions from detergent micelles to reconstitute into a lipid bilayer. We previously reported that the retinitis pigmentosa-linked F220C opsin mutant fails to dimerize in vitro, reconstituting as a monomer. Using fluorescence-based assays and molecular dynamics simulations we now report that whereas wild-type and F220C opsin display distinct dimerization propensities in vitro as previously shown, they both dimerize in the plasma membrane of HEK293 cells. Unexpectedly, molecular dynamics simulations show that F220C opsin forms an energetically favored dimer in the membrane when compared with the wild-type protein. The conformation of the F220C dimer is unique, with transmembrane helices 5 and 6 splayed apart, promoting widening of the intracellular vestibule of each protomer and influx of water into the protein interior. FRET experiments with SNAP-tagged wild-type and F220C opsin expressed in HEK293 cells are consistent with this conformational difference. We speculate that the unusual mode of dimerization of F220C opsin in the membrane may have physiological consequences.

Overlapping Properties of the Short Membrane-Active Peptide BP100 With (i) Polycationic TAT and (ii) α-helical Magainin Family Peptides

BP100 is a short, designer-made membrane-active peptide with multiple functionalities: antimicrobial, cell-penetrating, and fusogenic. Consisting of five lysines and 6 hydrophobic residues, BP100 was shown to bind to lipid bilayers as an amphipathic α-helix, but its mechanism of action remains unclear. With these features, BP100 embodies the characteristics of two distinctly different classes of membrane-active peptides, which have been studied in detail and where the mechanism of action is better understood. On the one hand, its amphiphilic helical structure is similar to the pore forming magainin family of antimicrobial peptides, though BP100 is much too short to span the membrane. On the other hand, its length and high charge density are reminiscent of the HIV-TAT family of cell penetrating peptides, for which inverted micelles have been postulated as translocation intermediates, amongst other mechanisms. Assays were performed to test the antimicrobial and hemolytic activity, the induced leakage and fusion of lipid vesicles, and cell uptake. From these results the functional profiles of BP100, HIV-TAT, and the magainin-like peptides magainin 2, PGLa, MSI-103, and MAP were determined and compared. It is observed that the activity of BP100 resembles most closely the much longer amphipathic α-helical magainin-like peptides, with high antimicrobial activity along with considerable fusogenic and hemolytic effects. In contrast, HIV-TAT shows almost no antimicrobial, fusogenic, or hemolytic effects. We conclude that the amphipathic helix of BP100 has a similar membrane-based activity as magainin-like peptides and may have a similar mechanism of action.

Tylopeptin B peptide antibiotic in lipid membranes at low concentrations: Self-assembling, mutual repulsion and localization

The medium-length peptide Tylopeptin B possesses activity against Gram-positive bacteria. It binds to bacterial membranes altering their mechanical properties and increasing their permeability. This action is commonly related with peptide self-assembling, resulting in the formation of membrane channels. Here, pulsed double electron-electron resonance (DEER) data for spin-labeled Tylopeptin B in palmitoyl-oleoyl-glycero-phosphocholine (POPC) model membrane reveal that peptide self-assembling starts at concentration as low as 0.1 mol%; above 0.2 mol% it attains a saturation-like dependence with a mean number of peptides in the cluster <n> = 3.3. Using the electron spin echo envelope modulation (ESEEM) technique, Tylopeptin B molecules are found to possess a planar orientation in the membrane. In the peptide concentration range between 0.1 and 0.2 mol%, DEER data show that the peptide clusters have tendency of mutual repulsion, with a circle of inaccessibility of radius around 20 nm. It may be proposed that within this radius the peptides destabilize the membrane, providing so the peptide antimicrobial activity. Exploiting spin-labeled stearic acids as a model for free fatty acids (FFA), we found that at concentrations of 0.1-0.2 mol% the peptide promotes formation of lipid-mediated FFA clusters; further increase in peptide concentration results in dissipation of these clusters.

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.

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.

Confinement in Nanodiscs Anisotropically Modifies Lipid Bilayer Elastic Properties

Lipid nanodiscs are small synthetic lipid bilayer structures that are stabilized in solution by special circumscribing (or scaffolding) proteins or polymers. Because they create native-like environments for transmembrane proteins, lipid nanodiscs have become a powerful tool for structural determination of this class of systems when combined with cryo-electron microscopy or nuclear magnetic resonance. The elastic properties of lipid bilayers determine how the lipid environment responds to membrane protein perturbations, and how the lipid in turn modifies the conformational state of the embedded protein. However, despite the abundant use of nanodiscs in determining membrane protein structure, the elastic material properties of even pure lipid nanodiscs (i.e., without embedded proteins) have not yet been quantitatively investigated. A major hurdle is due to the inherently nonlocal treatment of the elastic properties of lipid systems implemented by most existing methods, both experimental and computational. In addition, these methods are best suited for very large “infinite” size lipidic assemblies, or ones that contain periodicity, in the case of simulations. We have previously described a computational analysis of molecular dynamics simulations designed to overcome these limitations, so it allows quantification of the bending rigidity (K_C) and tilt modulus (κ_t) on a local scale even for finite, nonperiodic systems, such as lipid nanodiscs. Here we use this computational approach to extract values of K_C and κ_t for a set of lipid nanodisc systems that vary in size and lipid composition. We find that the material properties of lipid nanodiscs are different from those of infinite bilayers of corresponding lipid composition, highlighting the effect of nanodisc confinement. Nanodiscs tend to show higher stiffness than their corresponding macroscopic bilayers, and moreover, their material properties vary spatially within them. For small-size MSP1 nanodiscs, the stiffness decreases radially, from a value that is larger in their center than the moduli of the corresponding bilayers by a factor of ∼2–3. The larger nanodiscs (MSP1E3D1 and MSP2N2) show milder spatial changes of moduli that are composition dependent and can be maximal in the center or at some distance from it. These trends in moduli correlate with spatially varying structural properties, including the area per lipid and the nanodisc thickness. Finally, as has previously been reported, nanodiscs tend to show deformations from perfectly flat circular geometries to varying degrees, depending on size and lipid composition. The modulations of lipid elastic properties that we find should be carefully considered when making structural and functional inferences concerning embedded proteins.

Roles of Membrane Domains in Integrin-Mediated Cell Adhesion

The composition and organization of the plasma membrane play important functional and regulatory roles in integrin signaling, which direct many physiological and pathological processes, such as development, wound healing, immunity, thrombosis, and cancer metastasis. Membranes are comprised of regions that are thick or thin owing to spontaneous partitioning of long-chain saturated lipids from short-chain polyunsaturated lipids into domains defined as ordered and liquid-disorder domains, respectively. Liquid-ordered domains are typically 100 nm in diameter and sometimes referred to as lipid rafts. We posit that integrin β senses membrane thickness and that mechanical force on the membrane regulates integrin activation through membrane thinning. This review examines what we know about the nature and mechanism of the interaction of integrins with the plasma membrane and its effects on regulating integrins and its binding partners.

Role of the lipid bilayer in outer membrane protein folding in Gram-negative bacteria

β-Barrel outer membrane proteins (OMPs) represent the major proteinaceous component of the outer membrane (OM) of Gram-negative bacteria. These proteins perform key roles in cell structure and morphology, nutrient acquisition, colonization and invasion, and protection against external toxic threats such as antibiotics. To become functional, OMPs must fold and insert into a crowded and asymmetric OM that lacks much freely accessible lipid. This feat is accomplished in the absence of an external energy source and is thought to be driven by the high thermodynamic stability of folded OMPs in the OM. With such a stable fold, the challenge that bacteria face in assembling OMPs into the OM is how to overcome the initial energy barrier of membrane insertion. In this review, we highlight the roles of the lipid environment and the OM in modulating the OMP-folding landscape and discuss the factors that guide folding in vitro and in vivo We particularly focus on the composition, architecture, and physical properties of the OM and how an understanding of the folding properties of OMPs in vitro can help explain the challenges they encounter during folding in vivo Current models of OMP biogenesis in the cellular environment are still in flux, but the stakes for improving the accuracy of these models are high. OMP folding is an essential process in all Gram-negative bacteria, and considering the looming crisis of widespread microbial drug resistance it is an attractive target. To bring down this vital OMP-supported barrier to antibiotics, we must first understand how bacterial cells build it.

Interaction of the Mechanosensitive Channel, MscS, with the Membrane Bilayer through Lipid Intercalation into Grooves and Pockets

All membrane proteins have dynamic and intimate relationships with the lipids of the bilayer that may determine their activity. Mechanosensitive channels sense tension through their interaction with the lipids of the membrane. We have proposed a mechanism for the bacterial channel of small conductance, MscS, that envisages variable occupancy of pockets in the channel by lipid chains. Here, we analyze protein-lipid interactions for MscS by quenching of tryptophan fluorescence with brominated lipids. By this strategy, we define the limits of the bilayer for TM1, which is the most lipid exposed helix of this protein. In addition, we show that residues deep in the pockets, created by the oligomeric assembly, interact with lipid chains. On the cytoplasmic side, lipids penetrate as far as the pore-lining helices and lipid molecules can align along TM3b perpendicular to lipids in the bilayer. Cardiolipin, free fatty acids, and branched lipids can access the pockets where the latter have a distinct effect on function. Cholesterol is excluded from the pockets. We demonstrate that introduction of hydrophilic residues into TM3b severely impairs channel function and that even "conservative" hydrophobic substitutions can modulate the stability of the open pore. The data provide important insights into the interactions between phospholipids and MscS and are discussed in the light of recent developments in the study of Piezo1 and TrpV4.

The desmosome is a mesoscale lipid raft-like membrane domain

Desmogleins (Dsgs) are cadherin family adhesion molecules essential for epidermal integrity. Previous studies have shown that desmogleins associate with lipid rafts, but the significance of this association was not clear. Here, we report that the desmoglein transmembrane domain (TMD) is the primary determinant of raft association. Further, we identify a novel mutation in the DSG1 TMD (G562R) that causes severe dermatitis, multiple allergies, and metabolic wasting syndrome. Molecular modeling predicts that this G-to-R mutation shortens the DSG1 TMD, and experiments directly demonstrate that this mutation compromises both lipid raft association and desmosome incorporation. Finally, cryo-electron tomography indicates that the lipid bilayer within the desmosome is ∼10% thicker than adjacent regions of the plasma membrane. These findings suggest that differences in bilayer thickness influence the organization of adhesion molecules within the epithelial plasma membrane, with cadherin TMDs recruited to the desmosome via the establishment of a specialized mesoscale lipid raft-like membrane domain.

Perturbations of Native Membrane Protein Structure in Alkyl Phosphocholine Detergents: A Critical Assessment of NMR and Biophysical Studies

Membrane proteins perform a host of vital cellular functions. Deciphering the molecular mechanisms whereby they fulfill these functions requires detailed biophysical and structural investigations. Detergents have proven pivotal to extract the protein from its native surroundings. Yet, they provide a milieu that departs significantly from that of the biological membrane, to the extent that the structure, the dynamics, and the interactions of membrane proteins in detergents may considerably vary, as compared to the native environment. Understanding the impact of detergents on membrane proteins is, therefore, crucial to assess the biological relevance of results obtained in detergents. Here, we review the strengths and weaknesses of alkyl phosphocholines (or foscholines), the most widely used detergent in solution-NMR studies of membrane proteins. While this class of detergents is often successful for membrane protein solubilization, a growing list of examples points to destabilizing and denaturing properties, in particular for α-helical membrane proteins. Our comprehensive analysis stresses the importance of stringent controls when working with this class of detergents and when analyzing the structure and dynamics of membrane proteins in alkyl phosphocholine detergents.

Emerging roles for sphingolipids in cellular aging

Aging is a gradual loss of physiological functions as organisms' progress in age. Although aging in multicellular organisms is complex, some fundamental mechanisms and pathways may be shared from the single cellular yeast to human. Budding yeast Saccharomyces cerevisiae has been established model system for aging studies. A yeast cell divides asymmetrically to produce two cells that differ in size and age. The one that is smaller coming from bud is a newborn cell that with a full replicative potential head irrespective of the replicative age of its mother-the larger cell from which the bud grows out before division. The age asymmetry between daughter and mother is thought to be dependent on asymmetric segregation of certain factors such as protein aggregates, extrachromosomal DNA (ERCs) and dysfunctional organelles during successive cell divisions of the yeast replicative lifespan (RLS). It is also thought that certain plasma membrane proteins, in particular multidrug-resistant (MDR) proteins, asymmetrically partition between the mother and the bud based on the age of the polypeptides. Functional decline associated with the molecular aging of those proteins contributes to the fitness decline at advance age. In our recent study, we showed that sphingolipids facilitate the age-dependent segregation of MDRs between daughter and mother cell. In this review, we highlight and discuss the potential mechanisms by which sphingolipids regulate the aging process in yeast and cells of vertebrate animals including human.

An engineered opsin monomer scrambles phospholipids

The G protein-coupled receptor opsin is a phospholipid scramblase that facilitates rapid transbilayer phospholipid exchange in liposomes. The mechanism by which opsin scrambles lipids is unknown. It has been proposed that lipid translocation may occur at protein-protein interfaces of opsin dimers. To test this possibility, we rationally engineered QUAD opsin by tryptophan substitution of four lipid-facing residues in transmembrane helix 4 (TM4) that is known to be important for dimerization. Atomistic molecular dynamics simulations of wild type and QUAD opsins combined with continuum modeling revealed that the tryptophan substitutions lower the energetically unfavorable residual hydrophobic mismatch between TM4 and the membrane, reducing the drive of QUAD opsin to dimerize. We purified thermostable wild type and QUAD opsins, with or without a SNAP tag for fluorescence labeling. Single molecule fluorescence measurements of purified SNAP-tagged constructs revealed that both proteins are monomers. Fluorescence-based activity assays indicated that QUAD opsin is a fully functional scramblase. However, unlike wild type opsin which dimerizes en route to insertion into phospholipid vesicles, QUAD opsin reconstitutes as a monomer. We conclude that an engineered opsin monomer can scramble phospholipids, and that the lipid-exposed face of TM4 is unlikely to contribute to transbilayer phospholipid exchange.

Complex biomembrane mimetics on the sub-nanometer scale

Biomimetic lipid vesicles are indispensable tools for gaining insight into the biophysics of cell physiology on the molecular level. The level of complexity of these model systems has steadily increased, and now spans from domain-forming lipid mixtures to asymmetric lipid bilayers. Here, we review recent progress in the development and application of elastic neutron and X-ray scattering techniques for studying these systems in situ and under physiologically relevant conditions on the nanometer to sub-nanometer length scales. In particular, we focus on: (1) structural details of coexisting liquid-ordered and liquid-disordered domains, including their thickness and lipid packing mismatch as a function of a size transition from nanoscopic to macroscopic domains; (2) membrane-mediated protein partitioning into lipid domains; (3) the role of the aqueous medium in tuning interactions between membranes and domains; and (4) leaflet-specific structure in asymmetric bilayers and passive lipid flip-flop.

Influence of the Length and Charge on the Activity of α-Helical Amphipathic Antimicrobial Peptides

Hydrophobic mismatch is important for pore-forming amphipathic antimicrobial peptides, as demonstrated recently [Grau-Campistany, A., et al. (2015) Sci. Rep. 5, 9388]. A series of different length peptides have been generated with the heptameric repeat sequence KIAGKIA, called KIA peptides, and it was found that only those helices sufficiently long to span the hydrophobic thickness of the membrane could induce leakage in lipid vesicles; there was also a clear length dependence of the antimicrobial and hemolytic activities. For the original KIA sequences, the cationic charge increased with peptide length. The goal of this work is to examine whether the charge also has an effect on activity; hence, we constructed two further series of peptides with a sequence similar to those of the KIA peptides, but with a constant charge of +7 for all lengths from 14 to 28 amino acids. For both of these new series, a clear length dependence similar to that of KIA peptides was observed, indicating that charge has only a minor influence. Both series also showed a distinct threshold length for peptides to be active, which correlates directly with the thickness of the membrane. Among the longer peptides, the new series showed activities only slightly lower than those of the original KIA peptides of the same length that had a higher charge. Shorter peptides, in which Gly was replaced with Lys, showed activities similar to those of KIA peptides of the same length, but peptides in which Ile was replaced with Lys lost their helicity and were less active.

Lipoprotein insertion into membranes of various complexity: lipid sorting, interfacial adsorption and protein clustering

In a combined chemical-biological and biophysical approach we explored the membrane partitioning of the lipidated signaling proteins N-Ras and K-Ras4B into membrane systems of different complexity, ranging from one-component lipid bilayers and anionic binary and ternary heterogeneous membrane systems even up to partitioning studies on protein-free and protein-containing giant plasma membrane vesicles (GPMVs). To yield a pictorial view of the localization process, imaging using confocal laser scanning and atomic force microscopy was performed. The results reveal pronounced isoform-specific differences regarding the lateral distribution and formation of protein-rich membrane domains. Line tension is one of the key parameters controlling not only the size and dynamic properties of segregated lipid domains but also the partitioning process of N-Ras that acts as a lineactant. The formation of N-Ras protein clusters is even recorded for almost vanishing hydrophobic mismatch. Conversely, for K-Ras4B, selective localization and clustering are electrostatically mediated by its polybasic farnesylated C-terminus. The formation of K-Ras4B clusters is also observed for the multi-component GPMV membrane, i.e., it seems to be a general phenomenon, largely independent of the details of the membrane composition, including the anionic charge density of lipid headgroups. Our data indicate that unspecific and entropy-driven membrane-mediated interactions play a major role in the partitioning behavior, thus relaxing the need for a multitude of fine-tuned interactions. Such a scenario seems also to be reasonable recalling the high dynamic nature of cellular membranes. Finally, we note that even relatively simple models of heterogeneous membranes are able to reproduce many of the properties of much more complex biological membranes.

Interplay of G Protein-Coupled Receptors with the Membrane: Insights from Supra-Atomic Coarse Grain Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are central to many fundamental cellular signaling pathways. They transduce signals from the outside to the inside of cells in physiological processes ranging from vision to immune response. It is extremely challenging to look at them individually using conventional experimental techniques. Recently, a pseudo atomistic molecular model has emerged as a valuable tool to access information on GPCRs, more specifically on their interactions with their environment in their native cell membrane and the consequences on their supramolecular organization. This approach uses the Martini coarse grain (CG) model to describe the receptors, lipids, and solvent in molecular dynamics (MD) simulations and in enough detail to allow conserving the chemical specificity of the different molecules. The elimination of unnecessary degrees of freedom has opened up large-scale simulations of the lipid-mediated supramolecular organization of GPCRs. Here, after introducing the Martini CGMD method, we review these studies carried out on various members of the GPCR family, including rhodopsin (visual receptor), opioid receptors, adrenergic receptors, adenosine receptors, dopamine receptor, and sphingosine 1-phosphate receptor. These studies have brought to light an interesting set of novel biophysical principles. The insights range from revealing localized and heterogeneous deformations of the membrane bilayer at the surface of the protein, specific interactions of lipid molecules with individual GPCRs, to the effect of the membrane matrix on global GPCR self-assembly. The review ends with an overview of the lessons learned from the use of the CGMD method, the biophysical-chemical findings on lipid-protein interplay.

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.

Hydrophobic Mismatch Drives the Interaction of E5 with the Transmembrane Segment of PDGF Receptor

The oncogenic E5 protein from bovine papillomavirus is a short (44 amino acids long) integral membrane protein that forms homodimers. It activates platelet-derived growth factor receptor (PDGFR) β in a ligand-independent manner by transmembrane helix-helix interactions. The nature of this recognition event remains elusive, as numerous mutations are tolerated in the E5 transmembrane segment, with the exception of one hydrogen-bonding residue. Here, we examined the conformation, stability, and alignment of the E5 protein in fluid lipid membranes of substantially varying bilayer thickness, in both the absence and presence of the PDGFR transmembrane segment. Quantitative synchrotron radiation circular dichroism analysis revealed a very long transmembrane helix for E5 of ∼26 amino acids. Oriented circular dichroism and solid-state ¹⁵N-NMR showed that the alignment and stability of this unusually long segment depend critically on the membrane thickness. When reconstituted alone in exceptionally thick DNPC lipid bilayers, the E5 helix was found to be inserted almost upright. In moderately thick bilayers (DErPC and DEiPC), it started to tilt and became slightly deformed, and finally it became aggregated in conventional DOPC, POPC, and DMPC membranes due to hydrophobic mismatch. On the other hand, when E5 was co-reconstituted with the transmembrane segment of PDGFR, it was able to tolerate even the most pronounced mismatch and was stabilized by binding to the receptor, which has the same hydrophobic length. As E5 is known to activate PDGFR within the thin membranes of the Golgi compartment, we suggest that the intrinsic hydrophobic mismatch of these two interaction partners drives them together. They seem to recognize each other by forming a closely packed bundle of mutually aligned transmembrane helices, which is further stabilized by a specific pair of hydrogen-bonding residues.

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.

On scattered waves and lipid domains: detecting membrane rafts with X-rays and neutrons

In order to understand the biological role of lipids in cell membranes, it is necessary to determine the mesoscopic structure of well-defined model membrane systems. Neutron and X-ray scattering are non-invasive, probe-free techniques that have been used extensively in such systems to probe length scales ranging from angstroms to microns, and dynamics occurring over picosecond to millisecond time scales. Recent developments in the area of phase separated lipid systems mimicking membrane rafts will be presented, and the underlying concepts of the different scattering techniques used to study them will be discussed in detail.

Assessment of the functionality and stability of detergent purified nAChR from Torpedo using lipidic matrixes and macroscopic electrophysiology

In our previous study we examined the functionality and stability of nicotinic acetylcholine receptor (nAChR)-detergent complexes (nAChR-DCs) from affinity-purified Torpedo californica (Tc) using fluorescence recovery after photobleaching (FRAP) in Lipidic Cubic Phase (LCP) and planar lipid bilayer (PLB) recordings for phospholipid and cholesterol like detergents. In the present study we enhanced the functional characterization of nAChR-DCs by recording macroscopic ion channel currents in Xenopus oocytes using the two electrode voltage clamp (TEVC). The use of TEVC allows for the recording of macroscopic currents elicited by agonist activation of nAChR-DCs that assemble in the oocyte plasma membrane. Furthermore, we examined the stability of nAChR-DCs, which is obligatory for the nAChR crystallization, using a 30 day FRAP assay in LCP for each detergent. The present results indicate a marked difference in the fractional fluorescence recovery (ΔFFR) within the same detergent family during the 30 day period assayed. Within the cholesterol analog family, sodium cholate and CHAPSO displayed a minimum ΔFFR and a mobile fraction (MF) over 80%. In contrast, CHAPS and BigCHAP showed a marked decay in both the mobile fraction and diffusion coefficient. nAChR-DCs containing phospholipid analog detergents with an alkylphosphocholine (FC) and lysofoscholine (LFC) of 16 carbon chains (FC-16, LFC-16) were more effective in maintaining a mobile fraction of over 80% compared to their counterparts with shorter acyl chain (C12, C14). The significant differences in macroscopic current amplitudes, activation and desensitization rates among the different nAChR-DCs evaluated in the present study allow to dissect which detergent preserves both, agonist activation and ion channel function. Functionality assays using TEVC demonstrated that LFC16, LFC14, and cholate were the most effective detergents in preserving macroscopic ion channel function, however, the nAChR-cholate complex display a significant delay in the ACh-induce channel activation. In summary, these results suggest that the physical properties of the lipid analog detergents (headgroup and acyl chain length) are the most effective in maintaining both the stability and functionality of the nAChR in the detergent solubilized complex.

Hydrophobic mismatch demonstrated for membranolytic peptides, and their use as molecular rulers to measure bilayer thickness in native cells

Hydrophobic mismatch is a well-recognized principle in the interaction of transmembrane proteins with lipid bilayers. This concept was extended here to amphipathic membranolytic α-helices. Nine peptides with lengths between 14 and 28 amino acids were designed from repeated KIAGKIA motifs, and their helical nature was confirmed by circular dichroism spectroscopy. Biological assays for antimicrobial activity and hemolysis, as well as fluorescence vesicle leakage and solid-state NMR spectroscopy, were used to correlate peptide length with membranolytic activity. These data show that the formation of transmembrane pores is only possible under the condition of hydrophobic matching: the peptides have to be long enough to span the hydrophobic bilayer core to be able to induce vesicle leakage, kill bacteria, and cause hemolysis. By correlating the threshold lengths for biological activity with the biophysical results on model vesicles, the peptides could be utilized as molecular rulers to measure the membrane thickness in different cells.

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.

AMPs and OMPs: Is the folding and bilayer insertion of β-stranded outer membrane proteins governed by the same biophysical principles as for α-helical antimicrobial peptides?

The folding and function of membrane proteins is controlled not only by specific but also by unspecific interactions with the constituent lipids. In this review, we focus on the influence of the spontaneous lipid curvature on the folding and insertion of peptides and proteins in membranes. Amphiphilic α-helical peptides, as represented by various antimicrobial sequences, are compared with β-barrel proteins, which are found in the outer membrane of Gram-negative bacteria. It has been shown that cationic amphiphilic peptides are always surface-bound in lipids with a negative spontaneous curvature like POPC, i.e. they are oriented parallel to the membrane plane. On the other hand, in lipids like DMPC with a positive curvature, these peptides can get tilted or completely inserted in a transmembrane state. Remarkably, the folding and spontaneous membrane insertion of β-barrel outer membrane proteins also proceeds more easily in lipids with a positive intrinsic curvature, while it is hampered by negative curvature. We therefore propose that a positive spontaneous curvature of the lipids promotes the ability of a surface-bound molecule to insert more deeply into the bilayer core, irrespective of the conformation, size, or shape of the peptide, protein, or folding intermediate. This article is part of a Special Issue entitled: Lipid-protein interactions.

Life at the border: adaptation of proteins to anisotropic membrane environment

This review discusses main features of transmembrane (TM) proteins which distinguish them from water-soluble proteins and allow their adaptation to the anisotropic membrane environment. We overview the structural limitations on membrane protein architecture, spatial arrangement of proteins in membranes and their intrinsic hydrophobic thickness, co-translational and post-translational folding and insertion into lipid bilayers, topogenesis, high propensity to form oligomers, and large-scale conformational transitions during membrane insertion and transport function. Special attention is paid to the polarity of TM protein surfaces described by profiles of dipolarity/polarizability and hydrogen-bonding capacity parameters that match polarity of the lipid environment. Analysis of distributions of Trp resides on surfaces of TM proteins from different biological membranes indicates that interfacial membrane regions with preferential accumulation of Trp indole rings correspond to the outer part of the lipid acyl chain region-between double bonds and carbonyl groups of lipids. These "midpolar" regions are not always symmetric in proteins from natural membranes. We also examined the hydrophobic effect that drives insertion of proteins into lipid bilayer and different free energy contributions to TM protein stability, including attractive van der Waals forces and hydrogen bonds, side-chain conformational entropy, the hydrophobic mismatch, membrane deformations, and specific protein-lipid binding.

Not just an oil slick: how the energetics of protein-membrane interactions impacts the function and organization of transmembrane proteins

The membrane environment, its composition, dynamics, and remodeling, have been shown to participate in the function and organization of a wide variety of transmembrane (TM) proteins, making it necessary to study the molecular mechanisms of such proteins in the context of their membrane settings. We review some recent conceptual advances enabling such studies, and corresponding computational models and tools designed to facilitate the concerted experimental and computational investigation of protein-membrane interactions. To connect productively with the high resolution achieved by cognate experimental approaches, the computational methods must offer quantitative data at an atomistically detailed level. We show how such a quantitative method illuminated the mechanistic importance of a structural characteristic of multihelical TM proteins, that is, the likely presence of adjacent polar and hydrophobic residues at the protein-membrane interface. Such adjacency can preclude the complete alleviation of the well-known hydrophobic mismatch between TM proteins and the surrounding membrane, giving rise to an energy cost of residual hydrophobic mismatch. The energy cost and biophysical formulation of hydrophobic mismatch and residual hydrophobic mismatch are reviewed in the context of their mechanistic role in the function of prototypical members of multihelical TM protein families: 1), LeuT, a bacterial homolog of mammalian neurotransmitter sodium symporters; and 2), rhodopsin and the β1- and β2-adrenergic receptors from the G-protein coupled receptor family. The type of computational analysis provided by these examples is poised to translate the rapidly growing structural data for the many TM protein families that are of great importance to cell function into ever more incisive insights into mechanisms driven by protein-ligand and protein-protein interactions in the membrane environment.

Protein and lipid interactions driving molecular mechanisms of in meso crystallization

The recent advances in the in meso crystallization technique for the structural characterization of G-protein coupled receptor (GPCR) proteins have established the usefulness of the lipidic-cubic phases (LCPs) in the field of crystallography of membrane proteins. It is surprising that despite the success of the approach, the molecular mechanisms of the in meso method are still not well understood. Therefore, the approach must rely on extensive screening for a suitable protein construct, for host and additive lipids, and for the appropriate precipitants and temperature. To shed light on the in meso crystallization mechanisms, we used extensive coarse-grained molecular dynamics simulations to study, in molecular detail, LCPs under different conditions (compositions and temperatures relevant to crystallogenesis) and their interactions with different types of GPCR constructs. The results presented show how the modulation of the lattice constant of the LCP (triggered by the addition of precipitant during the in meso assay), or of the host lipid type, can destabilize monomeric proteins in the bilayer of the LCP and thus drive their aggregation into the stacked lamellae, where the residual hydrophobic mismatch between the protein and the membrane can drive the formation of lateral contacts leading to nucleation and crystal growth. Moreover, we demonstrate how particular protein designs (such as transmembrane proteins engineered to contain large polar regions) can promote protein stacking interactions in the third, out-of-plane, dimension. The insights provided by the new aspects of the specific molecular mechanisms responsible for protein-protein interactions inside the cubic phase presented here should be helpful in guiding the rational design of future in meso trials with successful outcomes.

How the dynamic properties and functional mechanisms of GPCRs are modulated by their coupling to the membrane environment

Experimental observations of the dependence of function and organization of G protein-coupled receptors (GPCRs) on their lipid environment have stimulated new quantitative studies of the coupling between the proteins and the membrane. It is important to develop such a quantitative understanding at the molecular level because the effects of the coupling are seen to be physiologically and clinically significant. Here we review findings that offer insight into how membrane-GPCR coupling is connected to the structural characteristics of the GPCR, from sequence to 3D structural detail, and how this coupling is involved in the actions of ligands on the receptor. The application of a recently developed computational approach designed for quantitative evaluation of membrane remodeling and the energetics of membrane-protein interactions brings to light the importance of the radial asymmetry of the membrane-facing surface of GPCRs in their interaction with the surrounding membrane. As the radial asymmetry creates adjacencies of hydrophobic and polar residues at specific sites of the GPCR, the ability of membrane remodeling to achieve complete hydrophobic matching is limited, and the residual mismatch carries a significant energy cost. The adjacencies are shown to be affected by ligand-induced conformational changes. Thus, functionally important organization of GPCRs in the cell membrane can depend both on ligand-determined properties and on the lipid composition of various membrane regions with different remodeling capacities. That this functionally important reorganization can be driven by oligomerization patterns that reduce the energy cost of the residual mismatch, suggests a new perspective on GPCR dimerization and ligand-GPCR interactions. The relation between the modulatory effects on GPCRs from the binding of specific cell-membrane components, e.g., cholesterol, and those produced by the non-local energetics of hydrophobic mismatch are discussed in this context.

Membrane driven spatial organization of GPCRs

Spatial organization of G-protein coupled receptors (GPCRs) into dimers and higher order oligomers has been demonstrated in vitro and in vivo. The pharmacological readout was shown to depend on the specific interfaces, but why particular regions of the GPCR structure are involved, and how ligand-determined states change them remains unknown. Here we show why protein-membrane hydrophobic matching is attained upon oligomerization at specific interfaces from an analysis of coarse-grained molecular dynamics simulations of the spontaneous diffusion-interaction of the prototypical beta2-adrenergic (β₂AR) receptors in a POPC lipid bilayer. The energy penalty from mismatch is significantly reduced in the spontaneously emerging oligomeric arrays, making the spatial organization of the GPCRs dependent on the pattern of mismatch in the monomer. This mismatch pattern is very different for β₂AR compared to the highly homologous and structurally similar β₁AR, consonant with experimentally observed oligomerization patterns of β₂AR and β₁AR. The results provide a mechanistic understanding of the structural context of oligomerization.

Helical membrane protein conformations and their environment

Evidence that membrane proteins respond conformationally and functionally to their environment is growing. Structural models, by necessity, have been characterized in preparations where the protein has been removed from its native environment. Different structural methods have used various membrane mimetics that have recently included lipid bilayers as a more native-like environment. Structural tools applied to lipid bilayer-embedded integral proteins are informing us about important generic characteristics of how membrane proteins respond to the lipid environment as compared with their response to other nonlipid environments. Here, we review the current status of the field, with specific reference to observations of some well-studied α-helical membrane proteins, as a starting point to aid the development of possible generic principles for model refinement.

Canonical azimuthal rotations and flanking residues constrain the orientation of transmembrane helices

In biological membranes the alignment of embedded proteins provides crucial structural information. The transmembrane (TM) parts have well-defined secondary structures, in most cases α-helices and their orientation is given by a tilt angle and an azimuthal rotation angle around the main axis. The tilt angle is readily visualized and has been found to be functionally relevant. However, there exist no general concepts on the corresponding azimuthal rotation. Here, we show that TM helices prefer discrete rotation angles. They arise from a combination of intrinsic properties of the helix geometry plus the influence of the position and type of flanking residues at both ends of the hydrophobic core. The helical geometry gives rise to canonical azimuthal angles for which the side chains of residues from the two ends of the TM helix tend to have maximum or minimum immersion within the membrane. This affects the preferential position of residues that fall near hydrophobic/polar interfaces of the membrane, depending on their hydrophobicity and capacity to form specific anchoring interactions. On this basis, we can explain the orientation and dynamics of TM helices and make accurate predictions, which correspond well to the experimental values of several model peptides (including dimers), and TM segments of polytopic membrane proteins.

The cost of living in the membrane: a case study of hydrophobic mismatch for the multi-segment protein LeuT

Many observations of the role of the membrane in the function and organization of transmembrane (TM) proteins have been explained in terms of hydrophobic mismatch between the membrane and the inserted protein. For a quantitative investigation of this mechanism in the lipid-protein interactions of functionally relevant conformations adopted by a multi-TM segment protein, the bacterial leucine transporter (LeuT), we employed a novel method, Continuum-Molecular Dynamics (CTMD), that quantifies the energetics of hydrophobic mismatch by combining the elastic continuum theory of membrane deformations with an atomistic level description of the radially asymmetric membrane-protein interface from MD simulations. LeuT has been serving as a model for structure-function studies of the mammalian neurotransmitter:sodium symporters (NSSs), such as the dopamine and serotonin transporters, which are the subject of intense research in the field of neurotransmission. The membrane models in which LeuT was embedded for these studies were composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid, or 3:1 mixture of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG) lipids. The results show that deformation of the host membrane alone is not sufficient to alleviate the hydrophobic mismatch at specific residues of LeuT. The calculations reveal significant membrane thinning and water penetration due to the specific local polar environment produced by the charged K288 of TM7 in LeuT, that is membrane-facing deep inside the hydrophobic milieu of the membrane. This significant perturbation is shown to result in unfavorable polar-hydrophobic interactions at neighboring hydrophobic residues in TM1a and TM7. We show that all the effects attributed to the K288 residue (membrane thinning, water penetration, and the unfavorable polar-hydrophobic interactions at TM1a and TM7), are abolished in calculations with the K288A mutant. The involvement of hydrophobic mismatch is somewhat different in the functionally distinct conformations (outward-open, occluded, inward-open) of LeuT, and the differences are shown to connect to structural elements (e.g., TM1a) known to play key roles in transport. This finding suggests a mechanistic hypothesis for the enhanced transport activity observed for the K288A mutant, suggesting that the unfavorable hydrophobic-hydrophilic interactions hinder the motion of TM1a in the functionally relevant conformational transition to the inward-open state. Various extents of such unfavorable interactions, involving exposure to the lipid environment of adjacent hydrophobic and polar residues, are common in multi-segment transmembrane proteins, and must be considered to affect functionally relevant conformational transitions.