Ganglioside-Lipid and Ganglioside-Protein Interactions Revealed by Coarse-Grained and Atomistic Molecular Dynamics Simulations

Viral fusion proteins of classes II and III recognize and reorganize complex biological membranes

Viral infection requires stable binding of viral fusion proteins to host membranes, which contain hundreds of lipid species. The mechanisms by which fusion proteins utilize specific host lipids to drive virus-host membrane fusion remains elusive. We conducted molecular simulations of classes I, II, and III fusion proteins interacting with membranes of diverse lipid compositions. Free energy calculations reveal that class I fusion proteins generally exhibit stronger membrane binding compared to classes II and III - a trend consistent across 74 fusion proteins from 13 viral families as suggested by sequence analysis. Class II fusion proteins utilize a lipid binding pocket formed by fusion protein monomers, stabilizing the initial binding of monomers to the host membrane prior to assembling into fusogenic trimers. In contrast, class III fusion proteins form a lipid binding pocket at the monomer-monomer interface through a unique fusion loop crossover. The distinct lipid binding modes correlate with the differing maturation pathways of classes II and III proteins. Binding affinity was predominantly controlled by cholesterol and gangliosides as well as via local enrichment of polyunsaturated lipids, thereby locally enhancing membrane disorder. Our study reveals energetics and atomic details underlying lipid recognition and reorganization by different viral fusion protein classes, offering insights into their specialized membrane fusion pathways.

Computational Methods for Modeling Lipid-Mediated Active Pharmaceutical Ingredient Delivery

Lipid-mediated delivery of active pharmaceutical ingredients (API) opened new possibilities in advanced therapies. By encapsulating an API into a lipid nanocarrier (LNC), one can safely deliver APIs not soluble in water, those with otherwise strong adverse effects, or very fragile ones such as nucleic acids. However, for the rational design of LNCs, a detailed understanding of the composition-structure-function relationships is missing. This review presents currently available computational methods for LNC investigation, screening, and design. The state-of-the-art physics-based approaches are described, with the focus on molecular dynamics simulations in all-atom and coarse-grained resolution. Their strengths and weaknesses are discussed, highlighting the aspects necessary for obtaining reliable results in the simulations. Furthermore, a machine learning, i.e., data-based learning, approach to the design of lipid-mediated API delivery is introduced. The data produced by the experimental and theoretical approaches provide valuable insights. Processing these data can help optimize the design of LNCs for better performance. In the final section of this Review, state-of-the-art of computer simulations of LNCs are reviewed, specifically addressing the compatibility of experimental and computational insights.

Gangliosides and Cholesterol: Dual Regulators of Neuronal Membrane Framework in Autism Spectrum Disorder

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with heterogeneous clinical presentation. Diagnosing ASD is complex, and the criteria for diagnosis, as well as the term ASD, have changed during the last decades. Diagnosis is made based on observation and accomplishment of specific diagnostic criteria, while a particular biomarker of ASD does not yet exist. However, studies universally report a disequilibrium in membrane lipid content, pointing to a unique neurolipid signature of ASD. This review sheds light on the possible role of cholesterol and gangliosides, complex membrane glycosphingolipids, in the development of ASD. In addition to maintaining membrane integrity, neuronal signaling, and synaptic plasticity, these lipids play a role in neurotransmitter release and calcium signaling. Evidence linking ASD to lipidome changes includes low cholesterol levels, unusual ganglioside levels, and unique metabolic profiles. ASD symptoms may be mitigated with therapeutic interventions targeting the lipid composition of membranes. However, restoring membrane equilibrium in the central nervous system remains a challenge. This review underscores the need for comprehensive research into lipid metabolism to uncover practical insights into ASD etiology and treatment as lipidomics emerges as a major area in ASD research.

Assessing the Martini 3 protein model: A review of its path and potential

Coarse-grained (CG) protein models have become indispensable tools for studying many biological protein details, from conformational dynamics to the organization of protein macro-complexes, and even the interaction of proteins with other molecules. The Martini force field is one of the most widely used CG models for bio-molecular simulations, partly because of the enormous success of its protein model. With the recent release of a new and improved version of the Martini force field - Martini 3 - a new iteration of its protein model was also made available. The Martini 3 protein force field is an evolution of its Martini 2 counterpart, aimed at improving many of the shortcomings that had been previously identified. In this mini-review, we first provide a general overview of the model and then focus on the successful advances made in the short time since its release, many of which would not have been possible before. Furthermore, we discuss reported limitations, potential directions for model improvement and comment on what the likely future development and application avenues are.

Effect of GM1 concentration change on plasma membrane: molecular dynamics simulation and analysis

Ganglioside GM1 is a class of glycolipids predominantly located in the nervous system. Comprising a ceramide anchor and an oligosaccharide chain containing sialic acid, GM1 plays a pivotal role in various cellular processes, including signal transduction, cell adhesion, and membrane organization. Moreover, GM1 has been implicated in the pathogenesis of several neurological disorders, such as Parkinson's disease, Alzheimer's disease, and stroke. In this study, by creating a neural cell model membrane simulation system and employing rigorous molecular models, we utilize a coarse-grained molecular dynamics approach to explore the structural and dynamic characteristics of multi-component neuronal plasma membranes at varying GM1 ganglioside concentrations. The simulation results reveal that as GM1 concentration increases, a greater number of hydrogen bonds form between GM1 molecules, resulting in the formation of larger clusters, which leads to reduced membrane fluidity, increased lipid ordering, decreased membrane thickness and surface area and higher levels of GM1 dissociation. Through a meticulous analysis, while considering GM1's structural attributes, we offer valuable insights into the structural and dynamic traits of the cell membrane. This study provides a robust methodology for exploring membrane characteristics and enhances our comprehension of GM1 molecules, serving as a resource for both experimental and computational researchers in this field.

Improved Protein Dynamics and Hydration in the Martini3 Coarse-Grain Model

