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

Direct effect of membrane environment on the activation of mGluR2 revealed by single-molecule FRET

The microenvironment of membrane receptors controls their mobility, structure, interactions, and dynamics, but a systematic understanding of how it modulates receptor function is often lacking. Using single-molecule Förster resonance energy transfer (smFRET), we characterized how detergents and cholesterol modulate the conformational dynamics of metabotropic glutamate receptor 2 (mGluR2), a class C G protein-coupled receptor (GPCR). We found that, within the resolution of our measurements, all tested detergents stabilize the same overall active and inactive structure of different domains of mGluR2. However, the degree of stabilization and the equilibrium between active and inactive conformations depended on the detergent. Detergents with a single hydrophobic tail increased the active state occupancy compared to those with long, branched tails. Adding cholesterol to micelles with branched hydrophobic tails shifted the equilibrium toward the inactive state. Mutagenesis identified residues potentially involved in cholesterol interaction with mGluR2. Targeting the cholesterol-binding site with synthetic molecules could be a viable therapeutic approach.

Oxysterols Modulate Protein-Sterol Interactions to Impair CXCR4 Signaling in Aging Cells

Organismal aging is accompanied by the accumulation of senescent cells in the body, which drives tissue dysfunction. Senescent cells have a distinctive profile, including proliferation arrest, resistance to apoptosis, altered gene expression, and high inflammation. Despite global signaling and metabolic dysregulation during senescence, the underlying reasons for changes in signaling remain unclear. GPCRs are pivotal in cellular signaling, dynamically mediating the complex interplay between cells and their surrounding environment to maintain cellular homeostasis. The chemokine receptor CXCR4 plays a crucial role in modulating immune responses and inflammation. It has been shown that the expression of CXCR4 increases in cells undergoing senescence, which enhances inflammation postactivation. Here, we examine CXCR4 signaling in deeply senescent cells (aged cells), where cholesterol and its oxidized derivatives, oxysterols, affect receptor function. We report elevated oxysterol levels in senescent cells, which altered classical CXCL12-mediated CXCR4 signaling. Tail-oxidized sterols disrupted signaling more than ring-oxidized counterparts. Molecular dynamics simulations revealed that 27-hydroxycholesterol displaces cholesterol and binds strongly to alter the conformation of critical signaling residues, modifying the sterol-CXCR4 interaction landscape. Our study provides a molecular view of the observed mitigated GPCR signaling in the presence of oxysterols, which switched G-protein signaling from Gαi/o to Gαs class. Overall, we present an altered paradigm of GPCR signaling, where cholesterol oxidation alters the signaling outcome in aged cells.

Structure and dynamics determine G protein coupling specificity at a class A GPCR

G protein-coupled receptors (GPCRs) exhibit varying degrees of selectivity for different G protein isoforms. Despite the abundant structures of GPCR-G protein complexes, little is known about the mechanism of G protein coupling specificity. The β2-adrenergic receptor is an example of GPCR with high selectivity for Gαs, the stimulatory G protein for adenylyl cyclase, and much weaker for the Gαi family of G proteins inhibiting adenylyl cyclase. By developing a Gαi-biased agonist (LM189), we provide structural and biophysical evidence supporting that distinct conformations at ICL2 and TM6 are required for coupling of the different G protein subtypes Gαs and Gαi. These results deepen our understanding of G protein specificity and bias and can accelerate the design of ligands that select for preferred signaling pathways.

Cholesterol inhibits assembly and oncogenic 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, we have a poor understanding 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 reveal that a reduction of plasma membrane cholesterol strongly promotes EphA2 self-assembly. Indeed, low cholesterol levels cause a similar effect to the EphA2 ligand ephrinA1-Fc. These results indicate that cholesterol inhibits EphA2 assembly. Phosphorylation studies in different cell lines reveal that low cholesterol increased phospho-serine levels in EphA2, 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 it also suggests that cholesterol acts as a molecular safeguard mechanism that prevents uncontrolled self-assembly and activation of EphA2.

Structure and dynamics determine G protein coupling specificity at a class A GPCR

G protein coupled receptors (GPCRs) exhibit varying degrees of selectivity for different G protein isoforms. Despite the abundant structures of GPCR-G protein complexes, little is known about the mechanism of G protein coupling specificity. The β2-adrenergic receptor is an example of GPCR with high selectivity for Gαs, the stimulatory G protein for adenylyl cyclase, and much weaker for the Gαi family of G proteins inhibiting adenylyl cyclase. By developing a new Gαi-biased agonist (LM189), we provide structural and biophysical evidence supporting that distinct conformations at ICL2 and TM6 are required for coupling of the different G protein subtypes Gαs and Gαi. These results deepen our understanding of G protein specificity and bias and can accelerate the design of ligands that select for preferred signaling pathways.

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.

Anisotropic interactions for continuum modeling of protein-membrane systems

In this work, a model for anisotropic interactions between proteins and cellular membranes is proposed for large-scale continuum simulations. The framework of the model is based on dynamic density functional theory, which provides a formalism to describe the lipid densities within the membrane as continuum fields while still maintaining the fidelity of the underlying molecular interactions. Within this framework, we extend recent results to include the anisotropic effects of protein-lipid interactions. As applications, we consider two membrane proteins of biological interest: a RAS-RAF complex tethered to the membrane and a membrane embedded G protein-coupled receptor. A strong qualitative and quantitative agreement is found between the numerical results and the corresponding molecular dynamics simulations. Combining the scope of continuum level simulations with the details from molecular level particle simulations enables research into protein-membrane behaviors at a more biologically relevant scale, which crucially can also be accessed via experiment.

