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

Tracking Membrane Protein Dynamics in Real Time

Abstract

Membrane proteins govern critical cellular processes and are central to human health and associated disease. Understanding of membrane protein function is obscured by the vast ranges of structural dynamics—both in the spatial and time regime—displayed in the protein and surrounding membrane. The membrane lipids have emerged as allosteric modulators of membrane protein function, which further adds to the complexity. In this review, we discuss several examples of membrane dependency. A particular focus is on how molecular dynamics (MD) simulation have aided to map membrane protein dynamics and how enhanced sampling methods can enable observing the otherwise inaccessible biological time scale. Also, time-resolved X-ray scattering in solution is highlighted as a powerful tool to track membrane protein dynamics, in particular when combined with MD simulation to identify transient intermediate states. Finally, we discuss future directions of how to further develop this promising approach to determine structural dynamics of both the protein and the surrounding lipids.

Graphic Abstract

Ligand based conformational space studies of the μ-opioid receptor

BACKGROUND

G protein-coupled receptors (GPCRs) comprise a family of membrane proteins that can be activated by a variety of external factors. The μ-opioid receptor (MOR), a class A GPCR, is the main target of morphine. Recently, enhanced sampling molecular dynamics simulations of a constitutively active mutant of MOR in its apo form allowed us to capture the novel intermediate states of activation, as well as the active state. This prompted us to apply the same techniques to wild type MOR in complex with ligands, in order to explore their contributions to the receptor conformational changes in the activation process.

METHODS

MOR was modeled in complex with agonists (morphine, BU72), a partial agonist (naloxone benzoylhydrazone) and an antagonist (naloxone). Replica exchange with solute tempering (REST2) molecular dynamics simulations were carried out for all systems. Trajectory frames were clustered, and the activation state of each cluster was assessed by two different methods.

RESULTS

Cluster sizes and activation indices show that while agonists stabilized structures in a higher activation state, the antagonist behaved oppositely. Morphine tends to drive the receptor towards increasing R165-T279 distances, while naloxone tends to increase the NPxxYA motif conformational change.

CONCLUSIONS

Despite not observing a full transition between inactive and active states, an important conformational change of transmembrane helix 5 was observed and associated with a ligand-driven step of the process.

GENERAL SIGNIFICANCE

The activation process of GPCRs is widely studied but still not fully understood. Here we carried out a step forward in the direction of gaining more details of this process.

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.

Cholesterol Localization around the Metabotropic Glutamate Receptor 2

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

Mitochondria Regulate Inflammatory Paracrine Signalling in Neurodegeneration

Mitochondrial dysfunction occurs in most neurodegenerative diseases, contributing to both their onset and progression. A recent breakthrough unveiled that propagation of the inflammatory response and subsequent neuronal injury are also mediated extracellularly by damaged mitochondria, which are released from active microglial cells into the brain milieu. These extracellular fragmented mitochondria can therefore generate sufficient toxicity to trigger neuronal death and widespread brain damage through activation of naïve astrocytes. Besides suggesting potential new pharmacological strategies of therapeutic intervention in neurodegeneration, this original work indicates that mitochondria might act as bioactive ligands exerting paracrine functions. This is an interesting, novel and impactful concept that deserves consideration by the scientific community, as the attention should now be focused on the identification of the specific receptors through which mitochondria mediate such an important extracellular signalling mechanism in neurological conditions.

The Glycosphingolipid GM3 Modulates Conformational Dynamics of the Glucagon Receptor

The extracellular domain (ECD) of class B1 G-protein-coupled receptors (GPCRs) plays a central role in signal transduction and is uniquely positioned to sense both the extracellular and membrane environments. Although recent studies suggest a role for membrane lipids in the modulation of class A and class F GPCR signaling properties, little is known about the effect of lipids on class B1 receptors. In this study, we employed multiscale molecular dynamics simulations to access the dynamics of the glucagon receptor (GCGR) ECD in the presence of native-like membrane bilayers. Simulations showed that the ECD could move about a hinge region formed by residues Q122–E126 to adopt both closed and open conformations relative to the transmembrane domain. ECD movements were modulated by binding of the glycosphingolipid GM3. These large-scale fluctuations in ECD conformation may affect the ligand binding and receptor activation properties. We also identify a unique phosphatidylinositol (4,5)-bisphosphate (PIP₂) interaction profile near intracellular loop (ICL) 2/TM3 at the G-protein-coupling interface, suggesting a mechanism of engaging G-proteins that may have a distinct dependence on PIP₂ compared with class A GPCRs. Given the structural conservation of class B1 GPCRs, the modulatory effects of GM3 and PIP₂ on GCGR may be conserved across these receptors, offering new insights into potential therapeutic targeting.

