Are specific nonannular cholesterol binding sites present in G-protein coupled receptors?

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.

Improved Highly Mobile Membrane Mimetic Model for Investigating Protein-Cholesterol Interactions

Cholesterol (CHL) plays an integral role in modulating the function and activity of various mammalian membrane proteins. Due to the slow dynamics of lipids, conventional computational studies of protein-CHL interactions rely on either long-time scale atomistic simulations or coarse-grained approximations to sample the process. A highly mobile membrane mimetic (HMMM) has been developed to enhance lipid diffusion and thus used to facilitate the investigation of lipid interactions with peripheral membrane proteins and, with customized in silico solvents to replace phospholipid tails, with integral membrane proteins. Here, we report an updated HMMM model that is able to include CHL, a nonphospholipid component of the membrane, henceforth called HMMM-CHL. To this end, we had to optimize the effect of the customized solvents on CHL behavior in the membrane. Furthermore, the new solvent is compatible with simulations using force-based switching protocols. In the HMMM-CHL, both improved CHL dynamics and accelerated lipid diffusion are integrated. To test the updated model, we have applied it to the characterization of protein-CHL interactions in two membrane protein systems, the human β2-adrenergic receptor (β2AR) and the mitochondrial voltage-dependent anion channel 1 (VDAC-1). Our HMMM-CHL simulations successfully identified CHL binding sites and captured detailed CHL interactions in excellent consistency with experimental data as well as other simulation results, indicating the utility of the improved model in applications where an enhanced sampling of protein-CHL interactions is desired.

Illuminating the structural basis of human neurokinin 1 receptor (NK1R) antagonism through classical all-atoms molecular dynamics simulations

Advances in structural biology have bestowed insights into the pleiotropic effects of neurokinin 1 receptors (NK1R) in diverse patho-physiological processes, thereby highlighting the potential therapeutic value of antagonists directed against NK1R. Herein, we investigate the mode of antagonist recognition to discern the obscure atomic facets germane for the function and molecular determinants of NK1R. To commence discernment of potent antagonists and the conformational changes in NK1R, induced upon antagonist binding, state-of-the-art classical all-atoms molecular dynamics (MD) simulations in lipid mimetic bilayers have been utilized. MD simulations of structural ensembles reveals the involvement of TM5 and TM6 in tight anchoring of antagonists through a network of interhelical hydrogen-bonds, while, the extracellular loop 2 (ECL2) governs the overall size and nature of the pocket, thereby modulating NK1R. Consistent comparison between experiments and MD simulation results discerns the predominant role of TM3, TM4, and TM6 in lipid-NK1R interaction. Correlation between hydrophobic index and helicity of TM domains elucidates their importance in maintaining the structural stability in addition to regulating NK1R antagonism. Taken together, we anticipate that our computational study marks a comprehensive structural basis of NK1R antagonism in lipid bilayers, which may facilitate designing of new therapeutics against associated diseases targeting human neurokinin receptors.

Physiological changes in bilayer thickness induced by cholesterol control GPCR rhodopsin function

We monitored the effect on function of the G-protein-coupled receptor (GPCR) rhodopsin from small, stepwise changes in bilayer thickness induced by cholesterol. Over a range of phosphatidylcholine bilayers with hydrophobic thickness from ≈21 Å to 38 Å, the metarhodopsin-I (MI)/metarhodopsin-II (MII) equilibrium was monitored with UV-visible spectroscopy while ordering of hydrocarbon chains was probed by 2H-NMR. Addition of cholesterol shifted equilibrium toward MII for bilayers thinner than the average length of hydrophobic transmembrane helices (27 Å) and to MI for thicker bilayers, while small bilayer thickness changes within the range of the protein hydrophobic thickness drastically up- or downregulated MII formation. The cholesterol-induced shifts toward MII for thinner membranes correlated with the cholesterol-induced increase of bilayer hydrophobic thickness measured by NMR, consistent with continuum elastic modeling. The energetic penalty of adding cholesterol to thick bilayers caused rhodopsin oligomerization and a shift toward MI. In membranes of physiological thickness, changes in bilayer mechanical properties induced by cholesterol potentiated the interplay between bilayer and protein thickness resulting in large swings of the MI-MII equilibrium. In membrane containing cholesterol, elastic deformations near the protein are a dominant energetic contribution to the functional equilibrium of the model GPCR rhodopsin.

Insights of ligand binding in modeled h5-HT1A receptor: homology modeling, docking, MM-GBSA, screening and molecular dynamics

The pharmacologically characterized receptor subtype of the serotonin family, the 5HT1A receptor is implicated in the pathophysiology and treatment of depression and anxiety-related disorders. Being the most extensively targeted receptor for developing novel antidepressants and anxiolytics, a near-ideal theoretical model can aid in high-throughput screening of promising drug candidates. However, the design of potential drug candidates suffers owing to a lack of complete structural information. In this work, homology models of 5-HT1A receptor are generated using two distinct alignments (CW and PSTA) and model building methods (KB and EB). The developed models are validated for virtual screening using a ligand dataset of agonists and antagonists. The best-suited model was efficient in discriminating agonist/antagonist binding. Correlation plots between pKi and docking (R²_agonist≥ 0.6, R²_antagonist≥ 0.7) and MM-GBSA dG bind values (R²_agonist≥ 0.5, R²_antagonist≥ 0.7) revealed optimum corroboration between in vitro and in silico outcomes, which further suggested the usefulness of the developed model for the design of high-affinity probes for the neurological disorders.Communicated by Ramaswamy H. Sarma.

