Interfacial Binding Sites for Cholesterol on G Protein-Coupled Receptors

Bayesian Nonparametric Analysis of Residence Times for Protein-Lipid Interactions in Molecular Dynamics Simulations

Molecular Dynamics (MD) simulations are a versatile tool to investigate the interactions of proteins within their environments, in particular, of membrane proteins with the surrounding lipids. However, quantitative analysis of lipid-protein binding kinetics has remained challenging due to considerable noise and low frequency of long binding events, even in hundreds of microseconds of simulation data. Here, we apply Bayesian nonparametrics to compute residue-resolved residence time distributions from MD trajectories. Such an analysis characterizes binding processes at different time scales (quantified by their kinetic off-rate) and assigns to each trajectory frame a probability of belonging to a specific process. In this way, we classify trajectory frames in an unsupervised manner and obtain, for example, different binding poses or molecular densities based on the time scale of the process. We demonstrate our approach by characterizing interactions of cholesterol with six different G-protein-coupled receptors (A2AAR, β2AR, CB1R, CB2R, CCK1R, and CCK2R) simulated with coarse-grained MD simulations with the MARTINI model. The nonparametric Bayesian analysis allows us to connect the coarse binding time series data to the underlying molecular picture and thus not only infers accurate binding kinetics with error distributions from MD simulations but also describes molecular events responsible for the broad range of kinetic rates.

Bayesian nonparametric analysis of residence times for protein-lipid interactions in Molecular Dynamics simulations

Molecular Dynamics (MD) simulations are a versatile tool to investigate the interactions of proteins within their environments, in particular of membrane proteins with the surrounding lipids. However, quantitative analysis of lipid-protein binding kinetics has remained challenging due to considerable noise and low frequency of long binding events, even in hundreds of microseconds of simulation data. Here we apply Bayesian nonparametrics to compute residue-resolved residence time distributions from MD trajectories. Such an analysis characterizes binding processes at different timescales (quantified by their kinetic off-rate) and assigns to each trajectory frame a probability of belonging to a specific process. In this way, we classify trajectory frames in an unsupervised manner and obtain, for example, different binding poses or molecular densities based on the timescale of the process. We demonstrate our approach by characterizing interactions of cholesterol with six different G-protein coupled receptors (A 2 A AR ,β 2 AR ,CB 1 R ,CB 2 R ,CCK 1 R ,CCK 2 R ) simulated with coarse-grained MD simulations with the MARTINI model. The nonparametric Bayesian analysis allows us to connect the coarse binding time series data to the underlying molecular picture and, thus, not only infers accurate binding kinetics with error distributions from MD simulations but also describes molecular events responsible for the broad range of kinetic rates.

The Role of Cholesterol in M2 Clustering and Viral Budding Explained

The influenza A M2 homotetrameric channel consists of four transmembrane (TM) and four amphipathic helices (AHs). This viral proton channel is suggested to form clusters in the catenoid budding neck areas in raft-like domains of the plasma membrane, resulting in cell membrane scission and viral release. The channel clustering environment is rich in cholesterol. Previous experiments have shown that cholesterol significantly contributes to lipid bilayer undulations in viral buds. However, a clear explanation of membrane curvature from the distribution of cholesterol around the M2TM-AH clusters is lacking. Using coarse-grained molecular dynamics simulations of M2TM-AH in bilayers, we observed that M2 channels form specific, C2-symmetric, clusters with conical shapes driven by the attraction of their AHs. We showed that cholesterol stabilized the formation of M2 channel clusters by filling and bridging the conical gap between M2 channels at specific sites in the N-termini of adjacent channels or via the C-terminal region of TM and AHs, with the latter sites displaying a longer interaction time and higher stability. The potential of mean force calculations showed that when cholesterols occupy the identified interfacial binding sites between two M2 channels, the dimer is stabilized by 11 kJ/mol. This translates to the cholesterol-bound dimer being populated by almost 2 orders of magnitude compared to a dimer lacking cholesterol. We demonstrated that the cholesterol-bridged M2 channels can exert a lateral force on the surrounding membrane to induce the necessary negative Gaussian curvature profile, which permits spontaneous scission of the catenoid membrane neck and leads to viral buds and scission.

The pleomorphic cholesterol sensing motifs of transmembrane proteins

