Emerging Diversity in Lipid-Protein Interactions

Lipid distributions and transleaflet cholesterol migration near heterogeneous surfaces in asymmetric bilayers

Specific and nonspecific protein-lipid interactions in cell membranes have important roles in an abundance of biological functions. We have used coarse-grained (CG) molecular dynamics (MD) simulations to assess lipid distributions and cholesterol flipping dynamics around surfaces in a model asymmetric plasma membrane containing one of six structurally distinct entities: aquaporin-1 (AQP1), the bacterial β-barrel outer membrane proteins OmpF and OmpX, the KcsA potassium channel, the WALP23 peptide and a carbon nanotube (CNT). Our findings revealed varied lipid partitioning and cholesterol flipping times around the different solutes and putative cholesterol binding sites in AQP1 and KcsA. The results suggest that protein-lipid interactions can be highly variable, and that surface-dependent lipid profiles are effectively manifested in CG simulations with the Martini force field.

Sphingomyelinase Disables Inactivation in Endogenous PIEZO1 Channels

Endogenous PIEZO1 channels of native endothelium lack the hallmark inactivation often seen when these channels are overexpressed in cell lines. Because prior work showed that the force of shear stress activates sphingomyelinase in endothelium, we considered if sphingomyelinase is relevant to endogenous PIEZO1. Patch clamping was used to quantify PIEZO1-mediated signals in freshly isolated murine endothelium exposed to the mechanical forces caused by shear stress and membrane stretch. Neutral sphingomyelinase inhibitors and genetic disruption of sphingomyelin phosphodiesterase 3 (SMPD3) cause PIEZO1 to switch to profoundly inactivating behavior. Ceramide (a key product of SMPD3) rescues non-inactivating channel behavior. Its co-product, phosphoryl choline, has no effect. In contrast to ceramide, sphingomyelin (the SMPD3 substrate) does not affect inactivation but alters channel force sensitivity. The data suggest that sphingomyelinase activity, ceramide, and sphingomyelin are determinants of native PIEZO gating that enable sustained activity.

•

SMPD3 sphingomyelinase enables long-lasting PIEZO1 activity in response to force

•

Ceramide, a key lipid product of SMPD3, promotes long-lasting activity

•

Sphingomyelin, the SMPD3 substrate, does not affect the duration of activity

•

Sphingomyelin alters PIEZO1 force sensitivity

Ultrafast energy transfer between lipid-linked chromophores and plant light-harvesting complex II

Light-Harvesting Complex II (LHCII) is a membrane protein found in plant chloroplasts that has the crucial role of absorbing solar energy and subsequently performing excitation energy transfer to the reaction centre subunits of Photosystem II. LHCII provides strong absorption of blue and red light, however, it has minimal absorption in the green spectral region where solar irradiance is maximal. In a recent proof-of-principle study, we enhanced the absorption in this spectral range by developing a biohybrid system where LHCII proteins together with lipid-linked Texas Red (TR) chromophores were assembled into lipid membrane vesicles. The utility of these systems was limited by significant LHCII quenching due to protein-protein interactions and heterogeneous lipid structures. Here, we organise TR and LHCII into a lipid nanodisc, which provides a homogeneous, well-controlled platform to study the interactions between TR molecules and single LHCII complexes. Fluorescence spectroscopy determined that TR-to-LHCII energy transfer has an efficiency of at least 60%, resulting in a 262% enhancement of LHCII fluorescence in the 525-625 nm range, two-fold greater than in the previous system. Ultrafast transient absorption spectroscopy revealed two time constants of 3.7 and 128 ps for TR-to-LHCII energy transfer. Structural modelling and theoretical calculations indicate that these timescales correspond to TR-lipids that are loosely- or tightly-associated with the protein, respectively, with estimated TR-to-LHCII separations of ∼3.5 nm and ∼1 nm. Overall, we demonstrate that a nanodisc-based biohybrid system provides an idealised platform to explore the photophysical interactions between extrinsic chromophores and membrane proteins with potential applications in understanding more complex natural or artificial photosynthetic systems.

Structure and Mechanism of Respiratory III-IV Supercomplexes in Bioenergetic Membranes

In the final steps of energy conservation in aerobic organisms, free energy from electron transfer through the respiratory chain is transduced into a proton electrochemical gradient across a membrane. In mitochondria and many bacteria, reduction of the dioxygen electron acceptor is catalyzed by cytochrome c oxidase (complex IV), which receives electrons from cytochrome bc1 (complex III), via membrane-bound or water-soluble cytochrome c. These complexes function independently, but in many organisms they associate to form supercomplexes. Here, we review the structural features and the functional significance of the nonobligate III2IV1/2 Saccharomyces cerevisiae mitochondrial supercomplex as well as the obligate III2IV2 supercomplex from actinobacteria. The analysis is centered around the Q-cycle of complex III, proton uptake by CytcO, as well as mechanistic and structural solutions to the electronic link between complexes III and IV.

The Role of the Membrane in Transporter Folding and Activity

The synthesis, folding, and function of membrane transport proteins are critical factors for defining cellular physiology. Since the stability of these proteins evolved amidst the lipid bilayer, it is no surprise that we are finding that many of these membrane proteins demonstrate coupling of their structure or activity in some way to the membrane. More and more transporter structures are being determined with some information about the surrounding membrane, and computational modeling is providing further molecular details about these solvation structures. Thus, the field is moving towards identifying which molecular mechanisms - lipid interactions, membrane perturbations, differential solvation, and bulk membrane effects - are involved in linking membrane energetics to transporter stability and function. In this review, we present an overview of these mechanisms and the growing evidence that the lipid bilayer is a major determinant of the fold, form, and function of membrane transport proteins in membranes.

Impact of Increased Membrane Realism on Conformational Sampling of Proteins

The realism and accuracy of lipid bilayer simulations through molecular dynamics (MD) are heavily dependent on the lipid composition. While the field is pushing toward implementing more heterogeneous and realistic membrane compositions, a lack of high-resolution lipidomic data prevents some membrane protein systems from being modeled with the highest level of realism. Given the additional diversity of real-world cellular membranes and protein-lipid interactions, it is still not fully understood how altering membrane complexity affects modeled membrane protein functions or if it matters over long-timescale simulations. This is especially true for organisms whose membrane environments have little to no computational study, such as the plant plasma membrane. Tackling these issues in tandem, a generalized, realistic, and asymmetric plant plasma membrane with more than 10 different lipid species is constructed herein. Classical MD simulations of pure membrane constructs were performed to evaluate how altering the compositional complexity of the membrane impacted the plant membrane properties. The apo form of a plant sugar transporter, OsSWEET2b, was inserted into membrane models where lipid diversity was calculated in either a size-dependent or size-independent manner. An adaptive sampling simulation regime validated by Markov-state models was performed to capture the gating dynamics of OsSWEET2b in each of these membrane constructs. In comparison to previous OsSWEET2b simulations performed in a pure POPC bilayer, we confirm that simulations performed within a native-like membrane composition alter the stabilization of apo OsSWEET2b conformational states by ∼1 kcal/mol. The free-energy barriers of intermediate conformational states decrease when realistic membrane complexity is simplified, albeit roughly within sampling error, suggesting that protein-specific responses to membranes differ due to altered packing caused by compositional fluctuations. This work serves as a case study where a more realistic bilayer composition makes unbiased conformational sampling easier to achieve than with simplified bilayers.

