Emerging Diversity in Lipid-Protein Interactions

Targeting thermoTRP ion channels: in silico preclinical approaches and opportunities

INTRODUCTION

A myriad of cellular pathophysiological responses are mediated by polymodal ion channels that respond to chemical and physical stimuli such as thermoTRP channels. Intriguingly, these channels are pivotal therapeutic targets with limited clinical pharmacology. In silico methods offer an unprecedented opportunity for discovering new lead compounds targeting thermoTRP channels with improved pharmacological activity and therapeutic index.

AREAS COVERED

This article reviews the progress on thermoTRP channel pharmacology because of (i) advances in solving their atomic structure using cryo-electron microscopy and, (ii) progress on computational techniques including homology modeling, molecular docking, virtual screening, molecular dynamics, ADME/Tox and artificial intelligence. Together, they have increased the number of lead compounds with clinical potential to treat a variety of pathologies. We used original and review articles from Pubmed (1997-2020), as well as the clinicaltrials.gov database, containing the terms thermoTRP, artificial intelligence, docking, and molecular dynamics.

EXPERT OPINION

The atomic structure of thermoTRP channels along with computational methods constitute a realistic first line strategy for designing drug candidates with improved pharmacology and clinical translation. In silico approaches can also help predict potential side-effects that can limit clinical development of drug candidates. Together, they should provide drug candidates with upgraded therapeutic properties.

ExPortal and the LiaFSR Regulatory System Coordinate the Response to Cell Membrane Stress in Streptococcus pyogenes

Bacterial two-component systems sense and induce transcriptional changes in response to environmental stressors, including antimicrobials and human antimicrobial peptides. Since the stresses imposed by the host’s defensive responses may act as markers of specific temporal stages of disease progression or host compartments, pathogens often coordinately regulate stress response programs with virulence factor expression. The mechanism by which bacteria recognize these stresses and subsequently induce transcriptional responses remains not well understood. In this study, we showed that LiaFSR senses cell envelope stress through colocalization of LiaF and LiaS with the group A Streptococcus (GAS) ExPortal and is activated in direct response to ExPortal disruption by antimicrobials or human antimicrobial peptides. Our studies shed new light on the sensing of cell envelope stress in Gram-positive bacteria and may contribute to the development of therapies targeting these processes.

LiaFSR is a gene regulatory system important for response to cell membrane stress in Gram-positive bacteria but is minimally studied in the important human pathogen group A Streptococcus (GAS). Using immunofluorescence and immunogold electron microscopy, we discovered that LiaF (a membrane-bound repressor protein) and LiaS (a sensor kinase) reside within the GAS membrane microdomain (ExPortal). Cell envelope stress induced by antimicrobials resulted in ExPortal disruption and activation of the LiaFSR system. The only human antimicrobial peptide whose presence resulted in ExPortal disruption and LiaFSR activation was the alpha-defensin human neutrophil peptide 1 (hNP-1). Elimination of membrane cardiolipin through targeted gene deletion resulted in loss of LiaS colocalization with the GAS ExPortal and activation of LiaFSR, whereas LiaF membrane localization was unaffected. Isogenic mutants lacking either LiaF or LiaS revealed a critical role of LiaF in ExPortal integrity. Thus, LiaF and LiaS colocalize with the GAS ExPortal by distinct mechanisms, further supporting codependence. These are the first data identifying a multicomponent signal system within the ExPortal, thereby providing new insight into bacterial intramembrane signaling in GAS that may serve as a paradigm for Gram-positive bacteria.

Preparation of Lipid Nanodiscs with Lipid Mixtures

Lipid nanodiscs provide a native-like lipid environment for membrane proteins, and they have become a valuable platform for the study of membrane biophysics. A range of biophysical and biochemical analyses are enabled when membrane proteins are captured in lipid nanodiscs. Two parameters that can be controlled when capturing membrane proteins in lipid nanodiscs are the radius, and hence the surface area of the lipid surface, and the composition of the lipid bilayer. Despite their emergence as a versatile tool, most studies with lipid nanodiscs in the literature have focused on nanodiscs of a single radius with a single lipid. In light of the complexity of biological membranes, it is likely that nanodiscs with multiple membrane components would be more sophisticated models for membrane research. It is possible to prepare nanodiscs with more complex lipid mixtures to probe the effects of lipid composition on several aspects of membrane biochemistry. Detailed protocols are described here for the preparation of nanodiscs with mixtures of phospholipids, incorporation of cholesterol, and incorporation of a spectroscopic lipid probe. These protocols provide starting points for the construction of nanodiscs with more physiological membrane compositions or with useful biophysical probes. © 2019 by John Wiley & Sons, Inc. Basic Protocol 1: Assembly of mixed lipid nanodiscs Basic Protocol 2: Assembly of nanodiscs with cholesterol Basic Protocol 3: Incorporation of laurdan into nanodiscs for membrane fluidity measurements.

Membrane dynamics of γ-secretase with the anterior pharynx-defective 1B subunit

The four-subunit protease complex γ-secretase cleaves many single-pass transmembrane (TM) substrates, including Notch and β-amyloid precursor protein to generate amyloid-β (Aβ), central to Alzheimer's disease. Two of the subunits anterior pharynx-defective 1 (APH-1) and presenilin (PS) exist in two homologous forms APH1-A and APH1-B, and PS1 and PS2. The consequences of these variations are poorly understood and could affect Aβ production and γ-secretase medicine. Here, we developed the first complete structural model of the APH-1B subunit using the published cryo-electron microscopy (cryo-EM) structures of APH1-A (Protein Data Bank: 5FN2, 5A63, and 6IYC). We then performed all-atom molecular dynamics simulations at 303 K in a realistic bilayer system to understand both APH-1B alone and in γ-secretase without and with substrate C83-bound. We show that APH-1B adopts a 7TM topology with a water channel topology similar to APH-1A. We demonstrate direct transport of water through this channel, mainly via Glu84, Arg87, His170, and His196. The apo and holo states closely resemble the experimental cryo-EM structures with APH-1A, however with subtle differences: The substrate-bound APH-1B γ-secretase was quite stable, but some TM helices of PS1 and APH-1B rearranged in the membrane consistent with the disorder seen in the cryo-EM data. This produces different accessibility of water molecules for the catalytic aspartates of PS1, critical for Aβ production. In particular, we find that the typical distance between the catalytic aspartates of PS1 and the C83 cleavage sites are shorter in APH-1B, that is, it represents a more closed state, due to interactions with the C-terminal fragment of PS1. Our structural-dynamic model of APH-1B alone and in γ-secretase suggests generally similar topology but some notable differences in water accessibility which may be relevant to the protein's existence in two forms and their specific function and location.

