The Glycosphingolipid GM3 Modulates Conformational Dynamics of the Glucagon Receptor

Cellular lipids regulate the conformational ensembles of the disordered intracellular loop 3 in β2-adrenergic receptor

The intracellular loops of G protein-coupled receptors (GPCRs) have been shown to play a key role in G protein coupling and selectivity. We recently showed that the intrinsically disordered third intracellular loop (ICL3) of β2-adrenergic receptor is dynamic and equilibrates between open and closed conformations to regulate the G protein coupling. In this study, using the extensive molecular dynamics simulations in multi-lipid bilayer models, we show that the lipid phosphatidylinositol 4,5-bisphosphate (PIP2) stabilizes the active state of β2-adrenergic receptor by keeping ICL3 in an open conformation. This stabilization results in a tilt of the receptor within the membrane. Additionally, the ganglioside lipid, GM3 interacts with extracellular loops, impacting the ligand binding site allosterically. This demonstrates the active role of the chemistry of lipids in stabilizing specific GPCR conformations.

Cellular Lipids Regulate the Conformational Ensembles of the Disordered Intracellular Loop 3 in β2 Adrenergic Receptor

The structurally disordered intracellular loops (ICLs) of G protein-coupled receptors (GPCRs) play a critical role in G protein coupling. In our previous work, we used a combination of FRET-based and computational methodologies to show that the third intracellular loop (ICL3) modulates the activity and G protein coupling selectivity in GPCRs. In the current study, we have uncovered the role of several lipid components in modulating the conformational ensemble of ICL3 of the β2-adrenergic receptor (β2AR). Our findings indicate that phosphatidylinositol 4,5-bisphosphate (PIP2) in the inner leaflet of the membrane bilayer acts as a stabilizing anchor for ICL3, opening the intracellular cavity to facilitate G protein coupling. This interaction between PIP2 and ICL3 causes tilting of β2AR within the cellular membrane. Notably, this tilting of the receptor is supported by ganglioside GM3 stabilizing the extracellular loops on the outer leaflet of the bilayer, thereby exerting an allosteric effect on the orthosteric ligand binding pocket. Our results underscore the significance of lipids in modulating GPCR activity, proposing an allosteric mechanism that occurs through the receptor's orientation within the membrane.

Cardiolipin occupancy profiles of YidC paralogs reveal the significance of respective TM2 helix residues in determining paralog-specific phenotypes

YidC belongs to an evolutionarily conserved family of insertases, YidC/Oxa1/Alb3, in bacteria, mitochondria, and chloroplasts, respectively. Unlike Gram-negative bacteria, Gram-positives including Streptococcus mutans harbor two paralogs of YidC. The mechanism for paralog-specific phenotypes of bacterial YidC1 versus YidC2 has been partially attributed to the differences in their cytoplasmic domains. However, we previously identified a W138R gain-of-function mutation in the YidC1 transmembrane helix 2. YidC1W138R mostly phenocopied YidC2, yet the mechanism remained unknown. Primary sequence comparison of streptococcal YidCs led us to identify and mutate the YidC1W138 analog, YidC2S152 to W/A, which resulted in a loss of YidC2- and acquisition of YidC1-like phenotype. The predicted lipid-facing side chains of YidC1W138/YidC2S152 led us to propose a role for membrane phospholipids in specific-residue dependent phenotypes of S. mutans YidC paralogs. Cardiolipin (CL), a prevalent phospholipid in the S. mutans cytoplasmic membrane during acid stress, is encoded by a single gene, cls. We show a concerted mechanism for cardiolipin and YidC2 under acid stress based on similarly increased promoter activities and similar elimination phenotypes. Using coarse grain molecular dynamics simulations with the Martini2.2 Forcefield, YidC1 and YidC2 wild-type and mutant interactions with CL were assessed in silico. We observed substantially increased CL interaction in dimeric versus monomeric proteins, and variable CL occupancy in YidC1 and YidC2 mutant constructs that mimicked characteristics of the other wild-type paralog. Hence, paralog-specific amino acid- CL interactions contribute to YidC1 and YidC2-associated phenotypes that can be exchanged by point mutation at positions 138 or 152, respectively.

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.

Cholesterol Biases the Conformational Landscape of the Chemokine Receptor CCR3: A MAS SSNMR-Filtered Molecular Dynamics Study

