Role of cholesterol-mediated effects in GPCR heterodimers

Probing the energy landscape of the lipid interactions of the serotonin1A receptor

G protein-coupled receptors (GPCRs) are lipid-regulated transmembrane proteins that play a central role in cell signaling and pharmacology. Although the role of membrane lipids in GPCR function is well established, the underlying GPCR-lipid interactions have not been thermodynamically characterized due to the complexity of these interactions. In this work, we estimate the energetics and dynamics of lipid association from coarse-grain simulations of the serotonin1A receptor embedded in a complex membrane. We show that lipids bind to the receptor with varying energetics of 1-4 kT, and timescales of 1-10 μs. The most favorable energetics and longest residence times are observed for cholesterol, glycosphingolipid GM1, phosphatidylethanolamine (PE) and phosphatidylserine (PS) lipids. Multi-exponential fitting of the contact probability suggests distinct dynamic regimes, corresponding to ps, ns and μs timescales, that we correlate with the annular, intermediate and non-annular lipid sites. The timescales of lipid binding correspond to high barrier heights, despite their relatively weaker energetics. Our results highlight that GPCR-lipid interactions are driven by both thermodynamic interactions and the dynamical features of lipid binding.

Relevance of G protein-coupled receptor (GPCR) dynamics for receptor activation, signalling bias and allosteric modulation

G protein-coupled receptors (GPCRs) are one of the major drug targets. In recent years, computational drug design for GPCRs has mainly focused on static structures obtained through X-ray crystallography, cryogenic electron microscopy (cryo-EM) or in silico modelling as a starting point for virtual screening campaigns. However, GPCRs are highly flexible entities with the ability to adopt different conformational states that elicit different physiological responses. Including this knowledge in the drug discovery pipeline can help to tailor novel conformation-specific drugs with an improved therapeutic profile. In this review, we outline our current knowledge about GPCR dynamics that is relevant for receptor activation, signalling bias and allosteric modulation. Ultimately, we highlight new technological implementations such as time-resolved X-ray crystallography and cryo-EM as well as computational algorithms that can contribute to a more comprehensive understanding of receptor dynamics and its relevance for GPCR functionality.

Cholesterol Dependence of the Conformational Changes in Metabotropic Glutamate Receptor 1

Metabotropic glutamate receptors (mGluRs) are class C G protein-coupled receptors that function as obligate dimers in regulating neurotransmission and synaptic plasticity in the central nervous system. The mGluR1 subtype has been shown to be modulated by the membrane lipid environment, particularly cholesterol, though the molecular mechanisms remain elusive. In this study, we employed all-atom molecular dynamics simulations to investigate the effects of cholesterol on the conformational dynamics of the mGluR1 seven-transmembrane (7TM) domain in an inactive state model. Simulations were performed with three different cholesterol concentrations (0%, 10%, and 25%) in a palmitoyl-oleoyl phosphatidylcholine (POPC) lipid bilayer system. Our results demonstrate that cholesterol induces conformational changes in the mGluR1 dimer more significantly than in the individual protomers. Notably, cholesterol modulates the dynamics and conformations of the TM1 and TM2 helices at the dimer interface. Interestingly, an intermediate cholesterol concentration of 10% elicits more pronounced conformational changes compared to both cholesterol-depleted (0%) and cholesterol-enriched (25%) systems. Specific electrostatic interaction unique to the 10% cholesterol system further corroborate these conformational differences. Given the high sequence conservation of the 7TM domains across mGluR subtypes, the cholesterol-dependent effects observed in mGluR1 are likely applicable to other members of this receptor family. Our findings provide atomistic insights into how cholesterol modulates the conformational landscape of mGluRs, which could impact their function and signaling mechanisms.

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.

Palmitoylation of the Glucagon-like Peptide-1 Receptor Modulates Cholesterol Interactions at the Receptor-Lipid Microenvironment