Regulating Lipid Composition Rationalizes Acyl Tail Saturation Homeostasis in Ectotherms

Cell membranes mainly consist of lipid bilayers with an actively regulated composition. The underlying processes are still poorly understood, in particular, how the hundreds of components are controlled. Cholesterol has been found to correlate with phospholipid saturation for reasons that remain unclear. To better understand the link between cell membrane regulation and chemical composition, we establish a computational framework based on chemical reaction networks, resulting in multiple semigrand canonical ensembles. By running computer simulations, we show that regulating the chemical potential of lipid species is sufficient to reproduce the experimentally observed increase in acyl tail saturation with added cholesterol. Our model proposes a different picture of lipid regulation in which components can be regulated passively instead of actively. In this picture, phospholipid acyl tail composition naturally adapts to added molecules such as cholesterol or proteins. A comparison between regulated membranes with commonly studied ternary model membranes shows stark differences: for instance, correlation lengths and viscosities observed are independent of lipid chemical affinity.

The Fluidity of the Bacterial Outer Membrane Is Species Specific: Bacterial Lifestyles and the Emergence of a Fluid Outer Membrane

The outer membrane (OM) is an essential barrier that guards Gram-negative bacteria from diverse environmental insults. Besides functioning as a chemical gatekeeper, the OM also contributes towards the strength and stiffness of cells and allows them to sustain mechanical stress. Largely influenced by studies of Escherichia coli, the OM is viewed as a rigid barrier where OM proteins and lipopolysaccharides display restricted mobility. Here the discussion is extended to other bacterial species, with a focus on Myxococcus xanthus. In contrast to the rigid OM paradigm, myxobacteria possess a relatively fluid OM. It is concluded that the fluidity of the OM varies across environmental species, which is likely linked to their evolution and adaptation to specific ecological niches. Importantly, a fluid OM can endow bacteria with distinct functions for cell-cell and cell-environment interactions.

Protein-ligand binding with the coarse-grained Martini model

The detailed understanding of the binding of small molecules to proteins is the key for the development of novel drugs or to increase the acceptance of substrates by enzymes. Nowadays, computer-aided design of protein–ligand binding is an important tool to accomplish this task. Current approaches typically rely on high-throughput docking essays or computationally expensive atomistic molecular dynamics simulations. Here, we present an approach to use the recently re-parametrized coarse-grained Martini model to perform unbiased millisecond sampling of protein–ligand interactions of small drug-like molecules. Remarkably, we achieve high accuracy without the need of any a priori knowledge of binding pockets or pathways. Our approach is applied to a range of systems from the well-characterized T4 lysozyme over members of the GPCR family and nuclear receptors to a variety of enzymes. The presented results open the way to high-throughput screening of ligand libraries or protein mutations using the coarse-grained Martini model.

Computer-aided design of protein-ligand binding is important for the development of novel drugs. Here authors present an approach to use the recently re-parametrized coarse-grained Martini model to perform unbiased millisecond sampling of protein-ligand binding interactions of small drug-like molecules.

Interfacial Binding Sites for Cholesterol on Kir, Kv, K2P, and Related Potassium Channels

Inwardly rectifying, voltage-gated, two-pore domain, and related K+ channels are located in eukaryotic membranes rich in cholesterol. Here, molecular docking is used to detect specific binding sites ("hot spots") for cholesterol on K+ channels with characteristics that match those of known cholesterol binding sites. The transmembrane surfaces of all available high-resolution structures for K+ channels were swept for potential binding sites. Cholesterol poses were found to be located largely in hollows between protein ridges. A comparison between cholesterol poses and resolved phospholipids suggests that not all cholesterol molecules binding to the transmembrane surface of a K+ channel will result in displacement of a phospholipid molecule from the surface. Competition between cholesterol binding and binding of anionic phospholipids essential for activity could explain some of the effects of cholesterol on channel function.

The Glycosphingolipid GM3 Modulates Conformational Dynamics of the Glucagon Receptor

The extracellular domain (ECD) of class B1 G-protein-coupled receptors (GPCRs) plays a central role in signal transduction and is uniquely positioned to sense both the extracellular and membrane environments. Although recent studies suggest a role for membrane lipids in the modulation of class A and class F GPCR signaling properties, little is known about the effect of lipids on class B1 receptors. In this study, we employed multiscale molecular dynamics simulations to access the dynamics of the glucagon receptor (GCGR) ECD in the presence of native-like membrane bilayers. Simulations showed that the ECD could move about a hinge region formed by residues Q122–E126 to adopt both closed and open conformations relative to the transmembrane domain. ECD movements were modulated by binding of the glycosphingolipid GM3. These large-scale fluctuations in ECD conformation may affect the ligand binding and receptor activation properties. We also identify a unique phosphatidylinositol (4,5)-bisphosphate (PIP₂) interaction profile near intracellular loop (ICL) 2/TM3 at the G-protein-coupling interface, suggesting a mechanism of engaging G-proteins that may have a distinct dependence on PIP₂ compared with class A GPCRs. Given the structural conservation of class B1 GPCRs, the modulatory effects of GM3 and PIP₂ on GCGR may be conserved across these receptors, offering new insights into potential therapeutic targeting.

Coarse-Grained Molecular Model for the Glycosylphosphatidylinositol Anchor with and without Protein

Glycosylphosphatidylinositol (GPI) anchors are a unique class of complex glycolipids that anchor a great variety of proteins to the extracellular leaflet of plasma membranes of eukaryotic cells. These anchors can exist either with or without an attached protein called GPI-anchored protein (GPI-AP) both in vitro and in vivo. Although GPIs are known to participate in a broad range of cellular functions, it is to a large extent unknown how these are related to GPI structure and composition. Their conformational flexibility and microheterogeneity make it difficult to study them experimentally. Simplified atomistic models are amenable to all-atom computer simulations in small lipid bilayer patches but not suitable for studying their partitioning and trafficking in complex and heterogeneous membranes. Here, we present a coarse-grained model of the GPI anchor constructed with a modified version of the MARTINI force field that is suited for modeling carbohydrates, proteins, and lipids in an aqueous environment using MARTINI's polarizable water. The nonbonded interactions for sugars were reparametrized by calculating their partitioning free energies between polar and apolar phases. In addition, sugar-sugar interactions were optimized by adjusting the second virial coefficients of osmotic pressures for solutions of glucose, sucrose, and trehalose to match with experimental data. With respect to the conformational dynamics of GPI-anchored green fluorescent protein, the accessible time scales are now at least an order of magnitude larger than for the all-atom system. This is particularly important for fine-tuning the mutual interactions of lipids, carbohydrates, and amino acids when comparing to experimental results. We discuss the prospective use of the coarse-grained GPI model for studying protein-sorting and trafficking in membrane models.