Millions of years of phylogenetic evolution have shaped the crosstalk between sterols and membrane-embedded proteins. This lengthy process, which began before the appearance of eukaryotic cells, has sculpted the two types of molecules to cover a wide spectrum of structural interconnectedness, ranging from rapid touch-and-go hits of low-affinity between surfaces to stronger lock-and-key type structural contacts. The former usually involve relatively loose contacts between linear amino acid sequences on the membrane-exposed transmembrane domains of the protein, readily accessible to the sterols as they briefly visit clefts between adjacent transmembrane segments while in rapid exchange with the bulk lipid bilayer. This operational mode is probably the most ancestral one, since it was already present in primitive bacteria interacting with hopanoid lipids. At the other end of this spectrum are more complex cholesterol binding sites that have required the acquisition of complex 3D non-sequential segments of the membrane protein to establish stereochemically elaborate 3D designs complementary to the rough and smooth surfaces of the eukaryotic neutral lipid, cholesterol. This short review explores cholesterol-membrane protein interactions using membrane protein paradigms having in common their participation in intercellular communications neurotransmission, hormone signalling, amino acid/neurotransmitter transport- and in cancer.

In Silico Tools to Score and Predict Cholesterol-Protein Interactions

Cholesterol is structurally distinct from other lipids, which confers it with singular roles in membrane organization and protein function. As a signaling molecule, cholesterol engages in discrete interactions with transmembrane, peripheral, and certain soluble proteins to control cellular responses. Accordingly, the cholesterol-protein interface is central to cholesterol-related diseases and is an essential consideration in drug design. However, cholesterol's hydrophobic, un-drug-like nature presents a unique challenge to traditional in silico analyses. In this Perspective, we survey a collection of tools designed to predict and evaluate cholesterol binding sites in proteins, including classical sequence motifs, molecular docking, template-based strategies, molecular dynamics simulations, and recent artificial intelligence approaches. We then comment on contemporary tools to evaluate ligand-protein interactions, their applicability to cholesterol, and the yet-untapped potential of cholesterol-protein interactions in human health and disease.

Residues of TRPM8 at the Lipid-Water-Interface have Coevolved with Cholesterol Interaction and are Relevant for Diverse Health Disorders

TRPM8 is a non-selective cation channel that is expressed in several tissues and cells and also has a unique property to be activated by low-temperature. In this work, we have analyzed the conservation of amino acids that are present in the lipid-water-interface (LWI) region of TRPM8, the region which experiences a microenvironment near the membrane surface. We demonstrate that the amino acids present in the LWI region are more conserved than the transmembrane or even full-length TRPM8, suggesting strong selection pressure in these residues. TRPM8 also has several conserved cholesterol-binding motifs where cholesterol can bind in different modes and energies. We suggest that mutations and/or physiological conditions can potentially alter these TRPM8-cholesterol complexes and can lead to physiological disorders or even apparently irreversible diseases such as cancer and neurodegeneration.

Differential Behavior of Conformational Dynamics in Active and Inactive States of Cannabinoid Receptor 1

Cannabinoid receptor 1 (CB1) is a G protein-coupled receptor that regulates critical physiological processes including pain, appetite, and cognition. Understanding the conformational dynamics of CB1 associated with transitions between inactive and active signaling states is imperative for developing targeted modulators. Using microsecond-level all-atom molecular dynamics simulations, we identified marked differences in the conformational ensembles of inactive and active CB1 in apo. The inactive state exhibited substantially increased structural heterogeneity and plasticity compared to the more rigidified active state in the absence of stabilizing ligands. Transmembrane helices TM3 and TM7 were identified as distinguishing factors modulating the state-dependent dynamics. TM7 displayed amplified fluctuations selectively in the inactive state simulations attributed to disruption of conserved electrostatic contacts anchoring it to surrounding helices in the active state. Additionally, we identified significant reorganizations in key salt bridge and hydrogen bond networks contributing to the CB1 activation/inactivation. For instance, D213-Y224 hydrogen bond and D184-K192 salt bridge showed marked rearrangements between the states. Collectively, these findings reveal the specialized role of TM7 in directing state-dependent CB1 dynamics through electrostatic switch mechanisms. By elucidating the intrinsic enhanced flexibility of inactive CB1, this study provides valuable insights into the conformational landscape enabling functional transitions. Our perspective advances understanding of CB1 activation mechanisms and offers opportunities for structure-based drug discovery targeting the state-specific conformational dynamics of this receptor.

