Two classes of cholesterol binding sites for the β2AR revealed by thermostability and NMR

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.

Cholesterol inhibits assembly and activation of the EphA2 receptor

The receptor tyrosine kinase EphA2 drives cancer malignancy by facilitating metastasis. EphA2 can be found in different self-assembly states: as a monomer, dimer, and oligomer. However, our understanding remains limited regarding which EphA2 state is responsible for driving pro-metastatic signaling. To address this limitation, we have developed SiMPull-POP, a single-molecule method for accurate quantification of membrane protein self-assembly. Our experiments revealed that a reduction of plasma membrane cholesterol strongly promoted EphA2 self-assembly. Indeed, low cholesterol caused a similar effect to the EphA2 ligand ephrinA1-Fc. These results indicate that cholesterol inhibits EphA2 assembly. Phosphorylation studies in different cell lines revealed that low cholesterol increased phospho-serine levels, the signature of oncogenic signaling. Investigation of the mechanism that cholesterol uses to inhibit the assembly and activity of EphA2 indicate an in-trans effect, where EphA2 is phosphorylated by protein kinase A downstream of beta-adrenergic receptor activity, which cholesterol also inhibits. Our study not only provides new mechanistic insights on EphA2 oncogenic function, but also suggests that cholesterol acts as a molecular safeguard mechanism that prevents uncontrolled self-assembly and activation of EphA2.

The mechanism of 25-hydroxycholesterol-mediated suppression of atrial β1-adrenergic responses

25-Hydroxycholesterol (25HC) is a biologically active oxysterol, whose production greatly increases during inflammation by macrophages and dendritic cells. The inflammatory reactions are frequently accompanied by changes in heart regulation, such as blunting of the cardiac β-adrenergic receptor (AR) signaling. Here, the mechanism of 25HC-dependent modulation of responses to β-AR activation was studied in the atria of mice. 25HC at the submicromolar levels decreased the β-AR-mediated positive inotropic effect and enhancement of the Ca2+ transient amplitude, without changing NO production. Positive inotropic responses to β1-AR (but not β2-AR) activation were markedly attenuated by 25HC. The depressant action of 25HC on the β1-AR-mediated responses was prevented by selective β3-AR antagonists as well as inhibitors of Gi protein, Gβγ, G protein-coupled receptor kinase 2/3, or β-arrestin. Simultaneously, blockers of protein kinase D and C as well as a phosphodiesterase inhibitor did not preclude the negative action of 25HC on the inotropic response to β-AR activation. Thus, 25HC can suppress the β1-AR-dependent effects via engaging β3-AR, Gi protein, Gβγ, G protein-coupled receptor kinase, and β-arrestin. This 25HC-dependent mechanism can contribute to the inflammatory-related alterations in the atrial β-adrenergic signaling.

Sterol-activated amyloid beta fibril formation

The metabolic processes that link Alzheimer's disease (AD) to elevated cholesterol levels in the brain are not fully defined. Amyloid beta (Aβ) plaque accumulation is believed to begin decades prior to symptoms and to contribute significantly to the disease. Cholesterol and its metabolites accelerate plaque formation through as-yet-undefined mechanisms. Here, the mechanism of cholesterol (CH) and cholesterol 3-sulfate (CS) induced acceleration of Aβ42 fibril formation is examined in quantitative ligand binding, Aβ42 fibril polymerization, and molecular dynamics studies. Equilibrium and pre-steady-state binding studies reveal that monomeric Aβ42•ligand complexes form and dissociate rapidly relative to oligomerization, that the ligand/peptide stoichiometry is 1-to-1, and that the peptide is likely saturated in vivo. Analysis of Aβ42 polymerization progress curves demonstrates that ligands accelerate polymer synthesis by catalyzing the conversion of peptide monomers into dimers that nucleate the polymerization reaction. Nucleation is accelerated ∼49-fold by CH, and ∼13,000-fold by CS - a minor CH metabolite. Polymerization kinetic models predict that at presumed disease-relevant CS and CH concentrations, approximately half of the polymerization nuclei will contain CS, small oligomers of neurotoxic dimensions (∼12-mers) will contain substantial CS, and fibril-formation lag times will decrease 13-fold relative to unliganded Aβ42. Molecular dynamics models, which quantitatively predict all experimental findings, indicate that the acceleration mechanism is rooted in ligand-induced stabilization of the peptide in non-helical conformations that readily form polymerization nuclei.