Lipid-protein interactions modulate the conformational equilibrium of a potassium channel

Cell membranes actively participate in the regulation of protein structure and function. In this work, we conduct molecular dynamics simulations to investigate how different membrane environments affect protein structure and function in the case of MthK, a potassium channel. We observe different ion permeation rates of MthK in membranes with different properties, and ascribe them to a shift of the conformational equilibrium between two states of the channel that differ according to whether a transmembrane helix has a kink. Further investigations indicate that two key residues in the kink region mediate a crosstalk between two gates at the selectivity filter and the central cavity, respectively. Opening of one gate eventually leads to closure of the other. Our simulations provide an atomistic model of how lipid-protein interactions affect the conformational equilibrium of a membrane protein. The gating mechanism revealed for MthK may also apply to other potassium channels.

Calcium-Lipid Interactions Observed with Isotope-Edited Infrared Spectroscopy

Calcium ions bind to lipid membranes containing anionic lipids; however, characterizing the specific ion-lipid interactions in multicomponent membranes has remained challenging because it requires nonperturbative lipid-specific probes. Here, using a combination of isotope-edited infrared spectroscopy and molecular dynamics simulations, we characterize the effects of a physiologically relevant (2 mM) Ca²⁺ concentration on zwitterionic phosphatidylcholine and anionic phosphatidylserine lipids in mixed lipid membranes. We show that Ca²⁺ alters hydrogen bonding between water and lipid headgroups by forming a coordination complex involving the lipid headgroups and water. These interactions distort interfacial water orientations and prevent hydrogen bonding with lipid ester carbonyls. We demonstrate, experimentally, that these effects are more pronounced for the anionic phosphatidylserine lipids than for zwitterionic phosphatidylcholine lipids in the same membrane.

Rotational Diffusion of Membrane Proteins in Crowded Membranes

Membrane proteins travel along cellular membranes and reorient themselves to form functional oligomers and protein-lipid complexes. Following the Saffman-Delbrück model, protein radius sets the rate of this diffusive motion. However, it is unclear how this model, derived for ideal and dilute membranes, performs under crowded conditions of cellular membranes. Here, we study the rotational motion of membrane proteins using molecular dynamics simulations of coarse-grained membranes and 2-dimensional Lennard-Jones fluids with varying levels of crowding. We find that the Saffman-Delbrück model captures the size-dependency of rotational diffusion under dilute conditions where protein-protein interactions are negligible, whereas stronger scaling laws arise under crowding. Together with our recent work on lateral diffusion, our results reshape the description of protein dynamics in native membrane environments: The translational and rotational motions of proteins with small transmembrane domains are rapid, whereas larger proteins or protein complexes display substantially slower dynamics.

Membrane mediated toppling mechanism of the folate energy coupling factor transporter

Energy coupling factor (ECF) transporters are responsible for the uptake of micronutrients in bacteria and archaea. They consist of an integral membrane unit, the S-component, and a tripartite ECF module. It has been proposed that the S-component mediates the substrate transport by toppling over in the membrane when docking onto an ECF module. Here, we present multi-scale molecular dynamics simulations and in vitro experiments to study the molecular toppling mechanism of the S-component of a folate-specific ECF transporter. Simulations reveal a strong bending of the membrane around the ECF module that provides a driving force for toppling of the S-component. The stability of the toppled state depends on the presence of non-bilayer forming lipids, as confirmed by folate transport activity measurements. Together, our data provide evidence for a lipid-dependent toppling-based mechanism for the folate-specific ECF transporter, a mechanism that might apply to other ECF transporters.

Measuring Remodeling of the Lipid Environment Surrounding Membrane Proteins with Lipid Exchange and Native Mass Spectrometry

Due to their crucial biochemical roles, membrane proteins are important drug targets. Although it is clear that lipids can influence membrane protein function, the chemistry of lipid binding remains difficult to study because protein-lipid interactions are polydisperse, competitive, and transient. Furthermore, detergents, which are often used to solubilize membrane proteins in micelles, may disrupt lipid interactions that occur in bilayers. Here, we present two new approaches to quantify protein-lipid interactions in bilayers and understand how membrane proteins remodel their surrounding lipid environment. First, we used mass spectrometry (MS) to measure the exchange of lipids between lipoprotein nanodiscs with and without an embedded membrane protein. Shifts in the lipid distribution toward the membrane protein nanodiscs revealed lipid binding, and titrations allowed measurement of the optimal lipid composition for the membrane protein. Second, we used native or nondenaturing MS to ionize membrane protein nanodiscs with heterogeneous lipids. Ejecting the membrane protein complex with bound lipids in the mass spectrometer revealed enrichment of specific lipids around the membrane protein. Both new approaches showed that the E. coli ammonium transporter AmtB prefers phosphatidylglycerol lipids overall but has a minor affinity for phosphatidylcholine lipids.

Understanding and modelling the interactions of peptides with membranes: from partitioning to self-assembly

Atomic detail simulations are starting to reveal how flexible polypeptides interact with fluid lipid bilayers. These insights are transforming our understanding of one of the fundamental processes in biology: membrane protein folding and assembly. Advanced molecular dynamics (MD) simulation techniques enable accurate prediction of protein structure, folding pathways and assembly in microsecond-timescales. Such simulations show how membrane-active peptides self-assemble in cell membranes, revealing their binding, folding, insertion, and aggregation, while at the same time providing atomic resolution details of peptide-lipid interactions. Essential to the impact of simulations are experimental approaches that enable calibration and validation of the computational models and techniques. In this review, we summarize the current development of applying unbiased atomic detail MD simulations and the relation to experimental techniques, to study peptide folding and provide our perspective of the field.