Impact of aging on the GABAB receptor-mediated connectome

GABA B receptors (GABABRs) are heterodimeric seven-transmembrane receptors that interact with a range of proteins and form large protein complexes on cholesterol-rich membrane microdomains. As the brain ages, membrane cholesterol levels exhibit alterations, although it remains unclear how these changes impact protein-protein interactions and downstream signaling. Herein, we studied the structural bases for the interaction between GABABR and the KCC2 transporter, including their protein expression and distribution, and we compared data between young and aged rat cerebella. Also, we analyzed lipid profiles for both groups, and we used molecular dynamics simulations on three plasma membrane systems with different cholesterol concentrations, to further explore the GABABR-transporter interaction. Based on our results, we report that a significant decrease in GABAB2 subunit expression occurs in the aged rat cerebella. After performing a comparative co-immunoprecipitation analysis, we confirm that GABABR and KCC2 form a protein complex in adult and aged rat cerebella, although their interaction levels are reduced substantially as the cerebellum ages. On the other hand, our lipid analyses reveal a significant increase in cholesterol and sphingomyelin levels of the aged cerebella. Finally, we used the Martini coarse-grained model to conduct molecular dynamics simulations, from which we observed that membrane cholesterol concentrations can dictate whether the GABABR tail domains physically establish G protein-independent contacts with a transporter, and the timing when those associations eventually occur. Taken together, our findings illustrate how age-related alterations in membrane cholesterol levels affect protein-protein interactions, and how they could play a crucial role in regulating GABABR's interactome-mediated signaling.

Significance Statement

This study elucidates age-related changes in cerebellar GABAB receptors (GABABRs), KCC2, and plasma membrane lipids, shedding light on mechanisms underlying neurological decline. Molecular dynamics simulations reveal how membrane lipids influence protein-protein interactions, offering insights into age-related neurodegeneration. The findings underscore the broader impact of cerebellar aging on motor functions, cognition, and emotional processing in the elderly. By elucidating plasma membrane regulation and GABAergic dynamics, this research lays the groundwork for understanding aging-related neurological disorders and inspires further investigation into therapeutic interventions.

Differential Behavior of Conformational Dynamics in Active and Inactive States of Cannabinoid Receptor 1

The cannabinoid receptor CB1 is a G protein-coupled receptor that regulates critical physiological processes including pain, appetite, and cognition. Understanding the conformational dynamics of CB1 associated with transitions between inactive and active signaling states is imperative for developing targeted modulators. Using microsecond-level all-atom molecular dynamics simulations, we identified marked differences in the conformational ensembles of inactive and active CB1 states in apo conditions. The inactive state exhibited substantially increased structural heterogeneity and plasticity compared to the more rigidified active state in the absence of stabilizing ligands. Transmembrane helices TM3 and TM7 were identified as distinguishing factors modulating the state-dependent dynamics. TM7 displayed amplified fluctuations selectively in the inactive state simulations attributed to disruption of conserved electrostatic contacts anchoring it to surrounding helices in the active state. Additionally, we identified significant reorganization of key salt bridge and hydrogen bond networks known to control CB1 activation between states. For instance, a conserved D213-Y224 hydrogen bond and D184-K192 salt bridge interactions showed marked rearrangements between the states. Collectively, these findings reveal the specialized role of TM7 in directing state-dependent CB1 dynamics through electrostatic switch mechanisms. By elucidating the intrinsic enhanced flexibility of inactive CB1, this study provides valuable insights into the conformational landscape enabling functional transitions. Our perspective advances understanding of CB1 activation mechanisms and offers opportunities for structure-based drug discovery targeting the state-specific conformational dynamics of this receptor.

LipIDens: simulation assisted interpretation of lipid densities in cryo-EM structures of membrane proteins