Alchemical Free Energy Calculations on Membrane-Associated Proteins

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Illuminating GPCR signaling mechanisms by NMR spectroscopy with stable-isotope labeled receptors

G protein-coupled receptors (GPCRs) exhibit remarkable structural plasticity, which underlies their capacity to recognize a wide range of extracellular molecules and interact with intracellular partner proteins. Nuclear magnetic resonance (NMR) spectroscopy is uniquely well-suited to investigate GPCR structural plasticity, enabled by stable-isotope "probes" incorporated into receptors that inform on structure and dynamics. Progress with stable-isotope labeling methods in Eukaryotic expression systems has enabled production of native or nearly-native human receptors with varied and complementary distributions of NMR probes. These advances have opened up new avenues for investigating the roles of conformational dynamics in signaling processes, including by mapping allosteric communication networks, understanding the specificity of GPCR interactions with partner proteins and exploring the impact of membrane environments on GPCR function.

Dual mechanisms of cholesterol-GPCR interactions that depend on membrane phospholipid composition

Cholesterol is a critical component of mammalian cell membranes and an allosteric modulator of G protein-coupled receptors (GPCRs), but divergent views exist on the mechanisms by which cholesterol influences receptor functions. Leveraging the benefits of lipid nanodiscs, i.e., quantitative control of lipid composition, we observe distinct impacts of cholesterol in the presence and absence of anionic phospholipids on the function-related conformational dynamics of the human A2A adenosine receptor (A2AAR). Direct receptor-cholesterol interactions drive activation of agonist-bound A2AAR in membranes containing zwitterionic phospholipids. Intriguingly, the presence of anionic lipids attenuates cholesterol's impact through direct interactions with the receptor, highlighting a more complex role for cholesterol that depends on membrane phospholipid composition. Targeted amino acid replacements at two frequently predicted cholesterol interaction sites showed distinct impacts of cholesterol at different receptor locations, demonstrating the ability to delineate different roles of cholesterol in modulating receptor signaling and maintaining receptor structural integrity.

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.

Comparative Study of Receptor-, Receptor State-, and Membrane-Dependent Cholesterol Binding Sites in A2A and A1 Adenosine Receptors Using Coarse-Grained Molecular Dynamics Simulations

We used coarse-grained molecular dynamics (CG MD) simulations to study protein-cholesterol interactions for different activation states of the A_2A adenosine receptor (A_2AR) and the A₁ adenosine receptor (A₁R) and predict new cholesterol binding sites indicating amino acid residues with a high residence time in three biologically relevant membranes. Compared to 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)-cholesterol and POPC-phosphatidylinositol-bisphosphate (PIP₂)-cholesterol, the plasma mimetic membrane best described the cholesterol binding sites previously detected for the inactive state of A_2AR and revealed the binding sites with long-lasting amino acid residues. We observed that using the plasma mimetic membrane and plotting residues with cholesterol residence time ≥2 μs, our CG MD simulations captured most obviously the cholesterol-protein interactions. For the inactive A_2AR, we identified one more binding site in which cholesterol is bound to residues with a long residence time compared to the previously detected, for the active A₁R, three binding sites, and for the inactive A₁R, two binding sites. We calculated that for the active states, cholesterol binds to residues with a much longer residence time compared to the inactive state for both A_2AR and A₁R. The stability of the identified binding sites to A₁R or A_2AR with CG MD simulations was additionally investigated with potential of mean force calculations using umbrella sampling. We observed that the binding sites with residues to which cholesterol has a long residence time in A_2AR have shallow binding free energy minima compared to the related binding sites in A₁R, suggesting a stronger binding for cholesterol to A₁R. The differences in binding sites in which cholesterol is stabilized and interacts with residues with a long residence time between active and inactive states of A₁R and A_2AR can be important for differences in functional activity and orthosteric agonist or antagonist affinity and can be used for the design of allosteric modulators, which can bind through lipid pathways. We observed a stronger binding for cholesterol to A₁R (i.e., generally higher association rates) compared to A_2AR, which remains to be demonstrated. For the active states, cholesterol binds to residues with much longer residence times compared to the inactive state for both A_2AR and A₁R. Taken together, binding sites of active A₁R may be considered as promising allosteric targets.