The G protein-coupled receptor (GPCR) superfamily of cell surface receptors has been shown to be functionally modulated by post-translational modifications. The glucagon-like peptide receptor-1 (GLP-1R), which is a drug target in diabetes and obesity, undergoes agonist-dependent palmitoyl tail conjugation. The palmitoylation in the C-terminal domain of GLP-1R has been suggested to modulate the receptor-lipid microenvironment. In this work, we have performed coarse-grain molecular dynamics simulations of palmitoylated and nonpalmitoylated GLP-1R to analyze the differential receptor-lipid interactions. Interestingly, the placement and dynamics of the C-terminal domain of GLP-1R are found to be directly dependent on the palmitoyl tail. We observe that both cholesterol and phospholipids interact with the receptor but display differential interactions in the presence and absence of the palmitoyl tail. We characterize important cholesterol-binding sites and validate sites that have been previously reported in experimentally resolved structures of the receptor. We show that the receptor acts like a conduit for cholesterol flip-flop by stabilizing cholesterol in the membrane core. Taken together, our work represents an important step in understanding the molecular effects of lipid modifications in GPCRs.

Cholesterol in Class C GPCRs: Role, Relevance, and Localization

G-protein coupled receptors (GPCRs), one of the largest superfamilies of cell-surface receptors, are heptahelical integral membrane proteins that play critical roles in virtually every organ system. G-protein-coupled receptors operate in membranes rich in cholesterol, with an imbalance in cholesterol level within the vicinity of GPCR transmembrane domains affecting the structure and/or function of many GPCRs, a phenomenon that has been linked to several diseases. These effects of cholesterol could result in indirect changes by altering the mechanical properties of the lipid environment or direct changes by binding to specific sites on the protein. There are a number of studies and reviews on how cholesterol modulates class A GPCRs; however, this area of study is yet to be explored for class C GPCRs, which are characterized by a large extracellular region and often form constitutive dimers. This review highlights specific sites of interaction, functions, and structural dynamics involved in the cholesterol recognition of the class C GPCRs. We summarize recent data from some typical family members to explain the effects of membrane cholesterol on the structural features and functions of class C GPCRs and speculate on their corresponding therapeutic potential.

Cofilin-Membrane Interactions: Electrostatic Effects in Phosphoinositide Lipid Binding

The actin cytoskeleton interacts with the cell membrane primarily through the indirect interactions of actin-binding proteins such as cofilin-1. The molecular mechanisms underlying the specific interactions of cofilin-1 with membrane lipids are still unclear. Here, we performed coarse-grain molecular dynamics simulations of cofilin-1 with complex lipid bilayers to analyze the specificity of protein-lipid interactions. We observed the maximal interactions with phosphoinositide (PIP) lipids, especially PIP₂ and PIP₃ lipids. A good match was observed between the residues predicted to interact and previous experimental studies. The clustering of PIP lipids around the membrane bound protein leads to an overall lipid demixing and gives rise to persistent membrane curvature. Further, through a series of control simulations, we observe that both electrostatics and geometry are critical for specificity of lipid binding. Our current study is a step towards understanding the physico-chemical basis of cofilin-PIP lipid interactions.

Synergistic and Competitive Lipid Interactions in the Serotonin1A Receptor Microenvironment

The interaction of lipids with G-protein-coupled receptors (GPCRs) has been shown to modulate and dictate several aspects of GPCR organization and function. Diverse lipid interaction sites have been identified from structural biology, bioinformatics, and molecular dynamics studies. For example, multiple cholesterol interaction sites have been identified in the serotonin_1A receptor, along with distinct and overlapping sphingolipid interaction sites. How these lipids interact with each other and what is the resultant effect on the receptor is still not clear. In this work, we have analyzed lipid-lipid crosstalk at the receptor of the serotonin_1A receptor embedded in a membrane bilayer that mimics the neuronal membrane composition by long coarse-grain simulations. Using a set of similarity coefficients, we classified lipids that bind at the receptor together as synergistic cobinding, and those that bind individually as competitive. Our results show that certain lipids interact with the serotonin_1A receptor in synergy with each other. Not surprisingly, the ganglioside GM1 and cholesterol show a synergistic cobinding, along with the relatively uncommon GM1-phosphatidylethanolamine (PE) and cholesterol-PE synergy. In contrast, certain lipid pairs such as cholesterol and sphingomyelin appear to be in competition at several sites, despite their coexistence in lipid nanodomains. In addition, we observed intralipid competition between two lipid tails, with the receptor exhibiting increased interactions with the unsaturated lipid tails. We believe our work represents an important step in understanding the diversity of GPCR-lipid interactions and exploring synergistic cobinding and competition in natural membranes.