Molecular evolutionary insights from PRLR in mammals

Prolactin (PRL) is a pleiotropic neurohormone secreted by the mammalian pituitary gland into the blood, thus reaching many tissues and organs beyond the brain. PRL binds to its receptor, PRLR, eliciting a molecular signaling cascade. This system modulates essential mammalian behaviors and promotes notable modifications in the reproductive female tissues and organs. Here, we explore how the intracellular domain of PRLR (PRLR-ICD) modulates the expression of the PRLR gene. Despite differences in the reproductive strategies between eutherian and metatherian mammals, there is no clear distinction between PRLR-ICD functional motifs. However, we found selection signatures that showed differences between groups, with many conserved functional elements strongly maintained through purifying selection across the class Mammalia. We observed a few residues under relaxed selection, the levels of which were more pronounced in Eutheria and particularly striking in primates (Simiiformes), which could represent a pre-adaptive genetic element protected from purifying selection. Alternative, new motifs, such as YLDP (318-321) and others with residues Y283 and Y290, may already be functional. These motifs would have been co-opted in primates as part of a complex genetic repertoire related to some derived adaptive phenotypes, but these changes would have no impact on the primordial functions that characterize the mammals as a whole and that are related to the PRL-PRLR system.

Identification and assessment of cardiolipin interactions with E. coli inner membrane proteins

Integral membrane proteins are localized and/or regulated by lipids present in the surrounding bilayer. While bacteria have relatively simple membranes, there is ample evidence that many bacterial proteins bind to specific lipids, especially the anionic lipid cardiolipin. Here, we apply molecular dynamics simulations to assess lipid binding to 42 different Escherichia coli inner membrane proteins. Our data reveal an asymmetry between the membrane leaflets, with increased anionic lipid binding to the inner leaflet regions of the proteins, particularly for cardiolipin. From our simulations, we identify >700 independent cardiolipin binding sites, allowing us to identify the molecular basis of a prototypical cardiolipin binding site, which we validate against structures of bacterial proteins bound to cardiolipin. This allows us to construct a set of metrics for defining a high-affinity cardiolipin binding site on bacterial membrane proteins, paving the way for a heuristic approach to defining other protein-lipid interactions.

Aggregation and mobility of membrane proteins interplay with local lipid order in the plasma membrane of T cells

To disentangle the elusive lipid-protein interactions in T-cell activation, we investigate how externally imposed variations in mobility of key membrane proteins (T-cell receptor [TCR], kinase Lck, and phosphatase CD45) affect the local lipid order and protein colocalisation. Using spectral imaging with polarity-sensitive membrane probes in model membranes and live Jurkat T cells, we find that partial immobilisation of proteins (including TCR) by aggregation or ligand binding changes their preference towards a more ordered lipid environment, which can recruit Lck. Our data suggest that the cellular membrane is poised to modulate the frequency of protein encounters upon alterations of their mobility, for example in ligand binding, which offers new mechanistic insight into the involvement of lipid-mediated interactions in membrane-hosted signalling events.

Modulation of Phospholipid Bilayer Properties by Simvastatin

Simvastatin (Zocor) is one of the most prescribed drugs for reducing high cholesterol. Although simvastatin is ingested in its inactive lactone form, it is converted to its active dihydroxyheptanoate form by carboxylesterases in the liver. The dihydroxyheptanoate form can also be converted back to its original lactone form. Unfortunately, some of the side effects associated with the intake of simvastatin and other lipophilic statins at higher doses include statin-associated myopathy (SAM) and, in more severe cases, kidney failure. While the cause of SAM is unknown, it is hypothesized that these side effects are dependent on the localization of statins in lipid bilayers and their impact on bilayer properties. In this work, we carry out all-atom molecular dynamics simulations on both the lactone and dihydroxyheptanoate forms of simvastatin (termed "SN" and "SA", respectively) with a pure 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer and a POPC/cholesterol (30 mol %) binary mixture as membrane models. Additional simulations were carried out with multiple simvastatin molecules to mimic in vitro conditions that produced pleiotropic effects. Both SN and SA spontaneously diffused into the lipid bilayer, and a longer simulation time of 4 μs was needed for the complete incorporation of multiple SAs into the bilayer. By constructing potential mean force and electron density profiles, we find that SN localizes deeper within the hydrophobic interior of the bilayer and that SA has a greater tendency to form hydrogen-bonding interactions with neighboring water molecules and lipid headgroups. For the pure POPC bilayer, both SN and SA increase membrane order, while membrane fluidity increases for the POPC/cholesterol bilayer.

A Data-Driven Dimensionality Reduction Approach to Compare and Classify Lipid Force Fields

Molecular dynamics simulations of all-atom and coarse-grained lipid bilayer models are increasingly used to obtain useful insights for understanding the structural dynamics of these assemblies. In this context, one crucial point concerns the comparison of the performance and accuracy of classical force fields (FFs), which sometimes remains elusive. To date, the assessments performed on different classical potentials are mostly based on the comparison with experimental observables, which typically regard average properties. However, local differences of the structure and dynamics, which are poorly captured by average measurements, can make a difference, but these are nontrivial to catch. Here, we propose an agnostic way to compare different FFs at different resolutions (atomistic, united-atom, and coarse-grained), by means of a high-dimensional similarity metrics built on the framework of Smooth Overlap of Atomic Position (SOAP). We compare and classify a set of 13 FFs, modeling 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers. Our SOAP kernel-based metrics allows us to compare, discriminate, and correlate different FFs at different model resolutions in an unbiased, high-dimensional way. This also captures differences between FFs in modeling nonaverage events (originating from local transitions), for example, the liquid-to-gel phase transition in dipalmitoylphosphatidylcholine (DPPC) bilayers, for which our metrics allows us to identify nucleation centers for the phase transition, highlighting some intrinsic resolution limitations in implicit versus explicit solvent FFs.

Aggregation of Lysozyme in the Presence of a Mixed Bilayer of POPC and POPG

Understanding the molecular mechanisms by which amyloidogenic proteins interact with membranes is a challenging task. Amyloid accumulates from many human diseases have been observed to contain membrane lipids. In this work, coarse-grained molecular dynamics simulations have been used to inspect hen egg white lysozyme (HEWL) aggregation and membrane association in the presence of a pure POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) bilayer and a POPC and POPG (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylglycerol) mixed bilayer. It was observed that, in both cases, two HEWLs formed aggregates. In the presence of a mixed bilayer, after aggregation, the aggregated system started to interact with the membrane. It has been found that one of the lysozymes which came closer to the mixed bilayer unfolded more. The process of the initial insertion of an aggregated system in the mixed bilayer has been analyzed. The structural rearrangements of the protein and lipids were analyzed as well along the course of the simulation. Although with a pure POPC bilayer, aggregation was observed, the aggregated system moved away from the membrane. We believe that our study will provide considerable insights into lysozyme aggregation in the presence of a membrane environment.

