Long-chain GM1 gangliosides alter transmembrane domain registration through interdigitation

Comparison and Possible Binding Orientations of SARS-CoV-2 Spike N-Terminal Domain for Gangliosides GM3 and GM1

SARS-CoV-2 spike glycoprotein is anchored by gangliosides. The sialic acid in the ganglioside headgroup is responsible for virus attachment and entry into host cells. We used coarse-grained (CG) molecular dynamics simulations to expand on our previous study of GM1 interaction with two different orientations of the SARS-CoV-2 S1 subunit N-terminal domain (NTD) and to confirm the role of sialic acid receptors in driving the viral receptor; GM3 was used as another ganglioside on the membrane. Because of the smaller headgroup, sialic acid is crucial in GM3 interactions, whereas GM1 interacts with NTD via both the sialic acid and external galactose. In line with our previous findings for NTD orientations in GM1 binding, we identified two orientations, "compact" and "distributed", comprising sugar receptor-interacting residues in GM3-embedded lipid bilayers. Gangliosides in closer proximity to the compact NTD orientation might cause relatively greater restrictions to penetrate the bilayer. However, the attachment of a distributed NTD orientation with more negative interaction energies appears to facilitate GM1/GM3 to move quickly across the membrane. Our findings likely shed some light on the orientations that the NTD receptor acquires during the early phases of interaction with GM1 and GM3 in a membrane environment.

Conformations of Three Types of Ultra-Long-Chain Fatty Acids in Multicomponent Lipid Bilayers

Ultra-long-chain fatty acids (ULCFAs) are biosynthesized in certain types of tissues, but their biological roles remain unknown. Here, we report how the conformation of ULCFAs depends on the length and unsaturated-bond ratio of the ultra-long chains and the composition of the host bilayer membrane using molecular dynamics simulations. The ultra-long chain of ULCFAs flips between the two leaflets and fluctuates among three conformations: elongated, L-shaped, and turned. Furthermore, we found that the saturated ultra-long chain exhibited an elongated conformation more frequently than the unsaturated chain. In addition, the truncation of the ultra-long chain at C26 had little effect on the remaining ULCFAs. ULCFAs respond to lipid-density differences in the two leaflets, and the ratio of the elongated and turned conformations changed to reduce this difference. However, in cholesterol-containing membranes, ULCFAs exhibit no density difference after the flip-flop of cholesterol removes the difference.

Computational Approaches to Explore Bacterial Toxin Entry into the Host Cell

Many bacteria secrete toxic protein complexes that modify and disrupt essential processes in the infected cell that can lead to cell death. To conduct their action, these toxins often need to cross the cell membrane and reach a specific substrate inside the cell. The investigation of these protein complexes is essential not only for understanding their biological functions but also for the rational design of targeted drug delivery vehicles that must navigate across the cell membrane to deliver their therapeutic payload. Despite the immense advances in experimental techniques, the investigations of the toxin entry mechanism have remained challenging. Computer simulations are robust complementary tools that allow for the exploration of biological processes in exceptional detail. In this review, we first highlight the strength of computational methods, with a special focus on all-atom molecular dynamics, coarse-grained, and mesoscopic models, for exploring different stages of the toxin protein entry mechanism. We then summarize recent developments that are significantly advancing our understanding, notably of the glycolipid–lectin (GL-Lect) endocytosis of bacterial Shiga and cholera toxins. The methods discussed here are also applicable to the design of membrane-penetrating nanoparticles and the study of the phenomenon of protein phase separation at the surface of the membrane. Finally, we discuss other likely routes for future development.

Glycolipids: Linchpins in the Organization and Function of Membrane Microdomains