Lipid-Protein Interactions Are a Unique Property and Defining Feature of G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are membrane-bound proteins that depend on their lipid environment to carry out their physiological function. Combined efforts from many theoretical and experimental studies on the lipid-protein interaction profile of several GPCRs hint at an intricate relationship of these receptors with their surrounding membrane environment, with several lipids emerging as particularly important. Using coarse-grained molecular dynamics simulations, we explore the lipid-protein interaction profiles of 28 different GPCRs, spanning different levels of classification and conformational states and totaling to 1 ms of simulation time. We find a close relationship with lipids for all GPCRs simulated, in particular, cholesterol and phosphatidylinositol phosphate (PIP) lipids, but the number, location, and estimated strength of these interactions is dependent on the specific GPCR as well as its conformational state. Although both cholesterol and PIP lipids bind specifically to GPCRs, they utilize distinct mechanisms. Interactions with PIP lipids are mediated by charge-charge interactions with intracellular loop residues and stabilized by one or both of the transmembrane helices linked by the loop. Interactions with cholesterol, on the other hand, are mediated by a hydrophobic environment, usually made up of residues from more than one helix, capable of accommodating its ring structure and stabilized by interactions with aromatic and charged/polar residues. Cholesterol binding to GPCRs occurs in a small number of sites, some of which (like the binding site on the extracellular side of transmembrane 6/7) are shared among many class A GPCRs. Combined with a thorough investigation of the local membrane structure, our results provide a detailed picture of GPCR-lipid interactions. Additionally, we provide an accompanying website to interactively explore the lipid-protein interaction profile of all GPCRs simulated to facilitate analysis and comparison of our data.

Capturing Choline-Aromatics Cation-π Interactions in the MARTINI Force Field

Cation−π interactions play an important role in biomolecular recognition, including interactions between membrane phosphatidylcholine lipids and aromatic amino acids of peripheral proteins. While molecular mechanics coarse grain (CG) force fields are particularly well suited to simulate membrane proteins in general, they are not parameterized to explicitly reproduce cation−π interactions. We here propose a modification of the polarizable MARTINI coarse grain (CG) model enabling it to model membrane binding events of peripheral proteins whose aromatic amino acid interactions with choline headgroups are crucial for their membrane binding. For this purpose, we first collected and curated a dataset of eight peripheral proteins from different families. We find that the MARTINI CG model expectedly underestimates aromatics–choline interactions and is unable to reproduce membrane binding of the peripheral proteins in our dataset. Adjustments of the relevant interactions in the polarizable MARTINI force field yield significant improvements in the observed binding events. The orientation of each membrane-bound protein is comparable to reference data from all-atom simulations and experimental binding data. We also use negative controls to ensure that choline–aromatics interactions are not overestimated. We finally check that membrane properties, transmembrane proteins, and membrane translocation potential of mean force (PMF) of aromatic amino acid side-chain analogues are not affected by the new parameter set. This new version “MARTINI 2.3P” is a significant improvement over its predecessors and is suitable for modeling membrane proteins including peripheral membrane binding of peptides and proteins.

Coarse-grained MD simulations reveal beta-amyloid fibrils of various sizes bind to interfacial liquid-ordered and liquid-disordered regions in phase separated lipid rafts with diverse membrane-bound conformational states

The membrane binding behaviors of beta-amyloid fibrils, dimers to pentamers, from solution to lipid raft surfaces, were investigated using coarse-grained (CG) MD simulations. Our CG rafts contain phospholipid, cholesterol (with or without tail- or headgroup modifications), and with or without asymmetrically distributed monosialotetrahexosylganglioside (GM1). All rafts exhibited liquid-ordered (Lo), liquid-disordered (Ld), and interfacial Lo/Ld (Lod) domains, with domain sizes depending on cholesterol structure. For rafts without GM1, all fibrils bound to the Lod domains. Specifically, dimer fibrils bound exclusively via the C-terminal, while larger fibrils could bind via other protein regions. Interestingly, a membrane-inserted state was detected for a trimer fibril in a raft with tail-group modified cholesterol. For rafts containing GM1, fibrils bound either to the GM1-clusters, with numerous membrane-bound conformations, or to the non-GM1-containing-Lod domains via the C-terminal. Our results indicate beta-amyloid fibrils bind to Lod domains or GM1, with diversified membrane-bound conformations, in structurally heterogeneous lipid membranes.

Distinct structural modulation of photosystem I and lipid environment stabilizes its tetrameric assembly

Photosystem I (PSI) is able to form different oligomeric states across various species. To reveal the structural basis for PSI dimerization and tetramerization, we structurally investigated PSI from the cyanobacterium Anabaena. This revealed a disrupted trimerization domain due to lack of the terminal residues of PsaL in the lumen, which resulted in PSI dimers with loose connections between monomers and weaker energy-coupled chlorophylls than in the trimer. At the dimer surface, specific phospholipids, cofactors and interactions in combination facilitated recruitment of another dimer to form a tetramer. Taken together, the relaxed luminal connections and lipid specificity at the dimer interface account for membrane curvature. PSI tetramer assembly appears to increase the surface area of the thylakoid membrane, which would contribute to PSI crowding.

Non-oxidative band-3 clustering agents cause the externalization of phosphatidylserine on erythrocyte surfaces by a calcium-independent mechanism

Aging of red blood cells (RBCs) is associated with alteration in a wide range of RBC features, occurring each on its own timescale. A number of these changes are interrelated and initiate a cascade of biochemical and structural transformations, including band-3 clustering and phosphatidylserine (PS) externalization. Using specific band-3 clustering agents (acridine orange (AO) and ZnCl₂), we examined whether treatment of RBCs with these agents may affects PS externalization and whether this process is Ca²⁺-dependent. RBCs were isolated from the blood of eight healthy donors upon obtaining their informed consent. The suspension was supplemented with increasing concentrations of AO or ZnCl₂ (from 0.5 to 2.0 mM) and incubated at 25 °C for 60 min. To detect PS at the RBC surface, we used allophycocyanin-conjugated recombinant human Annexin V. We demonstrated, that treatment of RBCs with both clustering agents caused an elevation in the percent of cells positively labeled by Annexin-V (RBC^PS), and that this value was not dependent on the presence of calcium in the buffer: RBCs treated with AO in the presence of either EDTA, EGTA or calcium exhibited similar percentage of RBC^PS. Moreover, the active influx of Zn²⁺ into RBCs induced by their co-incubation with both ZnCl₂ and A23187 did not increase the percent of RBC^PS as compared to RBCs incubated with ZnCl₂ alone. Taken together, these results demonstrate that the band-3 clustering agents (AO or ZnCl₂) induce PS externalization in a Ca²⁺ independent manner, and we hereby suggest a possible scenario for this phenomenon.