Allosteric Modulation of Muscarinic Receptors by Cholesterol, Neurosteroids and Neuroactive Steroids

Muscarinic acetylcholine receptors are membrane receptors involved in many physiological processes. Malfunction of muscarinic signaling is a cause of various internal diseases, as well as psychiatric and neurologic conditions. Cholesterol, neurosteroids, neuroactive steroids, and steroid hormones are molecules of steroid origin that, besides having well-known genomic effects, also modulate membrane proteins including muscarinic acetylcholine receptors. Here, we review current knowledge on the allosteric modulation of muscarinic receptors by these steroids. We give a perspective on the research on the non-genomic effects of steroidal compounds on muscarinic receptors and drug development, with an aim to ultimately exploit such knowledge.

Detergent-Free Isolation of Membrane Proteins and Strategies to Study Them in a Near-Native Membrane Environment

Atomic-resolution structural studies of membrane-associated proteins and peptides in a membrane environment are important to fully understand their biological function and the roles played by them in the pathology of many diseases. However, the complexity of the cell membrane has severely limited the application of commonly used biophysical and biochemical techniques. Recent advancements in NMR spectroscopy and cryoEM approaches and the development of novel membrane mimetics have overcome some of the major challenges in this area. For example, the development of a variety of lipid-nanodiscs has enabled stable reconstitution and structural and functional studies of membrane proteins. In particular, the ability of synthetic amphipathic polymers to isolate membrane proteins directly from the cell membrane, along with the associated membrane components such as lipids, without the use of a detergent, has opened new avenues to study the structure and function of membrane proteins using a variety of biophysical and biological approaches. This review article is focused on covering the various polymers and approaches developed and their applications for the functional reconstitution and structural investigation of membrane proteins. The unique advantages and limitations of the use of synthetic polymers are also discussed.

Cancer-Related Somatic Mutations in Transmembrane Helices Alter Adenosine A1 Receptor Pharmacology

Overexpression of the adenosine A₁ receptor (A₁AR) has been detected in various cancer cell lines. However, the role of A₁AR in tumor development is still unclear. Thirteen A₁AR mutations were identified in the Cancer Genome Atlas from cancer patient samples. We have investigated the pharmacology of the mutations located at the 7-transmembrane domain using a yeast system. Concentration-growth curves were obtained with the full agonist CPA and compared to the wild type hA₁AR. H78L^3.23 and S246T^6.47 showed increased constitutive activity, while only the constitutive activity of S246T^6.47 could be reduced to wild type levels by the inverse agonist DPCPX. Decreased constitutive activity was observed on five mutant receptors, among which A52V^2.47 and W188C^5.46 showed a diminished potency for CPA. Lastly, a complete loss of activation was observed in five mutant receptors. A selection of mutations was also investigated in a mammalian system, showing comparable effects on receptor activation as in the yeast system, except for residues pointing toward the membrane. Taken together, this study will enrich the view of the receptor structure and function of A₁AR, enlightening the consequences of these mutations in cancer. Ultimately, this may provide an opportunity for precision medicine for cancer patients with pathological phenotypes involving these mutations.

A cholesterol analog stabilizes the human β2-adrenergic receptor nonlinearly with temperature

In cell membranes, G protein-coupled receptors (GPCRs) interact with cholesterol, which modulates their assembly, stability, and conformation. Previous studies have shown how cholesterol modulates the structural properties of GPCRs at ambient temperature. Here, we characterized the mechanical, kinetic, and energetic properties of the human β2-adrenergic receptor (β2AR) in the presence and absence of the cholesterol analog cholesteryl hemisuccinate (CHS) at room temperature (25°C), at physiological temperature (37°C), and at high temperature (42°C). We found that CHS stabilized various structural regions of β2AR differentially, which changed nonlinearly with temperature. Thereby, the strongest effects were observed for structural regions that are important for receptor signaling. Moreover, at 37°C, but not at 25° or 42°C, CHS caused β2AR to increase and stabilize conformational substates to adopt to basal activity. These findings indicate that the nonlinear, temperature-dependent action of CHS in modulating the structural and functional properties of this GPCR is optimized for 37°C.

Impact of lipid rafts on transient receptor potential channel activities

Members of the transient receptor potential (TRP) superfamily are cation channels that are expressed in nearly every mammalian cell type and respond as cellular sensors to various environmental stimuli. Light, pressure, osmolarity, temperature, and other stimuli can induce TRP calcium conductivity and correspondingly trigger many signaling processes in cells. Disruption of TRP channel activity, as a rule, harms cellular function. Despite numerous studies, the mechanisms of TRP channel regulation are not yet sufficiently clear, in part, because TRP channels are regulated by a broad set of ligands having diverse physical and chemical features. It is now known that some TRP members are located in membrane microdomains termed lipid rafts. Moreover, interaction between specific raft-associated lipids with channels may be a key regulation mechanism. This review examines recent findings related to the roles of lipid rafts in regulation of TRP channel activity. The mechanistic events of channel interactions with the main lipid raft constituent, cholesterol, are being clarified. Better understanding of mechanisms behind such interactions would help establish the key elements of TRP channel regulation and hence allow control of cellular responses to environmental stimuli.

Molecular dynamics study of water transport through AQP5-R188C mutant causing palmoplantar keratoderma (PPK) using the gating mechanism concept