β2-Adrenergic Regulation of the Neuromuscular Transmission and Its Lipid-Dependent Switch

β2-Adrenoceptors (β2-ARs) are the most abundant subtype of adrenergic receptors in skeletal muscles. Their activation via a stabilization of postsynaptic architecture has beneficial effects in certain models of neuromuscular disorders. However, the ability of β2-ARs to regulate neuromuscular transmission at the presynaptic level is poorly understood. Using electrophysiological recordings and fluorescent FM dyes, we found that β2-AR activation with fenoterol enhanced an involvement of synaptic vesicles in exocytosis and neurotransmitter release during intense activity at the neuromuscular junctions of mouse diaphragm. This was accompanied by an improvement of contractile responses to phrenic nerve stimulation (but not direct stimulation of the muscle fibers) at moderate-to-high frequencies. β2-ARs mainly reside in lipid microdomains enriched with cholesterol and sphingomyelin. The latter is hydrolyzed by sphingomyelinases, whose upregulation occurs in many conditions characterized by muscle atrophy and sympathetic nerve hyperactivity. Sphingomyelinase treatment reversed the effects of β2-AR agonist on the neurotransmitter release and synaptic vesicle recruitment to the exocytosis during intense activity. Inhibition of Gi protein with pertussis toxin completely prevented the sphingomyelinase-mediated inversion in the β2-AR agonist action. Note that lipid raft disrupting enzyme cholesterol oxidase had the same effect on β2-AR agonist-mediated changes in neurotransmission as sphingomyelinase. Thus, β2-AR agonist fenoterol augmented recruitment and release of synaptic vesicles during intense activity in the diaphragm neuromuscular junctions. Sphingomyelin hydrolysis inversed the effects of β2-AR agonist on neurotransmission probably via switching to Gi protein-dependent signaling. This phenomenon may reflect a dependence of the β2-AR signaling on lipid raft integrity in the neuromuscular junctions.

Targeting CCR3 with antagonist SB 328437 sensitizes 5‑fluorouracil‑resistant gastric cancer cells: Experimental evidence and computational insights

Gastric cancer (GC) ranks fifth globally in cancer diagnoses and third for cancer-related deaths. Chemotherapy with 5-fluorouracil (5-FU), a primary treatment, faces challenges due to the development of chemoresistance. Tumor microenvironment factors, including C-C motif chemokine receptor 3 (CCR3), can contribute to chemoresistance. The present study evaluated the effect of CCR3 receptor inhibition using the antagonist SB 328437 and the molecular dynamics of this interaction on resistance to 5-FU in gastric cancer cells. The 5-FU-resistant AGS cell line (AGS R-5FU) demonstrated notable tolerance to higher concentrations of 5-FU, with a 2.6-fold increase compared with the parental AGS cell line. Furthermore, the mRNA expression levels of thymidylate synthase (TS), a molecular marker for 5-FU resistance, were significantly elevated in AGS R-5FU cells. CCR3 was shown to be expressed at significantly higher levels in these resistant cells. Combining SB 328437 with 5-FU resulted in a significant decrease in cell viability, particularly at higher concentrations of 5-FU. Furthermore, when SB 328437 was combined with 5-FU at a high concentration, the relative mRNA expression levels of CCR3 and TS decreased significantly. Computational analysis of CCR3 demonstrated dynamic conformational changes, especially in extracellular loop 2 region, which indicated potential alterations in ligand recognition. Docking simulations demonstrated that SB 328437 bound to the allosteric site of CCR3, inducing a conformational change in ECL2 and hindering ligand recognition. The present study provides comprehensive information on the molecular and structural aspects of 5-FU resistance and CCR3 modulation, highlighting the potential for therapeutic application of SB 328437 in GC treatment.

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

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

Cholesterol Dependence of the Conformational Changes in Metabotropic Glutamate Receptor 1

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

The density-threshold affinity: Calculating lipid binding affinities from unbiased coarse-grained molecular dynamics simulations

Many membrane proteins are sensitive to their local lipid environment. As structural methods for membrane proteins have improved, there is growing evidence of direct, specific binding of lipids to protein surfaces. Unfortunately the workhorse of understanding protein-small molecule interactions, the binding affinity for a given site, is experimentally inaccessible for these systems. Coarse-grained molecular dynamics simulations can be used to bridge this gap, and are relatively straightforward to learn. Such simulations allow users to observe spontaneous binding of lipids to membrane proteins and quantify localized densities of individual lipids or lipid fragments. In this chapter we outline a protocol for extracting binding affinities from these localized distributions, known as the "density threshold affinity." The density threshold affinity uses an adaptive and flexible definition of site occupancy that alleviates the need to distinguish between "bound'' lipids and bulk lipids that are simply diffusing through the site. Furthermore, the method allows "bead-level" resolution that is suitable for the case where lipids share binding sites, and circumvents ambiguities about a relevant reference state. This approach provides a convenient and straightforward method for comparing affinities of a single lipid species for multiple sites, multiple lipids for a single site, and/or a single lipid species modeled using multiple forcefields.