The energetics of protein-lipid interactions as viewed by molecular simulations

Membranes are formed from a bilayer containing diverse lipid species with which membrane proteins interact. Integral, membrane proteins are embedded in this bilayer, where they interact with lipids from their surroundings, whilst peripheral membrane proteins bind to lipids at the surface of membranes. Lipid interactions can influence the function of membrane proteins, either directly or allosterically. Both experimental (structural) and computational approaches can reveal lipid binding sites on membrane proteins. It is, therefore, important to understand the free energies of these interactions. This affords a more complete view of the engagement of a particular protein with the biological membrane surrounding it. Here, we describe many computational approaches currently in use for this purpose, including recent advances using both free energy and unbiased simulation methods. In particular, we focus on interactions of integral membrane proteins with cholesterol, and with anionic lipids such as phosphatidylinositol 4,5-bis-phosphate and cardiolipin. Peripheral membrane proteins are exemplified via interactions of PH domains with phosphoinositide-containing membranes. We summarise the current state of the field and provide an outlook on likely future directions of investigation.

Prebiotic Peptides: Molecular Hubs in the Origin of Life

The fundamental roles that peptides and proteins play in today's biology makes it almost indisputable that peptides were key players in the origin of life. Insofar as it is appropriate to extrapolate back from extant biology to the prebiotic world, one must acknowledge the critical importance that interconnected molecular networks, likely with peptides as key components, would have played in life's origin. In this review, we summarize chemical processes involving peptides that could have contributed to early chemical evolution, with an emphasis on molecular interactions between peptides and other classes of organic molecules. We first summarize mechanisms by which amino acids and similar building blocks could have been produced and elaborated into proto-peptides. Next, non-covalent interactions of peptides with other peptides as well as with nucleic acids, lipids, carbohydrates, metal ions, and aromatic molecules are discussed in relation to the possible roles of such interactions in chemical evolution of structure and function. Finally, we describe research involving structural alternatives to peptides and covalent adducts between amino acids/peptides and other classes of molecules. We propose that ample future breakthroughs in origin-of-life chemistry will stem from investigations of interconnected chemical systems in which synergistic interactions between different classes of molecules emerge.

Microsecond-timescale simulations suggest 5-HT-mediated preactivation of the 5-HT3A serotonin receptor

Aided by efforts to improve their speed and efficiency, molecular dynamics (MD) simulations provide an increasingly powerful tool to study the structure-function relationship of pentameric ligand-gated ion channels (pLGICs). However, accurate reporting of the channel state and observation of allosteric regulation by agonist binding with MD remains difficult due to the timescales necessary to equilibrate pLGICs from their artificial and crystalized conformation to a more native, membrane-bound conformation in silico. Here, we perform multiple all-atom MD simulations of the homomeric 5-hydroxytryptamine 3A (5-HT3A) serotonin receptor for 15 to 20 μs to demonstrate that such timescales are critical to observe the equilibration of a pLGIC from its crystalized conformation to a membrane-bound conformation. These timescales, which are an order of magnitude longer than any previous simulation of 5-HT3A, allow us to observe the dynamic binding and unbinding of 5-hydroxytryptamine (5-HT) (i.e., serotonin) to the binding pocket located on the extracellular domain (ECD) and allosteric regulation of the transmembrane domain (TMD) from synergistic 5-HT binding. While these timescales are not long enough to observe complete activation of 5-HT3A, the allosteric regulation of ion gating elements by 5-HT binding is indicative of a preactive state, which provides insight into molecular mechanisms that regulate channel activation from a resting state. This mechanistic insight, enabled by microsecond-timescale MD simulations, will allow a careful examination of the regulation of pLGICs at a molecular level, expanding our understanding of their function and elucidating key structural motifs that can be targeted for therapeutic regulation.

The Fats of Life: Using Computational Chemistry to Characterise the Eukaryotic Cell Membrane

Our current knowledge of the structural dynamics and complexity of lipid bilayers is still developing. Computational techniques, especially molecular dynamics simulations, have increased our understanding significantly as they allow us to model functions that cannot currently be experimentally resolved. Here we review available computational tools and techniques, the role of the major lipid species, insights gained into lipid bilayer structure and function from molecular dynamics simulations, and recent progress towards the computational modelling of the physiological complexity of eukaryotic lipid bilayers.

More than one way to bind to cholesterol: atypical variants of membrane-binding domain of perfringolysin O selected by ribosome display

Herein, we report a high-throughput approach for the selection of peripheral protein domains that bind specifically to cholesterol in lipid membranes. We discovered variants of perfringolysin O, with non-conserved amino acid substitutions at regions crucial for cholesterol recognition, demonstrating an unprecedented amino acid sequence variability with binding ability for cholesterol. The developed approach provides an effective platform for a comprehensive study of protein lipid interactions.

By using developed ribosomal display, we discovered variants of perfringolysin O, a pore forming toxin from bacteria Clostridium perfringens, with non-conserved amino acid substitutions at regions crucial for cholesterol recognition.

Integrated regulation of the phosphatidylinositol cycle and phosphoinositide-driven lipid transport at ER-PM contact sites

Among the structural phospholipids that form the bulk of eukaryotic cell membranes, phosphatidylinositol (PtdIns) is unique in that it also serves as the common precursor for low-abundance regulatory lipids, collectively referred to as polyphosphoinositides (PPIn). The metabolic turnover of PPIn species has received immense attention because of the essential functions of these lipids as universal regulators of membrane biology and their dysregulation in numerous human pathologies. The diverse functions of PPIn lipids occur, in part, by orchestrating the spatial organization and conformational dynamics of peripheral or integral membrane proteins within defined subcellular compartments. The emerging role of stable contact sites between adjacent membranes as specialized platforms for the coordinate control of ion exchange, cytoskeletal dynamics, and lipid transport has also revealed important new roles for PPIn species. In this review, we highlight the importance of membrane contact sites formed between the endoplasmic reticulum (ER) and plasma membrane (PM) for the integrated regulation of PPIn metabolism within the PM. Special emphasis will be placed on non-vesicular lipid transport during control of the PtdIns biosynthetic cycle as well as toward balancing the turnover of the signaling PPIn species that define PM identity.