It is widely known that any disruption to the water regulation in aquaporins (AQPs) leads to numerous important diseases. However, studies of dynamics and energetics of disease-causing mutations in the aquaporins on the molecular level are still limited. In the present work, the effects of a skin disease-causing mutant, R188C, on the structure of AQP5 and water transport mechanism within this mutated aquaporin are investigated using the concept of gating mechanism. Our results have revealed that the R188C mutation causes a remarkable increase in the pore radius inside the selectivity filter (SF) region facilitating the passage of water molecules. This observation is supported by plotting the free energy profiles of water molecules transport and calculating permeability values through AQP5-R188C, such that the energy barrier in the SF region of the pores was substantially reduced by this mutation, and therefore, the translocation of water molecules was improved. The total averaged osmotic permeability for R188C has been computed as about 11-fold of the wild-type permeability. However, a comparison between the osmotic permeability values related to the open conformation of CE revealed that this coefficient for AQP5-R188C is about 6.5 times larger than that of wt-AQP5, which can be a more accurate value according to the gating mechanism associated with the constriction region of the aquaporin.

Cholesterol in GPCR Structures: Prevalence and Relevance

Bound cholesterol molecules are emerging as important hallmarks of GPCR structures. In this commentary, we analyze their statistical prevalence and biological relevance.

Cholesterol footprint in high-resolution structures of serotonin receptors: Where are we now and what does it mean?

An emerging feature of several high-resolution GPCR structures is the presence of closely bound cholesterol molecules. In this Perspective, we share the excitement of the recent advancements in GPCR structural biology. We further highlight our laboratory's journey in comprehensively elucidating functional sensitivity of GPCRs (using the serotonin_1A receptor as a representative neurotransmitter GPCR) to membrane cholesterol and validation using a variety of assays and molecular dynamics simulations. Although high-resolution structures of many GPCRs have been reported in the last few years, the structure of the serotoin_1A receptor proved to be elusive for a long time. Very recently the cryo-EM structure of the serotoin_1A receptor displaying 10 bound cholesterol molecules has been reported. We conclude by providing a critical analysis of caveats involved in GPCR structure determination.

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.

Increasing complexity in small-angle X-ray and neutron scattering experiments: from biological membrane mimics to live cells

Small-angle X-ray and neutron scattering are well-established, non-invasive experimental techniques to interrogate global structural properties of biological membrane mimicking systems under physiologically relevant conditions. Recent developments, both in bottom-up sample preparation techniques for increasingly complex model systems, and in data analysis techniques have opened the path toward addressing long standing issues of biological membrane remodelling processes. These efforts also include emerging quantitative scattering studies on live cells, thus enabling a bridging of molecular to cellular length scales. Here, we review recent progress in devising compositional models for joint small-angle X-ray and neutron scattering studies on diverse membrane mimics - with a specific focus on membrane structural coupling to amphiphatic peptides and integral proteins - and live Escherichia coli. In particular, we outline the present state-of-the-art in small-angle scattering methods applied to complex membrane systems, highlighting how increasing system complexity must be followed by an advance in compositional modelling and data-analysis tools.

Structure and Dynamics of GPCRs in Lipid Membranes: Physical Principles and Experimental Approaches

Over the past decade, the vast amount of information generated through structural and biophysical studies of GPCRs has provided unprecedented mechanistic insight into the complex signalling behaviour of these receptors. With this recent information surge, it has also become increasingly apparent that in order to reproduce the various effects that lipids and membranes exert on the biological function for these allosteric receptors, in vitro studies of GPCRs need to be conducted under conditions that adequately approximate the native lipid bilayer environment. In the first part of this review, we assess some of the more general effects that a membrane environment exerts on lipid bilayer-embedded proteins such as GPCRs. This is then followed by the consideration of more specific effects, including stoichiometric interactions with specific lipid subtypes. In the final section, we survey a range of different membrane mimetics that are currently used for in vitro studies, with a focus on NMR applications.

Lipid-dependent conformational landscape of the ErbB2 growth factor receptor dimers

Altered lipid metabolism has been linked to cancer development and progression. Several roles have been attributed to the increased saturation and length of lipid acyl tails observed in tumors, but its effect on signaling receptors is still emerging. In this work, we have analyzed the lipid dependence of the ErbB2 growth factor receptor dimerization that plays an important role in the pathogenesis of breast cancer. We have performed coarse-grain ensemble molecular dynamics simulations to comprehensively sample the ErbB2 monomer-dimer association. Our results indicate a dynamic dimer state with a complex conformational landscape that is modulated with increasing lipid tail length. We resolve the native N-terminal "active" and C-terminal "inactive" conformations in all membrane compositions. However, the relative population of the N-terminal and C-terminal conformers is dependent on length of the saturated lipid tails. In short-tail membranes, additional non-specific dimers are observed which are reduced or absent in long-tailed bilayers. Our results indicate that the relative population as well as the structure of the dimer state is modulated by membrane composition. We have correlated these differences to local perturbations of the membrane around the receptor. Our work is an important step in characterizing ErbB dimers in healthy and diseased states and emphasize the importance of sampling lipid dynamics in understanding receptor association.

Expression and purification of recombinant G protein-coupled receptors: A review