Comparative docking analysis of cholesterol analogs to ion channels to discriminate between stereospecific binding vs. stereospecific response

Cholesterol is a major component of the membrane and a key regulator of many ion channels. Multiple studies showed that cholesterol regulates ion channels in a stereospecific manner, with cholesterol but not its chiral isomers having a functional effect. This stereospecificity has been universally attributed to the specificity of cholesterol binding, with the assumption that only native cholesterol binds to the channels whereas its isomers do not. In this study, we challenge this paradigm by docking analyses of cholesterol and its chiral isomers to five ion channels whose response to cholesterol was shown to be stereospecific, Kir2.2, KirBac1.1, TRPV1, GABA_A and BK. The analysis is performed using AutoDock Vina to predict the binding poses and energies of the sterols to the channels and identify amino acids interacting with the sterol molecules. We found that for every ion channel tested herein all three sterols showed similar binding poses and significant overlap in the set of the amino acids that comprise the predicted binding sites, along with similar energetic favorability to these overlapping sites. We also found, however, that specific orientations of the three sterols within the binding sites of the channels are distinct, so that a subset of the interacting amino acids is unique to each sterol. We propose therefore, that contrary to previous thought, stereospecific effects of cholesterol should be attributed not to the lack of binding of the stereoisomers but to specific, unique interactions between the cholesterol molecule and the residues within the binding sites of the channels.

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

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

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

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

Computer simulations of protein-membrane systems

The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.

Natural allosteric modulators and their biological targets: molecular signatures and mechanisms

Covering: 2008 to 2018Over the last decade more than two hundred single natural products were confirmed as natural allosteric modulators (alloNPs) of proteins. The compounds are presented and discussed with the support of a chemical space, constructed using a principal component analysis (PCA) of molecular descriptors from chemical compounds of distinct databases. This analysis showed that alloNPs are dispersed throughout the majority of the chemical space defined by natural products in general. Moreover, a cluster of alloNPs was shown to occupy a region almost devoid of allosteric modulators retrieved from a dataset composed mainly of synthetic compounds, further highlighting the importance to explore the entire natural chemical space for probing allosteric mechanisms. The protein targets which alloNPs bind to comprised 81 different proteins, which were classified into 5 major groups, with enzymes, in particular hydrolases, being the main representative group. The review also brings a critical interpretation on the mechanisms by which alloNPs display their molecular action on proteins. In the latter analysis, alloNPs were classified according to their final effect on the target protein, resulting in 3 major categories: (i) local alteration of the orthosteric site; (ii) global alteration in protein dynamics that change function; and (iii) oligomer stabilisation or protein complex destabilisation via protein-protein interaction in sites distant from the orthosteric site. G-protein coupled receptors (GPCRs), which use a combination of the three types of allosteric regulation found, were also probed by natural products. In summary, the natural allosteric modulators reviewed herein emphasise their importance for exploring alternative chemotherapeutic strategies, potentially pushing the boundaries of the druggable space of pharmacologically relevant drug targets.

Role of cholesterol-mediated effects in GPCR heterodimers

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

Cholesterol impacts chemokine CCR5 receptor ligand-binding activity