The Martini coarse-grain force-field has emerged as an important framework to probe cellular processes at experimentally relevant time- and length-scales. However, the recently developed version, the Martini3 force-field with the implemented Go̅ model (Martini3Go̅), as well as previous variants of the Martini model have not been benchmarked and rigorously tested for globular proteins. In this study, we consider three globular proteins, ubiquitin, lysozyme, and cofilin, and compare protein dynamics and hydration with observables from experiments and all-atom simulations. We show that the Martini3Go̅ model is able to accurately model the structural and dynamic features of small globular proteins. Overall, the structural integrity of the proteins is maintained, as validated by contact maps, radii of gyration (Rg), and SAXS profiles. The chemical shifts predicted from the ensemble sampled in the simulations are consistent with the experimental data. Further, a good match is observed in the protein-water interaction energetics, and the hydration levels of the residues are similar to atomistic simulations. However, the protein-water interaction dynamics is not accurately represented and appears to depend on the protein structural complexity, residue specificity, and water dynamics. Our work is a step toward testing and assessing the Martini3Go̅ model and provides insights into future efforts to refine Martini models with improved solvation effects and better correspondence to the underlying all-atom systems.

Advances in Mass Spectrometry of Gangliosides Expressed in Brain Cancers

Gangliosides are highly abundant in the human brain where they are involved in major biological events. In brain cancers, alterations of ganglioside pattern occur, some of which being correlated with neoplastic transformation, while others with tumor proliferation. Of all techniques, mass spectrometry (MS) has proven to be one of the most effective in gangliosidomics, due to its ability to characterize heterogeneous mixtures and discover species with biomarker value. This review highlights the most significant achievements of MS in the analysis of gangliosides in human brain cancers. The first part presents the latest state of MS development in the discovery of ganglioside markers in primary brain tumors, with a particular emphasis on the ion mobility separation (IMS) MS and its contribution to the elucidation of the gangliosidome associated with aggressive tumors. The second part is focused on MS of gangliosides in brain metastases, highlighting the ability of matrix-assisted laser desorption/ionization (MALDI)-MS, microfluidics-MS and tandem MS to decipher and structurally characterize species involved in the metastatic process. In the end, several conclusions and perspectives are presented, among which the need for development of reliable software and a user-friendly structural database as a search platform in brain tumor diagnostics.

Dynamic lipid interactions in the plasma membrane Na+,K+-ATPase

The function of ion-transporting Na+,K+-ATPases depends on the surrounding lipid environment in biological membranes. Two established lipid-interaction sites A and B within the transmembrane domain have been observed to induce protein activation and stabilization, respectively. In addition, lipid-mediated inhibition has been assigned to a site C, but with the exact location not experimentally confirmed. Also, possible effects on lipid interactions by disease mutants dwelling in the membrane-protein interface remain relatively uncharacterized. We simulated human Na+,K+-ATPase α1β1FXYD homology models in E1 and E2 states in an asymmetric, multicomponent plasma membrane to determine both wild-type and disease mutant lipid-protein interactions. The simulated wild-type lipid interactions at the established sites A and B were in agreement with experimental results thereby confirming the membrane-protein model system. The less well-characterized, proposed inhibitory site C was dominated by lipids lacking inhibitory properties. Instead, two sites hosting inhibitory lipids were identified at the extracellular side and also a cytoplasmic CHL-binding site that provide putative alternative locations of Na+,K+-ATPase inhibition. Three disease mutations, Leu302Arg, Glu840Arg and Met859Arg resided in the lipid-protein interface and caused drastic changes in the lipid interactions. The simulation results show that lipid interactions to the human Na+,K+-ATPase α1β1FXYD protein in the plasma membrane are highly state-dependent and can be disturbed by disease mutations located in the lipid interface, which can open up for new venues to understand genetic disorders.

Which Moiety Drives Gangliosides to Form Nanodomains?

Gangliosides are important glycosphingolipids involved in a multitude of physiological functions. From a physicochemical standpoint, this is related to their ability to self-organize into nanoscopic domains, even at molar concentrations of one per 1000 lipid molecules. Despite recent experimental and theoretical efforts suggesting that a hydrogen bonding network is crucial for nanodomain stability, the specific ganglioside moiety decisive for the development of these nanodomains has not yet been identified. Here, we combine an experimental technique achieving nanometer resolution (Förster resonance energy transfer analyzed by Monte Carlo simulations) with atomistic molecular dynamic simulations to demonstrate that the sialic acid (Sia) residue(s) at the oligosaccharide headgroup dominates the hydrogen bonding network between gangliosides, driving the formation of nanodomains even in the absence of cholesterol or sphingomyelin. Consequently, the clustering pattern of asialoGM₁, a Sia-depleted glycosphingolipid bearing three glyco moieties, is more similar to that of structurally distant sphingomyelin than that of the closely related gangliosides GM₁ and GD_1a with one and two Sia groups, respectively.

Blood pH Analysis in Combination with Molecular Medical Tools in Relation to COVID-19 Symptoms

The global outbreak of SARS-CoV-2/COVID-19 provided the stage to accumulate an enormous biomedical data set and an opportunity as well as a challenge to test new concepts and strategies to combat the pandemic. New research and molecular medical protocols may be deployed in different scientific fields, e.g., glycobiology, nanopharmacology, or nanomedicine. We correlated clinical biomedical data derived from patients in intensive care units with structural biology and biophysical data from NMR and/or CAMM (computer-aided molecular modeling). Consequently, new diagnostic and therapeutic approaches against SARS-CoV-2 were evaluated. Specifically, we tested the suitability of incretin mimetics with one or two pH-sensitive amino acid residues as potential drugs to prevent or cure long-COVID symptoms. Blood pH values in correlation with temperature alterations in patient bodies were of clinical importance. The effects of biophysical parameters such as temperature and pH value variation in relation to physical-chemical membrane properties (e.g., glycosylation state, affinity of certain amino acid sequences to sialic acids as well as other carbohydrate residues and lipid structures) provided helpful hints in identifying a potential Achilles heel against long COVID. In silico CAMM methods and in vitro NMR experiments (including ³¹P NMR measurements) were applied to analyze the structural behavior of incretin mimetics and SARS-CoV fusion peptides interacting with dodecylphosphocholine (DPC) micelles. These supramolecular complexes were analyzed under physiological conditions by ¹H and ³¹P NMR techniques. We were able to observe characteristic interaction states of incretin mimetics, SARS-CoV fusion peptides and DPC membranes. Novel interaction profiles (indicated, e.g., by ³¹P NMR signal splitting) were detected. Furthermore, we evaluated GM1 gangliosides and sialic acid-coated silica nanoparticles in complex with DPC micelles in order to create a simple virus host cell membrane model. This is a first step in exploring the structure-function relationship between the SARS-CoV-2 spike protein and incretin mimetics with conserved pH-sensitive histidine residues in their carbohydrate recognition domains as found in galectins. The applied methods were effective in identifying peptide sequences as well as certain carbohydrate moieties with the potential to protect the blood-brain barrier (BBB). These clinically relevant observations on low blood pH values in fatal COVID-19 cases open routes for new therapeutic approaches, especially against long-COVID symptoms.