Given their extensive role in cell signalling, GPCRs are significant drug targets; despite this, many of these receptors have limited or no available prophylaxis. Novel drug design and discovery significantly rely on structure determination, of which GPCRs are typically elusive. Progress has been made thus far to produce sufficient quantity and quality of protein for downstream analysis. As such, this review highlights the systems available for recombinant GPCR expression, with consideration of their advantages and disadvantages, as well as examples of receptors successfully expressed in these systems. Additionally, an overview is given on the use of detergents and the styrene maleic acid (SMA) co-polymer for membrane solubilisation, as well as purification techniques.

Cholesterol interaction motifs in G protein-coupled receptors: Slippery hot spots?

G protein-coupled receptors (GPCRs) are cell membrane associated signaling hubs that orchestrate a multitude of cellular functions upon binding to a diverse variety of extracellular ligands. Since GPCRs are integral membrane proteins with seven-transmembrane domain architecture, their function, organization and dynamics are intimately regulated by membrane lipids, such as cholesterol. Cholesterol is an extensively studied lipids in terms of its effects on GPCR structure and function. One of the possible mechanisms underlying modulation of GPCR function by cholesterol is via specific interaction of GPCRs with membrane cholesterol. These interactions of GPCRs with membrane cholesterol are often attributed to structural features of GPCRs that could facilitate their preferential association with cholesterol. In this backdrop, cholesterol interaction motifs represent putative interaction sites on GPCRs that could facilitate cholesterol-sensitive function of these receptors. In this review, we provide an overview of cholesterol interaction motifs found in GPCRs, which have been identified through a combination of crystallography, bioinformatics analysis, and functional studies. In addition, we will highlight, using specific examples, why mere presence of a cholesterol interaction motif at a given site may not directly implicate its role in interaction with membrane cholesterol. We therefore believe that experimental approaches, followed by functional analysis of cholesterol sensitivity of GPCRs, would provide a better understanding of the role played by these motifs in cholesterol-sensitive function. We envision that a comprehensive knowledge of cholesterol interaction sites in GPCRs would allow us to develop a better understanding of GPCR structure-function paradigm, and could be useful in future therapeutics. This article is categorized under: Models of Systems Properties and Processes > Mechanistic Models Analytical and Computational Methods > Computational Methods Laboratory Methods and Technologies > Macromolecular Interactions, Methods.

Role of cholesterol-mediated effects in GPCR heterodimers

G protein-coupled receptors (GPCRs) are transmembrane receptors that mediate a large number of cellular responses. The organization of GPCRs into dimers and higher-order oligomers is known to allow a larger repertoire of downstream signaling events. In this context, a crosstalk between the adenosine and dopamine receptors has been reported, indicating the presence of heterodimers that are functionally relevant. In this paper, we explored the effect of membrane cholesterol on the adenosine_2A (A_2A) and dopamine D₃ (D₃) receptors using coarse-grain molecular dynamics simulations. We analyzed cholesterol interaction sites on the A_2A receptor and were able to reproduce the sites indicated by crystallography and previous atomistic simulations. We predict novel cholesterol interaction sites on the D₃ receptor that could be important in the reported cholesterol sensitivity in receptor function. Further, we analyzed the formation of heterodimers between the two receptors. Our results suggest that membrane cholesterol modulates the relative population of several co-existing heterodimer conformations. Both direct receptor-cholesterol interaction and indirect membrane effects contribute toward the modulation of heterodimer conformations. These results constitute one of the first examples of modulation of GPCR hetero-dimerization by membrane cholesterol, and could prove to be useful in designing better therapeutic strategies.

Regulation of Membrane Calcium Transport Proteins by the Surrounding Lipid Environment

Calcium ions (Ca2+) are major messengers in cell signaling, impacting nearly every aspect of cellular life. Those signals are generated within a wide spatial and temporal range through a large variety of Ca2+ channels, pumps, and exchangers. More and more evidences suggest that Ca2+ exchanges are regulated by their surrounding lipid environment. In this review, we point out the technical challenges that are currently being overcome and those that still need to be defeated to analyze the Ca2+ transport protein-lipid interactions. We then provide evidences for the modulation of Ca2+ transport proteins by lipids, including cholesterol, acidic phospholipids, sphingolipids, and their metabolites. We also integrate documented mechanisms involved in the regulation of Ca2+ transport proteins by the lipid environment. Those include: (i) Direct interaction inside the protein with non-annular lipids; (ii) close interaction with the first shell of annular lipids; (iii) regulation of membrane biophysical properties (e.g., membrane lipid packing, thickness, and curvature) directly around the protein through annular lipids; and (iv) gathering and downstream signaling of several proteins inside lipid domains. We finally discuss recent reports supporting the related alteration of Ca2+ and lipids in different pathophysiological events and the possibility to target lipids in Ca2+-related diseases.

Changes in the function of angiotensin II type 1 receptor due to cholesterol depletion from cell membrane

Blockers of G-protein coupled receptors (GPCRs), angiotensin II (Ang II) type 1 (AT₁) receptor and β₁-adrenergic (Ad) receptor, have been shown to improve the prognosis of cardiovascular disease. Cholesterol molecules in the cell membrane are needed to stabilize GPCRs as well as the cell membrane itself. We determined whether the functions of AT₁ and β₁-Ad receptors were changed by cholesterol depletion from cardiovascular cell membranes. Ang II-induced inositol phosphate production through AT₁ receptor was suppressed by cholesterol depletion from cell membranes using rosuvastatin or methyl-β-cyclodextrin (MβCD), whereas isoproterenol-induced cyclic AMP production through β₁-Ad receptor did not change after cholesterol depletion. In addition, the binding affinities of Ang II and AT₁ receptor blocker after cholesterol depletion were significantly lower than those before depletion. Although AT₁ receptor expression levels did not change after cholesterol depletion, the expression levels of AT₁ receptor that could bind to Ang II significantly decreased after depletion. The changes in the structure of AT₁ receptor due to depletion were confirmed by substituted-cysteine accessibility mapping. In conclusion, Ang II-induced activation of AT₁ receptor is reduced without affecting the function of β₁-Ad receptor after cholesterol depletion from cardiovascular cell membranes.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