Labelfree mapping and profiling of altered lipid homeostasis in the rat hippocampus after traumatic stress: Role of oxidative homeostasis

Oxidative and lipid homeostasis are altered by stress and trauma and post-traumatic stress disorder (PTSD) is associated with alterations to lipid species in plasma. Stress-induced alterations to lipid oxidative and homeostasis may exacerbate PTSD pathology, but few preclinical investigations of stress-induced lipidomic changes in the brain exist. Currently available techniques for the quantification of lipid species in biological samples require tissue extraction and are limited in their ability to retrieve spatial information. Raman imaging can overcome this limitation through the quantification of lipid species in situ in minimally processed tissue slices. Here, we utilized a predator exposure and psychosocial stress (PE/PSS) model of traumatic stress to standardize Raman imaging of lipid species in the hippocampus using LC-MS based lipidomics and these data were confirmed with qRT-PCR measures of mRNA expression of relevant enzymes and transporters. Electron Paramagnetic Resonance Spectroscopy (EPR) was used to measure free radical production and an MDA assay to measure oxidized polyunsaturated fatty acids. We observed that PE/PSS is associated with increased cholesterol, altered lipid concentrations, increased free radical production and reduced oxidized polyunsaturated fats (PUFAs) in the hippocampus (HPC), indicating shifts in lipid and oxidative homeostasis in the HPC after traumatic stress.

Molecular Mechanisms Underlying Caveolin-1 Mediated Membrane Curvature

Caveolin-1 is one of the main protein components of caveolae that acts as a mechanosensor at the cell membrane. The interactions of caveolin-1 with membranes have been shown to lead to complex effects such as curvature and the clustering of specific lipids. Here, we review the emerging concepts on the molecular interactions of caveolin-1, with a focus on insights from coarse-grain molecular dynamics simulations. Consensus structural models of caveolin-1 report a helix-turn-helix core motif with flanking domains of higher disorder that could be membrane composition dependent. Caveolin-1 appears to be mainly surface-bound and does not embed very deep in the membrane to which it is bound. The most interesting aspect of caveolin-1 membrane binding is the interplay of cholesterol clustering and membrane curvature. Although cholesterol has been reported to cluster in the vicinity of caveolin-1 by several approaches, simulations show that the clustering is maximal in membrane leaflet opposing the surface-bound caveolin-1. The intrinsic negative curvature of cholesterol appears to stabilize the negative curvature in the opposing leaflet. In fact, the simulations show that blocking cholesterol clustering (through artificial position restraints) blocks membrane curvature, and vice versa. Concomitant with cholesterol clustering is sphingomyelin clustering, again in the opposing leaflet, but in a concentration-dependent manner. The differential stress due to caveolin-1 binding and the inherent asymmetry of the membrane leaflets could be the determinant for membrane curvature and needs to be further probed. The review is an important step to reconcile the molecular level details emerging from simulations with the mesoscopic details provided by state of the art experimental approaches.

Molecular determinants of GPCR pharmacogenetics: Deconstructing the population variants in β2-adrenergic receptor

G protein-coupled receptors (GPCRs) are membrane proteins that play a central role in cell signaling and constitute one of the largest classes of drug targets. The molecular mechanisms underlying GPCR function have been characterized by several experimental and computational methods and provide an understanding of their role in physiology and disease. Population variants arising from nsSNPs affect the native function of GPCRs and have been implicated in differential drug response. In this chapter, we provide an overview on GPCR structure and activation, with a special focus on the β₂-adrenergic receptor (β₂-AR). First, we discuss the current understanding of the structural and dynamic features of the wildtype receptor. Subsequently, the population variants identified in this receptor from clinical and large-scale genomic studies are described. We show how computational approaches such as bioinformatics tools and molecular dynamics simulations can be used to characterize the variant receptors in comparison to the wildtype receptor. In particular, we discuss three examples of clinically important variants and discuss how the structure and function of these variants differ from the wildtype receptor at a molecular level. Overall, the chapter provides an overview of structure and function of GPCR variants and is a step towards the study of inter-individual differences and personalized medicine.

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.

Conformational plasticity and dynamic interactions of the N-terminal domain of the chemokine receptor CXCR1