A solution NMR view of Lipidic Cubic Phases: Structure, dynamics, and beyond

Nuclear magnetic resonance (NMR) spectroscopy is well-established nowadays for the elucidation of the 3D structures of proteins and protein complexes, the evaluation of biomolecular dynamics with atomistic resolution across a range of time scales, the screening of drug candidates with site specificity, and for the quantitation of molecular translational diffusion. Lyotropic lipidic cubic phases (LCPs) are lipid bilayer-based materials with a complex geometry, formed via the spontaneous self-assembly of certain lipids in an aqueous environment at specific temperature ranges. LCPs have been successfully applied to the in meso crystallization of membrane proteins for structural studies by X-ray crystallography, and have also shown promising potential for serving as matrices for drug and nutrient delivery/release in vivo. The characterization of the structural and dynamics properties of LCPs is of significant interest for the application of these materials. Here we present a systematic review detailing the characterization of LCPs by solution NMR. Using LCPs formed by monoolein (MO) as an example, various aspects of LCPs readily accessible by solution NMR are covered, including spectral perturbation in the presence of additives, quantification of hydration levels, ¹³C relaxation-based measurements for studying atom-specific dynamics along the MO hydrocarbon chain, PGSE NMR measurement of translational diffusion and its correlation with release profiles, and the encapsulation of soluble proteins in LCPs. A brief discussion of future perspectives for the characterization of LCPs by solution NMR is also presented.

Dimerization of β2-adrenergic receptor is responsible for the constitutive activity subjected to inverse agonism

Dimerization of beta 2-adrenergic receptor (β₂-AR) has been observed across various physiologies. However, the function of dimeric β₂-AR is still elusive. Here, we revealed that dimerization of β₂-AR is responsible for the constitutive activity of β₂-AR generating inverse agonism. Using a co-immunoimmobilization assay, we found that transient β₂-AR dimers exist in a resting state, and the dimer was disrupted by the inverse agonists. A Gα_s preferentially interacts with dimeric β₂-AR, but not monomeric β₂-AR, in a resting state, resulting in the production of a resting cAMP level. The formation of β₂-AR dimers requires cholesterol on the plasma membrane. The cholesterol did not interfere with the agonist-induced activation of monomeric β₂-AR, unlike the inverse agonists, implying that the cholesterol is a specific factor regulating the dimerization of β₂-AR. Our model not only shows the function of dimeric β₂-AR but also provides a molecular insight into the mechanism of the inverse agonism of β₂-AR.

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.

Thermostability-based binding assays reveal complex interplay of cation, substrate and lipid binding in the bacterial DASS transporter, VcINDY

The divalent anion sodium symporter (DASS) family of transporters (SLC13 family in humans) are key regulators of metabolic homeostasis, disruption of which results in protection from diabetes and obesity, and inhibition of liver cancer cell proliferation. Thus, DASS transporter inhibitors are attractive targets in the treatment of chronic, age-related metabolic diseases. The characterisation of several DASS transporters has revealed variation in the substrate selectivity and flexibility in the coupling ion used to power transport. Here, using the model DASS co-transporter, VcINDY from Vibrio cholerae, we have examined the interplay of the three major interactions that occur during transport: the coupling ion, the substrate, and the lipid environment. Using a series of high-throughput thermostability-based interaction assays, we have shown that substrate binding is Na⁺-dependent; a requirement that is orchestrated through a combination of electrostatic attraction and Na⁺-induced priming of the binding site architecture. We have identified novel DASS ligands and revealed that ligand binding is dominated by the requirement of two carboxylate groups in the ligand that are precisely distanced to satisfy carboxylate interaction regions of the substrate-binding site. We have also identified a complex relationship between substrate and lipid interactions, which suggests a dynamic, regulatory role for lipids in VcINDY's transport cycle.

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.