Role of cholesterol-mediated effects in GPCR heterodimers

G protein-coupled receptors (GPCRs) are transmembrane receptors that mediate a large number of cellular responses. The organization of GPCRs into dimers and higher-order oligomers is known to allow a larger repertoire of downstream signaling events. In this context, a crosstalk between the adenosine and dopamine receptors has been reported, indicating the presence of heterodimers that are functionally relevant. In this paper, we explored the effect of membrane cholesterol on the adenosine_2A (A_2A) and dopamine D₃ (D₃) receptors using coarse-grain molecular dynamics simulations. We analyzed cholesterol interaction sites on the A_2A receptor and were able to reproduce the sites indicated by crystallography and previous atomistic simulations. We predict novel cholesterol interaction sites on the D₃ receptor that could be important in the reported cholesterol sensitivity in receptor function. Further, we analyzed the formation of heterodimers between the two receptors. Our results suggest that membrane cholesterol modulates the relative population of several co-existing heterodimer conformations. Both direct receptor-cholesterol interaction and indirect membrane effects contribute toward the modulation of heterodimer conformations. These results constitute one of the first examples of modulation of GPCR hetero-dimerization by membrane cholesterol, and could prove to be useful in designing better therapeutic strategies.

Structure and Dynamics of the Acyl Chains in the Membrane Trafficking and Enzymatic Processing of Lipids

The regulatory chemical mechanisms of lipid trafficking and degradation are involved in many pathophysiological processes, being implicated in severe pain, inflammation, and cancer. In addition, the processing of lipids is also relevant for industrial and environmental applications. However, there is poor understanding of the chemical features that control lipid membrane trafficking and allow lipid-degrading enzymes to efficiently select and hydrolyze specific fatty acids from a complex cellular milieu of bioactive lipids. This is particularly true for lipid acyl chains, which have diverse structures that can critically affect the many complex reactions needed to elongate, desaturate, or transport fatty acids. Building upon our own contributions in this field, we will discuss how molecular simulations, integrated with experimental evidence, have revealed that the structure and dynamics of the lipid tail are actively involved in modulating membrane trafficking at cellular organelles, and enzymatic reactions at cell membranes. Further evidence comes from recent crystal structures of lipid receptors and remodeling enzymes. Taken together, these recent works have identified those structural features of the lipid acyl chain that are crucial for the regioselectivity and stereospecificity of essential desaturation reactions. In this context, we will first illustrate how atomistic and coarse-grained simulations have elucidated the structure-function relationships between the chemical composition of the lipid's acyl chains and the molecular properties of lipid bilayers. Particular emphasis will be given to the prominent chemical role of the number of double carbon-carbon bonds along the lipid acyl chain, that is, discriminating between saturated, monounsaturated, and polyunsaturated lipids. Different levels of saturation in fatty acid molecules dramatically influence the biophysical properties of lipid assemblies and their interaction with proteins. We will then discuss the processing of lipids by membrane-bound enzymes. Our focus will be on lipids such as anandamide and 2-arachidonoylglycerol. These are the main molecules that act as neurotransmitters in the endocannabinoid system. Specifically, recent findings indicate a crucial interplay between the level of saturation of the lipid tail, its energetically and sterically favored conformations, and the hydrophobic accessory cavities in lipid-degrading enzymes, which help form catalytically active conformations of the selected substrate. This Account will emphasize how the specific chemical structure of acyl chains affects the molecular mechanisms for modulating membrane trafficking and selective hydrolysis. The results examined here show that, by using molecular simulations to investigate lipid plasticity and substrate flexibility, researchers can enrich their interpretation of experimental results about the structure-function relationships of lipids. This could positively impact chemical and biological studies in the field and ultimately support protein engineering studies and structure-based drug discovery to target lipid-processing enzymes.

Atoms to Phenotypes: Molecular Design Principles of Cellular Energy Metabolism

We report a 100-million atom-scale model of an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium, that reveals the cascade of energy conversion steps culminating in the generation of ATP from sunlight. Molecular dynamics simulations of this vesicle elucidate how the integral membrane complexes influence local curvature to tune photoexcitation of pigments. Brownian dynamics of small molecules within the chromatophore probe the mechanisms of directional charge transport under various pH and salinity conditions. Reproducing phenotypic properties from atomistic details, a kinetic model evinces that low-light adaptations of the bacterium emerge as a spontaneous outcome of optimizing the balance between the chromatophore's structural integrity and robust energy conversion. Parallels are drawn with the more universal mitochondrial bioenergetic machinery, from whence molecular-scale insights into the mechanism of cellular aging are inferred. Together, our integrative method and spectroscopic experiments pave the way to first-principles modeling of whole living cells.

The Effect of Transmembrane Protein Shape on Surrounding Lipid Domain Formation by Wetting

Signal transduction through cellular membranes requires the highly specific and coordinated work of specialized proteins. Proper functioning of these proteins is provided by an interplay between them and the lipid environment. Liquid-ordered lipid domains are believed to be important players here, however, it is still unclear whether conditions for a phase separation required for lipid domain formation exist in cellular membranes. Moreover, membrane leaflets are compositionally asymmetric, that could be an obstacle for the formation of symmetric domains spanning the lipid bilayer. We theoretically show that the presence of protein in the membrane leads to the formation of a stable liquid-ordered lipid phase around it by the mechanism of protein wetting by lipids, even in the absence of conditions necessary for the global phase separation in the membrane. Moreover, we show that protein shape plays a crucial role in this process, and protein conformational rearrangement can lead to changes in the size and characteristics of surrounding lipid domains.

Molecular model of the ferroportin intracellular gate and implications for the human iron transport cycle and hemochromatosis type 4A