A collage of cholesterol interaction motifs in the serotonin1A receptor: An evolutionary implication for differential cholesterol interaction

The serotonin_1A receptor is a representative member of the G protein-coupled receptor (GPCR) superfamily and acts as an important drug target. In our previous work, we comprehensively demonstrated that membrane cholesterol is necessary in the organization, dynamics and function of the serotonin_1A receptor. In this context, analysis of high-resolution GPCR crystal structures in general and in silico studies of the serotonin_1A receptor in particular, have suggested the presence of cholesterol interaction sites (hotspots) in various regions of the receptor. In this work, we have identified an evolutionarily conserved collage of four categories of cholesterol interaction motifs associated with transmembrane helix V and the adjacent intracellular loop 3 fragment of the vertebrate serotonin_1A receptor. This collage of motifs represents a total of twenty diverse context-dependent cholesterol interaction configurations. We envision that the gamut of cholesterol interaction sites, characterized by sequence plasticity in cholesterol interaction, could be relevant in receptor-cholesterol interaction in membranes of varying cholesterol content and organization, as found in diverse cell types. We conclude that an evolutionarily conserved mechanism of GPCR-cholesterol interaction allows the serotonin_1A receptor to adapt to diverse membrane cholesterol levels during natural evolution.

Cholesterol-Recognition Motifs in Membrane Proteins

The impact of cholesterol on the structure and function of membrane proteins was recognized several decades ago, but the molecular mechanisms underlying these effects have remained elusive. There appear to be multiple mechanisms by which cholesterol interacts with proteins. A complete understanding of cholesterol-sensing motifs is still undergoing refinement. Initially, cholesterol was thought to exert only non-specific effects on membrane fluidity. It was later shown that this lipid could specifically interact with membrane proteins and affect both their structure and function. In this article, we have summarized and critically analyzed our evolving understanding of the affinity, specificity and stereoselectivity of the interactions of cholesterol with membrane proteins. We review the different computational approaches that are currently used to identify cholesterol binding sites in membrane proteins and the biochemical logic that governs each type of site, including CRAC, CARC, SSD and amphipathic helix motifs. There are physiological implications of these cholesterol-recognition motifs for G-protein coupled receptors (GPCR) and ion channels, in membrane trafficking and membrane fusion (SNARE) proteins. There are also pathological implications of cholesterol binding to proteins involved in neurological disorders (Alzheimer, Parkinson, Creutzfeldt-Jakob) and HIV fusion. In each case, our discussion is focused on the key molecular aspects of the cholesterol and amino acid motifs in membrane-embedded regions of membrane proteins that define the physiologically relevant crosstalk between the two. Our understanding of the factors that determine if these motifs are functional in cholesterol binding will allow us enhanced predictive capabilities.

A Critical Analysis of Molecular Mechanisms Underlying Membrane Cholesterol Sensitivity of GPCRs

G protein-coupled receptors (GPCRs) are the largest and a diverse family of proteins involved in signal transduction across biological membranes. GPCRs mediate a wide range of physiological processes and have emerged as major targets for the development of novel drug candidates in all clinical areas. Since GPCRs are integral membrane proteins, regulation of their organization, dynamics, and function by membrane lipids, in particular membrane cholesterol, has emerged as an exciting area of research. Cholesterol sensitivity of GPCRs could be due to direct interaction of cholesterol with the receptor (specific effect). Alternately, GPCR function could be influenced by the effect of cholesterol on membrane physical properties (general effect). In this review, we critically analyze the specific and general mechanisms of the modulation of GPCR function by membrane cholesterol, taking examples from representative GPCRs. While evidence for both the proposed mechanisms exists, there appears to be no clear-cut distinction between these two mechanisms, and a combination of these mechanisms cannot be ruled out in many cases. We conclude that classifying the mechanism underlying cholesterol sensitivity of GPCR function merely into these two mutually exclusive classes could be somewhat arbitrary. A more holistic approach could be suitable for analyzing GPCR-cholesterol interaction.

Dexamethasone and Fludrocortisone Inhibit Hedgehog Signaling in Embryonic Cells

The hedgehog (Hh) pathway plays a central role in the development and repair of our bodies. Therefore, dysregulation of the Hh pathway is responsible for many developmental diseases and cancers. Basal cell carcinoma and medulloblastoma have well-established links to the Hh pathway, as well as many other cancers with Hh-dysregulated subtypes. A smoothened (SMO) receptor plays a central role in regulating the Hh signaling in the cells. However, the complexities of the receptor structural mechanism of action and other pathway members make it difficult to find Hh pathway inhibitors efficient in a wide range. Recent crystal structure of SMO with cholesterol indicates that it may be a natural ligand for SMO activation. Structural similarity of fluorinated corticosterone derivatives to cholesterol motivated us to study the effect of dexamethasone, fludrocortisone, and corticosterone on the Hh pathway activity. We identified an inhibitory effect of these three drugs on the Hh pathway using a functional assay in NIH3T3 glioma response element cells. Studies using BODIPY-cyclopamine and 20(S)-hydroxy cholesterol [20(S)-OHC] as competitors for the transmembrane (TM) and extracellular cysteine-rich domain (CRD) binding sites showed a non-competitive effect and suggested an alternative or allosteric binding site for the three drugs. Furthermore, the three steroids showed an additive effect on Hh pathway inhibition when tested in combination with cyclopamine. Our study reports the antagonistic effect of dexamethasone, fludrocortisone, and corticosterone on the Hh pathway using functional assay and confirmed that they do not bind to the CRD or adjacent TM binding cavities of SMO. The study also suggests that dexamethasone could be additionally beneficial as the adjuvant therapy for cancer patients with an established link to the dysregulated Hh pathway.