Siglec-6 mediates the uptake of extracellular vesicles through a noncanonical glycolipid binding pocket

Immunomodulatory Siglecs are controlled by their glycoprotein and glycolipid ligands. Siglec-glycolipid interactions are often studied outside the context of a lipid bilayer, missing the complex behaviors of glycolipids in a membrane. Through optimizing a liposomal formulation to dissect Siglec-glycolipid interactions, it is shown that Siglec-6 can recognize glycolipids independent of its canonical binding pocket, suggesting that Siglec-6 possesses a secondary binding pocket tailored for recognizing glycolipids in a bilayer. A panel of synthetic neoglycolipids is used to probe the specificity of this glycolipid binding pocket on Siglec-6, leading to the development of a neoglycolipid with higher avidity for Siglec-6 compared to natural glycolipids. This neoglycolipid facilitates the delivery of liposomes to Siglec-6 on human mast cells, memory B-cells and placental syncytiotrophoblasts. A physiological relevance for glycolipid recognition by Siglec-6 is revealed for the binding and internalization of extracellular vesicles. These results demonstrate a unique and physiologically relevant ability of Siglec-6 to recognize glycolipids in a membrane.

Identification of possible binding modes of SARS-CoV-2 spike N-terminal domain for ganglioside GM1

Coarse-grained molecular dynamics simulations of the lipid bilayer mixture of POPC and cholesterol were carried out in the presence and absence of ganglioside monosialo 1 (GM1) with N - terminal domain (NTD) of SARS-CoV-2 spike glycoprotein. The interactions of GM1 with two different NTD orientations were compared. NTD orientation I compactly bind GM1 predominantly through the sialic acid and the external galactose moieties providing more restriction to GM1 mobility whereas orientation II is more distributed on the lipid surface and due to the relaxed mobility of GM1 there, presumably, the NTD receptor penetrates more through the membrane.

Martini 3 Coarse-Grained Force Field for Carbohydrates

The Martini 3 force field is a full reparametrization of the Martini coarse-grained model for biomolecular simulations. Due to the improved interaction balance, it allows for a more accurate description of condensed phase systems. In the present work, we develop a consistent strategy to parametrize carbohydrate molecules accurately within the framework of Martini 3. In particular, we develop a canonical mapping scheme which decomposes arbitrarily large carbohydrates into a limited number of fragments. Bead types for these fragments have been assigned by matching physicochemical properties of mono- and disaccharides. In addition, guidelines for assigning bonds, angles, and dihedrals were developed. These guidelines enable a more accurate description of carbohydrate conformations than in the Martini 2 force field. We show that models obtained with this approach are able to accurately reproduce osmotic pressures of carbohydrate water solutions. Furthermore, we provide evidence that the model differentiates correctly the solubility of the polyglucoses dextran (water-soluble) and cellulose (water insoluble but soluble in ionic liquids). Finally, we demonstrate that the new building blocks can be applied to glycolipids. We show they are able to reproduce membrane properties and induce binding of peripheral membrane proteins. These test cases demonstrate the validity and transferability of our approach.

Molecular dynamics simulations and solid-state nuclear magnetic resonance spectroscopy measurements of C-H bond order parameters and effective correlation times in a POPC-GM3 bilayer

Glycolipids such as gangliosides affect the properties of lipid membranes and in extension the interactions between membranes and other biomolecules like proteins. To better understand how the properties of individual lipid molecules can contribute to shape the functional aspects of a membrane, the spatial restriction and dynamics of C-H bond segments can be measured using nuclear magnetic resonance (NMR) spectroscopy. We combine solid-state NMR spectroscopy with all-atom molecular dynamics (MD) simulations to investigate how ganglioside GM3 affects the bilayer structure and dynamics of C-H bond segments. These two methods yield reorientational correlation functions, molecular profiles of C-H bond order parameters |S_CH| and effective correlation times τ_e, which we compare for lipids in POPC bilayers with and without 30 mol% GM3. Our results revealed that all C-H segments of POPC reorient slower in the presence of GM3 and that the defining features of the GM3-POPC bilayer lie in the GM3 headgroup; it gives the bilayer an extended headgroup layer with high order (|S_CH| up to 0.3-0.4) and slow dynamics (τ_e up to 100 ns), a character that may be mechanistically important in ganglioside interactions with other biomolecules.

A cell free biomembrane platform for multimodal study of influenza virus hemagglutinin and for evaluation of entry-inhibitors against hemagglutinin

Seasonal periodic pandemics and epidemics caused by Influenza A viruses (IAVs) are associated with high morbidity and mortality worldwide. They are frequent and unpredictable in severity so there is a need for biophysical platforms that can be used to provide both mechanistic insights into influenza virulence and its potential treatment by anti-IAV agents. Host membrane viral association through the glycoprotein hemagglutinin (HA) of IAVs is one of the primary steps in infection. HA is thus a potential target for drug discovery and development against influenza. Deconvolution of the multivalent interactions of HA at the interfaces of the host cell membrane can help unravel therapeutic targets. In this contribution, we reported the effect of a multivalent HA glycoprotein association on various glycosphingolipid receptors (GD1a, GM3, GM1) doped asymmetrically into an artificial host membrane spanned across an aqueous filled microcavity array. The extent of HA association and its impact on membrane resistance, capacitance, and diffusivity was measured using highly sensitive electrochemical impedance spectroscopy (EIS) and fluorescence lifetime correlation spectroscopy (FLCS). Furthermore, we investigated the inhibition of the influenza HA glycoprotein association with the host mimetic surface by natural and synthetic sialic acid-based inhibitors (sialic acid, Siaα2,3-GalOMe, FB127, 3-sialyl lactose) using electrochemical impedance spectroscopy and observe that while all inhibit, they do not prevent host binding. Overall, the work demonstrates the platform provides a label-free screening platform for the biophysical evaluation of new inhibitors in the development of potential therapeutics for IAV infection prevention and treatment.