Cryo-electron microscopy (cryo-EM) enables the determination of membrane protein structures in native-like environments. Characterising how membrane proteins interact with the surrounding membrane lipid environment is assisted by resolution of lipid-like densities visible in cryo-EM maps. Nevertheless, establishing the molecular identity of putative lipid and/or detergent densities remains challenging. Here we present LipIDens, a pipeline for molecular dynamics (MD) simulation-assisted interpretation of lipid and lipid-like densities in cryo-EM structures. The pipeline integrates the implementation and analysis of multi-scale MD simulations for identification, ranking and refinement of lipid binding poses which superpose onto cryo-EM map densities. Thus, LipIDens enables direct integration of experimental and computational structural approaches to facilitate the interpretation of lipid-like cryo-EM densities and to reveal the molecular identities of protein-lipid interactions within a bilayer environment. We demonstrate this by application of our open-source LipIDens code to ten diverse membrane protein structures which exhibit lipid-like densities.

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.

Docking cholesterol to integral membrane proteins with Rosetta

Lipid molecules such as cholesterol interact with the surface of integral membrane proteins (IMP) in a mode different from drug-like molecules in a protein binding pocket. These differences are due to the lipid molecule's shape, the membrane's hydrophobic environment, and the lipid's orientation in the membrane. We can use the recent increase in experimental structures in complex with cholesterol to understand protein-cholesterol interactions. We developed the RosettaCholesterol protocol consisting of (1) a prediction phase using an energy grid to sample and score native-like binding poses and (2) a specificity filter to calculate the likelihood that a cholesterol interaction site may be specific. We used a multi-pronged benchmark (self-dock, flip-dock, cross-dock, and global-dock) of protein-cholesterol complexes to validate our method. RosettaCholesterol improved sampling and scoring of native poses over the standard RosettaLigand baseline method in 91% of cases and performs better regardless of benchmark complexity. On the β2AR, our method found one likely-specific site, which is described in the literature. The RosettaCholesterol protocol quantifies cholesterol binding site specificity. Our approach provides a starting point for high-throughput modeling and prediction of cholesterol binding sites for further experimental validation.

The role of cholesterol binding in the control of cholesterol by the Scap-Insig system

Scap and Insig, two proteins embedded in the membrane of the endoplasmic reticulum (ER), regulate the synthesis of cholesterol in animal cells by forming a dimer in the presence of high concentrations of cholesterol. Cryo-electron microscopic structures for the Scap-Insig dimer show a sterol-binding site at the dimer interface, but none of the structures include cholesterol itself. Here, a molecular docking approach developed to characterise cholesterol binding to the transmembrane (TM) regions of membrane proteins is used to characterise cholesterol binding to sites on the TM surface of the dimer and to the interfacial binding site. Binding of cholesterol is also observed at sites on the extra-membranous luminal domains of Scap, but the properties of these sites suggest that they will be unoccupied in vivo. Comparing the structure of Scap in the dimer with that predicted by AlphaFold for monomeric Scap suggests that dimer formation could result in relocation of TM helix 7 of Scap and of the loop between TM6 and 7, and that this could be the key change on Scap that signals that there is a high concentration of cholesterol in the ER.

Exploring TRPC3 Interaction with Cholesterol through Coarse-Grained Molecular Dynamics Simulations

Transient receptor potential canonical 3 (TRPC3) channel belongs to the superfamily of transient receptor potential (TRP) channels which mediate Ca2+ influx into the cell. These channels constitute essential elements of cellular signalling and have been implicated in a wide range of diseases. TRPC3 is primarily gated by lipids and its surface expression has been shown to be dependent on cholesterol, yet a comprehensive exploration of its interaction with this lipid has thus far not emerged. Here, through 80 µs of coarse-grained molecular dynamics simulations, we show that cholesterol interacts with multiple elements of the transmembrane machinery of TRPC3. Through our approach, we identify an annular binding site for cholesterol on the pre-S1 helix and a non-annular site at the interface between the voltage-sensor-like domain and pore domains. Here, cholesterol interacts with exposed polar residues and possibly acts to stabilise the domain interface.

Endothelial Cell Plasma Membrane Biomechanics Mediates Effects of Pro-Inflammatory Factors on Endothelial Mechanosensors: Vicious Circle Formation in Atherogenic Inflammation