Unravelling the structural complexity of protein-lipid interactions with neutron reflectometry

Neutron reflectometry (NR) is a large-facility technique used to examine structure at interfaces. In this brief review an introduction to the utilisation of NR in the study of protein-lipid interactions is given. Cold neutron beams penetrate matter deeply, have low energies, wavelengths in the Ångstrom regime and are sensitive to light elements. High differential hydrogen sensitivity (between protium and deuterium) enables solution and sample isotopic labelling to be utilised to enhance or diminish the scattering signal of individual components within complex biological structures. The combination of these effects means NR can probe buried structures such as those at the solid-liquid interface and encode molecular level structural information on interfacial protein-lipid complexes revealing the relative distribution of components as well as the overall structure. Model biological membrane sample systems can be structurally probed to examine phenomena such as antimicrobial mode of activity, as well as structural and mechanistic properties peripheral/integral proteins within membrane complexes. Here, the example of the antimicrobial protein α1-purothionin binding to a model Gram negative bacterial outer membrane is used to highlight the utilisation of this technique, detailing how changes in the protein/lipid distributions across the membrane before and after the protein interaction can be easily encoded using hydrogen isotope labelling.

Assessing the Role of Lipids in the Molecular Mechanism of Membrane Proteins

Membrane proteins have evolved to work optimally within the complex environment of the biological membrane. Consequently, interactions with surrounding lipids are part of their molecular mechanism. Yet, the identification of lipid-protein interactions and the assessment of their molecular role is an experimental challenge. Recently, biophysical approaches have emerged that are compatible with the study of membrane proteins in an environment closer to the biological membrane. These novel approaches revealed specific mechanisms of regulation of membrane protein function. Lipids have been shown to play a role in oligomerization, conformational transitions or allosteric coupling. In this review, we summarize the recent biophysical approaches, or combination thereof, that allow to decipher the role of lipid-protein interactions in the mechanism of membrane proteins.

ProLint: a web-based framework for the automated data analysis and visualization of lipid-protein interactions

The functional activity of membrane proteins is carried out in a complex lipid environment. Increasingly, it is becoming clear that lipids are an important player in regulating or generally modulating their activity. A routinely used method to gain insight into this interplay between lipids and proteins are Molecular Dynamics (MD) simulations, since they allow us to study interactions at atomic or near-atomic detail as a function of time. A major bottleneck, however, is analyzing and visualizing lipid-protein interactions, which, in practice, is a time-demanding task. Here, we present ProLint (www.prolint.ca), a webserver that completely automates analysis of MD generated files and visualization of lipid-protein interactions. Analysis is modular allowing users to select their preferred method, and visualization is entirely interactive through custom built applications that enable a detailed qualitative and quantitative exploration of lipid-protein interactions. ProLint also includes a database of published MD results that have been processed through the ProLint workflow and can be visualized by anyone regardless of their level of experience with MD. The automated analysis, feature-rich visualization, database integration, and open-source distribution with an easy to install process, will allow ProLint to become a routine workflow in lipid-protein interaction studies.

Multi-scale simulations of the T cell receptor reveal its lipid interactions, dynamics and the arrangement of its cytoplasmic region

The T cell receptor (TCR-CD3) initiates T cell activation by binding to peptides of Major Histocompatibility Complexes (pMHC). The TCR-CD3 topology is well understood but the arrangement and dynamics of its cytoplasmic tails remains unknown, limiting our grasp of the signalling mechanism. Here, we use molecular dynamics simulations and modelling to investigate the entire TCR-CD3 embedded in a model membrane. Our study demonstrates conformational changes in the extracellular and transmembrane domains, and the arrangement of the TCR-CD3 cytoplasmic tails. The cytoplasmic tails formed highly interlaced structures while some tyrosines within the immunoreceptor tyrosine-based activation motifs (ITAMs) penetrated the hydrophobic core of the membrane. Interactions between the cytoplasmic tails and phosphatidylinositol phosphate lipids in the inner membrane leaflet led to the formation of a distinct anionic lipid fingerprint around the TCR-CD3. These results increase our understanding of the TCR-CD3 dynamics and the importance of membrane lipids in regulating T cell activation.

Impact of Membrane Lipids on UapA and AzgA Transporter Subcellular Localization and Activity in Aspergillus nidulans

Recent biochemical and biophysical evidence have established that membrane lipids, namely phospholipids, sphingolipids and sterols, are critical for the function of eukaryotic plasma membrane transporters. Here, we study the effect of selected membrane lipid biosynthesis mutations and of the ergosterol-related antifungal itraconazole on the subcellular localization, stability and transport kinetics of two well-studied purine transporters, UapA and AzgA, in Aspergillus nidulans. We show that genetic reduction in biosynthesis of ergosterol, sphingolipids or phosphoinositides arrest A. nidulans growth after germling formation, but solely blocks in early steps of ergosterol (Erg11) or sphingolipid (BasA) synthesis have a negative effect on plasma membrane (PM) localization and stability of transporters before growth arrest. Surprisingly, the fraction of UapA or AzgA that reaches the PM in lipid biosynthesis mutants is shown to conserve normal apparent transport kinetics. We further show that turnover of UapA, which is the transporter mostly sensitive to membrane lipid content modification, occurs during its trafficking and by enhanced endocytosis, and is partly dependent on autophagy and Hect-type HulA^Rsp5 ubiquitination. Our results point out that the role of specific membrane lipids on transporter biogenesis and function in vivo is complex, combinatorial and transporter-dependent.

Effect of Membrane Environment on the Ligand-Binding Properties of the Terminal Oxidase Cytochrome bd-I from Escherichia coli

Cytochrome bd-I is a terminal oxidase of the Escherichia coli respiratory chain. This integral membrane protein contains three redox-active prosthetic groups (hemes b₅₅₈, b₅₉₅, and d) and couples the electron transfer from quinol to molecular oxygen to the generation of proton motive force, as one of its important physiological functions. The study was aimed at examining the effect of the membrane environment on the ligand-binding properties of cytochrome bd-I by absorption spectroscopy. The membrane environment was found to modulate the ligand-binding characteristics of the hemoprotein in both oxidized and reduced states. Absorption changes upon the addition of exogenous ligands, such as cyanide or carbon monoxide (CO), to the detergent-solubilized enzyme were much more significant and heterogeneous than those observed with the membrane-bound enzyme. In the native membranes, both cyanide and CO interacted mainly with heme d. An additional ligand-binding site (heme b₅₅₈) appeared in the isolated enzyme, as was evidenced by more pronounced changes in the absorption in the Soret band. This additional reactivity could also be detected after treatment of E. coli membranes with a detergent. The observed effect did not result from the enzyme denaturation, since reconstitution of the isolated enzyme into azolectin liposomes restored the ligand-binding pattern close to that observed for the intact membranes.