Recombination between 13C and 2H to Form Acetylide (13C22H-) Probes Nanoscale Interactions in Lipid Bilayers via Dynamic Secondary Ion Mass Spectrometry: Cholesterol and GM1 Clustering

Although it is thought that there is lateral heterogeneity of lipid and protein components within biological membranes, probing this heterogeneity has proven challenging. The difficulty in such experiments is due to both the small length scale over which such heterogeneity can occur, and the significant perturbation resulting from fluorescent or spin labeling on the delicate interactions within bilayers. Atomic recombination during dynamic nanoscale secondary ion imaging mass spectrometry (NanoSIMS) is a non-perturbative method for examining nanoscale bilayer interactions. Atomic recombination is a variation on conventional NanoSIMS imaging, whereby an isotope on one molecule combines with a different isotope on another molecule during the ionization process, forming an isotopically enriched polyatomic ion in a distance-dependent manner. We show that the recombinant ion, ¹³C₂²H^-, is formed in high yield from ¹³C- and ²H-labeled lipids. The low natural abundance of triply labeled acetylide also makes it an ideal ion to probe GM₁ clusters in model membranes and the effects of cholesterol on lipid-lipid interactions. We find evidence supporting the cholesterol condensation effect as well as the presence of nanoscale GM₁ clusters in model membranes.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Chemical control over Asialo-GM1: A dual ligand for pili and Lectin A that activates swarming motility and facilitates adherence of Pseudomonas aeruginosa

Glycolipid, ganglio-N-tetraosylceramide (asialo-GM1), on the mammalian cells are known to be recognized by type IV pili of Pseudomonas aeruginosa. In this work, we show that asialo-GM1 can also be recognized by Lectin A (LecA), another adhesin protein of the P. aeruginosa, by a fluorescent polarization assay, a label-free bacterial motility enabled binding assay, and bacterial mutant studies. On hydrated semi-solid gel surfaces, asialo-GM1 enables swarming and twitching motilities, while on solid surfaces facilitates the bacterial adherence of P. aeruginosa. These results indicate that asialo-GM1 can modulate bioactivities, adherence, and motilities, that are controlled by opposite signaling pathways. We demonstrate that when a solution of pilin monomers or LecA proteins are spread on hydrated gel surfaces, the asialo-GM1 mediated swarming motility is inhibited. Treatment of artificial liposomes containing asialo-GM1 as a component of lipid bilayer with pilin monomers or LecA proteins caused transient leakage of encapsulated dye from liposomes. These results suggest that pili and LecA proteins not only bind to asialo-GM1 but can also cause asialo-GM1 mediated leakage. We also show that both pili and LecA mutants of P. aeruginosa adhere to asialo-GM1 coated solid surfaces, and that a class of synthetic ligands for pili and LecA inhibits both pili and LecA-mediated adherence of P. aeruginosa on asialo-GM1-coated surfaces.

Start Me Up: How Can Surrounding Gangliosides Affect Sodium-Potassium ATPase Activity and Steer towards Pathological Ion Imbalance in Neurons?

Gangliosides, amphiphilic glycosphingolipids, tend to associate laterally with other membrane constituents and undergo extensive interactions with membrane proteins in cis or trans configurations. Studies of human diseases resulting from mutations in the ganglioside biosynthesis pathway and research on transgenic mice with the same mutations implicate gangliosides in the pathogenesis of epilepsy. Gangliosides are reported to affect the activity of the Na⁺/K⁺-ATPase, the ubiquitously expressed plasma membrane pump responsible for the stabilization of the resting membrane potential by hyperpolarization, firing up the action potential and ion homeostasis. Impaired Na⁺/K⁺-ATPase activity has also been hypothesized to cause seizures by several mechanisms. In this review we present different epileptic phenotypes that are caused by impaired activity of Na⁺/K⁺-ATPase or changed membrane ganglioside composition. We further discuss how gangliosides may influence Na⁺/K⁺-ATPase activity by acting as lipid sorting machinery providing the optimal stage for Na⁺/K⁺-ATPase function. By establishing a distinct lipid environment, together with other membrane lipids, gangliosides possibly modulate Na⁺/K⁺-ATPase activity and aid in “starting up” and “turning off” this vital pump. Therefore, structural changes of neuronal membranes caused by altered ganglioside composition can be a contributing factor leading to aberrant Na⁺/K⁺-ATPase activity and ion imbalance priming neurons for pathological firing.

Experimental and computational characterization of dynamic biomolecular interaction systems involving glycolipid glycans

On cell surfaces, carbohydrate chains that modify proteins and lipids mediate various biological functions, which are exerted not only through carbohydrate-protein interactions but also through carbohydrate-carbohydrate interactions. These glycans exhibit considerable degrees of conformational variability and often form clusters providing multiple binding sites. The integration of nuclear magnetic resonance spectroscopy and molecular dynamics simulation has made it possible to delineate the dynamical structures of carbohydrate chains. This approach has facilitated the remodeling of oligosaccharide conformational space in the prebound state to improve protein-binding affinity and has been applied to visualize dynamic carbohydrate-carbohydrate interactions that control glycoprotein-glycoprotein complex formation. Functional glycoclusters have been characterized by experimental and computational approaches applied to various model membranes and artificial self-assembling systems. This line of investigation has provided dynamic views of molecular assembling on glycoclusters, giving mechanistic insights into physiological and pathological molecular events on cell surfaces as well as clues for the design and creation of molecular systems exerting improved glycofunctions. Further development and accumulation of such studies will allow detailed understanding and artificial control of the "glycosynapse" foreseen by Dr. Sen-itiroh Hakomori.

Oligosaccharide Presentation Modulates the Molecular Recognition of Glycolipids by Galectins on Membrane Surfaces