Cholesterol directs the pathway of ligand-induced G protein-coupled receptor (GPCR) signal transduction. The GPCR C-C motif chemokine receptor 3 (CCR3) is the principal chemotactic receptor for eosinophils, with roles in cancer metastasis and autoinflammatory conditions. Recently, we discovered a direct correlation between bilayer cholesterol and increased agonist-triggered CCR3 signal transduction. However, the allosteric molecular mechanism escalating ligand affinity and G protein coupling is unknown. To study cholesterol-guided CCR3 conformational selection, we implement comparative, objective measurement of protein architectures by scoring shifts (COMPASS) to grade model structures from molecular dynamics simulations. In this workflow, we scored predicted chemical shifts against 2-dimensional solid-state NMR 13C-13C correlation spectra of U-15N,13C-CCR3 samples prepared with and without cholesterol. Our analysis of trajectory model structures uncovers that cholesterol induces site-specific conformational restraint of extracellular loop (ECL) 2 and conserved motion in transmembrane helices and ECL3 not observed in simulations of bilayers with only phosphatidylcholine lipids. PyLipID analysis implicates direct cholesterol agency in CCR3 conformational selection and dynamics. Residue-residue contact scoring shows that cholesterol biases the conformational selection of the orthosteric pocket involving Y411.39, Y1133.32, and E2877.39. Lastly, we observe contact remodeling in activation pathway residues centered on the initial transmission switch, Na+ pocket, and R3.50 in the DRY motif. Our observations have unique implications for understanding of CCR3 ligand recognition and specificity and provide mechanistic insight into how cholesterol functions as an allosteric regulator of CCR3 signal transduction.

Lipid Modulation of a Class B GPCR: Elucidating the Modulatory Role of PI(4,5)P2 Lipids

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂) lipids have been shown to stabilize an active conformation of class A G-protein coupled receptors (GPCRs) through a conserved binding site, not present in class B GPCRs. For class B GPCRs, previous molecular dynamics (MD) simulation studies have shown PI(4,5)P₂ interacting with the Glucagon receptor (GCGR), which constitutes an important target for diabetes and obesity therapeutics. In this work, we applied MD simulations supported by native mass spectrometry (nMS) to study lipid interactions with GCGR. We demonstrate how tail composition plays a role in modulating the binding of PI(4,5)P₂ lipids to GCGR. Specifically, we find the PI(4,5)P₂ lipids to have a higher affinity toward the inactive conformation of GCGR. Interestingly, we find that in contrast to class A GPCRs, PI(4,5)P₂ appear to stabilize the inactive conformation of GCGR through a binding site conserved across class B GPCRs but absent in class A GPCRs. This suggests differences in the regulatory function of PI(4,5)P₂ between class A and class B GPCRs.

Mutational Insight into Allosteric Regulation of Kir Channel Activity

Potassium (K⁺) channels are regulated in part by allosteric communication between the helical bundle crossing, or inner gate, and the selectivity filter, or outer gate. This network is triggered by gating stimuli. In concert, there is an allosteric network which is a conjugated set of interactions which correlate long-range structural rearrangements necessary for channel function. Inward-rectifier K⁺ (Kir) channels favor inward K⁺ conductance, are ligand-gated, and help establish resting membrane potentials. KirBac1.1 is a bacterial Kir (KirBac) channel homologous to human Kir (hKir) channels. Additionally, KirBac1.1 is gated by the anionic phospholipid ligand phosphatidylglycerol (PG). In this study, we use site-directed mutagenesis to investigate residues involved in the KirBac1.1 gating mechanism and allosteric network we previously proposed using detailed solid-state NMR (SSNMR) measurements. Using fluorescence-based K⁺ and sodium (Na⁺) flux assays, we identified channel mutants with impaired function that do not alter selectivity of the channel. In tandem, we performed coarse grain molecular dynamics simulations, observing changes in PG-KirBac1.1 interactions correlated with mutant channel activity and contacts between the two transmembrane helices and pore helix tied to this behavior. Lipid affinity is closely tied to the proximity of two tryptophan residues on neighboring subunits which lure anionic lipids to a cationic pocket formed by a cluster of arginine residues. Thus, these simulations establish a structural and functional basis for the role of each mutated site in the proposed allosteric network. The experimental and simulated data provide insight into key functional residues involved in gating and lipid allostery of K⁺ channels. Our findings also have direct implications on the physiology of hKir channels due to conservation of many of the residues identified in this work from KirBac1.1.