Modelling lipid rafts formation through chemo-mechanical interplay triggered by receptor-ligand binding

Cell membranes, mediator of many biological mechanisms from adhesion and metabolism up to mutation and infection, are highly dynamic and heterogeneous environments exhibiting a strong coupling between biochemical events and structural re-organisation. This involves conformational changes induced, at lower scales, by lipid order transitions and by the micro-mechanical interplay of lipids with transmembrane proteins and molecular diffusion. Particular attention is focused on lipid rafts, ordered lipid microdomains rich of signalling proteins, that co-localise to enhance substance trafficking and activate different intracellular biochemical pathways. In this framework, the theoretical modelling of the dynamic clustering of lipid rafts implies a full multiphysics coupling between the kinetics of phase changes and the mechanical work performed by transmembrane proteins on lipids, involving the bilayer elasticity. This mechanism produces complex interspecific dynamics in which membrane stresses and chemical potentials do compete by determining different morphological arrangements, alteration in diffusive walkways and coalescence phenomena, with a consequent influence on both signalling potential and intracellular processes. Therefore, after identifying the leading chemo-mechanical interactions, the present work investigates from a modelling perspective the spatio-temporal evolution of raft domains to theoretically explain co-localisation and synergy between proteins' activation and raft formation, by coupling diffusive and mechanical phenomena to observe different morphological patterns and clustering of ordered lipids. This could help to gain new insights into the remodelling of cell membranes and could potentially suggest mechanically based strategies to control their selectivity, by orienting intracellular functions and mechanotransduction.

Exploring Free Energies of Specific Protein Conformations Using the Martini Force Field

Coarse-grained (CG) level molecular dynamics simulations are routinely used to study various biomolecular processes. The Martini force field is currently the most widely adopted parameter set for such simulations. The functional form of this and several other CG force fields enforces secondary protein structure support by employing a variety of harmonic potentials or restraints that favor the protein's native conformation. We propose a straightforward method to calculate the energetic consequences of transitions between predefined conformational states in systems in which multiple factors can affect protein conformational equilibria. This method is designed for use within the Martini force field and involves imposing conformational transitions by linking a Martini-inherent elastic network to the coupling parameter λ. We demonstrate the applicability of our method using the example of five biomolecular systems that undergo experimentally characterized conformational transitions between well-defined structures (Staphylococcal nuclease, C-terminal segment of surfactant protein B, LAH4 peptide, and β2-adrenergic receptor) as well as between folded and unfolded states (GCN4 leucine zipper protein). The results show that the relative free energy changes associated with protein conformational transitions, which are affected by various factors, such as pH, mutations, solvent, and lipid membrane composition, are correctly reproduced. The proposed method may be a valuable tool for understanding how different conditions and modifications affect conformational equilibria in proteins.

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.

A cholesterol switch controls phospholipid scrambling by G protein-coupled receptors

Class A G protein-coupled receptors (GPCRs), a superfamily of cell membrane signaling receptors, moonlight as constitutively active phospholipid scramblases. The plasma membrane of metazoan cells is replete with GPCRs yet has a strong resting trans-bilayer phospholipid asymmetry, with the signaling lipid phosphatidylserine confined to the cytoplasmic leaflet. To account for the persistence of this lipid asymmetry in the presence of GPCR scramblases, we hypothesized that GPCR-mediated lipid scrambling is regulated by cholesterol, a major constituent of the plasma membrane. We now present a technique whereby synthetic vesicles reconstituted with GPCRs can be supplemented with cholesterol to a level similar to that of the plasma membrane and show that the scramblase activity of two prototypical GPCRs, opsin and the β1-adrenergic receptor, is impaired upon cholesterol loading. Our data suggest that cholesterol acts as a switch, inhibiting scrambling above a receptor-specific threshold concentration to disable GPCR scramblases at the plasma membrane.

Lipids in Psychiatric Disorders: Functional and Potential Diagnostic Role as Blood Biomarkers

Lipids are a crucial component of the human brain, serving important structural and functional roles. They are involved in cell function, myelination of neuronal projections, neurotransmission, neural plasticity, energy metabolism, and neuroinflammation. Despite their significance, the role of lipids in the development of mental disorders has not been well understood. This review focused on the potential use of lipids as blood biomarkers for common mental illnesses, such as major depressive disorder, anxiety disorders, bipolar disorder, and schizophrenia. This review also discussed the impact of commonly used psychiatric medications, such as neuroleptics and antidepressants, on lipid metabolism. The obtained data suggested that lipid biomarkers could be useful for diagnosing psychiatric diseases, but further research is needed to better understand the associations between blood lipids and mental disorders and to identify specific biomarker combinations for each disease.

A cholesterol switch controls phospholipid scrambling by G protein-coupled receptors

Class A G protein-coupled receptors (GPCRs), a superfamily of cell membrane signaling receptors, moonlight as constitutively active phospholipid scramblases. The plasma membrane of metazoan cells is replete with GPCRs, yet has a strong resting trans-bilayer phospholipid asymmetry, with the signaling lipid phosphatidylserine confined to the cytoplasmic leaflet. To account for the persistence of this lipid asymmetry in the presence of GPCR scramblases, we hypothesized that GPCR-mediated lipid scrambling is regulated by cholesterol, a major constituent of the plasma membrane. We now present a technique whereby synthetic vesicles reconstituted with GPCRs can be supplemented with cholesterol to a level similar to that of the plasma membrane and show that the scramblase activity of two prototypical GPCRs, opsin and the β1-adrenergic receptor, is impaired upon cholesterol loading. Our data suggest that cholesterol acts as a switch, inhibiting scrambling above a receptor-specific threshold concentration to disable GPCR scramblases at the plasma membrane.