Chronic low-grade vascular inflammation and endothelial dysfunction significantly contribute to the pathogenesis of cardiovascular diseases. In endothelial cells (ECs), anti-inflammatory or pro-inflammatory signaling can be induced by different patterns of the fluid shear stress (SS) exerted by blood flow on ECs. Laminar blood flow with high magnitude is anti-inflammatory, while disturbed flow and laminar flow with low magnitude is pro-inflammatory. Endothelial mechanosensors are the key upstream signaling proteins in SS-induced pro- and anti-inflammatory responses. Being transmembrane proteins, mechanosensors, not only experience fluid SS but also become regulated by the biomechanical properties of the lipid bilayer and the cytoskeleton. We review the apparent effects of pro-inflammatory factors (hypoxia, oxidative stress, hypercholesterolemia, and cytokines) on the biomechanics of the lipid bilayer and the cytoskeleton. An analysis of the available data suggests that the formation of a vicious circle may occur, in which pro-inflammatory cytokines enhance and attenuate SS-induced pro-inflammatory and anti-inflammatory signaling, respectively.

Insights into lipid-protein interactions from computer simulations

Lipid-protein interactions play an important direct role in the function of many membrane proteins. We argue they are key players in membrane structure, modulate membrane proteins in more subtle ways than direct binding, and are important for understanding the mechanism of classes of hydrophobic drugs. By directly comparing membrane proteins from different families in the same, complex lipid mixture, we found a unique lipid environment for every protein. Extending this work, we identified both differences and similarities in the lipid environment of GPCRs, dependent on which family they belong to and in some cases their conformational state, with particular emphasis on the distribution of cholesterol. More recently, we have been studying modes of coupling between protein conformation and local membrane properties using model proteins. In more applied approaches, we have used similar methods to investigate specific hypotheses on interactions of lipid and lipid-like molecules with ion channels. We conclude this perspective with some considerations for future work, including a new more sophisticated coarse-grained force field (Martini 3), an interactive visual exploration framework, and opportunities to improve sampling.

Recent Advances in Structure, Function, and Pharmacology of Class A Lipid GPCRs: Opportunities and Challenges for Drug Discovery

Great progress has been made over the past decade in understanding the structural, functional, and pharmacological diversity of lipid GPCRs. From the first determination of the crystal structure of bovine rhodopsin in 2000, much progress has been made in the field of GPCR structural biology. The extraordinary progress in structural biology and pharmacology of GPCRs, coupled with rapid advances in computational approaches to study receptor dynamics and receptor-ligand interactions, has broadened our comprehension of the structural and functional facets of the receptor family members and has helped usher in a modern age of structure-based drug design and development. First, we provide a primer on lipid mediators and lipid GPCRs and their role in physiology and diseases as well as their value as drug targets. Second, we summarize the current advancements in the understanding of structural features of lipid GPCRs, such as the structural variation of their extracellular domains, diversity of their orthosteric and allosteric ligand binding sites, and molecular mechanisms of ligand binding. Third, we close by collating the emerging paradigms and opportunities in targeting lipid GPCRs, including a brief discussion on current strategies, challenges, and the future outlook.

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.

Evidence that specific interactions play a role in the cholesterol sensitivity of G protein-coupled receptors

G protein-coupled receptors (GPCRs) are known to be modulated by membrane cholesterol levels, but whether or not the effects are caused by specific receptor-cholesterol interactions or cholesterol's general effects on the membrane is not well-understood. We performed coarse-grained molecular dynamics (CGMD) simulations coupled with structural bioinformatics approaches on the β2-adrenergic receptor (β2AR) and the cholecystokinin (CCK) receptor subfamily. The β2AR has been shown to be sensitive to membrane cholesterol and cholesterol molecules have been clearly resolved in numerous β2AR crystal structures. The two CCK receptors are highly homologous and preserve similar cholesterol recognition motifs but despite their homology, CCK1R shows functional sensitivity to membrane cholesterol while CCK2R does not. Our results offer new insights into how cholesterol modulates GPCR function by showing cholesterol interactions with β2AR that agree with previously published data; additionally, we observe differential and specific cholesterol binding in the CCK receptor subfamily while revealing a previously unreported Cholesterol Recognition Amino-acid Consensus (CRAC) sequence that is also conserved across 38% of class A GPCRs. A thermal denaturation assay (LCP-Tm) shows that mutation of a conserved CRAC sequence on TM7 of the β2AR affects cholesterol stabilization of the receptor in a lipid bilayer. The results of this study provide a better understanding of receptor-cholesterol interactions that can contribute to novel and improved therapeutics for a variety of diseases.

Structural insights into GIRK2 channel modulation by cholesterol and PIP2