The chemokine CCR5 receptor is target of maraviroc, a negative allosteric modulator of CCR5 that blocks the HIV protein gp120 from associating with the receptor, thereby inhibiting virus cellular entry. As noted with other G-protein-coupled receptor family members, the role of the lipid environment in CCR5 signaling remains obscure and very modestly investigated. Controversial literature on the impact of cholesterol (Chol) depletion in HIV infection and CCR5 signaling, including the hypothesis that Chol depletion could inhibit HIV infection, lead us to focus on the understanding of Chol impact in the first stages of receptor activation. To address this aim, the approach chosen was to employ reconstituted model lipid systems of controlled lipid composition containing CCR5 from two distinct expression systems: Pichia pastoris and cell-free expression. The characterization of receptor/ligand interaction in terms of total binding or competition binding assays was independently performed by plasmon waveguide resonance and fluorescence anisotropy, respectively. Maraviroc, a potent receptor antagonist, was the ligand investigated. Additionally, coarse-grained molecular dynamics simulation was employed to investigate Chol impact in the receptor-conformational flexibility and dynamics. Results obtained with receptor produced by different expression systems and using different biophysical approaches clearly demonstrate a considerable impact of Chol in the binding affinity of maraviroc to the receptor and receptor-conformational dynamics. Chol considerably decreases maraviroc binding affinity to the CCR5 receptor. The mechanisms by which this effect occurs seem to involve the adoption of distinct receptor-conformational states with restrained structural dynamics and helical motions in the presence of Chol.

Lipid Interactions of a Ciliary Membrane TRP Channel: Simulation and Structural Studies of Polycystin-2

Polycystin-2 (PC2) is a transient receptor potential (TRP) channel present in ciliary membranes of the kidney. PC2 shares a transmembrane fold with other TRP channels, in addition to an extracellular domain found in TRPP and TRPML channels. Using molecular dynamics (MD) simulations and cryoelectron microscopy we identify and characterize PIP₂ and cholesterol interactions with PC2. PC2 is revealed to have a PIP binding site close to the equivalent vanilloid/lipid binding site in the TRPV1 channel. A 3.0-Å structure reveals a binding site for cholesterol on PC2. Cholesterol interactions with the channel at this site are characterized by MD simulations. The two classes of lipid binding sites are compared with sites observed in other TRPs and in Kv channels. These findings suggest PC2, in common with other ion channels, may be modulated by both PIPs and cholesterol, and position PC2 within an emerging model of the roles of lipids in the regulation and organization of ciliary membranes.

•

Lipid interactions of PC2 channels have been explored by MD simulation and cryo-EM

•

PIP₂ binds to a site corresponding to the vanilloid/lipid binding site of TRPV1

•

Cholesterol binds between the S3 and S4 helices and S6 of the adjacent subunit

•

PC2, in common with other channels, may be modulated by PIPs and cholesterol

Integrating hydrogen-deuterium exchange mass spectrometry with molecular dynamics simulations to probe lipid-modulated conformational changes in membrane proteins

Biological membranes define the boundaries of cells and are composed primarily of phospholipids and membrane proteins. It has become increasingly evident that direct interactions of membrane proteins with their surrounding lipids play key roles in regulating both protein conformations and function. However, the exact nature and structural consequences of these interactions remain difficult to track at the molecular level. Here, we present a protocol that specifically addresses this challenge. First, hydrogen-deuterium exchange mass spectrometry (HDX-MS) of membrane proteins incorporated into nanodiscs of controlled lipid composition is used to obtain information on the lipid species that are involved in modulating the conformational changes in the membrane protein. Then molecular dynamics (MD) simulations in lipid bilayers are used to pinpoint likely lipid-protein interactions, which can be tested experimentally using HDX-MS. By bringing together the MD predictions with the conformational readouts from HDX-MS, we have uncovered key lipid-protein interactions implicated in stabilizing important functional conformations. This protocol can be applied to virtually any integral membrane protein amenable to classic biophysical studies and for which a near-atomic-resolution structure or homology model is available. This protocol takes ~4 d to complete, excluding the time for data analysis and MD simulations, which depends on the size of the protein under investigation.

A Perspective: Active Role of Lipids in Neurotransmitter Dynamics