The ligand-bound state of a G protein-coupled receptor stabilizes the interaction of functional cholesterol molecules

Cholesterol is a major component of mammalian plasma membranes that not only affects the physical properties of the lipid bilayer but also is the function of many membrane proteins including G protein-coupled receptors. The oxytocin receptor (OXTR) is involved in parturition and lactation of mammals and in their emotional and social behaviors. Cholesterol acts on OXTR as an allosteric modulator inducing a high-affinity state for orthosteric ligands through a molecular mechanism that has yet to be determined. Using the ion channel-coupled receptor technology, we developed a functional assay of cholesterol modulation of G protein-coupled receptors that is independent of intracellular signaling pathways and operational in living cells. Using this assay, we discovered a stable binding of cholesterol molecules to the receptor when it adopts an orthosteric ligand-bound state. This stable interaction preserves the cholesterol-dependent activity of the receptor in cholesterol-depleted membranes. This mechanism was confirmed using time-resolved FRET experiments on WT OXTR expressed in CHO cells. Consequently, a positive cross-regulation sequentially occurs in OXTR between cholesterol and orthosteric ligands.

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.

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.

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.

Identification of CB1 Receptor Allosteric Sites Using Force-Biased MMC Simulated Annealing and Validation by Structure-Activity Relationship Studies

Positive allosteric modulation of the cannabinoid 1 receptor (CB1R) has demonstrated distinct therapeutic advantages that address several limitations associated with orthosteric agonism and has opened a promising therapeutic avenue for further drug development. To advance the development of CB1R positive allosteric modulators, it is important to understand the molecular architecture of CB1R allosteric site(s). The goal of this work was to use Force-Biased MMC Simulated Annealing to identify binding sites for GAT228 (R), a partial allosteric agonist, and GAT229 (S), a positive allosteric modulator (PAM) at the CB1R. Our studies suggest that GAT228 binds in an intracellular (IC) TMH1-2-4 exosite that would allow this compound to act as a CB1 allosteric agonist as well as a CB1 PAM. In contrast, GAT229 binds at the extracellular (EC) ends of TMH2/3, just beneath the EC1 loop. At this site, this compound can act as CB1 PAM only. Finally, these results were successfully validated through the synthesis and biochemical evaluation of a focused library of compounds.

Elucidating ligand-bound structures of membrane proteins using solid-state NMR spectroscopy

Magic-angle-spinning (MAS) solid-state NMR spectroscopy is a versatile technique to elucidate functionally important protein-ligand interactions in lipid membranes. Here, we review recent solid-state NMR studies of membrane protein interactions with cholesterol, lipids, transported substrates, and peptide ligands. These studies are conducted in synthetic or native lipid bilayers to provide an accurate environment for ligand binding. The solid-state NMR approaches include multinuclear detection to gain comprehensive structural information, distance measurements to locate ligand-binding sites, and dynamic nuclear polarization and 1H detection to enhance spectral sensitivity. These studies provide novel insights into the mechanisms of virus budding, virus entry into cells, transmembrane signaling, substrate transport, antibacterial action, and many other biological processes.

NMR investigation of protein-ligand interactions for G-protein coupled receptors

In this review, we report NMR studies of ligand-GPCR interactions, including both ligand-observed and protein-observed NMR experiments. Published studies exemplify how NMR can be used as a powerful tool to design novel GPCR ligands and investigate the ligand-induced conformational changes of GPCRs. The strength of NMR also lies in its capability to explore the diverse signaling pathways and probe the allosteric modulation of these highly dynamic receptors. By offering unique opportunities for the identification, structural and functional characterization of GPCR ligands, NMR will likely play a major role for the generation of novel molecules both as new tools for the understanding of the GPCR function and as therapeutic compounds for a large diversity of pathologies.

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.