Anandamide Revisited: How Cholesterol and Ceramides Control Receptor-Dependent and Receptor-Independent Signal Transmission Pathways of a Lipid Neurotransmitter

Anandamide is a lipid neurotransmitter derived from arachidonic acid, a polyunsaturated fatty acid. The chemical differences between anandamide and arachidonic acid result in a slightly enhanced solubility in water and absence of an ionisable group for the neurotransmitter compared with the fatty acid. In this review, we first analyze the conformational flexibility of anandamide in aqueous and membrane phases. We next study the interaction of the neurotransmitter with membrane lipids and discuss the molecular basis of the unexpected selectivity of anandamide for cholesterol and ceramide from among other membrane lipids. We show that cholesterol behaves as a binding partner for anandamide, and that following an initial interaction mediated by the establishment of a hydrogen bond, anandamide is attracted towards the membrane interior, where it forms a molecular complex with cholesterol after a functional conformation adaptation to the apolar membrane milieu. The complex is then directed to the anandamide cannabinoid receptor (CB1) which displays a high affinity binding pocket for anandamide. We propose that cholesterol may regulate the entry and exit of anandamide in and out of CB1 by interacting with low affinity cholesterol recognition sites (CARC and CRAC) located in transmembrane helices. The mirror topology of cholesterol binding sites in the seventh transmembrane domain is consistent with the delivery, extraction and flip-flop of anandamide through a coordinated cholesterol-dependent mechanism. The binding of anandamide to ceramide illustrates another key function of membrane lipids which may occur independently of protein receptors. Interestingly, ceramide forms a tight complex with anandamide which blocks the degradation pathway of both lipids and could be exploited for anti-cancer therapies.

A direct interaction of cholesterol with the dopamine transporter prevents its out-to-inward transition

Monoamine transporters (MATs) carry out neurotransmitter reuptake from the synaptic cleft, a key step in neurotransmission, which is targeted in the treatment of neurological disorders. Cholesterol (CHOL), a major component of the synaptic plasma membrane, has been shown to exhibit a modulatory effect on MATs. Recent crystal structures of the dopamine transporter (DAT) revealed the presence of two conserved CHOL-like molecules, suggesting a functional protein-CHOL direct interaction. Here, we present extensive atomistic molecular dynamics (MD) simulations of DAT in an outward-facing conformation. In the absence of bound CHOL, DAT undergoes structural changes reflecting early events of dopamine transport: transition to an inward-facing conformation. In contrast, in the presence of bound CHOL, these conformational changes are inhibited, seemingly by an immobilization of the intracellular interface of transmembrane helix 1a and 5 by CHOL. We also provide evidence, from coarse grain MD simulations that the CHOL sites observed in the DAT crystal structures are preserved in all human monoamine transporters (dopamine, serotonin and norepinephrine), suggesting that our findings might extend to the entire family.

Cholecystokinin responsiveness varies across the population dependent on metabolic phenotype

Background: Cholecystokinin (CCK) is an important satiety factor, acting at type 1 receptors (CCK1Rs) on vagal afferent neurons; however, CCK agonists have failed clinical trials for obesity. We postulated that CCK1R function might be defective in such patients due to abnormal membrane composition, such as that observed in cholesterol gallstone disease.Objective: Due to the challenges in directly studying CCK1Rs relevant to appetite control, our goal was to develop and apply a method to determine the impact of a patient's own cellular environment on CCK stimulus-activity coupling and to determine whether CCK sensitivity correlated with the metabolic phenotype of a high-risk population.Design: Wild-type CCK1Rs were expressed on leukocytes from 112 Hispanic patients by using adenoviral transduction and 24-h culture, with quantitation of cholesterol composition and intracellular calcium responses to CCK. Results were correlated with clinical, biochemical, and morphometric characteristics.Results: Broad ranges of cellular cholesterol and CCK responsiveness were observed, with elevated cholesterol correlated with reduced CCK sensitivity. This was prominent with increasing degrees of obesity and the presence of diabetes, particularly when poorly controlled. No single standard clinical metric correlated directly with CCK responsiveness. Reduced CCK sensitivity best correlated with elevated serum triglycerides in normal-weight participants and with low HDL concentrations and elevated glycated hemoglobin in obese and diabetic patients.Conclusions: CCK responsiveness varies widely across the population, with reduced signaling in patients with obesity and diabetes. This could explain the failure of CCK agonists in previous clinical trials and supports the rationale to develop corrective modulators to reverse this defective servomechanism for appetite control. This trial was registered at www.clinicaltrials.gov as NCT03121755.

Complex biomembrane mimetics on the sub-nanometer scale