YPIBP: A repository for phosphoinositide-binding proteins in yeast

Phosphoinositides (PIs) are a family of eight lipids consisting of phosphatidylinositol (PtdIns) and its seven phosphorylated forms. PIs have important regulatory functions in the cell including lipid signaling, protein transport, and membrane trafficking. Yeast has been recognized as a eukaryotic model system to study lipid-protein interactions. Hundreds of yeast PI-binding proteins have been identified, but this research knowledge remains scattered. Besides, the complete PI-binding spectrum and potential PI-binding domains have not been interlinked. No comprehensive databases are available to support the lipid-protein interaction research on phosphoinositides. Here we constructed the first knowledgebase of Yeast Phosphoinositide-Binding Proteins (YPIBP), a repository consisting of 679 PI-binding proteins collected from high-throughput proteome-array and lipid-array studies, QuickGO, and a rigorous literature mining. The YPIBP also contains protein domain information in categories of lipid-binding domains, lipid-related domains and other domains. The YPIBP provides search and browse modes along with two enrichment analyses (PI-binding enrichment analysis and domain enrichment analysis). An interactive visualization is given to summarize the PI-domain-protein interactome. Finally, three case studies were given to demonstrate the utility of YPIBP. The YPIBP knowledgebase consolidates the present knowledge and provides new insights of the PI-binding proteins by bringing comprehensive and in-depth interaction network of the PI-binding proteins. YPIBP is available at http://cosbi7.ee.ncku.edu.tw/YPIBP/.

Physical Characterization of Triolein and Implications for Its Role in Lipid Droplet Biogenesis

Lipid droplets (LDs) are neutral lipid-storing organelles surrounded by a phospholipid (PL) monolayer. At present, how LDs are formed in the endoplasmic reticulum (ER) bilayer is poorly understood. In this study, we present a revised all-atom (AA) triolein (TG) model, the main constituent of the LD core, and characterize its properties in a bilayer membrane to demonstrate the implications of its behavior in LD biogenesis. In bilayer simulations, TG resides at the surface, adopting PL-like conformations (denoted in this work as SURF-TG). Free energy sampling simulation results estimate the barrier for TG relocating from the bilayer surface to the bilayer center to be ∼2 kcal/mol in the absence of an oil lens. SURF-TG is able to modulate membrane properties by increasing PL ordering, decreasing bending modulus, and creating local negative curvature. The other neutral lipid, dioleoyl-glycerol (DAG), also reduces the membrane bending modulus and populates negative curvature regions. A phenomenological coarse-grained (CG) model is also developed to observe larger-scale SURF-TG-mediated membrane deformation. CG simulations confirm that TG nucleates between the bilayer leaflets at a critical concentration when SURF-TG is evenly distributed. However, when one monolayer contains more SURF-TG, the membrane bends toward the other leaflet, followed by TG nucleation if a concentration is higher than the critical threshold. The central conclusion of this study is that SURF-TG is a negative curvature inducer, as well as a membrane modulator. To this end, a model is proposed in which the accumulation of SURF-TG in the luminal leaflet bends the ER bilayer toward the cytosolic side, followed by TG nucleation.

Binding of Ca2+-independent C2 domains to lipid membranes: A multi-scale molecular dynamics study

C2 domains facilitate protein interactions with lipid bilayers in either a Ca²⁺-dependent or -independent manner. We used molecular dynamics (MD) simulations to explore six Ca²⁺-independent C2 domains, from KIBRA, PI3KC2α, RIM2, PTEN, SHIP2, and Smurf2. In coarse-grained MD simulations these C2 domains formed transient interactions with zwitterionic bilayers, compared with longer-lived interactions with anionic bilayers containing phosphatidylinositol bisphosphate (PIP₂). Type I C2 domains bound non-canonically via the front, back, or side of the β sandwich, whereas type II C2 domains bound canonically, via the top loops. C2 domains interacted strongly with membranes containing PIP₂, causing bound anionic lipids to cluster around the protein. Binding modes were refined via atomistic simulations. For PTEN and SHIP2, CG simulations of their phosphatase plus C2 domains with PIP₂-containing bilayers were also performed, and the roles of the two domains in membrane localization compared. These studies establish a simulation protocol for membrane-recognition proteins.

•

Binding of Ca²⁺-independent C2 domains to membranes was explored by MD simulation

•

C2 domains from KIBRA, PI3KC2α, RIM2, PTEN, SHIP2, and Smurf2 were compared

•

C2 domains formed longer-lived interactions with lipid bilayers containing PIP₂

•

For PTEN and SHIP2, simulations of their phosphatase plus C2 domains were performed

Interfacial binding sites for cholesterol on GABAA receptors and competition with neurosteroids

γ-Aminobutyric acid type A (GABA_A) receptors in the brain are located in the outer membranes of brain cells where the concentration of cholesterol is high. Of the 25 available high-resolution structures available for GABA_A receptors, none were determined in the presence of cholesterol, but four include resolved molecules of cholesterol hemisuccinate (CHS). Here, a molecular docking procedure is used to sweep the transmembrane (TM) surfaces of the receptors for cholesterol binding sites. Cholesterol docking poses determined in this way match 89% of the resolved CHS when CHS molecules deemed unlikely to represent typical bound cholesterols are excluded. The receptors are pentameric, and their TM surfaces consist of a set of five facets, each including pairs of TM helices from two adjacent subunits. Each facet contains hydrophobic hollows running from one side of the membrane to the other, within which are six potential binding sites for cholesterol, three on each side of the membrane. High-resolution structures of GABA_A receptors with bound neurosteroids show that neurosteroids bind in these cholesterol binding sites, so the binding of neurosteroids and cholesterol will be competitive.

Computer simulations of a heterogeneous membrane with enhanced sampling techniques

Computational determination of the equilibrium state of heterogeneous phospholipid membranes is a significant challenge. We wish to explore the rich phase diagram of these multi-component systems. However, the diffusion and mixing times in membranes are long compared to typical time scales of computer simulations. Here, we evaluate the combination of the enhanced sampling techniques molecular dynamics with alchemical steps and Monte Carlo with molecular dynamics with a coarse-grained model of membranes (Martini) to reduce the number of steps and force evaluations that are needed to reach equilibrium. We illustrate a significant gain compared to straightforward molecular dynamics of the Martini model by factors between 3 and 10. The combination is a useful tool to enhance the study of phase separation and the formation of domains in biological membranes.

Molecular mechanisms of ion conduction and ion selectivity in TMEM16 lipid scramblases

TMEM16 lipid scramblases transport lipids and also operate as ion channels with highly variable ion selectivities and various physiological functions. However, their molecular mechanisms of ion conduction and selectivity remain largely unknown. Using computational electrophysiology simulations at atomistic resolution, we identified the main ion-conductive state of TMEM16 lipid scramblases, in which an ion permeation pathway is lined by lipid headgroups that directly interact with permeating ions in a voltage polarity-dependent manner. We found that lipid headgroups modulate the ion-permeability state and regulate ion selectivity to varying degrees in different scramblase isoforms, depending on the amino-acid composition of the pores. Our work has defined the structural basis of ion conduction and selectivity in TMEM16 lipid scramblases and uncovered the mechanisms responsible for the direct effects of membrane lipids on the conduction properties of ion channels.