Interfacial Binding Sites for Cholesterol on G Protein-Coupled Receptors

A docking procedure is described that allows the transmembrane surface of a G protein-coupled receptor (GPCR) to be swept rapidly for potential binding sites for cholesterol at the bilayer interfaces on the two sides of the membrane. The procedure matches 89% of the cholesterols resolved in published GPCR crystal structures, when cholesterols likely to be crystal packing artifacts are excluded. Docking poses are shown to form distinct clusters on the protein surface, the clusters corresponding to "greasy hollows" between protein ridges. Docking poses depend on the angle of tilt of the GPCR in the surrounding lipid bilayer. It is suggested that thermal motion could alter the optimal binding pose for a cholesterol molecule, with the range of binding poses within a cluster providing a guide to the range of thermal motions likely for a cholesterol within a binding site.

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.

Role of membrane cholesterol in differential sensitivity of muscarinic receptor subtypes to persistently bound xanomeline

Xanomeline (3-(Hexyloxy)-4-(1-methyl-1,2,5,6-tetrahydropyridin-3-yl)-1,2,5-thiadiazole) is a muscarinic agonist that is considered to be functionally selective for the M₁/M₄ receptor subtypes. Part of xanomeline binding is resistant to washing. Wash-resistant xanomeline activates muscarinic receptors persistently, except for the M₅ subtype. Mutation of leucine 6.46 to isoleucine at M₁ or M₄ receptors abolished persistent activation by wash-resistant xanomeline. Reciprocal mutation of isoleucine 6.46 to leucine at the M₅ receptor made it sensitive to activation by wash-resistant xanomeline. Lowering of membrane cholesterol made M₁ and M₄ mutants and M₅ wild type receptors sensitive to activation by wash-resistant xanomeline. Molecular docking revealed a cholesterol binding site in the groove between transmembrane helices 6 and 7. Molecular dynamics showed that interaction of cholesterol with this binding site attenuates receptor activation. We hypothesize that differences in cholesterol binding to this site between muscarinic receptor subtypes may constitute the basis for xanomeline apparent functional selectivity and may have notable therapeutic implications. Differences in receptor-membrane interactions, rather than in agonist-receptor interactions, represent a novel possibility to achieve pharmacological selectivity. Our findings may be applicable to other G protein coupled receptors.

A Streamlined, General Approach for Computing Ligand Binding Free Energies and Its Application to GPCR-Bound Cholesterol

The theory of receptor-ligand binding equilibria has long been well-established in biochemistry, and was primarily constructed to describe dilute aqueous solutions. Accordingly, few computational approaches have been developed for making quantitative predictions of binding probabilities in environments other than dilute isotropic solution. Existing techniques, ranging from simple automated docking procedures to sophisticated thermodynamics-based methods, have been developed with soluble proteins in mind. Biologically and pharmacologically relevant protein-ligand interactions often occur in complex environments, including lamellar phases like membranes and crowded, nondilute solutions. Here, we revisit the theoretical bases of ligand binding equilibria, avoiding overly specific assumptions that are nearly always made when describing receptor-ligand binding. Building on this formalism, we extend the asymptotically exact Alchemical Free Energy Perturbation technique to quantifying occupancies of sites on proteins in a complex bulk, including phase-separated, anisotropic, or nondilute solutions, using a thermodynamically consistent and easily generalized approach that resolves several ambiguities of current frameworks. To incorporate the complex bulk without overcomplicating the overall thermodynamic cycle, we simplify the common approach for ligand restraints by using a single distance-from-bound-configuration (DBC) ligand restraint during AFEP decoupling from protein. DBC restraints should be generalizable to binding modes of most small molecules, even those with strong orientational dependence. We apply this approach to compute the likelihood that membrane cholesterol binds to known crystallographic sites on three GPCRs (β₂-adrenergic, 5HT-2B, and μ-opioid) at a range of concentrations. Nonideality of cholesterol in a binary cholesterol:phosphatidylcholine (POPC) bilayer is characterized and consistently incorporated into the interpretation. We find that the three sites exhibit very different affinities for cholesterol: The site on the adrenergic receptor is predicted to be high affinity, with 50% occupancy for 1:10⁹ CHOL:POPC mixtures. The sites on the 5HT-2B and μ-opioid receptor are predicted to be lower affinity, with 50% occupancy for 1:10³ CHOL:POPC and 1:10² CHOL:POPC, respectively. These results could not have been predicted from the crystal structures alone.

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.