G-protein-gated inwardly rectifying potassium (GIRK) channels are important for determining neuronal excitability. In addition to G proteins, GIRK channels are potentiated by membrane cholesterol, which is elevated in the brains of people with neurodegenerative diseases such as Alzheimer’s dementia and Parkinson’s disease. The structural mechanism of cholesterol modulation of GIRK channels is not well understood. In this study, we present cryo-electron microscopy (cryoEM) structures of GIRK2 in the presence and absence of the cholesterol analog cholesteryl hemisuccinate (CHS) and phosphatidylinositol 4,5-bisphosphate (PIP₂). The structures reveal that CHS binds near PIP₂ in lipid-facing hydrophobic pockets of the transmembrane domain. Our structural analysis suggests that CHS stabilizes PIP₂ interaction with the channel and promotes engagement of the cytoplasmic domain onto the transmembrane region. Mutagenesis of one of the CHS binding pockets eliminates cholesterol-dependent potentiation of GIRK2. Elucidating the structural mechanisms underlying cholesterol modulation of GIRK2 channels could facilitate the development of therapeutics for treating neurological diseases.

Ion channels are important in determining neuronal excitability. Elevated cholesterol levels found in some neurodegenerative diseases can affect the function of ion channels. Mathiharan et al. take a structural and functional approach to identifying physical sites for cholesterol, and they provide details on how cholesterol potentiates ion channel activity.

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.

Interfacial binding sites for cholesterol on GABAA receptors and competition with neurosteroids

γ-Aminobutyric acid type A (GABA_A) receptors in the brain are located in the outer membranes of brain cells where the concentration of cholesterol is high. Of the 25 available high-resolution structures available for GABA_A receptors, none were determined in the presence of cholesterol, but four include resolved molecules of cholesterol hemisuccinate (CHS). Here, a molecular docking procedure is used to sweep the transmembrane (TM) surfaces of the receptors for cholesterol binding sites. Cholesterol docking poses determined in this way match 89% of the resolved CHS when CHS molecules deemed unlikely to represent typical bound cholesterols are excluded. The receptors are pentameric, and their TM surfaces consist of a set of five facets, each including pairs of TM helices from two adjacent subunits. Each facet contains hydrophobic hollows running from one side of the membrane to the other, within which are six potential binding sites for cholesterol, three on each side of the membrane. High-resolution structures of GABA_A receptors with bound neurosteroids show that neurosteroids bind in these cholesterol binding sites, so the binding of neurosteroids and cholesterol will be competitive.

Allosteric Modulation of GPCRs of Class A by Cholesterol

G-protein coupled receptors (GPCRs) are membrane proteins that convey extracellular signals to the cellular milieu. They represent a target for more than 30% of currently marketed drugs. Here we review the effects of membrane cholesterol on the function of GPCRs of Class A. We review both the specific effects of cholesterol mediated via its direct high-affinity binding to the receptor and non-specific effects mediated by cholesterol-induced changes in the properties of the membrane. Cholesterol binds to many GPCRs at both canonical and non-canonical binding sites. It allosterically affects ligand binding to and activation of GPCRs. Additionally, it changes the oligomerization state of GPCRs. In this review, we consider a perspective of the potential for the development of new therapies that are targeted at manipulating the level of membrane cholesterol or modulating cholesterol binding sites on to GPCRs.

Impact of Cholesterol Concentration and Lipid Phase on Structure and Fluctuation of Amyloid Precursor Protein