Novel Roles of SH2 and SH3 Domains in Lipid Binding

Signal transduction, the ability of cells to perceive information from the surroundings and alter behavior in response, is an essential property of life. Studies on tyrosine kinase action fundamentally changed our concept of cellular regulation. The induced assembly of subcellular hubs via the recognition of local protein or lipid modifications by modular protein interactions is now a central paradigm in signaling. Such molecular interactions are mediated by specific protein interaction domains. The first such domain identified was the SH2 domain, which was postulated to be a reader capable of finding and binding protein partners displaying phosphorylated tyrosine side chains. The SH3 domain was found to be involved in the formation of stable protein sub-complexes by constitutively attaching to proline-rich surfaces on its binding partners. The SH2 and SH3 domains have thus served as the prototypes for a diverse collection of interaction domains that recognize not only proteins but also lipids, nucleic acids, and small molecules. It has also been found that particular SH2 and SH3 domains themselves might also bind to and rely on lipids to modulate complex assembly. Some lipid-binding properties of SH2 and SH3 domains are reviewed here.

All-Atom Molecular Dynamics Elucidating Molecular Mechanisms of Single-Transmembrane Model Peptide Dimerization in a Lipid Bilayer

Protein-protein interactions between transmembrane helices are essential elements for membrane protein structures and functions. To understand the effects of peptide sequences and lipid compositions on these interactions, single-molecule experiments using model systems comprising artificial peptides and membranes have been extensively performed. However, their dynamic behavior at the atomic level remains largely unclear. In this study, we applied the all-atom molecular dynamics (MD) method to simulate the interactions of single-transmembrane helical peptide dimers in membrane environments, which has previously been analyzed by single-molecule experiments. The simulations were performed with two peptides (Ala- and Leu-based artificially designed peptides, termed "host peptide", and the host peptide added with the GXXXG motif, termed "GXXXG peptide"), two membranes (pure-POPC and POPC mixed with 30% cholesterols), and two dimer directions (parallel and antiparallel), consistent with those in the previous experiment. As a result, the MD simulations with parallel dimers reproduced the experimentally observed tendency that introducing cholesterols weakened the interactions in the GXXXG dimer and facilitated those in the host dimer. Our simulation suggested that the host dimer formed hydrogen bonds but the GXXXG dimer did not. However, some discrepancies were also observed between the experiments and simulations. Limitations in the space and time scales of simulations restrict the large-scale undulation and peristaltic motions of the membranes, resulting in differences in lateral pressure profiles. This effect could cause a discrepancy in the rotation angles of helices against the membrane normal.

Identification of multiple substrate binding sites in SLC4 transporters in the outward-facing conformation: Insights into the transport mechanism

Solute carrier family 4 (SLC4) transporters mediate the transmembrane transport of HCO₃^-, CO₃^2-, and Cl^- necessary for pH regulation, transepithelial H⁺/base transport, and ion homeostasis. Substrate transport with varying stoichiometry and specificity is achieved through an exchange mechanism and/or through coupling of the uptake of anionic substrates to typically co-transported Na⁺. Recently solved outward-facing structures of two SLC4 members (human anion exchanger 1 [hAE1] and human electrogenic sodium bicarbonate cotransporter 1 [hNBCe1]) with different transport modes (Cl^-/HCO₃^- exchange versus Na⁺-CO₃^2- symport) revealed highly conserved three-dimensional organization of their transmembrane domains. However, the exact location of the ion binding sites and their protein-ion coordination motifs are still unclear. In the present work, we combined site identification by ligand competitive saturation mapping and extensive molecular dynamics sampling with functional mutagenesis studies which led to the identification of two substrate binding sites (entry and central) in the outward-facing states of hAE1 and hNBCe1. Mutation of residues in the identified binding sites led to impaired transport in both proteins. We also showed that R730 in hAE1 is crucial for anion binding in both entry and central sites, whereas in hNBCe1, a Na⁺ acts as an anchor for CO₃^2- binding to the central site. Additionally, protonation of the central acidic residues (E681 in hAE1 and D754 in hNBCe1) alters the ion dynamics in the permeation cavity and may contribute to the transport mode differences in SLC4 proteins. These results provide a basis for understanding the functional differences between hAE1 and hNBCe1 and may facilitate potential drug development for diseases such as proximal and distal renal tubular acidosis.

Evaluating the Efficiency of the Martini Force Field to Study Protein Dimerization in Aqueous and Membrane Environments

Protein-protein complex assembly is one of the major drivers of biological response. Understanding the mechanisms of protein oligomerization/dimerization would allow one to elucidate how these complexes participate in biological activities and could ultimately lead to new approaches in designing novel therapeutic agents. However, determining the exact association pathways and structures of such complexes remains a challenge. Here, we use parallel tempering metadynamics simulations in the well-tempered ensemble to evaluate the performance of Martini 2.2P and Martini open-beta 3 (Martini 3) force fields in reproducing the structure and energetics of the dimerization process of membrane proteins and proteins in an aqueous solution in reasonable accuracy and throughput. We find that Martini 2.2P systematically overestimates the free energy of association by estimating large barriers in distinct areas, which likely leads to overaggregation when multiple monomers are present. In comparison, the less viscous Martini 3 results in a systematic underestimation of the free energy of association for proteins in solution, while it performs well in describing the association of membrane proteins. In all cases, the near-native dimer complexes are identified as minima in the free energy surface albeit not always as the lowest minima. In the case of Martini 3, we find that the spurious supramolecular protein aggregation present in Martini 2.2P multimer simulations is alleviated and thus this force field may be more suitable for the study of protein oligomerization. We propose that the use of enhanced sampling simulations with a refined coarse-grained force field and appropriately defined collective variables is a robust approach for studying the protein dimerization process, although one should be cautious of the ranking of energy minima.

Molecular interactions of the M and E integral membrane proteins of SARS-CoV-2

Specific lipid-protein interactions are key for cellular processes, and even more so for the replication of pathogens. The COVID-19 pandemic has drastically changed our lives and cause the death of nearly three million people worldwide, as of this writing. SARS-CoV-2 is the virus that causes the disease and has been at the center of scientific research over the past year. Most of the research on the virus is focused on key players during its initial attack and entry into the cellular host; namely the S protein, its glycan shield, and its interactions with the ACE2 receptors of human cells. As cases continue to raise around the globe, and new mutants are identified, there is an urgent need to understand the mechanisms of this virus during different stages of its life cycle. Here, we consider two integral membrane proteins of SARS-CoV-2 known to be important for viral assembly and infectivity. We have used microsecond-long all-atom molecular dynamics to examine the lipid-protein and protein-protein interactions of the membrane (M) and envelope (E) structural proteins of SARS-CoV-2 in a complex membrane model. We contrast the two proposed protein complexes for each of these proteins, and quantify their effect on their local lipid environment. This ongoing work also aims to provide molecular-level understanding of the mechanisms of action of this virus to possibly aid in the design of novel treatments.