Synaptic neurotransmission is generally considered as a function of membrane-embedded receptors and ion channels in response to the neurotransmitter (NT) release and binding. This perspective aims to widen the protein-centric view by including another vital component—the synaptic membrane—in the discussion. A vast set of atomistic molecular dynamics simulations and biophysical experiments indicate that NTs are divided into membrane-binding and membrane-nonbinding categories. The binary choice takes place at the water-membrane interface and follows closely the positioning of the receptors’ binding sites in relation to the membrane. Accordingly, when a lipophilic NT is on route to a membrane-buried binding site, it adheres on the membrane and, then, travels along its plane towards the receptor. In contrast, lipophobic NTs, which are destined to bind into receptors with extracellular binding sites, prefer the water phase. This membrane-based sorting splits the neurotransmission into membrane-independent and membrane-dependent mechanisms and should make the NT binding into the receptors more efficient than random diffusion would allow. The potential implications and notable exceptions to the mechanisms are discussed here. Importantly, maintaining specific membrane lipid compositions (MLCs) at the synapses, especially regarding anionic lipids, affect the level of NT-membrane association. These effects provide a plausible link between the MLC imbalances and neurological diseases such as depression or Parkinson’s disease. Moreover, the membrane plays a vital role in other phases of the NT life cycle, including storage and release from the synaptic vesicles, transport from the synaptic cleft, as well as their synthesis and degradation.

Lipid-Dependent Regulation of Ion Channels and G Protein-Coupled Receptors: Insights from Structures and Simulations

Ion channels and G protein-coupled receptors (GPCRs) are regulated by lipids in their membrane environment. Structural studies combined with biophysical and molecular simulation investigations reveal interaction sites for specific lipids on membrane protein structures. For K channels, PIP₂ plays a key role in regulating Kv and Kir channels. Likewise, several recent cryo-EM structures of TRP channels have revealed bound lipids, including PIP₂ and cholesterol. Among the pentameric ligand-gated ion channel family, structural and biophysical studies suggest the M4 TM helix may act as a lipid sensor, e.g., forming part of the binding sites for neurosteroids on the GABA_A receptor. Structures of GPCRs have revealed multiple cholesterol sites, which may modulate both receptor dynamics and receptor oligomerization. PIP₂ also interacts with GPCRs and may modulate their interactions with G proteins. Overall, it is evident that multiple lipid binding sites exist on channels and receptors that modulate their function allosterically and are potential druggable sites.

Agonist-induced membrane nanodomain clustering drives GLP-1 receptor responses in pancreatic beta cells

The glucagon-like peptide-1 receptor (GLP-1R), a key pharmacological target in type 2 diabetes (T2D) and obesity, undergoes rapid endocytosis after stimulation by endogenous and therapeutic agonists. We have previously highlighted the relevance of this process in fine-tuning GLP-1R responses in pancreatic beta cells to control insulin secretion. In the present study, we demonstrate an important role for the translocation of active GLP-1Rs into liquid-ordered plasma membrane nanodomains, which act as hotspots for optimal coordination of intracellular signaling and clathrin-mediated endocytosis. This process is dynamically regulated by agonist binding through palmitoylation of the GLP-1R at its carboxyl-terminal tail. Biased GLP-1R agonists and small molecule allosteric modulation both influence GLP-1R palmitoylation, clustering, nanodomain signaling, and internalization. Downstream effects on insulin secretion from pancreatic beta cells indicate that these processes are relevant to GLP-1R physiological actions and might be therapeutically targetable.

Nanodomain segregation and clustering of the glucagon-like peptide-1 receptor, a key target for type 2 diabetes therapy, is regulated by agonist binding, leading to compartmentalization of downstream signaling and clathrin-dependent internalization and impacting pancreatic beta cell responses.

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.

Reduced level of docosahexaenoic acid shifts GPCR neuroreceptors to less ordered membrane regions

G protein-coupled receptors (GPCRs) control cellular signaling and responses. Many of these GPCRs are modulated by cholesterol and polyunsaturated fatty acids (PUFAs) which have been shown to co-exist with saturated lipids in ordered membrane domains. However, the lipid compositions of such domains extracted from the brain cortex tissue of individuals suffering from GPCR-associated neurological disorders show drastically lowered levels of PUFAs. Here, using free energy techniques and multiscale simulations of numerous membrane proteins, we show that the presence of the PUFA DHA helps helical multi-pass proteins such as GPCRs partition into ordered membrane domains. The mechanism is based on hybrid lipids, whose PUFA chains coat the rough protein surface, while the saturated chains face the raft environment, thus minimizing perturbations therein. Our findings suggest that the reduction of GPCR partitioning to their native ordered environments due to PUFA depletion might affect the function of these receptors in numerous neurodegenerative diseases, where the membrane PUFA levels in the brain are decreased. We hope that this work inspires experimental studies on the connection between membrane PUFA levels and GPCR signaling.