Membrane microdomains, also called lipid rafts, are areas on membrane enriched in glycolipids, sphingolipids, and cholesterol. Although membrane microdomains are thought to play key roles in many cellular functions, their structures, properties, and biological functions remain obscure. Cellular membranes contain several types of glycoproteins, glycolipids, and other lipids, including cholesterol, glycerophospholipids, and sphingomyelin. Depending on their physicochemical properties, especially the characteristics of their glycolipids, various microdomains form on these cell membranes, providing structural or functional contextures thought to be essential for biological activities. For example, the plasma membranes of human neutrophils are enriched in lactosylceramide (LacCer) and phosphatidylglucoside (PtdGlc), each of which forms different membrane microdomains with different surrounding molecules and is involved in different functions of neutrophils. Specifically, LacCer forms Lyn-coupled lipid microdomains, which mediate neutrophil chemotaxis, phagocytosis, and superoxide generation, whereas PtdGlc-enriched microdomains mediate neutrophil differentiation and spontaneous apoptosis. However, the mechanisms by which these glycolipids form different nano/meso microdomains and mediate their specialized functions remain incompletely understood. This review describes current understanding of the roles of glycolipids and sphingolipids in their enriched contextures on cellular membranes, including their mechanisms of facilitation and regulation of intracellular signaling. This review also introduces new concepts about the roles of glycolipid and sphingolipid-dependent contextures in immunological functions.

Conformation of ultra-long-chain fatty acid in lipid bilayer: Molecular dynamics study

Ultra-long-chain fatty acids (ULCFAs) are biosynthesized in the restricted tissues such as retina, testis, and skin. The conformation of a single ULCFA, in which the sn-1 unsaturated chain has 32 carbons, in three types of phospholipid bilayers is studied by molecular dynamics simulations. It is found that the ultra-long tail of the ULCFA flips between two leaflets and fluctuates among an elongation into the opposite leaflet, lies between two leaflets, and turns back. As the number ratio of lipids in the opposite leaflet increases, the ratio of the elongated shape linearly decreases in all three cases. Thus, ULCFAs can sense the density differences between the two leaflets and respond to these changes.

Pivotal Role of Interdigitation in Interleaflet Interactions: Implications from Molecular Dynamics Simulations

The asymmetric lipid composition in plasma membranes within the inner leaflet is not typically suitable for domain formation. Thus elucidation of the likelihood of the formation or stability of a raft-like domain in the inner leaflet is necessary. Herein we investigated the phase behavior of asymmetric membranes using coarse-grained molecular dynamics simulations. The lipid leaflet comprising dioleoylphosphatidylcholine (DOPC) and cholesterol (Chol) does not typically show well-developed domains in symmetric bilayer membranes; however, it does separate into liquid ordered (L_o) and liquid disordered (L_d) phases when the opposing leaflet containing sphingomyelin (SM), DOPC, and Chol demonstrates domain formation. We determine that interdigitated acyl chains modulated the partitioning of Chol in the opposing leaflet, resulting in phase separation. Similarly, the acyl chain length of SM within the opposing leaflet affected the phase behavior of the leaflet. Our results reveal the crucial role of interdigitation in determining the phase status in asymmetric membranes.

Computer simulations of protein-membrane systems

The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.

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.

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.

Understanding the Role of Lipids in Signaling Through Atomistic and Multiscale Simulations of Cell Membranes

Cell signaling controls essentially all cellular processes. While it is often assumed that proteins are the key architects coordinating cell signaling, recent studies have shown more and more clearly that lipids are also involved in signaling processes in a number of ways. Lipids do, for instance, act as messengers, modulate membrane receptor conformation and dynamics, and control membrane receptor partitioning. Further, through structural modifications such as oxidation, the functions of lipids as part of signaling processes can be modified. In this context, in this article we discuss the understanding recently revealed by atomistic and coarse-grained computer simulations of nanoscale processes and underlying physicochemical principles related to lipids’ functions in cellular signaling.

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.

Cholesterol Flip-Flop Impacts Domain Registration in Plasma Membrane Models

The plasma membrane is a highly complex multicomponent system that is central to the functioning of cells. Cholesterol, a key lipid component of the plasma membrane, promotes the formation of nanodomains. These nanodomains are often correlated across the two membrane leaflets, but the underlying physical mechanism remains unclear. Using coarse-grained molecular dynamics simulations, we investigate the influence of cholesterol flip-flop on membrane properties, in particular, the interleaflet coupling of cholesterol-enriched domains. We show that the cholesterol density correlation between the leaflets of an average mammalian plasma membrane is significantly reduced by suppressing interleaflet cholesterol population. Our results suggest an amplifying role of cholesterol in signal transduction across the leaflets.

Lipid self-assembly and lectin-induced reorganization of the plasma membrane