Hepatocyte cholesterol content modulates glucagon receptor signalling

OBJECTIVE

To determine whether glucagon receptor (GCGR) actions are modulated by cellular cholesterol levels.

METHODS

We determined the effects of experimental cholesterol depletion and loading on glucagon-mediated cAMP production, ligand internalisation and glucose production in human hepatoma cells, mouse and human hepatocytes. GCGR interactions with lipid bilayers were explored using coarse-grained molecular dynamic simulations. Glucagon responsiveness was measured in mice fed a high cholesterol diet with or without simvastatin to modulate hepatocyte cholesterol content.

RESULTS

GCGR cAMP signalling was reduced by higher cholesterol levels across different cellular models. Ex vivo glucagon-induced glucose output from mouse hepatocytes was enhanced by simvastatin treatment. Mice fed a high cholesterol diet had increased hepatic cholesterol and a blunted hyperglycaemic response to glucagon, both of which were partially reversed by simvastatin. Simulations identified likely membrane-exposed cholesterol binding sites on the GCGR, including a site where cholesterol is a putative negative allosteric modulator.

CONCLUSIONS

Our results indicate that cellular cholesterol content influences glucagon sensitivity and indicate a potential molecular basis for this phenomenon. This could be relevant to the pathogenesis of non-alcoholic fatty liver disease, which is associated with both hepatic cholesterol accumulation and glucagon resistance.

Understanding the functional role of membrane confinements in TNF-mediated signaling by multiscale simulations

The interaction between TNFα and TNFR1 is essential in maintaining tissue development and immune responses. While TNFR1 is a cell surface receptor, TNFα exists in both soluble and membrane-bound forms. Interestingly, it was found that the activation of TNFR1-mediated signaling pathways is preferentially through the soluble form of TNFα, which can also induce the clustering of TNFR1 on plasma membrane of living cells. We developed a multiscale simulation framework to compare receptor clustering induced by soluble and membrane-bound ligands. Comparing with the freely diffusive soluble ligands, we hypothesize that the conformational dynamics of membrane-bound ligands are restricted, which affects the clustering of ligand-receptor complexes at cell-cell interfaces. Our simulation revealed that only small clusters can form if TNFα is bound on cell surface. In contrast, the clustering triggered by soluble TNFα is more dynamic, and the size of clusters is statistically larger. We therefore demonstrated the impact of membrane-bound ligand on dynamics of receptor clustering. Moreover, considering that larger TNFα-TNFR1 clusters is more likely to provide spatial platform for downstream signaling pathway, our studies offer new mechanistic insights about why the activation of TNFR1-mediated signaling pathways is not preferred by membrane-bound form of TNFα.

The Bak core dimer focuses triacylglycerides in the membrane

Apoptosis, the intrinsic programmed cell death process, is mediated by the Bcl-2 family members Bak and Bax. Activation via formation of symmetric core dimers and oligomerization on the mitochondrial outer membrane (MOM) leads to permeabilization and cell death. Although this process is linked to the MOM, the role of the membrane in facilitating such pores is poorly understood. We recently described Bak core domain dimers, revealing lipid binding sites and an initial role of lipids in oligomerization. Here we describe simulations that identified localized clustering and interaction of triacylglycerides (TAGs) with a minimized Bak dimer construct. Coalescence of TAGs occurred beneath this Bak dimer, mitigating dimer-induced local membrane thinning and curvature in representative coarse-grain MOM and model membrane systems. Furthermore, the effects observed as a result of coarse-grain TAG cluster formation was concentration dependent, scaling from low physiological MOM concentrations to those found in other organelles. We find that increasing the TAG concentration in liposomes mimicking the MOM decreased the ability of activated Bak to permeabilize these liposomes. These results suggest that the presence of TAGs within a Bak-lipid membrane preserves membrane integrity and is associated with reduced membrane stress, suggesting a possible role of TAGs in Bak-mediated apoptosis.

Relative Affinities of Protein-Cholesterol Interactions from Equilibrium Molecular Dynamics Simulations