Galectins are a family of glycan binding proteins that stand out for the wide range of biological phenomena in which they are involved. Most galectin functions are associated with their glycan binding capacities, which are generally well characterized at the oligosaccharide level, but not at the glycoprotein or glycolipid level. Glycolipids form the part of cell membranes where they can act as galectin cellular receptors. In this scenario, glycan presentation as well as the membrane chemical and structural features are expected to have a strong impact in these molecular association processes. Herein, liposomes were used as membrane mimicking scaffolds for the presentation of glycosphingolipids (GSLs) and to investigate their interaction with Galectin-3 and the N-domain of Galectin-8 (Gal8N). The binding towards GM3 and GM1 and their non-silaylated GSLs was compared to the binding to the free glycans, devoid of lipid. The analysis was carried out using a combination of NMR methods, membrane perturbation studies, and molecular modeling. Our results showed a different tendency of the two galectins in their binding capacities towards the glycans, depending on whether they were free oligosaccharides or as part of GSL inserted into a lipid bilayer, highlighting the significance of GSL glycan presentation on membranes in lectin binding.

Aquaporin ion conductance properties defined by membrane environment, protein structure, and cell physiology

Aquaporins (AQPs) are multifunctional transmembrane channel proteins permeable to water and an expanding array of solutes. AQP-mediated ion channel activity was first observed when purified AQP0 from bovine lens was incorporated into lipid bilayers. Electrophysiological properties of ion-conducting AQPs since discovered in plants, invertebrates, and mammals have been assessed using native, reconstituted, and heterologously expressed channels. Accumulating evidence is defining amino acid residues that govern differential solute permeability through intrasubunit and central pores of AQP tetramers. Rings of charged and hydrophobic residues around pores influence AQP selectivity, and are candidates for further work to define motifs that distinguish ion conduction capability, versus strict water and glycerol permeability. Similarities between AQP ion channels thus far include large single channel conductances and long open times, but differences in ionic selectivity, permeability to divalent cations, and mechanisms of gating (e.g., by voltage, pH, and cyclic nucleotides) are unique to subtypes. Effects of lipid environments in modulating parameters such as single channel amplitude could explain in part the variations in AQP ion channel properties observed across preparations. Physiological roles of the ion-conducting AQP classes span diverse processes including regulation of cell motility, organellar pH, neural development, signaling, and nutrient acquisition. Advances in computational methods can generate testable predictions of AQP structure-function relationships, which combined with innovative high-throughput assays could revolutionize the field in defining essential properties of ion-conducting AQPs, discovering new AQP ion channels, and understanding the effects of AQP interactions with proteins, signaling cascades, and membrane lipids.

Developing initial conditions for simulations of asymmetric membranes: a practical recommendation

It has been proposed that the surface tension difference between leaflets (or differential stress) in asymmetric bilayers is generally nonvanishing. This implies that there is no unique approach to generate initial conditions for simulations of asymmetric bilayers in the absence of experimentally derived constraints. Current generation methods include individual area per lipid (APL) based, leaflet surface area (SA) matching, and zero leaflet tension based (0-DS). This work adds a bilayer-based approach that aims for achieving partial chemical equilibrium by interleaflet switching of selected lipids via P2₁ periodic boundary conditions. Based on a recently proposed theoretical framework, we obtained expressions for tensions in asymmetric bilayers from both the bending and area strains. We also developed a quantitative measure for the energetic penalty from the differential stress. The impacts of APL-, SA-, and 0-DS-based approaches on mechanical properties are assessed for two different asymmetric bilayers. The lateral pressure profile and its moments differ significantly for each method, whereas the area compressibility modulus is relatively insensitive. Application of P2₁ periodic boundary conditions (APL/P2₁, SA/P2₁, and 0-DS/P2₁) results in better agreement in mechanical properties between asymmetric bilayers generated by APL-, SA-, and 0-DS-based approaches, in which changes are the smallest for bilayers from the SA-based method. The estimated differential stress from the theory shows good agreement with that from the simulations. These simulation results and the good agreement between the predicted and observed differential stress further support the theoretical framework in which bilayer mechanical properties are outcomes of the interplay between intrinsic bending and asymmetric lipid packing. Based on the simulation results and theoretical predictions, the SA/P2₁-based, or at least the SA-based (when the differential stress is small), approach is recommended as a practical method for developing initial conditions for asymmetric bilayer simulations.

Fusion pores with low conductance are cation selective

Many neurotransmitters are organic ions that carry a net charge, and their release from secretory vesicles is therefore an electrodiffusion process. The selectivity of early exocytotic fusion pores is investigated by combining electrodiffusion theory, measurements of amperometric foot signals from chromaffin cells with anion substitution, and molecular dynamics simulation. The results reveal that very narrow fusion pores are cation selective, but more dilated fusion pores become anion permeable. The transition occurs around a fusion pore conductance of ~300 pS. The cation selectivity of a narrow fusion pore accelerates the release of positively charged transmitters such as dopamine, noradrenaline, adrenaline, serotonin, and acetylcholine, while glutamate release may require a more dilated fusion pore.

For transmission, a fusion pore forms when vesicle and target membranes are brought together by SNARE proteins. Delacruz et al. demonstrate that selectivity of the pore accelerates release of positively charged transmitters such as dopamine, noradrenaline, adrenaline, serotonin, and acetylcholine, while glutamate release may require a more dilated fusion pore.

The role of plasmalogens, Forssman lipids, and sphingolipid hydroxylation in modulating the biophysical properties of the epithelial plasma membrane

A coarse-grain model of the epithelial plasma membrane was developed from high-resolution lipidomic data and simulated using the MARTINI force field to characterize its biophysical properties. Plasmalogen lipids, Forssman glycosphingolipids, and hydroxylated Forssman glycosphingolipids and sphingomyelin were systematically added to determine their structural effects. Plasmalogen lipids have a minimal effect on the overall biophysical properties of the epithelial plasma membrane. In line with the hypothesized role of Forssman lipids in the epithelial apical membrane, the introduction of Forssman lipids initiates the formation of glycosphingolipid-rich nanoscale lipid domains, which also include phosphatidylethanolamine (PE), sphingomyelin (SM), and cholesterol (CHOL). This decreases the lateral diffusion in the extracellular leaflet, as well as the area per lipid of domain forming lipids, most notably PE. Finally, hydroxylation of the Forssman glycosphingolipids and sphingomyelin further modulates the lateral organization of the membrane. Through comparison to the previously studied average and neuronal plasma membranes, the impact of membrane lipid composition on membrane properties was characterized. Overall, this study furthers our understanding of the biophysical properties of complex membranes and the impact of lipid diversity in modulating membrane properties.