Elevated levels of cellular cholesterol have been identified as one factor contributing to the onset of Alzheimer's disease (AD). Specific interaction between cholesterol and the amyloid precursor protein (APP), investigated via NMR experiments and computational studies, has been proposed to play a critical role in the processing of APP by secretases and the biogenesis of amyloid-β (Aβ) protein. We present all-atom molecular dynamics simulations of the 40-residue congener of the C-terminal domain of APP, C9916-55 (C99), in cholesterol-enriched DMPC lipid bilayers. We investigated the effect of cholesterol concentration on the conformational ensemble of wild-type C99 and C99-cholesterol associations at the low pH of endosomal environments, at which residues E22 and D23 are neutral. C99 was also characterized in liquid ordered domains for Dutch (E22Q) and Iowa (D23N) Familial AD mutants at low pH and for the wild-type sequence using protonation states characteristic of neutral pH. Our results reproduce the equilibrium constant of past NMR characterizations of the C99-cholesterol interaction but are not consistent with the C99-cholesterol binding hypothesis. We find that the lifetimes of both DMPC and cholesterol complexed with C99 display a power-law distribution of residence lifetimes. Longer-lived C99-DMPC and C99-cholesterol complexes are primarily stabilized by salt bridges and hydrogen bonds of lysine amines to phosphate and hydroxyl groups. Nevertheless, specific interfaces for C99-cholesterol association which are not present for DMPC can be identified. Changes to C99-cholesterol interfaces are found to depend on C99 tilt angle and orientation of the juxtamembrane domain of C99 containing residues E22 and D23. These observations support a more nuanced view of the C99-cholesterol interaction than has previously been suggested. We propose that cholesterol modulates the conformation and activity of C99 and other small transmembrane proteins indirectly through induction of the liquid ordered phase and directly through hydrogen bonding. This suggests a critical role for membrane heterogeneity introduced by cholesterol in modulating the structural ensemble of C99 and the production of Aβ.

Cholesterol Localization around the Metabotropic Glutamate Receptor 2

The metabotropic glutamate receptor (mGluR) 2 plays a key role in the central nervous system. mGluR2 has been shown to be regulated by its surrounding lipid environment, especially by cholesterol, by an unknown mechanism. Here, using a combination of biochemical approaches, photo-cross-linking experiments, and molecular dynamics simulations we show the interaction of cholesterol with at least two, but potentially five more, preferential sites on the mGluR2 transmembrane domain. Our simulations demonstrate that surface matching, rather than electrostatic interactions with specific amino acids, is the main factor defining cholesterol localization. Moreover, the cholesterol localization observed here is similar to the sterol-binding pattern previously described in silico for other members of the mGluR family. Biochemical assays suggest little influence of cholesterol on trafficking or dimerization of mGluR2. Nevertheless, simulations revealed a significant reduction of residue-residue contacts together with an alteration in the internal mechanical stress at the cytoplasmic side of the helical bundle when cholesterol was present in the membrane. These alterations may be related to destabilization of the basal state of mGluR2. Due to the high sequence conservation of the transmembrane domains of mGluRs, the molecular interaction of cholesterol and mGluR2 described here is also likely to be relevant for other members of the mGLuR family.

Interfacial Binding Sites for Cholesterol on Kir, Kv, K2P, and Related Potassium Channels

Inwardly rectifying, voltage-gated, two-pore domain, and related K+ channels are located in eukaryotic membranes rich in cholesterol. Here, molecular docking is used to detect specific binding sites ("hot spots") for cholesterol on K+ channels with characteristics that match those of known cholesterol binding sites. The transmembrane surfaces of all available high-resolution structures for K+ channels were swept for potential binding sites. Cholesterol poses were found to be located largely in hollows between protein ridges. A comparison between cholesterol poses and resolved phospholipids suggests that not all cholesterol molecules binding to the transmembrane surface of a K+ channel will result in displacement of a phospholipid molecule from the surface. Competition between cholesterol binding and binding of anionic phospholipids essential for activity could explain some of the effects of cholesterol on channel function.

Lipid-Protein Interactions Are a Unique Property and Defining Feature of G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are membrane-bound proteins that depend on their lipid environment to carry out their physiological function. Combined efforts from many theoretical and experimental studies on the lipid-protein interaction profile of several GPCRs hint at an intricate relationship of these receptors with their surrounding membrane environment, with several lipids emerging as particularly important. Using coarse-grained molecular dynamics simulations, we explore the lipid-protein interaction profiles of 28 different GPCRs, spanning different levels of classification and conformational states and totaling to 1 ms of simulation time. We find a close relationship with lipids for all GPCRs simulated, in particular, cholesterol and phosphatidylinositol phosphate (PIP) lipids, but the number, location, and estimated strength of these interactions is dependent on the specific GPCR as well as its conformational state. Although both cholesterol and PIP lipids bind specifically to GPCRs, they utilize distinct mechanisms. Interactions with PIP lipids are mediated by charge-charge interactions with intracellular loop residues and stabilized by one or both of the transmembrane helices linked by the loop. Interactions with cholesterol, on the other hand, are mediated by a hydrophobic environment, usually made up of residues from more than one helix, capable of accommodating its ring structure and stabilized by interactions with aromatic and charged/polar residues. Cholesterol binding to GPCRs occurs in a small number of sites, some of which (like the binding site on the extracellular side of transmembrane 6/7) are shared among many class A GPCRs. Combined with a thorough investigation of the local membrane structure, our results provide a detailed picture of GPCR-lipid interactions. Additionally, we provide an accompanying website to interactively explore the lipid-protein interaction profile of all GPCRs simulated to facilitate analysis and comparison of our data.