The dynamic interactions between G protein-coupled receptors (GPCRs) and their cognate protein partners are central to several cell signaling pathways. For example, the association of CXC chemokine receptor 1 (CXCR1) with its cognate chemokine, interleukin-8 (IL8 or CXCL8) initiates pathways leading to neutrophil-mediated immune responses. The N-terminal domain of chemokine receptors confers ligand selectivity, but unfortunately the conformational dynamics of this intrinsically disordered region remains unresolved. In this work, we have explored the interaction of CXCR1 with IL8 by microsecond time scale coarse-grain simulations, complemented by atomistic models and NMR chemical shift predictions. We show that the conformational plasticity of the apo-receptor N-terminal domain is restricted upon ligand binding, driving it to an open C-shaped conformation. Importantly, we corroborated the dynamic complex sampled in our simulations against chemical shift perturbations reported by previous NMR studies and show that the trends are similar. Our results indicate that chemical shift perturbation is often not a reporter of residue contacts in such dynamic associations. We believe our results represent a step forward in devising a strategy to understand intrinsically disordered regions in GPCRs and how they acquire functionally important conformational ensembles in dynamic protein-protein interfaces.

Bacterial Analogs to Cholesterol Affect Dimerization of Proteorhodopsin and Modulates Preferred Dimer Interface

Hopanoids, the bacterial analogues of sterols, are ubiquitous in bacteria and play a significant role in organismal survival under stressful environments. Unlike sterols, hopanoids have a high degree of variation in the size and chemical nature of the substituent attached to the ring moiety, leading to different effects on the structure and dynamics of biological membranes. While it is understood that hopanoids can indirectly tune membrane physical properties, little is known on the role that hopanoids may play in affecting the organization and behavior of bacterial membrane proteins. In this work we used coarse-grained molecular dynamics simulations to characterize the effects of two hopanoids, diploptene (DPT) and bacteriohopanetetrol (BHT), on the oligomerization of proteorhodopsin (PR) in a model membrane composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phophoethanolamine (POPE) and 1-palmitoyl-2-oleoyl-sn-3-phosphoglycerol (POPG). PR is a bacterial membrane protein that functions as a light-activated proton pump. We chose PR based on its ability to adopt a distribution of oligomeric states in different membrane environments. Furthermore, the efficiency of proton pumping in PR is intimately linked to its organization into oligomers. Our results reveal that both BHT and DPT indirectly affect dimerization by tuning membrane properties in a fashion that is concentration-dependent. Variation in their interaction with PR in the membrane-embedded and the cytoplasmic regions leads to distinctly different effects on the plasticity of the dimer interface. BHT has the ability to intercalate between monomers in the dimeric interface, whereas DPT shifts dimerization interactions via packing of the interleaflet region of the membrane. Our results show a direct relationship between hopanoid structure and lateral organization of PR, providing a first glimpse at how these bacterial analogues to eukaryotic sterols produce very similar biophysical effects within the cell membrane.

Role of Cholesterol in Transmembrane Dimerization of the ErbB2 Growth Factor Receptor

The association of ErbB2 growth factor receptors is critical for cell growth and potentiates tumor proliferation in several cancer types. An important aspect in ErbB2 association is the role of lipids such as cholesterol, especially since their metabolism is often reprogrammed in cancer cells. Here, we have coupled metadynamics with coarse-grain simulations to identify cholesterol effects in the transmembrane dimerization of ErbB2 receptors. Overall, cholesterol interactions are observed with the receptor that directly tunes the association energetics. Several dimer conformations are identified both in the presence and absence of cholesterol, although the dimer regime appears to be more favorable in the presence of cholesterol. We observe an overall modulation of the underlying energy profile and the symmetric active and inactive conformational states are not distinguished in the presence of cholesterol. We show that cholesterol binds to the receptor transmembrane domain at a site (CRAC motif) that overlaps with the dimer interface (SmXXXSm motif). The competition between the transmembrane interactions and cholesterol interactions decides the final conformational landscape. Our work is an important step toward characterizing cholesterol effects in ErbB2 membrane receptor function.

Guardians of the Cell: State-of-the-Art of Membrane Proteins from a Computational Point-of-View

Membrane proteins (MPs) encompass a large family of proteins with distinct cellular functions, and although representing over 50% of existing pharmaceutical drug targets, their structural and functional information is still very scarce. Over the last years, in silico analysis and algorithm development were essential to characterize MPs and overcome some limitations of experimental approaches. The optimization and improvement of these methods remain an ongoing process, with key advances in MPs' structure, folding, and interface prediction being continuously tackled. Herein, we discuss the latest trends in computational methods toward a deeper understanding of the atomistic and mechanistic details of MPs.