Multiscale Molecular Dynamics Studies Reveal Different Modes of Receptor Clustering by Gb3-Binding Lectins

The recognition of carbohydrate receptors on host cell membranes by pathogenic lectins is a crucial step in the microbial invasion. Two bacterial lectins, the B-subunit of Shiga toxin from Shigella dysenteria (StxB) and lectin I from Pseudomonas aeruginosa (LecA), are specific to the same galactolipid-globotriaosylceramide (Gb3). In this study we present a coarse-grained (cg) model of Gb3, which we further apply to unravel the molecular details of glycolipid binding by two lectins on the surface of a DOPC/cholesterol/Gb3 bilayer. In cg molecular dynamics simulations with time scales of dozens of microseconds, Gb3 was randomly distributed. The binding of both StxB or LecA is accompanied by Gb3 clustering in a cholesterol environment and with exclusion of DOPC in protein vicinity. StxB being bound by all 15 binding sites induced membrane bending, while LecA interacted with two out of four binding sites for most of the time causing a smaller inward curvature of the model membrane. Stable interactions occurred preferably when LecA was normal to the membrane surface. Furthermore, all-atom simulations revealed that LecA bound Gb3's headgroup at only one out of two possible conformations of the carbohydrate moiety observed at protein-free conditions. The results shed light on the mechanism of interactions between two lectins and Gb3 on the membrane surface and offer a coarse-grained model to study more complex systems at large spatiotemporal scales.

Out of Sight, Out of Mind: The Effect of the Equilibration Protocol on the Structural Ensembles of Charged Glycolipid Bilayers

Molecular dynamics (MD) simulations represent an essential tool in the toolbox of modern chemistry, enabling the prediction of experimental observables for a variety of chemical systems and processes and majorly impacting the study of biological membranes. However, the chemical diversity of complex lipids beyond phospholipids brings new challenges to well-established protocols used in MD simulations of soft matter and requires continuous assessment to ensure simulation reproducibility and minimize unphysical behavior. Lipopolysaccharides (LPS) are highly charged glycolipids whose aggregation in a lamellar arrangement requires the binding of numerous cations to oppositely charged groups deep inside the membrane. The delicate balance between the fully hydrated carbohydrate region and the smaller hydrophobic core makes LPS membranes very sensitive to the choice of equilibration protocol. In this work, we show that the protocol successfully used to equilibrate phospholipid bilayers when applied to complex lipopolysaccharide membranes occasionally leads to a small expansion of the simulation box very early in the equilibration phase. Although the use of a barostat algorithm controls the system dimension and particle distances according to the target pressure, fluctuation in the fleeting pressure occasionally enables a few water molecules to trickle into the hydrophobic region of the membrane, with spurious solvent buildup. We show that this effect stems from the initial steps of NPT equilibration, where initial pressure can be fairly high. This can be solved with the use of a stepwise-thermalization NVT/NPT protocol, as demonstrated for atomistic MD simulations of LPS/DPPE and lipid-A membranes in the presence of different salts using an extension of the GROMOS forcefield within the GROMACS software. This equilibration protocol should be standard procedure for the generation of consistent structural ensembles of charged glycolipids starting from atomic coordinates not previously pre-equilibrated. Although different ways to deal with this issue can be envisioned, we investigated one alternative that could be readily available in major MD engines with general users in mind.

Neoglycolipids as Glycosphingolipid Surrogates for Protein Binding Studies Using Nanodiscs and Native Mass Spectrometry

Interactions between glycan-binding proteins (GBPs) and glycosphingolipids (GSLs) in the membranes of cells are implicated in a wide variety of normal and pathophysiological processes. Despite the critical biological roles these interactions play, the GSL ligands of most GBPs have not yet been identified. The limited availability of purified GSLs represents a significant challenge to the discovery and characterization of biologically relevant GBP-GSL interactions. The present work investigates the use of neoglycolipids (NGLs) as surrogates for GSLs for catch-and-release-electrospray ionization mass spectrometry (CaR-ESI-MS)-based screening, implemented with nanodiscs, for the discovery of GSL ligands. Three pairs of NGLs based on the blood group type A and B trisaccharides, with three different lipid head groups but all with "ring-closed" monosaccharide residue at the reducing end, were synthesized. The incorporation efficiencies (into nanodiscs) of the NGLs and their affinities for a fragment of family 51 carbohydrate-binding module (CBM) identified an amide-linked 1,3-di-O-hexadecyl-glycerol moiety as the optimal lipid structure. Binding measurements performed on cholera toxin B subunit homopentamer (CTB₅) and nanodiscs containing an NGL consisting of the optimal lipid moiety and the GM1 ganglioside pentasaccharide yielded affinities similar, within a factor of 2, to those of native GM1. Finally, nanodiscs containing the optimal A and B trisaccharide NGLs, as well as the corresponding NGLs of lactose, A type 2 tetrasaccharide, and the GM1 and GD2 pentasaccharides were screened against the family 51 CBM, human galectin-7, and CTB₅ to illustrate the potential of NGLs to accelerate the discovery of GSL ligands of GBPs.

Amphiphilic gold nanoparticles perturb phase separation in multidomain lipid membranes

Amphiphilic gold nanoparticles with diameters in the 2-4 nm range are promising as theranostic agents thanks to their spontaneous translocation through cell membranes. This study addresses the effects that these nanoparticles may have on a distinct feature of plasma membranes: lipid lateral phase separation. Atomic force microscopy, quartz crystal microbalance, and molecular dynamics are combined to study the interaction between model neuronal membranes, which spontaneously form ordered and disordered lipid domains, and amphiphilic gold nanoparticles having negatively charged surface functionalization. Nanoparticles are found to interact with the bilayer and form bilayer-embedded ordered aggregates. Nanoparticles also suppress lipid phase separation, in a concentration-dependent fashion. A general, yet simple thermodynamic model is developed to show that the change of lipid-lipid enthalpy is the dominant driving force towards the nanoparticle-induced destabilization of phase separation.

Understanding the Link between Lipid Diversity and the Biophysical Properties of the Neuronal Plasma Membrane