The energetics of protein-lipid interactions as viewed by molecular simulations

Membranes are formed from a bilayer containing diverse lipid species with which membrane proteins interact. Integral, membrane proteins are embedded in this bilayer, where they interact with lipids from their surroundings, whilst peripheral membrane proteins bind to lipids at the surface of membranes. Lipid interactions can influence the function of membrane proteins, either directly or allosterically. Both experimental (structural) and computational approaches can reveal lipid binding sites on membrane proteins. It is, therefore, important to understand the free energies of these interactions. This affords a more complete view of the engagement of a particular protein with the biological membrane surrounding it. Here, we describe many computational approaches currently in use for this purpose, including recent advances using both free energy and unbiased simulation methods. In particular, we focus on interactions of integral membrane proteins with cholesterol, and with anionic lipids such as phosphatidylinositol 4,5-bis-phosphate and cardiolipin. Peripheral membrane proteins are exemplified via interactions of PH domains with phosphoinositide-containing membranes. We summarise the current state of the field and provide an outlook on likely future directions of investigation.

Current Advances in Allosteric Modulation of Muscarinic Receptors

Allosteric modulators are ligands that bind to a site on the receptor that is spatially separated from the orthosteric binding site for the endogenous neurotransmitter. Allosteric modulators modulate the binding affinity, potency, and efficacy of orthosteric ligands. Muscarinic acetylcholine receptors are prototypical allosterically-modulated G-protein-coupled receptors. They are a potential therapeutic target for the treatment of psychiatric, neurologic, and internal diseases like schizophrenia, Alzheimer’s disease, Huntington disease, type 2 diabetes, or chronic pulmonary obstruction. Here, we reviewed the progress made during the last decade in our understanding of their mechanisms of binding, allosteric modulation, and in vivo actions in order to understand the translational impact of studying this important class of pharmacological agents. We overviewed newly developed allosteric modulators of muscarinic receptors as well as new spin-off ideas like bitopic ligands combining allosteric and orthosteric moieties and photo-switchable ligands based on bitopic agents.

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.

Interfacial Binding Sites for Cholesterol on TRP Ion Channels

Transient receptor potential (TRP) channels are members of a large family of ion channels located in membranes rich in cholesterol, some of whose functions are affected by the cholesterol content of the membrane. Here, cholesterol binding to TRPs is studied using a docking procedure that allows the transmembrane surface of a TRP to be swept rapidly for potential binding sites at the interfaces on the two sides of the membrane. Cholesterol docking poses determined in this way match 89% of the cholesterol hemisuccinate molecules in published TRP structures when cholesterol hemisuccinate molecules unlikely to represent typical bound cholesterols are excluded. TRPs are tetrameric, with large clefts at the interfaces between subunits; cholesterol poses are located in hollows, largely within these clefts. Comparison of cholesterol poses with phospholipid binding sites suggests that binding of cholesterol to a TRP need not result in displacement of phospholipid molecules from the TRP surface.

Local Enrichment of Unsaturated Chains around the A2A Adenosine Receptor

Two 15 μs all-atom simulations of the A_2A adenosine receptor were obtained in a ternary mixture of cholesterol, saturated phosphatidylcholine lipids, and unsaturated phosphatidylcholine lipids. An analysis of local lipid solvation is reported on the basis of a Voronoi tessellation of the upper and lower leaflets, identifying first and second solvation shells. The local environments of both the inactive state and the partially active state of the receptor are significantly enriched with unsaturated chains but depleted of cholesterol and saturated chains, relative to the bulk membrane composition. In spite of the local depletion of cholesterol, the partially active receptor binds cholesterol at three locations during the entire simulation trajectory. These long-lived interactions represent the extreme of a very broad distribution of first-solvation shell lipid lifetimes, confounding sharp distinctions between lipid interactions. The broad distributions of lifetimes also make equilibrating the local lipid environment difficult, necessitating long simulation times.