Martini 3 Coarse-Grained Force Field for Cholesterol

Cholesterol plays a crucial role in biomembranes by regulating various properties, such as fluidity, rigidity, permeability, and organization of lipid bilayers. The latest version of the Martini model, Martini 3, offers significant improvements in interaction balance, molecular packing, and inclusion of new bead types and sizes. However, the release of the new model resulted in the need to reparameterize many core molecules, including cholesterol. Here, we describe the development and validation of a Martini 3 cholesterol model, addressing issues related to its bonded setup, shape, volume, and hydrophobicity. The proposed model mitigates some limitations of its Martini 2 predecessor while maintaining or improving the overall behavior.

Membrane Mimetic-Dependence of GPCR Energy Landscapes

Protein function strongly depends on temperature, which is related to temperature-dependent changes in the equilibria of protein conformational states. We leveraged variable-temperature 19F-NMR spectroscopy to interrogate the temperature dependence of the conformational landscape of the human A2A adenosine receptor (A2AAR), a class A GPCR. Temperature-induced changes in the conformational equilibria of A2AAR in lipid nanodiscs were markedly dependent on the efficacy of bound drugs. While antagonist complexes displayed only modest changes as the temperature rose, both full and partial agonist complexes exhibited substantial increases in the active state population. Importantly, the temperature-dependent response of complexes with both full and partial agonists exhibited a pronounced sensitivity to the specific membrane mimetic employed. In striking contrast to observations within lipid nanodiscs, in detergent micelles the active state population exhibited different behavior for A2AAR complexes with both full and partial agonists. This underscores the importance of the protein environment in understanding the thermodynamics of GPCR activation.

How Neuromembrane Lipids Modulate Membrane Proteins: Insights from G-Protein-Coupled Receptors (GPCRs) and Receptor Tyrosine Kinases (RTKs)

Lipids play a diverse and critical role in cellular processes in all tissues. The unique lipid composition of nerve membranes is particularly interesting because it contains, among other things, polyunsaturated lipids, such as docosahexaenoic acid, which the body only gets through the diet. The crucial role of lipids in neurological processes, especially in receptor-mediated cell signaling, is emphasized by the fact that in many neuropathological diseases there are significant deviations in the lipid composition of nerve membranes compared to healthy individuals. The lipid composition of neuromembranes can significantly affect the function of receptors by regulating the physical properties of the membrane or by affecting specific interactions between receptors and lipids. In addition, it is worth noting that the ligand-binding pocket of many receptors is located inside the cell membrane, due to which lipids can even modulate the binding of ligands to their receptors. These mechanisms highlight the importance of lipids in the regulation of membrane receptor activation and function. In this article, we focus on two major protein families: G-protein-coupled receptors (GPCRs) and receptor tyrosine kinases (RTKs) and discuss how lipids affect their function in neuronal membranes, elucidating the basic mechanisms underlying neuronal function and dysfunction.

Quantitative Comparison against Experiments Reveals Imperfections in Force Fields' Descriptions of POPC-Cholesterol Interactions

Cholesterol is a central building block in biomembranes, where it induces orientational order, slows diffusion, renders the membrane stiffer, and drives domain formation. Molecular dynamics (MD) simulations have played a crucial role in resolving these effects at the molecular level; yet, it has recently become evident that different MD force fields predict quantitatively different behavior. Although easily neglected, identifying such limitations is increasingly important as the field rapidly progresses toward simulations of complex membranes mimicking the in vivo conditions: pertinent multicomponent simulations must capture accurately the interactions between their fundamental building blocks, such as phospholipids and cholesterol. Here, we define quantitative quality measures for simulations of binary lipid mixtures in membranes against the C-H bond order parameters and lateral diffusion coefficients from NMR spectroscopy as well as the form factors from X-ray scattering. Based on these measures, we perform a systematic evaluation of the ability of commonly used force fields to describe the structure and dynamics of binary mixtures of palmitoyloleoylphosphatidylcholine (POPC) and cholesterol. None of the tested force fields clearly outperforms the others across the tested properties and conditions. Still, the Slipids parameters provide the best overall performance in our tests, especially when dynamic properties are included in the evaluation. The quality evaluation metrics introduced in this work will, particularly, foster future force field development and refinement for multicomponent membranes using automated approaches.

Pro-phagocytic function and structural basis of GPR84 signaling

GPR84 is a unique orphan G protein-coupled receptor (GPCR) that can be activated by endogenous medium-chain fatty acids (MCFAs). The signaling of GPR84 is largely pro-inflammatory, which can augment inflammatory response, and GPR84 also functions as a pro-phagocytic receptor to enhance phagocytic activities of macrophages. In this study, we show that the activation of GPR84 by the synthetic agonist 6-OAU can synergize with the blockade of CD47 on cancer cells to induce phagocytosis of cancer cells by macrophages. We also determine a high-resolution structure of the GPR84-Gi signaling complex with 6-OAU. This structure reveals an occluded binding pocket for 6-OAU, the molecular basis of receptor activation involving non-conserved structural motifs of GPR84, and an unusual Gi-coupling interface. Together with computational docking and simulations studies, this structure also suggests a mechanism for the high selectivity of GPR84 for MCFAs and a potential routes of ligand binding and dissociation. These results provide a framework for understanding GPR84 signaling and developing new drugs targeting GPR84.