Endoplasmic Reticulum-Plasma Membrane Contact Sites as an Organizing Principle for Compartmentalized Calcium and cAMP Signaling

In eukaryotic cells, ultimate specificity in activation and action-for example, by means of second messengers-of the myriad of signaling cascades is primordial. In fact, versatile and ubiquitous second messengers, such as calcium (Ca2+) and cyclic adenosine monophosphate (cAMP), regulate multiple-sometimes opposite-cellular functions in a specific spatiotemporal manner. Cells achieve this through segregation of the initiators and modulators to specific plasma membrane (PM) subdomains, such as lipid rafts and caveolae, as well as by dynamic close contacts between the endoplasmic reticulum (ER) membrane and other intracellular organelles, including the PM. Especially, these membrane contact sites (MCSs) are currently receiving a lot of attention as their large influence on cell signaling regulation and cell physiology is increasingly appreciated. Depletion of ER Ca2+ stores activates ER membrane STIM proteins, which activate PM-residing Orai and TRPC Ca2+ channels at ER-PM contact sites. Within the MCS, Ca2+ fluxes relay to cAMP signaling through highly interconnected networks. However, the precise mechanisms of MCS formation and the influence of their dynamic lipid environment on their functional maintenance are not completely understood. The current review aims to provide an overview of our current understanding and to identify open questions of the field.

Assembly of Model Membrane Nanodiscs for Native Mass Spectrometry

Native mass spectrometry (MS) with nanodiscs is a promising technique for characterizing membrane protein and peptide interactions in lipid bilayers. However, prior studies have used nanodiscs made of only one or two lipids, which lack the complexity of a natural lipid bilayer. To better model specific biological membranes, we developed model mammalian, bacterial, and mitochondrial nanodiscs with up to four different phospholipids. Careful selection of lipids with similar masses that balance the fluidity and curvature enabled these complex nanodiscs to be assembled and resolved with native MS. We then applied this approach to characterize the specificity and incorporation of LL-37, a human antimicrobial peptide, in single-lipid nanodiscs versus model bacterial nanodiscs. Overall, development of these model membrane nanodiscs reveals new insights into the assembly of complex nanodiscs and provides a useful toolkit for studying membrane protein, peptide, and lipid interactions in model biological membranes.

An ω-3, but Not an ω-6 Polyunsaturated Fatty Acid Decreases Membrane Dipole Potential and Stimulates Endo-Lysosomal Escape of Penetratin

Although the largely positive intramembrane dipole potential (DP) may substantially influence the function of transmembrane proteins, its investigation is deeply hampered by the lack of measurement techniques suitable for high-throughput examination of living cells. Here, we describe a novel emission ratiometric flow cytometry method based on F66, a 3-hydroxiflavon derivative, and demonstrate that 6-ketocholestanol, cholesterol and 7-dehydrocholesterol, saturated stearic acid (SA) and ω-6 γ-linolenic acid (GLA) increase, while ω-3 α-linolenic acid (ALA) decreases the DP. These changes do not correlate with alterations in cell viability or membrane fluidity. Pretreatment with ALA counteracts, while SA or GLA enhances cholesterol-induced DP elevations. Furthermore, ALA (but not SA or GLA) increases endo-lysosomal escape of penetratin, a cell-penetrating peptide. In summary, we have developed a novel method to measure DP in large quantities of individual living cells and propose ALA as a physiological DP lowering agent facilitating cytoplasmic entry of penetratin.

Perspectives on High-Throughput Ligand/Protein Docking With Martini MD Simulations

Molecular docking is central to rational drug design. Current docking techniques suffer, however, from limitations in protein flexibility and solvation models and by the use of simplified scoring functions. All-atom molecular dynamics simulations, on the other hand, feature a realistic representation of protein flexibility and solvent, but require knowledge of the binding site. Recently we showed that coarse-grained molecular dynamics simulations, based on the most recent version of the Martini force field, can be used to predict protein/ligand binding sites and pathways, without requiring any a priori information, and offer a level of accuracy approaching all-atom simulations. Given the excellent computational efficiency of Martini, this opens the way to high-throughput drug screening based on dynamic docking pipelines. In this opinion article, we sketch the roadmap to achieve this goal.

Lipid regulation of hERG1 channel function

The lipid regulation of mammalian ion channel function has emerged as a fundamental mechanism in the control of electrical signalling and transport specificity in various cell types. In this work, we combine molecular dynamics simulations, mutagenesis, and electrophysiology to provide mechanistic insights into how lipophilic molecules (ceramide-sphingolipid probe) alter gating kinetics and K+ currents of hERG1. We show that the sphingolipid probe induced a significant left shift of activation voltage, faster deactivation rates, and current blockade comparable to traditional hERG1 blockers. Microseconds-long MD simulations followed by experimental mutagenesis elucidated ceramide specific binding locations at the interface between the pore and voltage sensing domains. This region constitutes a unique crevice present in mammalian channels with a non-swapped topology. The combined experimental and simulation data provide evidence for ceramide-induced allosteric modulation of the channel by a conformational selection mechanism.

Membrane signalosome: where biophysics meets systems biology

We opine on the recent advances in experiments and modeling of modular signaling complexes assembled on mammalian cell membranes (membrane signalosomes) in the context of several applications including intracellular trafficking, cell migration, and immune response. Characterizing the individual components of the membrane assemblies at the nanoscale, ranging from protein-lipid and protein-protein interactions, to membrane morphology, and the energetics of emergent assemblies at the subcellular to cellular scales pose significant challenges. Overcoming these challenges through the iterative coupling of multiscale modeling and experiment can be transformative in terms of addressing the gaps between structural biology and super-resolution microscopy, as it holds the key to the discovery of fundamental mechanisms behind the emergence of function in the membrane signalosome.

Detection of ligand binding to purified HCN channels using fluorescence-based size exclusion chromatography

Biochemical measurements of ligand binding to eukaryotic membrane proteins are challenging because they can require large amounts of purified protein. For this reason, ligand binding is preferentially evaluated on soluble domains rather than on the full length proteins. In this chapter, we describe the use of fluorescence size exclusion chromatography-based thermostability (FSEC-TS) as an assay to monitor ligand binding to the full length mammalian ion channel HCN4. FSEC-TS monitors the effect of the ligand on the thermal denaturation curve of the protein by following the fluorescence of a fused GFP protein. Changes in the melting temperature (Tm) provide a quantitative value for measuring ligand-protein interaction. As a proof of concept, we describe here the protocol for monitoring the binding of the second messenger cAMP and of the known HCN drug Ivabradine to the purified GFP-HCN4 channel. cTMP, a known non-binder of HCN channels, is used as a control. Due to the small amount of protein required, the assay represents a high-throughput screening system for evaluating binding of small molecules to full length proteins.

Effects of Melatonin on Dairy Herd Improvement (DHI) of Holstein Cow with High SCS