Cell membranes contain incredible diversity in the chemical structures of their individual lipid species and the ratios in which these lipids are combined to make membranes. Nevertheless, our current understanding of how each of these components affects the properties of the cell membrane remains elusive, in part due to the difficulties in studying the dynamics of membranes at high spatiotemporal resolution. In this work, we use coarse-grained molecular dynamics simulations to investigate how individual lipid species contribute to the biophysical properties of the neuronal plasma membrane. We progress through eight membranes of increasing chemical complexity, ranging from a simple POPC/CHOL membrane to a previously published neuronal plasma membrane [Ingólfsson, H. I., et al. (2017) Biophys. J. 113 (10), 2271-2280] containing 49 distinct lipid species. Our results show how subtle chemical changes can affect the properties of the membrane and highlight the lipid species that give the neuronal plasma membrane its unique biophysical properties. This work has potential far-reaching implications for furthering our understanding of cell membranes.

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.

A Trifunctional Linker for Palmitoylation and Peptide and Protein Localization in Biological Membranes

Attachment of lipophilic groups is an important post-translational modification of proteins, which involves the coupling of one or more anchors such as fatty acids, isoprenoids, phospholipids, or glycosylphosphatidyl inositols. To study its impact on the membrane partitioning of hydrophobic peptides or proteins, we designed a tyrosine-based trifunctional linker. The linker allows the facile incorporation of two different functionalities at a cysteine residue in a single step. We determined the effect of the lipid modification on the membrane partitioning of the synthetic α-helical model peptide WALP with or without here and in all cases below; palmitoyl groups in giant unilamellar vesicles that contain a liquid-ordered (Lo ) and liquid-disordered (Ld ) phase. Introduction of two palmitoyl groups did not alter the localization of the membrane peptides, nor did the membrane thickness or lipid composition. In all cases, the peptide was retained in the Ld phase. These data demonstrate that the Lo domain in model membranes is highly unfavorable for a single membrane-spanning peptide.

Lipid-specific interactions determine the organization and dynamics of membrane-active peptide melittin

The cell membranes of different cells deviate significantly in lipid compositions and thus provide varying biological environments to modulate the diffusion, organization and the resultant function of biomacromolecules. However, the detailed modulation mechanism remains elusive especially in consideration of the current overuse of the simplified membrane models such as the pure phosphatidylcholine (PC) membrane. In this work, with the typical membrane-active peptide melittin, we demonstrated that a more complicated membrane environment, such as the bacterial (IME) or plasma membrane (PM), would significantly change the organization and dynamics of melittin, by using molecular dynamics simulations as a "computational microscope". It was found that in these membrane systems, adding melittin would cause a varying degree of reduction in the lateral diffusion of lipids due to the different assembly states of peptides. Melittin tended to aggregate to oligomers in the pure PC membrane, mostly as a tetramer or trimer, while in IME or PM, its degree of oligomerization was significantly reduced. More surprisingly, melittin displayed a strong affinity with ganglioside GM3 in PM, leading to the formation of melittin-GM3 nanoclusters, which hindered its diffusion and further oligomerization. Additionally, small changes in the residue sequence of melittin could modulate the degree or structure of the peptide oligomer. Our work provides a typical example of a study on the organization and dynamics of pore-forming peptides in specific membrane environments and has great significance on the optimization of peptide sequences and the design of helix bundles in the membrane for target biological function.

Permeation of nanoparticles across the intestinal lipid membrane: dependence on shape and surface chemistry studied through molecular simulations

Nanoparticles are being explored for topical and oral drug delivery applications as they can cross various biological barriers, for example, the intestinal epithelium. The ability of nanoparticles to cross barriers depends on their morphological and surface properties such as size, surface chemistry and shape, among others. The effect of nanoparticle size on their membrane permeability has been well studied both experimentally and theoretically. However, less attention has been given to understand the role of nanoparticle shape in their translocation across biological barrier membranes. Here, we report on the influence of the nanoparticle's shape, surface chemistry and concentration on their permeation across a human intestinal apical cell membrane model. A representative multicomponent lipid bilayer model of the human intestinal apical membrane was built. The free energy of permeation of nanoparticles across the model lipid bilayer was calculated using multiple umbrella sampling simulations. The interaction of these nanoparticles with the model lipid bilayer was captured using extensive microsecond unrestrained molecular dynamics simulations. We observed that: (a) irrespective of the surface chemistry, the efficacy of nanoparticle penetration across the lipid layer was in the order of rod > disc > sphere; (b) irrespective of the shape, apolar and nonpolar nanoparticles were found to locate in the interior of the lipid bilayer, whereas charged and polar nanoparticles were either adsorbed on the lipid headgroups or remained in the water layer; (c) apolar and nonpolar disc shaped nanoparticles had higher efficacy in permeation across the lipid bilayer as compared to disc and sphere shaped nanoparticles; and (d) at a higher concentration of nanoparticles, sphere and disc shaped nanoparticles exhibited more agglomeration as compared to rod shaped nanoparticles. Based on these outcomes, a few nanoparticles were designed which penetrated readily into the lipid layer and these nanoparticles were also able to co-deliver a therapeutic protein inside the lipid layer. The apical model lipid membrane and protocols used in this study can thus be utilized for the in silico design of nanoparticles for the oral delivery of therapeutics.

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.

Microcavity array supported lipid bilayer models of ganglioside – influenza hemagglutinin1 binding

The binding of influenza receptor (HA₁) to membranes containing different glycosphingolipid receptors was investigated at Microcavity Supported Lipid Bilayers (MSLBs).

Gangliosides Destabilize Lipid Phase Separation in Multicomponent Membranes

Gangliosides (GMs) form an important class of lipids found in the outer leaflet of the plasma membrane. Typically, they colocalize with cholesterol and sphingomyelin in ordered membrane domains. However, detailed understanding of the lateral organization of GM-rich membranes is still lacking. To gain molecular insight, we performed molecular dynamics simulations of GMs in model membranes composed of coexisting liquid-ordered and liquid-disordered domains. We found that GMs indeed have a preference to partition into the ordered domains. At higher concentrations (>10 mol %), we observed a destabilizing effect of GMs on the phase coexistence. Further simulations with modified GMs show that the structure of the GM headgroup affects the phase separation, whereas the nature of the tail determines the preferential location. Together, our findings provide a molecular basis to understand the lateral organization of GM-rich membranes.