Lipid packing is disrupted in copolymeric nanodiscs compared with intact membranes

Discoidal lipid-protein nanoparticles known as nanodiscs are widely used tools in structural and membrane biology. Amphipathic, synthetic copolymers have recently become an attractive alternative to membrane scaffold proteins for the formation of nanodiscs. Such copolymers can directly intercalate into, and form nanodiscs from, intact membranes without detergents. Although these copolymer nanodiscs can extract native membrane lipids, it remains unclear whether native membrane properties are also retained. To determine the extent to which bilayer lipid packing is retained in nanodiscs, we measured the behavior of packing-sensitive fluorescent dyes in various nanodisc preparations compared with intact lipid bilayers. We analyzed styrene-maleic acid (SMA), diisobutylene-maleic acid (DIBMA), and polymethacrylate (PMA) as nanodisc scaffolds at various copolymer-to-lipid ratios and temperatures. Measurements of Laurdan spectral shifts revealed that dimyristoyl-phosphatidylcholine (DMPC) nanodiscs had increased lipid headgroup packing compared with large unilamellar vesicles (LUVs) above the lipid melting temperature for all three copolymers. Similar effects were observed for DMPC nanodiscs stabilized by membrane scaffolding protein MSP1E1. Increased lipid headgroup packing was also observed when comparing nanodiscs with intact membranes composed of binary mixtures of 1-palmitoyl-2-oleoyl-phosphocholine (POPC) and di-palmitoyl-phosphocholine (DPPC), which show fluid-gel-phase coexistence. Similarly, Laurdan reported increased headgroup packing in nanodiscs for biomimetic mixtures containing cholesterol, most notable for relatively disordered membranes. The magnitudes of these ordering effects were not identical for the various copolymers, with SMA being the most and DIBMA being the least perturbing. Finally, nanodiscs derived from mammalian cell membranes showed similarly increased lipid headgroup packing. We conclude that nanodiscs generally do not completely retain the physical properties of intact membranes.

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

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

Anionic phospholipids control mechanisms of GPCR-G protein recognition

G protein-coupled receptors (GPCRs) are embedded in phospholipids that strongly influence drug-stimulated signaling. Anionic lipids are particularly important for GPCR signaling complex formation, but a mechanism for this role is not understood. Using NMR spectroscopy, we explore the impact of anionic lipids on the function-related conformational equilibria of the human A2A adenosine receptor (A2AAR) in bilayers containing defined mixtures of zwitterionic and anionic phospholipids. Anionic lipids prime the receptor to form complexes with G proteins through a conformational selection process. Without anionic lipids, signaling complex formation proceeds through a less favorable induced fit mechanism. In computational models, anionic lipids mimic interactions between a G protein and positively charged residues in A2AAR at the receptor intracellular surface, stabilizing a pre-activated receptor conformation. Replacing these residues strikingly alters the receptor response to anionic lipids in experiments. High sequence conservation of the same residues among all GPCRs supports a general role for lipid-receptor charge complementarity in signaling.

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.

Regulation of membrane protein structure and function by their lipid nano-environment

Transmembrane proteins comprise ~30% of the mammalian proteome, mediating metabolism, signalling, transport and many other functions required for cellular life. The microenvironment of integral membrane proteins (IMPs) is intrinsically different from that of cytoplasmic proteins, with IMPs solvated by a compositionally and biophysically complex lipid matrix. These solvating lipids affect protein structure and function in a variety of ways, from stereospecific, high-affinity protein-lipid interactions to modulation by bulk membrane properties. Specific examples of functional modulation of IMPs by their solvating membranes have been reported for various transporters, channels and signal receptors; however, generalizable mechanistic principles governing IMP regulation by lipid environments are neither widely appreciated nor completely understood. Here, we review recent insights into the inter-relationships between complex lipidomes of mammalian membranes, the membrane physicochemical properties resulting from such lipid collectives, and the regulation of IMPs by either or both. The recent proliferation of high-resolution methods to study such lipid-protein interactions has led to generalizable insights, which we synthesize into a general framework termed the 'functional paralipidome' to understand the mutual regulation between membrane proteins and their surrounding lipid microenvironments.

Anionic Phospholipids Control Mechanisms of GPCR-G Protein Recognition

G protein-coupled receptors (GPCRs) are embedded in phospholipids that strongly influence drug-stimulated signaling. Anionic lipids are particularly important for GPCR signaling complex formation, but a mechanism for this role is not understood. Using NMR spectroscopy, we visualized the impact of anionic lipids on the function-related conformational equilibria of the human A 2A adenosine receptor (A 2A AR) in bilayers containing defined mixtures of zwitterionic and anionic phospholipids. Anionic lipids primed the receptor to form complexes with G proteins through a conformational selection process. Without anionic lipids, signaling complex formation proceeded through a less favorable induced fit mechanism. In computational models, anionic lipids mimicked interactions between a G protein and positively charged residues in A 2A AR at the receptor intracellular surface, stabilizing a pre-activated receptor conformation. Replacing these residues strikingly altered the receptor response to anionic lipids in experiments. High sequence conservation of the same residues among all GPCRs supports a general role for lipid-receptor charge complementarity in signaling.