Challenges and approaches to understand cholesterol-binding impact on membrane protein function: an NMR view

Experimental evidence for a direct role of lipids in determining the structure, dynamics, and function of membrane proteins leads to the term 'functional lipids'. In particular, the sterol molecule cholesterol modulates the activity of many membrane proteins. The precise nature of cholesterol-binding sites and the consequences of modulation of local membrane micro-viscosity by cholesterol, however, is often unknown. Here, we review the current knowledge of the interaction of cholesterol with transmembrane proteins, with a special focus on structural aspects of the interaction derived from nuclear magnetic resonance approaches. We highlight examples of the importance of cholesterol modulation of membrane protein function, discuss the specificity of cholesterol binding, and review the proposed binding motifs from a molecular perspective. We conclude with a short perspective on what could be future trends in research efforts targeted towards a better understanding of cholesterol/membrane protein interactions.

Enabling STD-NMR fragment screening using stabilized native GPCR: A case study of adenosine receptor

Structural studies of integral membrane proteins have been limited by the intrinsic conformational flexibility and the need to stabilize the proteins in solution. Stabilization by mutagenesis was very successful for structural biology of G protein-coupled receptors (GPCRs). However, it requires heavy protein engineering and may introduce structural deviations. Here we describe the use of specific calixarenes-based detergents for native GPCR stabilization. Wild type, full length human adenosine A_2A receptor was used to exemplify the approach. We could stabilize native, glycosylated, non-aggregated and homogenous A_2AR that maintained its ligand binding capacity. The benefit of the preparation for fragment screening, using the Saturation-Transfer Difference nuclear magnetic resonance (STD-NMR) experiment is reported. The binding of the agonist adenosine and the antagonist caffeine were observed and competition experiments with CGS-21680 and ZM241385 were performed, demonstrating the feasibility of the STD-based fragment screening on the native A_2A receptor. Interestingly, adenosine was shown to bind a second binding site in the presence of the agonist CGS-21680 which corroborates published results obtained with molecular dynamics simulation. Fragment-like compounds identified using STD-NMR showed antagonistic effects on A_2AR in the cAMP cellular assay. Taken together, our study shows that stabilization of native GPCRs represents an attractive approach for STD-based fragment screening and drug design.

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.

Determinants of Ligand Subtype-Selectivity at α1A-Adrenoceptor Revealed Using Saturation Transfer Difference (STD) NMR

α_1A- and α_1B-adrenoceptors (α_1A-AR and α_1B-AR) are closely related G protein-coupled receptors (GPCRs) that modulate the cardiovascular and nervous systems in response to binding epinephrine and norepinephrine. The GPCR gene superfamily is made up of numerous subfamilies that, like α_1A-AR and α_1B-AR, are activated by the same endogenous agonists but may modulate different physiological processes. A major challenge in GPCR research and drug discovery is determining how compounds interact with receptors at the molecular level, especially to assist in the optimization of drug leads. Nuclear magnetic resonance spectroscopy (NMR) can provide great insight into ligand-binding epitopes, modes, and kinetics. Ideally, ligand-based NMR methods require purified, well-behaved protein samples. The instability of GPCRs upon purification in detergents, however, makes the application of NMR to study ligand binding challenging. Here, stabilized α_1A-AR and α_1B-AR variants were engineered using Cellular High-throughput Encapsulation, Solubilization, and Screening (CHESS), allowing the analysis of ligand binding with Saturation Transfer Difference NMR (STD NMR). STD NMR was used to map the binding epitopes of epinephrine and A-61603 to both receptors, revealing the molecular determinants for the selectivity of A-61603 for α_1A-AR over α_1B-AR. The use of stabilized GPCRs for ligand-observed NMR experiments will lead to a deeper understanding of binding processes and assist structure-based drug design.