Our current picture of cellular membranes depicts them as laterally heterogeneous sheets of lipids crowded with membrane proteins. These proteins often require a specific lipid environment to efficiently perform their functions. Certain neuroreceptor proteins are regulated by membrane cholesterol that is considered to be enriched in ordered membrane domains. In the brain, these very same domains also contain a fair amount of polyunsaturated fatty acids (PUFAs) that have also been discovered to interact favorably with many receptor proteins. However, certain neurological diseases—associated with the inadequate functioning of the neuroreceptors—seem to result in the decrease of brain PUFA levels. We hypothesized that this decrease in PUFA levels somehow inhibits receptor partitioning to cholesterol-rich domains, which could further compromise their function. We verified our hypothesis by an extensive set of computer simulations. They demonstrated that the PUFA–receptor interaction indeed leads to favorable partitioning of the receptors in the cholesterol-rich ordered domains. Moreover, the underlying mechanism based on the shielding of the rough protein surface by the PUFAs seems to be exclusive for multi-helical protein structures, of which neuroreceptors are a prime example.

Folding and Misfolding of Human Membrane Proteins in Health and Disease: From Single Molecules to Cellular Proteostasis

Advances over the past 25 years have revealed much about how the structural properties of membranes and associated proteins are linked to the thermodynamics and kinetics of membrane protein (MP) folding. At the same time biochemical progress has outlined how cellular proteostasis networks mediate MP folding and manage misfolding in the cell. When combined with results from genomic sequencing, these studies have established paradigms for how MP folding and misfolding are linked to the molecular etiologies of a variety of diseases. This emerging framework has paved the way for the development of a new class of small molecule "pharmacological chaperones" that bind to and stabilize misfolded MP variants, some of which are now in clinical use. In this review, we comprehensively outline current perspectives on the folding and misfolding of integral MPs as well as the mechanisms of cellular MP quality control. Based on these perspectives, we highlight new opportunities for innovations that bridge our molecular understanding of the energetics of MP folding with the nuanced complexity of biological systems. Given the many linkages between MP misfolding and human disease, we also examine some of the exciting opportunities to leverage these advances to address emerging challenges in the development of therapeutics and precision medicine.

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.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

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

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

Understanding the Role of Lipids in Signaling Through Atomistic and Multiscale Simulations of Cell Membranes

Cell signaling controls essentially all cellular processes. While it is often assumed that proteins are the key architects coordinating cell signaling, recent studies have shown more and more clearly that lipids are also involved in signaling processes in a number of ways. Lipids do, for instance, act as messengers, modulate membrane receptor conformation and dynamics, and control membrane receptor partitioning. Further, through structural modifications such as oxidation, the functions of lipids as part of signaling processes can be modified. In this context, in this article we discuss the understanding recently revealed by atomistic and coarse-grained computer simulations of nanoscale processes and underlying physicochemical principles related to lipids’ functions in cellular signaling.

Structural insight into Class F receptors - What have we learnt regarding agonist-induced activation?

Class F receptors, including the ten Frizzleds (FZD_1-10 ) and SMO, mediate the effects of WNTs and hedgehog proteins and belong to the superfamily of G protein-coupled receptors (GPCRs). While the recent, high-resolution insight into mechanisms of GPCR activation provides a better understanding of receptor activation in Class A, B and C GPCRs, it remains unclear how Class F receptors bind their ligands, how ligand binding is translated to receptor activation and how signal initiation and specification are achieved. Here, we summarize recent efforts in elucidating Class F receptor structure and activation mechanisms and critically discuss the progress made in this area. A better understanding of the activation mechanisms of Class F receptors is required to engage in mechanism-based and structure-guided drug discovery to exploit the large therapeutic potential of targeting these receptors pharmacologically.

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.

State-dependent Lipid Interactions with the A2a Receptor Revealed by MD Simulations Using In Vivo-Mimetic Membranes