PI(4,5)P2 and Cholesterol: Synthesis, Regulation, and Functions

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂) is the most abundant membrane phosphoinositide and cholesterol is an essential component of the plasma membrane (PM). Both lipids play key roles in a variety of cellular functions including as signaling molecules and major regulators of protein function. This chapter provides an overview of these two important lipids. Starting from a brief description of their structure, synthesis, and regulation, the chapter continues to describe the primary functions and signaling processes in which PI(4,5)P₂ and cholesterol are involved. While PI(4,5)P₂ and cholesterol can act independently, they often act in concert or affect each other's impact. The chapters in this volume on "Cholesterol and PI(4,5)P₂ in Vital Biological Functions: From Coexistence to Crosstalk" focus on the emerging relationship between cholesterol and PI(4,5)P₂ in a variety of biological systems and processes. In this chapter, the next section provides examples from the ion channel field demonstrating that PI(4,5)P₂ and cholesterol can act via common mechanisms. The chapter ends with a discussion of future directions.

Synergistic and Competitive Lipid Interactions in the Serotonin1A Receptor Microenvironment

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

Effects of Gly Residue and Cholesterol on the GXXXG-Mediated Parallel Association of Transmembrane Helices: A Single-Pair FRET Study

Small residue-mediated interhelical packing is ubiquitous in helical membrane proteins: however, the lipid dependence of its stability remains unclear. We previously demonstrated that the introduction of a GXXXG sequence in the middle of de novo-designed (AALALAA)₃ helices (AALALAA AGLALGA AALALAA) facilitated their dimerization, which was abolished by cholesterol. Here single-pair FRET measurements revealed that a longer GXXXGXXXG segment (AALALAA A GLALGA AAGALAA) promoted helix dimerization in POPC/cholesterol bilayers, but not without cholesterol. The predicted dimer structures and degrees of helix packing suggested that helix dimers with small (∼10°) and large (∼55°) crossing angles were only stabilized in POPC and POPC/cholesterol membranes, respectively. A steric hindrance in the dimer interface and the large flexibility of helices prevented the formation of stable dimers. Therefore, amino acid sequences and lipid compositions distinctively constrain stable dimer structures in membranes.

Lysine 101 in the CRAC Motif in Transmembrane Helix 2 Confers Cholesterol-Induced Thermal Stability to the Serotonin1A Receptor

G protein-coupled receptors (GPCRs) constitute the largest class of membrane proteins that transduce signals across the plasma membrane and orchestrate a multitude of physiological processes within cells. The serotonin_1A receptor is a crucial neurotransmitter receptor in the GPCR family involved in a multitude of neurological, behavioral and cognitive functions. We have previously shown, using a combination of experimental and simulation approaches, that membrane cholesterol acts as a key regulator of organization, dynamics, signaling and endocytosis of the serotonin_1A receptor. In addition, we showed that membrane cholesterol stabilizes the serotonin_1A receptor against thermal deactivation. In the present work, we explored the molecular basis of cholesterol-induced thermal stability of the serotonin_1A receptor. For this, we explored the possible role of the K101 residue in a cholesterol recognition/interaction amino acid consensus (CRAC) motif in transmembrane helix 2 in conferring the thermal stability of the serotonin_1A receptor. Our results show that a mutation in the K101 residue leads to loss in thermal stability of the serotonin_1A receptor imparted by cholesterol, independent of membrane cholesterol content. We envision that our results could have potential implications in structural biological advancements of GPCRs and design of thermally stabilized receptors for drug development.

L-DOPA and Droxidopa: From Force Field Development to Molecular Docking into Human β2-Adrenergic Receptor

The increasing interest in the molecular mechanism of the binding of different agonists and antagonists to β₂-adrenergic receptor (β₂AR) inactive and active states has led us to investigate protein–ligand interactions using molecular docking calculations. To perform this study, the 3.2 Å X-ray crystal structure of the active conformation of human β₂AR in the complex with the endogenous agonist adrenaline has been used as a template for investigating the binding of two exogenous catecholamines to this adrenergic receptor. Here, we show the derivation of L-DOPA and Droxidopa OPLS all atom (AA) force field (FF) parameters via quantum mechanical (QM) calculations, molecular dynamics (MD) simulations in aqueous solutions of the two catecholamines and the molecular docking of both ligands into rigid and flexible β₂AR models. We observe that both ligands share with adrenaline similar experimentally observed binding anchor sites, which are constituted by Asp113/Asn312 and Ser203/Ser204/Ser207 side chains. Moreover, both L-DOPA and Droxidopa molecules exhibit binding affinities comparable to that predicted for adrenaline, which is in good agreement with previous experimental and computational results. L-DOPA and Droxidopa OPLS AA FFs have also been tested by performing MD simulations of these ligands docked into β₂AR proteins embedded in lipid membranes. Both hydrogen bonds and hydrophobic interaction networks observed over the 1 μs MD simulation are comparable with those derived from molecular docking calculations and MD simulations performed with the CHARMM FF.

Hepatocyte cholesterol content modulates glucagon receptor signalling

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Filling of a water-free void explains the allosteric regulation of the β1-adrenergic receptor by cholesterol