Ferroportin 1 (FPN1) is a major facilitator superfamily transporter that is essential for proper maintenance of human iron homeostasis at the systemic and cellular level. FPN1 dysfunction leads to the progressive accumulation of iron in reticuloendothelial cells, causing hemochromatosis type 4A (or ferroportin disease), an autosomal dominant disorder that displays large phenotypic heterogeneity. Although crystal structures have unveiled the outward- and inward-facing conformations of the bacterial homolog Bdellovibrio bacteriovorus Fpn (or Bd2019) and calcium has recently been identified as an essential cofactor, our molecular understanding of the iron transport mechanism remains incomplete. Here, we used a combination of molecular modeling, molecular dynamics simulations, and Ala site-directed mutagenesis, followed by complementary in vitro functional analyses, to explore the structural architecture of the human FPN1 intracellular gate. We reveal an interdomain network that involves 5 key amino acids and is likely very important for stability of the iron exporter facing the extracellular milieu. We also identify inter- and intradomain interactions that rely on the 2 Asp84 and Asn174 critical residues and do not exist in the bacterial homolog. These interactions are thought to play an important role in the modulation of conformational changes during the transport cycle. We interpret these results in the context of hemochromatosis type 4A, reinforcing the idea that different categories of loss-of-function mutations exist. Our findings provide an unprecedented view of the human FPN1 outward-facing structure and the particular function of the so-called "gating residues" in the mechanism of iron export.-Guellec, J., Elbahnsi, A., Le Tertre, M., Uguen, K., Gourlaouen, I., Férec, C., Ka, C., Callebaut, I., Le Gac, G. Molecular model of the ferroportin intracellular gate and implications for the human iron transport cycle and hemochromatosis type 4A.

Interfacial Binding Sites for Cholesterol on TRP Ion Channels

Transient receptor potential (TRP) channels are members of a large family of ion channels located in membranes rich in cholesterol, some of whose functions are affected by the cholesterol content of the membrane. Here, cholesterol binding to TRPs is studied using a docking procedure that allows the transmembrane surface of a TRP to be swept rapidly for potential binding sites at the interfaces on the two sides of the membrane. Cholesterol docking poses determined in this way match 89% of the cholesterol hemisuccinate molecules in published TRP structures when cholesterol hemisuccinate molecules unlikely to represent typical bound cholesterols are excluded. TRPs are tetrameric, with large clefts at the interfaces between subunits; cholesterol poses are located in hollows, largely within these clefts. Comparison of cholesterol poses with phospholipid binding sites suggests that binding of cholesterol to a TRP need not result in displacement of phospholipid molecules from the TRP surface.

Serine Phosphorylation of L-Selectin Regulates ERM Binding, Clustering, and Monocyte Protrusion in Transendothelial Migration

The migration of circulating leukocytes toward damaged tissue is absolutely fundamental to the inflammatory response, and transendothelial migration (TEM) describes the first cellular barrier that is breached in this process. Human CD14⁺ inflammatory monocytes express L-selectin, bestowing a non-canonical role in invasion during TEM. In vivo evidence supports a role for L-selectin in regulating TEM and chemotaxis, but the intracellular mechanism is poorly understood. The ezrin-radixin-moesin (ERM) proteins anchor transmembrane proteins to the cortical actin-based cytoskeleton and additionally act as signaling adaptors. During TEM, the L-selectin tail within transmigrating pseudopods interacts first with ezrin to transduce signals for protrusion, followed by moesin to drive ectodomain shedding of L-selectin to limit protrusion. Collectively, interaction of L-selectin with ezrin and moesin fine-tunes monocyte protrusive behavior in TEM. Using FLIM/FRET approaches, we show that ERM binding is absolutely required for outside-in L-selectin clustering. The cytoplasmic tail of human L-selectin contains two serine (S) residues at positions 364 and 367, and here we show that they play divergent roles in regulating ERM binding. Phospho-S364 blocks direct interaction with ERM, whereas molecular modeling suggests phospho-S367 likely drives desorption of the L-selectin tail from the inner leaflet of the plasma membrane to potentiate ERM binding. Serine-to-alanine mutagenesis of S367, but not S364, significantly reduced monocyte protrusive behavior in TEM under flow conditions. Our data propose a model whereby L-selectin tail desorption from the inner leaflet of the plasma membrane and ERM binding are two separable steps that collectively regulate protrusive behavior in TEM.

Membrane Remodeling by the Lytic Fragment of SticholysinII: Implications for the Toroidal Pore Model

Sticholysins are pore-forming toxins of biomedical interest and represent a prototype of proteins acting through the formation of protein-lipid or toroidal pores. Peptides spanning the N-terminus of sticholysins can mimic their permeabilizing activity and, together with the full-length toxins, have been used as a tool to understand the mechanism of pore formation in membranes. However, the lytic mechanism of these peptides and the lipid shape modulating their activity are not completely clear. In this article, we combine molecular dynamics simulations and experimental biophysical tools to dissect different aspects of the pore-forming mechanism of StII_1-30, a peptide derived from the N-terminus of sticholysin II (StII). With this combined approach, membrane curvature induction and flip-flop movement of the lipids were identified as two important membrane remodeling steps mediated by StII_1-30. Pore formation by this peptide was enhanced by the presence of the negatively curved lipid phosphatidylethanolamine in membranes. This lipid emerged not only as a facilitator of membrane interactions but also as a structural element of the StII_1-30 pore that is recruited to the ring upon its assembly. Collectively, these, to our knowledge, new findings support a toroidal model for the architecture of the pore formed by StII_1-30 and provide new molecular insight into the role of phosphatidylethanolamine as a membrane component that can easily integrate into the ring of toroidal pores, thus probably aiding in their stabilization. This study contributes to a better understanding of the molecular mechanism underlying the permeabilizing activity of StII_1-30 and peptides or proteins acting via a toroidal pore mechanism and offers an informative framework for the optimization of the biomedical application of this and similar molecules.

Local Enrichment of Unsaturated Chains around the A2A Adenosine Receptor

Two 15 μs all-atom simulations of the A_2A adenosine receptor were obtained in a ternary mixture of cholesterol, saturated phosphatidylcholine lipids, and unsaturated phosphatidylcholine lipids. An analysis of local lipid solvation is reported on the basis of a Voronoi tessellation of the upper and lower leaflets, identifying first and second solvation shells. The local environments of both the inactive state and the partially active state of the receptor are significantly enriched with unsaturated chains but depleted of cholesterol and saturated chains, relative to the bulk membrane composition. In spite of the local depletion of cholesterol, the partially active receptor binds cholesterol at three locations during the entire simulation trajectory. These long-lived interactions represent the extreme of a very broad distribution of first-solvation shell lipid lifetimes, confounding sharp distinctions between lipid interactions. The broad distributions of lifetimes also make equilibrating the local lipid environment difficult, necessitating long simulation times.

Insights into Membrane Protein-Lipid Interactions from Free Energy Calculations