Membranes are known to have modulatory effects on G protein-coupled receptors (GPCRs) via specific lipid interactions. However, the mechanisms of such modulations in physiological conditions and how they influence GPCR functions remain unclear. Here we report coarse-grained molecular dynamics simulations on the Adenosine A2a receptor in different conformational states embedded in an in vivo-mimetic membrane model. Nine lipid interaction sites were revealed. The strength of lipid interactions with these sites showed a degree of dependence on the conformational states of the receptor, suggesting that these lipids may regulate the conformational dynamics of the receptor. In particular, we revealed a dual role of PIP₂ on A2aR activation that involves both stabilization of the characteristic outward tilt of TM6 and enhancement of A2aR-mini-Gs association. Our results demonstrated that the bound lipids allosterically regulate the functional properties of GPCRs. These protein-lipid interactions provide a springboard for design of allosteric modulators of GPCRs.

Rapid diffusion of cholesterol along polyunsaturated membranes via deep dives

Cholesterol regulates the function of membrane proteins either via direct or membrane-mediated effects. Therefore, its ready availability is crucial for many protein-governed cellular processes. Recent studies suggest that cholesterol can partition to the core of polyunsaturated membranes, where cholesterol binding sites of many membrane proteins are also located. This core region is characterized by a lower viscosity. Therefore, we hypothesized that cholesterol partitioning into the membrane interior increases the rate of its diffusion in polyunsaturated membrane environments. We studied the behavior of cholesterol in membranes with increasing level of lipid chain unsaturation using a combination of atomistic and coarse-grained molecular dynamics simulations. Our simulations suggest a strong correlation between entropy-driven enhanced cholesterol partitioning to the membrane core and its faster lateral diffusion, which indicates that the less viscous membrane core indeed provides an efficient means for cholesterol movement in polyunsaturated membrane environments.

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.

An engineered thermal-shift screen reveals specific lipid preferences of eukaryotic and prokaryotic membrane proteins

Membrane bilayers are made up of a myriad of different lipids that regulate the functional activity, stability, and oligomerization of many membrane proteins. Despite their importance, screening the structural and functional impact of lipid–protein interactions to identify specific lipid requirements remains a major challenge. Here, we use the FSEC-TS assay to show cardiolipin-dependent stabilization of the dimeric sodium/proton antiporter NhaA, demonstrating its ability to detect specific protein-lipid interactions. Based on the principle of FSEC-TS, we then engineer a simple thermal-shift assay (GFP-TS), which facilitates the high-throughput screening of lipid- and ligand- interactions with membrane proteins. By comparing the thermostability of medically relevant eukaryotic membrane proteins and a selection of bacterial counterparts, we reveal that eukaryotic proteins appear to have evolved to be more dependent to the presence of specific lipids.

Membrane bilayers are made up of a myriad of different lipids that affect membrane proteins, but identifying those specific lipid requirements remains a challenge. Here authors present an engineered thermal-shift screen which reveals specific lipid preferences of eukaryotic and prokaryotic membrane proteins.

Microscopic view of lipids and their diverse biological functions

Biological membranes and their diverse lipid constituents play key roles in a broad spectrum of cellular and physiological processes. Characterization of membrane-associated phenomena at a microscopic level is therefore essential to our fundamental understanding of such processes. Due to the semi-fluid and dynamic nature of lipid bilayers, and their complex compositions, detailed characterization of biological membranes at an atomic scale has been refractory to experimental approaches. Computational modeling and simulation offer a highly complementary toolset with sufficient spatial and temporal resolutions to fill this gap. Here, we review recent molecular dynamics studies focusing on the diversity of lipid composition of biological membranes, or aiming at the characterization of lipid-protein interaction, with the overall goal of dissecting how lipids impact biological roles of the cellular membranes.

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.

Structural insights into positive and negative allosteric regulation of a G protein-coupled receptor through protein-lipid interactions