Recent high-pressure NMR results indicate that the preactive conformation of the β₁-adrenergic receptor (β₁AR) harbours completely empty cavities of ~100 ų volume, which disappear in the active conformation of the receptor. Here we have localized these cavities using X-ray crystallography of xenon-derivatized β₁AR crystals. One of the cavities is in direct contact with the cholesterol-binding pocket. Solution NMR shows that addition of the cholesterol analogue cholesteryl hemisuccinate impedes the formation of the active conformation of detergent-solubilized β₁AR by blocking conserved G protein-coupled receptor microswitches, concomitant with an affinity reduction of both isoprenaline and G protein-mimicking nanobody Nb80 for β₁AR detected by isothermal titration calorimetry. This wedge-like action explains the function of cholesterol as a negative allosteric modulator of β₁AR. A detailed understanding of G protein-coupled receptor regulation by cholesterol by filling of a dry void and the easy scouting for such voids by xenon may provide new routes for the development of allosteric drugs.

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.

Structural Rearrangement of the Serotonin Transporter Intracellular Gate Induced by Thr276 Phosphorylation

The reuptake of the neurotransmitter serotonin from the synaptic cleft by the serotonin transporter, SERT, is essential for proper neurological signaling. Biochemical studies have shown that Thr276 of transmembrane helix 5 is a site of PKG-mediated SERT phosphorylation, which has been proposed to shift the SERT conformational equilibria to promote inward-facing states, thus enhancing 5-HT transport. Recent structural and simulation studies have provided insights into the conformation transitions during substrate transport but have not shed light on SERT regulation via post-translational modifications. Using molecular dynamics simulations and Markov state models, we investigate how Thr276 phosphorylation impacts the SERT mechanism and its role in enhancing transporter stability and function. Our simulations show that Thr276 phosphorylation alters the hydrogen-bonding network involving residues on transmembrane helix 5. This in turn decreases the free energy barriers for SERT to transition to the inward-facing state, thus facilitating 5-HT import. The results provide atomistic insights into in vivo SERT regulation and can be extended to other pharmacologically important transporters in the solute carrier family.

PyLipID: A Python Package for Analysis of Protein-Lipid Interactions from Molecular Dynamics Simulations

Lipids play important modulatory and structural roles for membrane proteins. Molecular dynamics simulations are frequently used to provide insights into the nature of these protein-lipid interactions. Systematic comparative analysis requires tools that provide algorithms for objective assessment of such interactions. We introduce PyLipID, a Python package for the identification and characterization of specific lipid interactions and binding sites on membrane proteins from molecular dynamics simulations. PyLipID uses a community analysis approach for binding site detection, calculating lipid residence times for both the individual protein residues and the detected binding sites. To assist structural analysis, PyLipID produces representative bound lipid poses from simulation data, using a density-based scoring function. To estimate residue contacts robustly, PyLipID uses a dual-cutoff scheme to differentiate between lipid conformational rearrangements while bound from full dissociation events. In addition to the characterization of protein-lipid interactions, PyLipID is applicable to analysis of the interactions of membrane proteins with other ligands. By combining automated analysis, efficient algorithms, and open-source distribution, PyLipID facilitates the systematic analysis of lipid interactions from large simulation data sets of multiple species of membrane proteins.

Examining the Effect of Charged Lipids on Mitochondrial Outer Membrane Dynamics Using Atomistic Simulations

The outer mitochondrial membrane (OMM) is involved in multiple cellular functions such as apoptosis, inflammation and signaling via its membrane-associated and -embedded proteins. Despite the central role of the OMM in these vital phenomena, the structure and dynamics of the membrane have regularly been investigated in silico using simple two-component models. Accordingly, the aim was to generate the realistic multi-component model of the OMM and inspect its properties using atomistic molecular dynamics (MD) simulations. All major lipid components, phosphatidylinositol (PI), phosphatidylcholine (PC), phosphatidylethanolamine (PE), and phosphatidylserine (PS), were included in the probed OMM models. Because increased levels of anionic PS lipids have potential effects on schizophrenia and, more specifically, on monoamine oxidase B enzyme activity, the effect of varying the PS concentration was explored. The MD simulations indicate that the complex membrane lipid composition (MLC) behavior is notably different from the two-component PC-PE model. The MLC changes caused relatively minor effects on the membrane structural properties such as membrane thickness or area per lipid; however, notable effects could be seen with the dynamical parameters at the water-membrane interface. Increase of PS levels appears to slow down lateral diffusion of all lipids and, in general, the presence of anionic lipids reduced hydration and slowed down the PE headgroup rotation. In addition, sodium ions could neutralize the membrane surface, when PI was the main anionic component; however, a similar effect was not seen for high PS levels. Based on these results, it is advisable for future studies on the OMM and its protein or ligand partners, especially when wanting to replicate the correct properties on the water-membrane interface, to use models that are sufficiently complex, containing anionic lipid types, PI in particular.

Lipids and Phosphorylation Conjointly Modulate Complex Formation of β2-Adrenergic Receptor and β-arrestin2