Biomimetic lipid vesicles are indispensable tools for gaining insight into the biophysics of cell physiology on the molecular level. The level of complexity of these model systems has steadily increased, and now spans from domain-forming lipid mixtures to asymmetric lipid bilayers. Here, we review recent progress in the development and application of elastic neutron and X-ray scattering techniques for studying these systems in situ and under physiologically relevant conditions on the nanometer to sub-nanometer length scales. In particular, we focus on: (1) structural details of coexisting liquid-ordered and liquid-disordered domains, including their thickness and lipid packing mismatch as a function of a size transition from nanoscopic to macroscopic domains; (2) membrane-mediated protein partitioning into lipid domains; (3) the role of the aqueous medium in tuning interactions between membranes and domains; and (4) leaflet-specific structure in asymmetric bilayers and passive lipid flip-flop.

Differential Membrane Dipolar Orientation Induced by Acute and Chronic Cholesterol Depletion

Cholesterol plays a crucial role in cell membrane organization, dynamics and function. Depletion of cholesterol represents a popular approach to explore cholesterol-sensitivity of membrane proteins. An emerging body of literature shows that the consequence of membrane cholesterol depletion often depends on the actual process (acute or chronic), although the molecular mechanism underlying the difference is not clear. Acute depletion, using cyclodextrin-type carriers, is faster relative to chronic depletion, in which inhibitors of cholesterol biosynthesis are used. With the overall goal of addressing molecular differences underlying these processes, we monitored membrane dipole potential under conditions of acute and chronic cholesterol depletion in CHO-K1 cells, using a voltage-sensitive fluorescent dye in dual wavelength ratiometric mode. Our results show that the observed membrane dipole potential exhibits difference under acute and chronic cholesterol depletion conditions, even when cholesterol content was identical. To the best of our knowledge, these results provide, for the first time, molecular insight highlighting differences in dipolar reorganization in these processes. A comprehensive understanding of processes in which membrane cholesterol gets modulated would provide novel insight in its interaction with membrane proteins and receptors, thereby allowing us to understand the role of cholesterol in cellular physiology associated with health and disease.

Structural Stringency of Cholesterol for Membrane Protein Function Utilizing Stereoisomers as Novel Tools: A Review

Cholesterol is an important lipid in the context of membrane protein function. The function of a number of membrane proteins, including G protein-coupled receptors (GPCRs) and ion channels, has been shown to be dependent on membrane cholesterol. However, the molecular mechanism underlying such regulation is still being explored. In some cases, specific interaction between cholesterol and the protein has been implicated. In other cases, the effect of cholesterol on the membrane properties has been attributed for the regulation of protein function. In this article, we have provided an overview of experimental approaches that are useful for determining the degree of structural stringency of cholesterol for membrane protein function. In the process, we have highlighted the role of immediate precursors in cholesterol biosynthetic pathway in the function of membrane proteins. Special emphasis has been given to the application of stereoisomers of cholesterol in deciphering the structural stringency required for regulation of membrane protein function. A comprehensive examination of these processes would help in understanding the molecular basis of cholesterol regulation of membrane proteins in subtle details.

Mapping Cholesterol Interaction Sites on Serotonin Transporter through Coarse-Grained Molecular Dynamics

Serotonin transporter (SERT) modulates serotonergic signaling via re-uptake of serotonin in pre-synaptic cells. The inclusion in cholesterol-enriched membrane domains is crucial for SERT activity, suggesting a cross-talk between the protein and the sterol. Here, we develop a protocol to identify potential cholesterol interaction sites coupling statistical analysis to multi-microsecond coarse-grained molecular dynamics simulations of SERT in a previously validated raft-like membrane model. Six putative sites were found, including a putative CRAC motif on TM4 and a CARC motif on TM10. Among them, four hot-spots near regions related to ion binding, transport, and inhibition were detected. Our results encourage prospective studies to unravel mechanistic features of the transporter and related drug discovery implications.

Interplay of G Protein-Coupled Receptors with the Membrane: Insights from Supra-Atomic Coarse Grain Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are central to many fundamental cellular signaling pathways. They transduce signals from the outside to the inside of cells in physiological processes ranging from vision to immune response. It is extremely challenging to look at them individually using conventional experimental techniques. Recently, a pseudo atomistic molecular model has emerged as a valuable tool to access information on GPCRs, more specifically on their interactions with their environment in their native cell membrane and the consequences on their supramolecular organization. This approach uses the Martini coarse grain (CG) model to describe the receptors, lipids, and solvent in molecular dynamics (MD) simulations and in enough detail to allow conserving the chemical specificity of the different molecules. The elimination of unnecessary degrees of freedom has opened up large-scale simulations of the lipid-mediated supramolecular organization of GPCRs. Here, after introducing the Martini CGMD method, we review these studies carried out on various members of the GPCR family, including rhodopsin (visual receptor), opioid receptors, adrenergic receptors, adenosine receptors, dopamine receptor, and sphingosine 1-phosphate receptor. These studies have brought to light an interesting set of novel biophysical principles. The insights range from revealing localized and heterogeneous deformations of the membrane bilayer at the surface of the protein, specific interactions of lipid molecules with individual GPCRs, to the effect of the membrane matrix on global GPCR self-assembly. The review ends with an overview of the lessons learned from the use of the CGMD method, the biophysical-chemical findings on lipid-protein interplay.

Cholesterol-dependent Conformational Plasticity in GPCR Dimers