Force Field Comparison of GM1 in a DOPC Bilayer Validated with AFM and FRET Experiments

The great physiological relevance of glycolipids is being increasingly recognized, and glycolipid interactions have been shown to be central to cell-cell recognition, neuronal plasticity, protein-ligand recognition, and other important processes. However, detailed molecular-level understanding of these processes remains to be fully resolved. Molecular dynamics simulations could reveal the details of the glycolipid interactions, but the results may be influenced by the choice of the employed force field. Here, we have compared the behavior and properties of GM1, a common, biologically important glycolipid, using the CHARMM36, OPLS, GROMOS, and Amber99-GLYCAM06 (in bilayers comprising SLIPIDS and LIPID14 lipids) force fields in bilayers comprising 1,2-dioleoyl-sn-glycero-3-phosphocholine lipids and compared the results to atomic force microscopy and fluorescence resonance energy transfer experiments. We found discrepancies within the GM1 behavior displayed between the investigated force fields. Based on a direct comparison with complementary experimental results derived from fluorescence and AFM measurements, we recommend using the Amber99-GLYCAM force field in bilayers comprising LIPID14 or SLIPIDS lipids followed by CHARMM36 and OPLS force fields in simulations. The GROMOS force field is not recommended for reproducing the properties of the GM1 head group.

Bilayer Membranes with Frequent Flip-Flops Have Tensionless Leaflets

Biomembranes are built up from lipid bilayers with two leaflets that typically differ in their lipid composition. Each lipid molecule stays within one leaflet of the bilayer before it undergoes a transition, or flip-flop, to the other leaflet. The corresponding flip-flop times are very different for different lipid species and vary over several orders of magnitude. Here, we use molecular dynamics simulations to elucidate the consequences of this separation of time scales for compositionally asymmetric bilayers. We first study bilayers with two lipid components that do not undergo flip-flops on the accessible time scales. In such a situation, one must distinguish a bilayer state in which both leaflets have the same preferred area from another state in which each leaflet is tensionless. However, when we add a third lipid component that undergoes frequent flip-flops, the bilayer relaxes toward the state with tensionless leaflets, not to the state with equal preferred leaflet areas. Furthermore, we show that bilayers with compositional asymmetry acquire a significant spontaneous curvature even if both leaflets are tensionless. Our results can be extended to lipid bilayers with a large number of lipid components provided at least one of these components undergoes frequent flip-flops. For cellular membranes containing lipid pumps, the leaflet tensions also depend on the rates of protein-induced flip-flops.

Two distinct anionic phospholipid-dependent events involved in SecA-mediated protein translocation

Protein translocation across the bacterial cytoplasmic membrane is an essential process catalyzed by the Sec translocase, which in its minimal form consists of the protein-conducting channel SecYEG, and the motor ATPase SecA. SecA binds via its positively charged N-terminus to membranes containing anionic phospholipids, leading to a lipid-bound intermediate. This interaction induces a conformational change in SecA, resulting in a high-affinity association with SecYEG, which initiates protein translocation. Here, we examined the effect of anionic lipids on the SecA-SecYEG interaction in more detail, and discovered a second, yet unknown, anionic lipid-dependent event that stimulates protein translocation. Based on molecular dynamics simulations we identified an anionic lipid-enriched region in vicinity of the lateral gate of SecY. Here, the anionic lipid headgroup accesses the lateral gate, thereby stabilizing the pre-open state of the channel. The simulations suggest flip-flop movement of phospholipid along the lateral gate. Electrostatic contribution of the anionic phospholipids at the lateral gate may directly stabilize positively charged residues of the signal sequence of an incoming preprotein. Such a mechanism allows for the correct positioning of the entrant peptide, thereby providing a long-sought explanation for the role of anionic lipids in signal sequence folding during protein translocation.

Recent Advances in Coarse-Grained Models for Biomolecules and Their Applications

Molecular dynamics simulations have emerged as a powerful tool to study biological systems at varied length and timescales. The conventional all-atom molecular dynamics simulations are being used by the wider scientific community in routine to capture the conformational dynamics and local motions. In addition, recent developments in coarse-grained models have opened the way to study the macromolecular complexes for time scales up to milliseconds. In this review, we have discussed the principle, applicability and recent development in coarse-grained models for biological systems. The potential of coarse-grained simulation has been reviewed through state-of-the-art examples of protein folding and structure prediction, self-assembly of complexes, membrane systems and carbohydrates fiber models. The multiscale simulation approaches have also been discussed in the context of their emerging role in unravelling hierarchical level information of biosystems. We conclude this review with the future scope of coarse-grained simulations as a constantly evolving tool to capture the dynamics of biosystems.

Lipid Bilayer-like Mixed Self-Assembled Monolayers with Strong Mobility and Clustering-Dependent Lectin Affinity

Glycans at the surface of cellular membranes modulate biological activity via multivalent association with extracellular messengers. The lack of tuneable simplified models mimicking this dynamic environment complicates basic studies of these phenomena. We here present a series of mixed reversible self-assembled monolayers (rSAMs) that addresses this deficiency. Mixed rSAMs were prepared in water by simple immersion of a negatively charged surface in a mixture of sialic acid- and hydroxy-terminated benzamidine amphiphiles. Surface compositions derived from infrared reflection-absorption spectroscopy (IRAS) and film thickness information (atomic force microscopy, ellipsometry) suggest the latter to be statistically incorporated in the monolayer. These surfaces' affinity for the lectin hemagglutinin revealed a strong dependence of the affinity on the presentation, density, and mobility of the sialic acid ligands. Hence, a spacer length of 4 ethylene glycol and a surface density of 15% resulted in a dissociation constant K_d,multi of 1.3 × 10^-13 M, on par with the best di- or tri-saccharide-based binders reported to date, whereas a density of 20% demonstrated complete resistance to hemagglutinin binding. These results correlated with ligand mobility measured by fluorescence recovery after photobleaching which showed a dramatic drop in the same interval. The results have a direct bearing on biological cell surface multivalent recognition involving lipid bilayers and may guide the design of model surfaces and sensors for both fundamental and applied studies.

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.