Mastitis is a common disease in cows breeding. The milk quality will be significantly reduced with increased milk somatic cells, which often occurs in cows with mastitis. In this study, the influence of seasonal changes, age and lactation stages in the Dairy Herd Improvement (DHI) of cows was investigated. Then, the Dairy Herd Improvement (DHI) of cows with high somatic cell score (SCS) after melatonin treatment was systemically investigated. The results showed that melatonin significantly suppressed the milk somatic cell score under all of the tested conditions. The melatonin treatment also improved the milk nutritional value by reducing its fat but increasing its lactose and protein contents. The application of melatonin significantly improved the DHI. The beneficial effects of melatonin on DHI are likely attributed to the antioxidant and anti-inflammatory activities of melatonin.

State-dependent protein-lipid interactions of a pentameric ligand-gated ion channel in a neuronal membrane

Pentameric ligand-gated ion channels (pLGICs) are receptor proteins that are sensitive to their membrane environment, but the mechanism for how lipids modulate function under physiological conditions in a state dependent manner is not known. The glycine receptor is a pLGIC whose structure has been resolved in different functional states. Using a realistic model of a neuronal membrane coupled with coarse-grained molecular dynamics simulations, we demonstrate that some key lipid-protein interactions are dependent on the receptor state, suggesting that lipids may regulate the receptor's conformational dynamics. Comparison with existing structural data confirms known lipid binding sites, but we also predict further protein-lipid interactions including a site at the communication interface between the extracellular and transmembrane domain. Moreover, in the active state, cholesterol can bind to the binding site of the positive allosteric modulator ivermectin. These protein-lipid interaction sites could in future be exploited for the rational design of lipid-like allosteric drugs.

Influence of Antiplatelet Agents on the Lipid Composition of Platelet Plasma Membrane: A Lipidomics Approach with Ticagrelor and Its Active Metabolite

Lipids contained in the plasma membrane of platelets play an important role in platelet function. Modifications in the lipid composition can fluidify or rigidify the environment around embedded receptors, in order to facilitate the access of the receptor by the drug. However, data concerning the lipid composition of platelet plasma membrane need to be updated. In addition, data on the impact of drugs on plasma membrane composition, in particular antiplatelet agents, remain sparse. After isolation of platelet plasma membrane, we assessed, using lipidomics, the effect of ticagrelor, a P2Y12 antagonist, and its active metabolite on the lipid composition of these plasma membranes. We describe the exact lipid composition of plasma membrane, including all sub-species. Ticagrelor and its active metabolite significantly increased cholesterol and phosphatidylcholine ether with short saturated acyl chains 16:0/16:0, and decreased phosphatidylcholine, suggesting overall rigidification of the membrane. Furthermore, ticagrelor and its active metabolite decreased some arachidonylated plasmalogens, suggesting a decrease in availability of arachidonic acid from the membrane phospholipids for synthesis of biologically active mediators. To conclude, ticagrelor and its active metabolite seem to influence the lipid environment of receptors embedded in the lipid bilayer and modify the behavior of the plasma membrane.

Wrapping Up Viruses at Multiscale Resolution: Optimizing PACKMOL and SIRAH Execution for Simulating the Zika Virus

Simulating huge biomolecular complexes of million atoms at relevant biological time scales is becoming accessible to the broad scientific community. That proves to be crucial for urgent responses against emergent diseases in real time. Yet, there are still issues to sort regarding the system setup so that molecular dynamics (MD) simulations can be run in a simple and standard way. Here, we introduce an optimized pipeline for building and simulating enveloped virus-like particles (VLP). First, the membrane packing problem is tackled with new features and optimized options in PACKMOL. This allows preparing accurate membrane models of thousands of lipids in the context of a VLP within a few hours using a single CPU. Then, the assembly of the VLP system is done within the multiscale framework of the coarse-grained SIRAH force field. Finally, the equilibration protocol provides a system ready for production MD simulations within a few days on broadly accessible GPU resources. The pipeline is applied to study the Zika virus as a test case for large biomolecular systems. The VLP stabilizes at approximately 0.5 μs of MD simulation, reproducing correlations greater than 0.90 against experimental density maps from cryo-electron microscopy. Detailed structural analysis of the protein envelope also shows very good agreement in root-mean-square deviations and B-factors with the experimental data. The level of details attained shows for the first time a possible role for anionic phospholipids in stabilizing the envelope. Combining an efficient and reliable setup procedure with an accurate coarse-grained force field provides a valuable pipeline for simulating arbitrary viral systems or subcellular compartments, paving the way toward whole-cell simulations.

Cholesterol content in the membrane promotes key lipid-protein interactions in a pentameric serotonin-gated ion channel

Pentameric ligand-gated ion channels (pLGICs), embedded in the lipid membranes of nerve cells, mediate fast synaptic transmission and are major pharmaceutical targets. Because of their complexity and the limited knowledge of their structure, their working mechanisms have still to be fully unraveled at the molecular level. Over the past few years, evidence that the lipid membrane may modulate the function of membrane proteins, including pLGICs, has emerged. Here, we investigate, by means of molecular dynamics simulations, the behavior of the lipid membrane at the interface with the 5-HT_3A receptor (5-HT_3AR), a representative pLGIC which is the target of nausea-suppressant drugs, in a nonconductive state. Three lipid compositions are studied, spanning different concentrations of the phospholipids, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine, and of cholesterol, hence a range of viscosities. A variety of lipid interactions and persistent binding events to different parts of the receptor are revealed in the investigated models, providing snapshots of the dynamical environment at the membrane-receptor interface. Some of these events result in lipid intercalation within the transmembrane domain, and others reach out to protein key sections for signal transmission and receptor activation, such as the Cys-loop and the M2-M3 loop. In particular, phospholipids, with their long hydrophobic tails, play an important role in these interactions, potentially providing a bridge between these two structures. A higher cholesterol content appears to promote lipid persistent binding to the receptor.

Investigating the lipid fingerprint of SLC6 neurotransmitter transporters: a comparison of dDAT, hDAT, hSERT, and GlyT2

The local lipid annulus, or "fingerprint", of four SLC6 transporters (dDAT, hDAT, hSERT, and GlyT2) embedded in a complex neuronal membrane were compared and characterised using molecular dynamics. Our analysis included the development of new tools to improve membrane leaflet detection and the analysis of leaflet-dependent properties. Overall, the lipid fingerprints of the four transporters are comprised of similar lipids when grouped by headgroup or tail saturation. The enrichment and depletion of specific lipids, including sites of cholesterol contacts, varies between transporters. The subtle differences in lipid fingerprints results in varying membrane biophysical properties near the transporter. Our results highlight that the lipid-fingerprint of SLC6 transporters in complex membranes is highly dependent on membrane composition. Our results further characterize how the presence and identity of membrane proteins affects the complex interplay of lipid-protein interactions, influencing the local lipid environment and membrane biophysical properties.

Oxygen Transport Parameter in Plasma Membrane of Eye Lens Fiber Cells by Saturation Recovery EPR