Allosteric binding sites at the receptor-lipid bilayer interface: novel targets for GPCR drug discovery

As a superfamily of membrane receptors, G-protein-coupled receptors (GPCRs) have significant roles in human physiological processes, including cell proliferation, metabolism, and neuromodulation. GPCRs are vital targets of therapeutic drugs, and their allosteric regulation represents a novel direction for drug discovery. Given the numerous breakthroughs in structural biology, diverse allosteric sites on GPCRs have been identified within the extracellular and intracellular loops, and the seven core transmembrane helices. However, a unique type of allosteric site has also been discovered at the interface of the receptor-lipid bilayer, similar to the β₂-adrenergic receptor. Here, we review recent identifications of these allosteric sites and the detailed modulator-target interactions within the interface for each modulator to highlight the role of lipids in GPCR allosteric drug discovery.

Active cholesterol 20 years on

This review considers the following hypotheses, some well-supported and some speculative. Almost all of the sterol molecules in plasma membranes are associated with bilayer phospholipids in complexes of varied strength and stoichiometry. These complexes underlie many of the material properties of the bilayer. The small fraction of cholesterol molecules exceeding the binding capacity of the phospholipids is thermodynamically active and serves diverse functions. It circulates briskly among the cell membranes, particularly through contact sites linking the organelles. Active cholesterol provides the upstream feedback signal to multiple mechanisms governing plasma membrane homeostasis, pegging the sterol level to a threshold set by its phospholipids. Active cholesterol could also be the cargo for various inter-organelle transporters and the form excreted from cells by reverse transport. Furthermore, it is integral to the function of caveolae; a mediator of Hedgehog regulation; and a ligand for the binding of cytolytic toxins to membranes. Active cholesterol modulates a variety of plasma membrane proteins-receptors, channels and transporters-at least in vitro.

New perspectives on GPCRs: GPCR heterodimer formation (melanocortin receptor) and GPCR on primary cilia (melanin concentrating hormone receptor)

GPCRs are the largest family of receptors accounting for about 30% of the current drug targets. However, it is difficult to fully elucidate the mechanisms regulating intracellular GPCR signal regulation. It is thus important to consider and investigate GPCRs with respect to endogenous situations. Our group has been investigating GPCRs involved in body color (teleost and amphibian) and eating (vertebrate). Here, I review two independent GPCR systems (heterodimer formation and primary ciliated GPCR) that can be breakthroughs in GPCR research. In teleosts, MCRs form heterodimers, which significantly reduce their affinity for acetylated ligands. In mammals, MCHR1 is localized in the ciliary membrane and shortens the length of the primary cilia through a unique signal from the ciliary membrane. Considering these two new GPCR concepts is expected to advance the overall view of the GPCR system.

Cholesterol-mediated oligomerization pathways of serotonin G-coupled receptor 5-HT2C

Serotonin (5-HT) receptors have been shown to homodimerize and heterodimerize with other G protein-coupled receptors (GPCRs), although the details of this process have not yet been elucidated. Here we use coarse-grained molecular dynamics on monomeric 5-HT2C receptors to predict the transmembrane (TM) helices involved in such associations. All these simulations were carried out both in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayers and in mixed composition POPC-Cholesterol ones, to show whether the presence of cholesterol could directly influence and drive the dimeric association. The goal is to get insights on the self-assembly pathway leading to GPCRs 5-HT2C oligomerization, which is supposed to be the basis of its constitutional activity. From the analysis of the molecular dynamics trajectories, we observed the formation of 5-HT2C oligomers through self-assembly and we identified the main domains involved in the receptor dimerization. In particular, dimers and oligomers from the two different environments show TM4-TM5 and TM1-TM7-H8 as the preferential dimerization interfaces. Nevertheless, substantial differences arise for oligomers in POPC and in POPC-Chol membranes: in POPC-Chol the variability of dimers interfaces is strictly limited to the TM1-TM7-H8 and TM4-TM5 interfaces and the dimorphism depends on cholesterol that directly participates in its formation. These results are in agreement with both experimental evidences and other computational studies conducted on other GPCRs oligomerization.