Lipids are becoming known as essential allosteric modulators of G protein-coupled receptor (GPCRs). However, how they exert their effects on GPCR conformation at the atomic level is still unclear. In light of recent experimental data, we have performed several long-timescale molecular dynamics (MD) simulations, totalling 24 μs, to rigorously map allosteric modulation and conformational changes in the β₂ adrenergic receptor (β2AR) that occur as a result of interactions with three different phospholipids. In particular, we identify different sequential mechanisms behind receptor activation and deactivation, respectively, mediated by specific lipid interactions with key receptor regions. We show that net negatively charged lipids stabilize an active-like state of β2AR that is able to dock G_sα protein. Clustering of anionic lipids around the receptor with local distortion of membrane thickness is also apparent. On the other hand, net-neutral zwitterionic lipids inactivate the receptor, generating either fully inactive or intermediate states, with kinetics depending on lipid headgroup charge distribution and hydrophobicity. These chemical differences alter membrane thickness and density, which differentially destabilize the β2AR active state through lateral compression effects.

Closely related, yet unique: Distinct homo- and heterodimerization patterns of G protein coupled chemokine receptors and their fine-tuning by cholesterol

Chemokine receptors, a subclass of G protein coupled receptors (GPCRs), play essential roles in the human immune system, they are involved in cancer metastasis as well as in HIV-infection. A plethora of studies show that homo- and heterodimers or even higher order oligomers of the chemokine receptors CXCR4, CCR5, and CCR2 modulate receptor function. In addition, membrane cholesterol affects chemokine receptor activity. However, structural information about homo- and heterodimers formed by chemokine receptors and their interplay with cholesterol is limited. Here, we report homo- and heterodimer configurations of the chemokine receptors CXCR4, CCR5, and CCR2 at atomistic detail, as obtained from thousands of molecular dynamics simulations. The observed homodimerization patterns were similar for the closely related CC chemokine receptors, yet they differed significantly between the CC receptors and CXCR4. Despite their high sequence identity, cholesterol modulated the CC homodimer interfaces in a subtype-specific manner. Chemokine receptor heterodimers display distinct dimerization patterns for CXCR4/CCR5 and CXCR4/CCR2. Furthermore, associations between CXCR4 and CCR5 reveal an increased cholesterol-sensitivity as compared to CXCR4/CCR2 heterodimerization patterns. This work provides a first comprehensive structural overview over the complex interaction network between chemokine receptors and indicates how heterodimerization and the interaction with the membrane environment diversifies the function of closely related GPCRs.

Lipids Alter Rhodopsin Function via Ligand-like and Solvent-like Interactions

Rhodopsin, a prototypical G protein-coupled receptor, is a membrane protein that can sense dim light. This highly effective photoreceptor is known to be sensitive to the composition of its lipidic environment, but the molecular mechanisms underlying this fine-tuned modulation of the receptor's function and structural stability are not fully understood. There are two competing hypotheses to explain how this occurs: 1) lipid modulation occurs via solvent-like interactions, where lipid composition controls membrane properties like hydrophobic thickness, which in turn modulate the protein's conformational equilibrium; or 2) protein-lipid interactions are ligand-like, with specific hot spots and long-lived binding events. By analyzing an ensemble of all-atom molecular dynamics simulations of five different states of rhodopsin, we show that a local ordering effect takes place in the membrane upon receptor activation. Likewise, docosahexaenoic acid acyl tails and phosphatidylethanolamine headgroups behave like weak ligands, preferentially binding to the receptor in inactive-like conformations and inducing subtle but significant structural changes.

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.

Identification of Two New Cholesterol Interaction Sites on the A2A Adenosine Receptor

By mole, cholesterol is the most abundant component of animal cell plasma membranes. Many membrane proteins have been shown to be functionally dependent on cholesterol, several of which have also been shown to bind cholesterol at well-defined locations on their membrane-facing surface. In this work, a combination of coarse-grained "Martini" and all-atom simulations are used to identify two, to our knowledge, new cholesterol-binding sites on the A_2A adenosine receptor, a G-protein-coupled receptor that is a target for the treatment of Parkinson's disease. One of the sites is also observed to bind cholesterol in several recent, high-resolution crystal structures of the protein, and in the simulations, interacts with cholesterol only when bound to the inverse agonist ZM241385. Cataloguing cholesterol-binding sites is a vital step in the effort to understand cholesterol-dependent function of membrane proteins. Given that cholesterol content in plasma membranes varies with cell type and on administration of widely prescribed pharmaceuticals, such as statins, understanding cholesterol-dependent function is an important step toward exploiting membrane compositional variation for therapeutic purposes.