Integral membrane proteins are regulated by specific interactions with lipids from the surrounding bilayer. The structures of protein–lipid complexes can be determined through a combination of experimental and computational approaches, but the energetic basis of these interactions is difficult to resolve. Molecular dynamics simulations provide the primary computational technique to estimate the free energies of these interactions. We demonstrate that the energetics of protein–lipid interactions may be reliably and reproducibly calculated using three simulation-based approaches: potential of mean force calculations, alchemical free energy perturbation, and well-tempered metadynamics. We employ these techniques within the framework of a coarse-grained force field and apply them to both bacterial and mammalian membrane protein–lipid systems. We demonstrate good agreement between the different techniques, providing a robust framework for their automated implementation within a pipeline for annotation of newly determined membrane protein structures.

Fat SIRAH: Coarse-Grained Phospholipids To Explore Membrane-Protein Dynamics

The capability to handle highly heterogeneous molecular assemblies in a consistent manner is among the greatest challenges faced when deriving simulation parameters. This is particularly the case for coarse-grained (CG) simulations in which chemical functional groups are lumped into effective interaction centers for which transferability between different chemical environments is not guaranteed. Here, we introduce the parametrization of a set of CG phospholipids compatible with the latest version of the SIRAH force field for proteins. The newly introduced lipid species include different acylic chain lengths and partial unsaturation, as well as polar and acidic head groups that show a very good reproduction of structural membrane determinants, such as areas per lipid, thickness, order parameter, etc., and their dependence with temperature. Simulation of membrane proteins showed unprecedented accuracy in the unbiased description of the thickness-dependent membrane-protein orientation in systems where this information is experimentally available (namely, the SarcoEndoplasmic Reticulum Calcium-SERCA-pump and its regulator Phospholamban). The interactions that lead to this faithful reproduction can be traced down to the single amino acid-lipid interaction level and show full agreement with biochemical data present in the literature. Finally, the present parametrization is implemented in the GROMACS and AMBER simulation packages facilitating its use by a wide portion of the biocomputing community.

A Molecular Perspective on Mitochondrial Membrane Fusion: From the Key Players to Oligomerization and Tethering of Mitofusin

Mitochondria are dynamic organelles characterized by an ultrastructural organization which is essential in maintaining their quality control and ensuring functional efficiency. The complex mitochondrial network is the result of the two ongoing forces of fusion and fission of inner and outer membranes. Understanding the functional details of mitochondrial dynamics is physiologically relevant as perturbations of this delicate equilibrium have critical consequences and involved in several neurological disorders. Molecular actors involved in this process are large GTPases from the dynamin-related protein family. They catalyze nucleotide-dependent membrane remodeling and are widely conserved from bacteria to higher eukaryotes. Although structural characterization of different family members has contributed in understanding molecular mechanisms of mitochondrial dynamics in more detail, the complete structure of some members as well as the precise assembly of functional oligomers remains largely unknown. As increasing structural data become available, the domain modularity across the dynamin superfamily emerged as a foundation for transfering the knowledge towards less characterized members. In this review, we will first provide an overview of the main actors involved in mitochondrial dynamics. We then discuss recent example of computational methodologies for the study of mitofusin oligomers, and present how the usage of integrative modeling in conjunction with biochemical data can be an asset in progressing the still challenging field of membrane dynamics.

Life During Wartime: A Personal Recollection of the Circa 1990 Prestegard Lab and Its Contributions to Membrane Biophysics

A subjective account is presented of challenges and excitement of being a postdoctoral trainee in the lab of James H. Prestegard at Yale University in New Haven, Connecticut from 1989 to 1991. This includes accounts of the early development of bicelles and of oriented sample NMR results that contributed to our modern understanding of the properties of the water-lipid interface of disordered phase biological membranes.

Molecular Simulations of Intact Anion Exchanger 1 Reveal Specific Domain and Lipid Interactions

Anion exchanger 1 (AE1) is responsible for the exchange of bicarbonate and chloride across the erythrocyte plasma membrane. Human AE1 consists of a cytoplasmic and a membrane domain joined by a 33-residue flexible linker. Crystal structures of the individual domains have been determined, but the intact AE1 structure remains elusive. In this study, we use molecular dynamics simulations and modeling to build intact AE1 structures in a complex lipid bilayer that resembles the native erythrocyte plasma membrane. AE1 models were evaluated using available experimental data to provide an atomistic view of the interaction and dynamics of the cytoplasmic domain, the membrane domain, and the connecting linker in a complete model of AE1 in a lipid bilayer. Anionic lipids were found to interact strongly with AE1 at specific amino acid residues that are linked to diseases and blood group antigens. Cholesterol was found in the dimeric interface of AE1, suggesting that it may regulate subunit interactions and anion transport.

Protein surface charge effect on 3D domain swapping in cells for c-type cytochromes

Many c-type cytochromes (cyts) can form domain-swapped oligomers. The positively charged Hydrogenobacter thermophilus (HT) cytochrome (cyt) c₅₅₂ forms domain-swapped oligomers during expression in the Escherichia coli (E. coli) expression system, but the factors influencing the oligomerization remain unrevealed. Here, we found that the dimer of the negatively charged Shewanella violacea (SV) cyt c₅ exhibits a domain-swapped structure, in which the N-terminal helix is exchanged between protomers, similar to the structures of the HT cyt c₅₅₂ and Pseudomonas aeruginosa (PA) cyt c₅₅₁ domain-swapped dimers. Positively charged horse cyt c and HT cyt c₅₅₂ domain swapped during expression in E. coli, whereas negatively charged PA cyt c₅₅₁ and SV cyt c₅ did not. Oligomers were formed during expression in E. coli for HT cyt c₅₅₂ attached to either a co- or post-translational signal peptide for transportation through the cytoplasm membrane, but not for PA cyt c₅₅₁ attached to either signal peptide. HT cyt c₅₅₂ formed oligomers in E. coli in the presence and absence of rare codons. More oligomers were obtained from the in vitro folding of horse cyt c and HT cyt c₅₅₂ by the addition of negatively charged liposomes during folding, whereas the amount of oligomers for the in vitro folding of PA cyt c₅₅₁ and SV cyt c₅ did not change significantly by the addition. These results indicate that the protein surface charge affects the oligomerization of c-type cyts in cells; positively charged c-type cyts assemble on a negatively charged membrane, inducing formation of domain-swapped oligomers during folding.