Specific interactions of lipids with membrane proteins contribute to protein stability and function. Multiple lipid interactions surrounding a membrane protein are often identified in molecular dynamics (MD) simulations and are, increasingly, resolved in cryo-electron microscopy (cryo-EM) densities. Determining the relative importance of specific interaction sites is aided by determination of lipid binding affinities using experimental or simulation methods. Here, we develop a method for determining protein-lipid binding affinities from equilibrium coarse-grained MD simulations using binding saturation curves, designed to mimic experimental protocols. We apply this method to directly obtain affinities for cholesterol binding to multiple sites on a range of membrane proteins and compare our results with free energies obtained from density-based equilibrium methods and with potential of mean force calculations, getting good agreement with respect to the ranking of affinities for different sites. Thus, our binding saturation method provides a robust, high-throughput alternative for determining the relative consequence of individual sites seen in, e.g., cryo-EM derived membrane protein structures surrounded by an array of ancillary lipid densities.

Multiplicity of Glycosphingolipid-Enriched Microdomain-Driven Immune Signaling

Glycosphingolipids (GSLs), together with cholesterol, sphingomyelin (SM), and glycosylphosphatidylinositol (GPI)-anchored and membrane-associated signal transduction molecules, form GSL-enriched microdomains. These specialized microdomains interact in a cis manner with various immune receptors, affecting immune receptor-mediated signaling. This, in turn, results in the regulation of a broad range of immunological functions, including phagocytosis, cytokine production, antigen presentation and apoptosis. In addition, GSLs alone can regulate immunological functions by acting as ligands for immune receptors, and exogenous GSLs can alter the organization of microdomains and microdomain-associated signaling. Many pathogens, including viruses, bacteria and fungi, enter host cells by binding to GSL-enriched microdomains. Intracellular pathogens survive inside phagocytes by manipulating intracellular microdomain-driven signaling and/or sphingolipid metabolism pathways. This review describes the mechanisms by which GSL-enriched microdomains regulate immune signaling.

Insights into the Role of Membrane Lipids in the Structure, Function and Regulation of Integral Membrane Proteins

Membrane proteins exist within the highly hydrophobic membranes surrounding cells and organelles, playing key roles in cellular function. It is becoming increasingly clear that the membrane does not just act as an appropriate environment for these proteins, but that the lipids that make up these membranes are essential for membrane protein structure and function. Recent technological advances in cryogenic electron microscopy and in advanced mass spectrometry methods, as well as the development of alternative membrane mimetic systems, have allowed experimental study of membrane protein–lipid complexes. These have been complemented by computational approaches, exploiting the ability of Molecular Dynamics simulations to allow exploration of membrane protein conformational changes in membranes with a defined lipid content. These studies have revealed the importance of lipids in stabilising the oligomeric forms of membrane proteins, mediating protein–protein interactions, maintaining a specific conformational state of a membrane protein and activity. Here we review some of the key recent advances in the field of membrane protein–lipid studies, with major emphasis on respiratory complexes, transporters, channels and G-protein coupled receptors.

In Silico Identification of Cholesterol Binding Motifs in the Chemokine Receptor CCR3

CC motif chemokine receptor 3 (CCR3) is a Class A G protein-coupled receptor (GPCR) mainly responsible for the cellular trafficking of eosinophils. As such, it plays key roles in inflammatory conditions, such as asthma and arthritis, and the metastasis of many deadly forms of cancer. However, little is known about how CCR3 functionally interacts with its bilayer environment. Here, we investigate cholesterol binding sites in silico through Coarse-Grained Molecular Dynamics (MD) and Pylipid analysis using an extensively validated homology model based on the crystal structure of CCR5. These simulations identified several cholesterol binding sites containing Cholesterol Recognition/Interaction Amino Acid Consensus motif (CRAC) and its inversion CARC motifs in CCR3. One such site, a CARC site in TM1, in conjunction with aliphatic residues in TM7, emerged as a candidate for future investigation based on the cholesterol residency time within the binding pocket. This site forms the core of a cholesterol binding site previously observed in computational studies of CCR2 and CCR5. Most importantly, these cholesterol binding sites are conserved in other chemokine receptors and may provide clues to cholesterol regulation mechanisms in this subfamily of Class A GPCRs.

Allosteric modulation of ghrelin receptor signaling by lipids