The plasma membrane represents an outstanding example of self-organization in biology. It plays a vital role in protecting the integrity of the cell interior and regulates meticulously the import and export of diverse substances. Its major building blocks are proteins and lipids, which self-assemble to a fluid lipid bilayer driven mainly by hydrophobic forces. Even if the plasma membrane appears-globally speaking-homogeneous at physiological temperatures, the existence of specialized nano- to micrometre-sized domains of raft-type character within cellular and synthetic membrane systems has been reported. It is hypothesized that these domains are the origin of a plethora of cellular processes, such as signalling or vesicular trafficking. This review intends to highlight the driving forces of lipid self-assembly into a bilayer membrane and the formation of small, transient domains within the plasma membrane. The mechanisms of self-assembly depend on several factors, such as the lipid composition of the membrane and the geometry of lipids. Moreover, the dynamics and organization of glycosphingolipids into nanometre-sized clusters will be discussed, also in the context of multivalent lectins, which cluster several glycosphingolipid receptor molecules and thus create an asymmetric stress between the two membrane leaflets, leading to tubular plasma membrane invaginations.This article is part of the theme issue 'Self-organization in cell biology'.

Glycan Chains of Gangliosides: Functional Ligands for Tissue Lectins (Siglecs/Galectins)

Molecular signals on the cell surface are responsible for adhesion and communication. Of relevance in this respect, their chemical properties endow carbohydrates with the capacity to store a maximum of information in a minimum of space. One way to present glycans on the cell surface is their covalent conjugation to a ceramide anchor. Among the resulting glycosphingolipids, gangliosides are special due to the presence of at least one sialic acid in the glycan chains. Their spatial accessibility and the dynamic regulation of their profile are factors that argue in favor of a role of glycans of gangliosides as ligands (counterreceptors) for carbohydrate-binding proteins (lectins). Indeed, as discovered first for a bacterial toxin, tissue lectins bind gangliosides and mediate contact formation (trans) and signaling (cis). While siglecs have a preference for higher sialylated glycans, certain galectins also target the monosialylated pentasaccharide of ganglioside GM1. Enzymatic interconversion of ganglioside glycans by sialidase action, relevant for neuroblastoma cell differentiation and growth control in vitro, for axonogenesis and axon regeneration, as well as for proper communication between effector and regulatory T cells, changes lectin-binding affinity profoundly. The GD1a-to-GM1 "editing" is recognized by such lectins, for example, myelin-associated glycoprotein (siglec-4) losing affinity and galectin-1 gaining reactivity, and then translated into postbinding signaling. Orchestrations of loss/gain of affinity, of ganglioside/lectin expression, and of lectin presence in a network offer ample opportunities for fine-tuning. Thus glycans of gangliosides such as GD1a and GM1 are functional counterreceptors by a pairing with tissue lectins, an emerging aspect of ganglioside and lectin functionality.

Effects of carotenoids on lipid bilayers

Carotenoids have been found to be important in improving the integrity of biomembranes in eukaryotes. However, the molecular details of how carotenoids modulate the physical properties of biomembranes are unknown. To this end, we have conducted a series of molecular dynamics simulations of different biologically-relevant membranes in the presence of carotenoids. The carotenoid effect on the membrane was found to be specific to the identity of the carotenoid and the composition of the membrane itself. Therefore, different classes of carotenoids produce a different effect on the membrane, and different membrane phases are affected differently by carotenoids. It is apparent from our data that carotenoids do trigger the bilayer to become thinner. The mechanism by which this occurs depends on two competing factors, the ability of the lipid tails of opposing monolayers to either (1) compress or (2) interdigitate as the bilayer condenses. Indeed, carotenoids directly influence the physical properties via these two mechanisms, thus compacting the bilayer. However, the degree to which these competing mechanisms are utilized depends on the bilayer phase and the carotenoid identity.

Gangliosides of the Nervous System

This review begins by attempting to recount some of the pioneering discoveries that first identified the presence of gangliosides in the nervous system, their structures and topography. This is presented as prelude to the current emphasis on physiological function, about which much has been learned but still remains to be elucidated. These areas include ganglioside roles in nervous system development including stem cell biology, membranes and organelles within neurons and glia, ion transport mechanisms, receptor modulation including neurotrophic factor receptors, and importantly the pathophysiological role of ganglioside aberrations in neurodegenerative disorders. This relates to their potential as therapeutic agents, especially in those conditions characterized by deficiency of one or more specific gangliosides. Finally we attempt to speculate on future directions ganglioside research is likely to take so as to capitalize on the impressive progress to date.