The organization and function of the serotonin_1A receptor, an important member of the GPCR family, have been shown to be cholesterol-dependent, although the molecular mechanism is not clear. We performed a comprehensive structural and dynamic analysis of dimerization of the serotonin_1A receptor by coarse-grain molecular dynamics simulations totaling 3.6 ms to explore the molecular details of its cholesterol-dependent association. A major finding is that the plasticity and flexibility of the receptor dimers increase with increased cholesterol concentration. In particular, a dimer interface formed by transmembrane helices I-I was found to be sensitive to cholesterol. The modulation of dimer interface appears to arise from a combination of direct cholesterol occupancy and indirect membrane effects. Interestingly, the presence of cholesterol at the dimer interface is correlated with increased dimer plasticity and flexibility. These results represent an important step in characterizing the molecular interactions in GPCR organization with potential relevance to therapeutic interventions.

Membrane lipid rafts and neurobiology: age-related changes in membrane lipids and loss of neuronal function

A better understanding of the cellular physiological role that plasma membrane lipids, fatty acids and sterols play in various cellular systems may yield more insight into how cellular and whole organ function is altered during the ageing process. Membrane lipid rafts (MLRs) within the plasma membrane of most cells serve as key organizers of intracellular signalling and tethering points of cytoskeletal components. MLRs are plasmalemmal microdomains enriched in sphingolipids, cholesterol and scaffolding proteins; they serve as a platform for signal transduction, cytoskeletal organization and vesicular trafficking. Within MLRs are the scaffolding and cholesterol binding proteins named caveolin (Cav). Cavs not only organize a multitude of receptors including neurotransmitter receptors (NMDA and AMPA receptors), signalling proteins that regulate the production of cAMP (G protein-coupled receptors, adenylyl cyclases, phosphodiesterases (PDEs)), and receptor tyrosine kinases involved in growth (Trk), but also interact with components that modulate actin and tubulin cytoskeletal dynamics (e.g. RhoGTPases and actin binding proteins). MLRs are essential for the regulation of the physiology of organs such as the brain, and age-related loss of cholesterol from the plasma membrane leads to loss of MLRs, decreased presynaptic vesicle fusion, and changes in neurotransmitter release, all of which contribute to different forms of neurodegeneration. Thus, MLRs provide an active membrane domain that tethers and reorganizes the cytoskeletal machinery necessary for membrane and cellular repair, and genetic interventions that restore MLRs to normal cellular levels may be exploited as potential therapeutic means to reverse the ageing and neurodegenerative processes.

Can Specific Protein-Lipid Interactions Stabilize an Active State of the Beta 2 Adrenergic Receptor?

G-protein-coupled receptors are eukaryotic membrane proteins with broad biological and pharmacological relevance. Like all membrane-embedded proteins, their location and orientation are influenced by lipids, which can also impact protein function via specific interactions. Extensive simulations totaling 0.25 ms reveal a process in which phospholipids from the membrane's cytosolic leaflet enter the empty G-protein binding site of an activated β2 adrenergic receptor and form salt-bridge interactions that inhibit ionic lock formation and prolong active-state residency. Simulations of the receptor embedded in an anionic membrane show increased lipid binding, providing a molecular mechanism for the experimental observation that anionic lipids can enhance receptor activity. Conservation of the arginine component of the ionic lock among Rhodopsin-like G-protein-coupled receptors suggests that intracellular lipid ingression between receptor helices H6 and H7 may be a general mechanism for active-state stabilization.

Lateral pressure-mediated protein partitioning into liquid-ordered/liquid-disordered domains

We have studied the contributions of stored elastic energies in liquid-ordered (Lo) and liquid-disordered (Ld) domains to transmembrane proteins using the lateral pressure concept. In particular we applied previously reported experimental data for the membrane thickness, intrinsic curvature and bending elasticities of coexisting Lo/Ld domains to calculate whether proteins of simple geometric shapes would preferentially diffuse into Lo or Ld domains and form oligomers of a certain size. For the studied lipid mixture we generally found that proteins with convex shapes prefer sorting to Ld phases and the formation of large clusters. Lo domains in turn would be enriched in monomers of concave shaped proteins. We further observed that proteins which are symmetric with respect to the bilayer center prefer symmetric Lo or Ld domains, while asymmetric proteins favor a location in domains with Lo/Ld asymmetry. In the latter case we additionally retrieved a strong dependence on protein directionality, thus providing a mechanism for transmembrane protein orientation.

Cholesterol under oxidative stress-How lipid membranes sense oxidation as cholesterol is being replaced by oxysterols

The behavior of oxysterols in phospholipid membranes and their effects on membrane properties were investigated by means of dynamic light scattering, fluorescence spectroscopy, NMR, and extensive atomistic simulations. Two families of oxysterols were scrutinized-tail-oxidized sterols, which are mostly produced by enzymatic processes, and ring-oxidized sterols, formed mostly via reactions with free radicals. The former family of sterols was found to behave similar to cholesterol in terms of molecular orientation, roughly parallel to the bilayer normal, leading to increasing membrane stiffness and suppression of its membrane permeability. In contrast, ring-oxidized sterols behave quantitatively differently from cholesterol. They acquire tilted orientations and therefore disrupt the bilayer structure with potential implications for signaling and other biochemical processes in the membranes.

Unfolding and refolding of a protein by cholesterol and cyclodextrin: a single molecule study

Cholesterol induced unfolding of a globular protein, human serum albumin (HSA), and β-cyclodextrin induced refolding of the unfolded protein is demonstrated in this study.