GPCRs: What Can We Learn from Molecular Dynamics Simulations?

Advances in the structural biology of G-protein Coupled Receptors have resulted in a significant step forward in our understanding of how this important class of drug targets function at the molecular level. However, it has also become apparent that they are very dynamic molecules, and moreover, that the underlying dynamics is crucial in shaping the response to different ligands. Molecular dynamics simulations can provide unique insight into the dynamic properties of GPCRs in a way that is complementary to many experimental approaches. In this chapter, we describe progress in three distinct areas that are particularly difficult to study with other techniques: atomic level investigation of the conformational changes that occur when moving between the various states that GPCRs can exist in, the pathways that ligands adopt during binding/unbinding events and finally, the influence of lipids on the conformational dynamics of GPCRs.

Membrane cholesterol effect on the 5-HT2A receptor: Insights into the lipid-induced modulation of an antipsychotic drug target

The serotonin 5-hydroxytryptamine 2A (5-HT_2A ) receptor is a G-protein-coupled receptor (GPCR) relevant for the treatment of CNS disorders. In this regard, neuronal membrane composition in the brain plays a crucial role in the modulation of the receptor functioning. Since cholesterol is an essential component of neuronal membranes, we have studied its effect on the 5-HT_2A receptor dynamics through all-atom MD simulations. We find that the presence of cholesterol in the membrane increases receptor conformational variability in most receptor segments. Importantly, detailed structural analysis indicates that conformational variability goes along with the destabilization of hydrogen bonding networks not only within the receptor but also between receptor and lipids. In addition to increased conformational variability, we also find receptor segments with reduced variability. Our analysis suggests that this increased stabilization is the result of stabilizing effects of tightly bound cholesterol molecules to the receptor surface. Our finding contributes to a better understanding of membrane-induced alterations of receptor dynamics and points to cholesterol-induced stabilizing and destabilizing effects on the conformational variability of GPCRs.

Na+ and solute diffusion in aqueous channels of Myverol bicontinuous cubic phase: PGSE NMR and computer modelling

The apparent diffusion coefficients of ²³ Na⁺ ions and the solute 2-fluoroethylamine present in the aqueous domain of a Myverol/water bulk bicontinuous cubic phase (BCP) were measured using pulsed field-gradient spin echo (PGSE) NMR spectroscopy. The measured values were dependent on the diffusion time interval, which is a characteristic of restricted diffusion. The translational motion of ²³ Na⁺ and water in the aqueous channels of a cubic phase were simulated using a Monte-Carlo random walk algorithm, and the simulation results were compared with those from real PGSE NMR experiments. The simulations indicated that diffusion of ²³ Na⁺ ions and water would appear to be restricted even on the shortest timescales available to PGSE NMR experiments. The micro-viscosity of the aqueous domain of the BCPs was estimated from the longitudinal relaxation times of ²³ Na⁺ and 2-fluoroethylamine; this was three times higher than in free solution and suggests one of (but not the only) likely impediments to the release of hydrophilic drugs from stabilised aqueous dispersions of BCPs (cubosomes) when they are used therapeutically in vivo. Monte Carlo simulations of diffusive efflux from cubosomes suggest that the principal impediment to drug release is presented by a surfactant or lipid barrier at the cubosome surface, which separates the BCP aqueous channels from the bulk solution. The dynamics inferred from these studies informs quantitative predictions of drug delivery from cubosomes. Copyright © 2016 John Wiley & Sons, Ltd.

Hypercholesterolemia impairs oxytocin-induced uterine contractility in late pregnant mouse