A probability distribution of rate constants contained within an exponential-like saturation recovery (SR) electron paramagnetic resonance signal can be constructed using stretched exponential function fitting parameters. Previously (Stein et al. Appl. Magn. Reson. 2019.), application of this method was limited to the case where only one relaxation process, namely spin-lattice relaxations due to the rotational diffusion of the spin labels in the intact eye-lens membranes, contributed to an exponential-like SR signal. These conditions were achieved for thoroughly deoxygenated samples. Here, the case is described where the second relaxation process, namely Heisenberg exchange between the spin label and molecular oxygen that occurs during bimolecular collisions, contributes to the decay of SR signals. We have further developed the theory for application of stretched exponential function to analyze SR signals involving these two processes. This new approach allows separation of stretched exponential parameters, namely characteristic stretched rates and heterogeneity parameters for both processes. Knowing these parameters allowed us to separately construct the probability distributions of spin-lattice relaxation rates determined by the rotational diffusion of spin labels and the distribution of relaxations induced strictly by collisions with molecular oxygen. The later distribution is determined by the distribution of oxygen diffusion concentration products within the membrane, which forms a sensitive new way to describe membrane fluidity and heterogeneity. This method was validated in silico and by fitting SR signals from spin-labeled intact nuclear fiber cell plasma membranes extracted from porcine eye lenses equilibrated with different fractions of air.

Function-related asymmetry of the specific cardiolipin binding sites on the mitochondrial ADP/ATP carrier

The ADP/ATP carrier (AAC) transports matrix ATP and cytosolic ADP across the inner mitochondrial membrane (IMM). It is well known that cardiolipin (CL) plays an important role in regulating the function of AAC, yet the underlying mechanism still remains elusive. AAC is composed of three homologous domains, and three specific CL binding sites are located at the domain-domain interfaces near the matrix side. Here we report an in-depth investigation on the dynamic properties of the bound CL within the three specific sites through all-atom molecular dynamics simulations of up to 13 μs in total. Our results highlight the importance of the basic and polar residues in CL binding. The basic residues from the linker helix and/or the [Y/W/F][K/R]G motif enable the bound CL to form an intra-domain binding mode, and the canonical inter-domain binding mode only forms when these basic residues are occupied by an additional phospholipid. Of special significance, differences in the basic and polar residues lead to remarkable asymmetry among the three specific CL binding sites. We found that the bound CL at the interface of domains 2 and 3 predominantly adopts inter-domain binding mode, while CLs at the other two sites have much more intra-domain populations. This is consistent with the asymmetric crystal structure of the matrix state (m-state) AAC which implies an asymmetric transport mechanism. The dynamic equilibrium between the inter-domain and intra-domain binding modes observed in our simulations could be highly important for the bound CLs to adapt to the movements during state transitions.

Minimizing isotropic and deviatoric membrane energy - An unifying formation mechanism of different cellular membrane nanovesicle types

Tiny membrane-enclosed cellular fragments that can mediate interactions between cells and organisms have recently become a subject of increasing attention. In this work the mechanism of formation of cell membrane nanovesicles (CNVs) was studied experimentally and theoretically. CNVs were isolated by centrifugation and washing of blood cells and observed by optical microscopy and scanning electron microscopy. The shape of the biological membrane in the budding process, as observed in phospholipid vesicles, in erythrocytes and in CNVs, was described by an unifying model. Taking the mean curvature h and the curvature deviator d of the membrane surface as the relevant parameters, the shape and the distribution of membrane constituents were determined theoretically by minimization of membrane free energy. Considering these results and previous results on vesiculation of red blood cells it was interpreted that the budding processes may lead to formation of different types of CNVs as regards the compartment (exo/endovesicles), shape (spherical/tubular/torocytic) and composition (enriched/depleted in particular kinds of molecules). It was concluded that the specificity of pinched off nanovesicles derives from the shape of the membrane constituents and not primarily from their chemical identity, which explains evidences on great heterogeneity of isolated extracellular vesicles with respect to composition.

One of the amazing properties of a biological membrane is the ability to undergo dramatic changes of its shape. It may exhibit very high curvature and thereby enclose nano-sized compartments that pinch off from the mother membrane and become freely moving cellular nanovesicles (CNVs). CNVs externalize the pieces of the cell and make them available to other cells within the same organism or other organisms. Therefore they have been acknowledged as mediators of communication between microorganisms, plants, animals and human. Furthernore, they dwell on the border between living and non-living things. Recent findings report on heterogeneity of the size and composition of CNVs found in isolates from different biological samples. As communication between cells is involved in many physiological and patophysiological processes, it is of importance to understand the mechanisms of CNVs formation and recognize the natural laws that mainly govern them. We point to an unifying mechanism that explains stability of differently shaped and composed CNVs by taking into account that the biological membrane tends to attain the minimum of its relevant energy. Conveniently, the procedure can be described by a mathematical model which allows for transparent comparison between experimentally induced shapes of membrane-enclosed vesicular structures and numerical calculations.

Advances in understanding and in multi-disciplinary methodology used to assess lipid regulation of signalling cascades from the cancer cell plasma membrane

The lipid bilayer is a functional component of cells, forming a stable platform for the initiation of key biological processes, including cell signalling. There are distinct changes in the lipid composition of cell membranes during oncogenic transformation resulting in aberrant activation and inactivation of signalling transduction pathways. Studying the role of the cell membrane in cell signalling is challenging, since techniques are often limited to by timescale, resolution, sensitivity, and averaging. To overcome these limitations, combining 'computational', 'wet-lab' and 'semi-dry' approaches offers the best opportunity to resolving complex biological processes involved in membrane organisation. In this review, we highlight analytical tools that have been applied for the study of cell signalling initiation from the cancer cell membranes through computational microscopy, biological assays, and membrane biophysics. The cancer therapeutic potential of extracellular membrane-modulating agents, such as cholesterol-reducing agents is also discussed, as is the need for future collaborative inter-disciplinary research for studying the role of the cell membrane and its components in cancer therapy.

Pore-forming proteins: From defense factors to endogenous executors of cell death

Pore-forming proteins (PFPs) and small antimicrobial peptides (AMPs) represent a large family of molecules with the common ability to punch holes in cell membranes to alter their permeability. They play a fundamental role as infectious bacteria's defensive tools against host's immune system and as executors of endogenous machineries of regulated cell death in eukaryotic cells. Despite being highly divergent in primary sequence and 3D structure, specific folds of pore-forming domains have been conserved. In fact, pore formation is considered an ancient mechanism that takes place through a general multistep process involving: membrane partitioning and insertion, oligomerization and pore formation. However, different PFPs and AMPs assemble and form pores following different mechanisms that could end up either in the formation of protein-lined or protein-lipid pores. In this review, we analyze the current findings in the mechanism of action of different PFPs and AMPs that support a wide role of membrane pore formation in nature. We also provide the newest insights into the development of state-of-art techniques that have facilitated the characterization of membrane pores. To understand the physiological role of these peptides/proteins or develop clinical applications, it is essential to uncover the molecular mechanism of how they perforate membranes.