The membrane is an integral component of the G protein-coupled receptor signaling machinery. Here we demonstrate that lipids regulate the signaling efficacy and selectivity of the ghrelin receptor GHSR through specific interactions and bulk effects. We find that PIP2 shifts the conformational equilibrium of GHSR away from its inactive state, favoring basal and agonist-induced G protein activation. This occurs because of a preferential binding of PIP2 to specific intracellular sites in the receptor active state. Another lipid, GM3, also binds GHSR and favors G protein activation, but mostly in a ghrelin-dependent manner. Finally, we find that not only selective interactions but also the thickness of the bilayer reshapes the conformational repertoire of GHSR, with direct consequences on G protein selectivity. Taken together, this data illuminates the multifaceted role of the membrane components as allosteric modulators of how ghrelin signal could be propagated.

Structures and Dynamics of Native-State Transmembrane Protein Targets and Bound Lipids

Membrane proteins work within asymmetric bilayers of lipid molecules that are critical for their biological structures, dynamics and interactions. These properties are lost when detergents dislodge lipids, ligands and subunits, but are maintained in native nanodiscs formed using styrene maleic acid (SMA) and diisobutylene maleic acid (DIBMA) copolymers. These amphipathic polymers allow extraction of multicomponent complexes of post-translationally modified membrane-bound proteins directly from organ homogenates or membranes from diverse types of cells and organelles. Here, we review the structures and mechanisms of transmembrane targets and their interactions with lipids including phosphoinositides (PIs), as resolved using nanodisc systems and methods including cryo-electron microscopy (cryo-EM) and X-ray diffraction (XRD). We focus on therapeutic targets including several G protein-coupled receptors (GPCRs), as well as ion channels and transporters that are driving the development of next-generation native nanodiscs. The design of new synthetic polymers and complementary biophysical tools bodes well for the future of drug discovery and structural biology of native membrane:protein assemblies (memteins).

The Interplay of Glucagon-Like Peptide-1 Receptor Trafficking and Signalling in Pancreatic Beta Cells

The glucagon-like peptide 1 receptor (GLP-1R) is a class B G protein-coupled receptor (GPCR) which mediates the effects of GLP-1, an incretin hormone secreted primarily from L-cells in the intestine and within the central nervous system. The GLP-1R, upon activation, exerts several metabolic effects including the release of insulin and suppression of appetite, and has, accordingly, become an important target for the treatment for type 2 diabetes (T2D). Recently, there has been heightened interest in how the activated GLP-1R is trafficked between different endomembrane compartments, controlling the spatial origin and duration of intracellular signals. The discovery of "biased" GLP-1R agonists that show altered trafficking profiles and selective engagement with different intracellular effectors has added to the tools available to study the mechanisms and physiological importance of these processes. In this review we survey early and recent work that has shed light on the interplay between GLP-1R signalling and trafficking, and how it might be therapeutically tractable for T2D and related diseases.

A multiscale study on the mechanisms of spatial organization in ligand-receptor interactions on cell surfaces

The binding of cell surface receptors with extracellular ligands triggers distinctive signaling pathways, leading into the corresponding phenotypic variation of cells. It has been found that in many systems, these ligand-receptor complexes can further oligomerize into higher-order structures. This ligand-induced oligomerization of receptors on cell surfaces plays an important role in regulating the functions of cell signaling. The underlying mechanism, however, is not well understood. One typical example is proteins that belong to the tumor necrosis factor receptor (TNFR) superfamily. Using a generic multiscale simulation platform that spans from atomic to subcellular levels, we compared the detailed physical process of ligand-receptor oligomerization for two specific members in the TNFR superfamily: the complex formed between ligand TNFα and receptor TNFR1 versus the complex formed between ligand TNFβ and receptor TNFR2. Interestingly, although these two systems share high similarity on the tertiary and quaternary structural levels, our results indicate that their oligomers are formed with very different dynamic properties and spatial patterns. We demonstrated that the changes of receptor’s conformational fluctuations due to the membrane confinements are closely related to such difference. Consistent to previous experiments, our simulations also showed that TNFR can preassemble into dimers prior to ligand binding, while the introduction of TNF ligands induced higher-order oligomerization due to a multivalent effect. This study, therefore, provides the molecular basis to TNFR oligomerization and reveals new insights to TNFR-mediated signal transduction. Moreover, our multiscale simulation framework serves as a prototype that paves the way to study higher-order assembly of cell surface receptors in many other bio-systems.