G protein-coupled receptors (GPCRs) are the largest class of human membrane proteins that bind extracellular ligands at their orthosteric binding pocket to transmit signals to the cell interior. Ligand binding evokes conformational changes in GPCRs that trigger the binding of intracellular interaction partners (G proteins, G protein kinases, and arrestins), which initiate diverse cellular responses. It has become increasingly evident that the preference of a GPCR for a certain intracellular interaction partner is modulated by a diverse range of factors, e.g., ligands or lipids embedding the transmembrane receptor. Here, by means of molecular dynamics simulations of the β₂-adrenergic receptor and β-arrestin2, we study how membrane lipids and receptor phosphorylation regulate GPCR-arrestin complex conformation and dynamics. We find that phosphorylation drives the receptor’s intracellular loop 3 (ICL3) away from a native negatively charged membrane surface to interact with arrestin. If the receptor is embedded in a neutral membrane, the phosphorylated ICL3 attaches to the membrane surface, which widely opens the receptor core. This opening, which is similar to the opening in the G protein-bound state, weakens the binding of arrestin. The loss of binding specificity is manifested by shallower arrestin insertion into the receptor core and higher dynamics of the receptor-arrestin complex. Our results show that receptor phosphorylation and the local membrane composition cooperatively fine-tune GPCR-mediated signal transduction. Moreover, the results suggest that deeper understanding of complex GPCR regulation mechanisms is necessary to discover novel pathways of pharmacological intervention.

Mutants of the white ABCG Transporter in Drosophila melanogaster Have Deficient Olfactory Learning and Cholesterol Homeostasis

Drosophila's white gene encodes an ATP-binding cassette G-subfamily (ABCG) half-transporter. White is closely related to mammalian ABCG family members that function in cholesterol efflux. Mutants of white have several behavioral phenotypes that are independent of visual defects. This study characterizes a novel defect of white mutants in the acquisition of olfactory memory using the aversive olfactory conditioning paradigm. The w¹¹¹⁸ mutants learned slower than wildtype controls, yet with additional training, they reached wildtype levels of performance. The w¹¹¹⁸ learning phenotype is also found in the w^apricot and w^coral alleles, is dominant, and is rescued by genomic white and mini-white transgenes. Reducing dietary cholesterol strongly impaired olfactory learning for wildtype controls, while w¹¹¹⁸ mutants were resistant to this deficit. The w¹¹¹⁸ mutants displayed higher levels of cholesterol and cholesterol esters than wildtype under this low-cholesterol diet. Increasing levels of serotonin, dopamine, or both in the white mutants significantly improved w¹¹¹⁸ learning. However, serotonin levels were not lower in the heads of the w¹¹¹⁸ mutants than in wildtype controls. There were also no significant differences found in synapse numbers within the w¹¹¹⁸ brain. We propose that the w¹¹¹⁸ learning defect may be due to inefficient biogenic amine signaling brought about by altered cholesterol homeostasis.

Glucosylceramide synthase deficiency in the heart compromises β1-adrenergic receptor trafficking

AIMS

Cardiac injury and remodelling are associated with the rearrangement of cardiac lipids. Glycosphingolipids are membrane lipids that are important for cellular structure and function, and cardiac dysfunction is a characteristic of rare monogenic diseases with defects in glycosphingolipid synthesis and turnover. However, it is not known how cardiac glycosphingolipids regulate cellular processes in the heart. The aim of this study is to determine the role of cardiac glycosphingolipids in heart function.

METHODS AND RESULTS

Using human myocardial biopsies, we showed that the glycosphingolipids glucosylceramide and lactosylceramide are present at very low levels in non-ischaemic human heart with normal function and are elevated during remodelling. Similar results were observed in mouse models of cardiac remodelling. We also generated mice with cardiomyocyte-specific deficiency in Ugcg, the gene encoding glucosylceramide synthase (hUgcg-/- mice). In 9- to 10-week-old hUgcg-/- mice, contractile capacity in response to dobutamine stress was reduced. Older hUgcg-/- mice developed severe heart failure and left ventricular dilatation even under baseline conditions and died prematurely. Using RNA-seq and cell culture models, we showed defective endolysosomal retrograde trafficking and autophagy in Ugcg-deficient cardiomyocytes. We also showed that responsiveness to β-adrenergic stimulation was reduced in cardiomyocytes from hUgcg-/- mice and that Ugcg knockdown suppressed the internalization and trafficking of β1-adrenergic receptors.

CONCLUSIONS

Our findings suggest that cardiac glycosphingolipids are required to maintain β-adrenergic signalling and contractile capacity in cardiomyocytes and to preserve normal heart function.

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.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

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.

Hydrodynamics of immiscible binary fluids with viscosity contrast: a multiparticle collision dynamics approach

We present a multiparticle collision dynamics (MPC) implementation of layered immiscible fluids A and B of different shear viscosities separated by planar interfaces. The simulated flow profile for imposed steady shear motion and the time-dependent shear stress functions are in excellent agreement with our continuum hydrodynamics results for the composite fluid. The wave-vector dependent transverse velocity auto-correlation functions (TVAF) in the bulk-fluid regions of the layers decay exponentially, and agree with those of single-phase isotropic MPC fluids. In addition, we determine the hydrodynamic mobilities of an embedded colloidal sphere moving steadily parallel or transverse to a fluid-fluid interface, as functions of the distance from the interface. The obtained mobilities are in good agreement with hydrodynamic force multipoles calculations, for a no-slip sphere moving under creeping flow conditions near a clean, ideally flat interface. The proposed MPC fluid-layer model can be straightforwardly implemented, and it is computationally very efficient. Yet, owing to the spatial discretization inherent to the MPC method, the model can not reproduce all hydrodynamic features of an ideally flat interface between immiscible fluids.

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

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