High cholesterol is known to negatively affect uterine contractility inex vivoconditions. The aim of the present study was to reveal the effect ofin vivohypercholesterolemia on spontaneous and oxytocin-induced uterine contractility in late pregnant mouse uterus. Female Swiss albino mice were fed with high cholesterol (HC) diet (0.5% sodium cholate, 1.25% cholesterol and 15% fat) for 6 weeks and then throughout the gestation period after mating. On day 19 of gestation, serum cholesterol level was increased more than 3-fold while triglycerides level was reduced in HC diet-fed animals as compared to control animals fed with a standard diet. In tension experiments, neither the mean integral tension of spontaneous contractility nor the response to CaCl2in high K+-depolarized tissues was altered, but the oxytocin-induced concentration-dependent contractile response in uterine strips was attenuated in hypercholesterolemic mice as compared to control. Similarly, hypercholesterolemia dampened concentration-dependent uterine contractions elicited by a GNAQ protein activator,Pasteurella multocidatoxin. However, it had no effect on endogenous oxytocin level either in plasma or in uterine tissue. It also did not affect the prostaglandin release in oxytocin-stimulated tissues. Western blot data showed a significant increase in caveolin-1 and GRK6 proteins but decline in oxytocin receptor, GNAQ and RHOA protein expressions in hypercholesterolemic mouse uterus. The results of the present study suggest that hypercholesterolemia may attenuate the uterotonic action of oxytocin in late pregnancy by causing downregulation of oxytocin receptors and suppressing the signaling efficacy through GNAQ and RHOA proteins.

Mechanism of allosteric regulation of β2-adrenergic receptor by cholesterol

There is evidence that lipids can be allosteric regulators of membrane protein structure and activation. However, there are no data showing how exactly the regulation emerges from specific lipid-protein interactions. Here we show in atomistic detail how the human β₂-adrenergic receptor (β₂AR) - a prototypical G protein-coupled receptor - is modulated by cholesterol in an allosteric fashion. Extensive atomistic simulations show that cholesterol regulates β₂AR by limiting its conformational variability. The mechanism of action is based on the binding of cholesterol at specific high-affinity sites located near the transmembrane helices 5-7 of the receptor. The alternative mechanism, where the β₂AR conformation would be modulated by membrane-mediated interactions, plays only a minor role. Cholesterol analogues also bind to cholesterol binding sites and impede the structural flexibility of β₂AR, however cholesterol generates the strongest effect. The results highlight the capacity of lipids to regulate the conformation of membrane receptors through specific interactions.

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.

Interaction of G protein coupled receptors and cholesterol

G protein coupled receptors (GPCRs) form the largest receptor superfamily in eukaryotic cells. Owing to their seven transmembrane helices, large parts of these proteins are embedded in the cholesterol-rich plasma membrane bilayer. Thus, GPCRs are always in proximity to cholesterol. Some of them are functionally dependent on the specific presence of cholesterol. Over the last years, enormous progress on receptor structures has been achieved. While lipophilic ligands other than cholesterol have been shown to bind either inside the helix bundle or at the receptor-lipid interface, the binding site of cholesterol was either a single transmembrane helix or a groove between two or more transmembrane helices. A clear preference for one of the two membrane leaflets has not been observed. Not surprisingly, many hydrophobic residues (primarily leucine and isoleucine) were found to be involved in cholesterol binding. In most cases, the rough β-face of cholesterol contacted the transmembrane helix bundle rather than the surrounding lipid matrix. The polar hydroxy group of cholesterol was localized near the water-membrane interface with potential hydrogen bonding to residues in receptor loop regions. Although a canonical motif, designated as CCM site, was detected as a specific cholesterol binding site in case of the β2AR, this site was not found to be occupied by cholesterol in other GPCRs possessing the same motif. Cholesterol-receptor interactions can increase the compactness of the receptor structure and are able to enhance the conformational stability towards active or inactive receptor states. Overall, all current data suggest a high plasticity of cholesterol interaction sites in GPCRs.