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

Activation Microswitches in Adenosine Receptor A2A Function as Rheostats in the Cell Membrane

Although multiple components of the cell membrane modulate the stability and activation of G protein-coupled receptors (GPCRs), insights into the dynamics of GPCR structures come from biophysical studies conducted in detergents. This is because of the challenges of studying activation in a multicomponent lipid bilayer. To understand the role of cellular membrane lipids and cations in GPCR activation, we performed multiscale molecular dynamics simulations (56 μs) on three different conformational states of adenosine receptor A2AR, in both the cell membrane-like lipid bilayer and in detergent micelles. Molecular dynamics (MD) simulations show that the phosphatidylinositol bisphosphate (PIP2) interacts with the basic residues in the intracellular regions of A2AR, thereby reducing the flexibility of the receptor in the inactive state and limiting the transition to the active-intermediate state. In the G protein-coupled fully active state, PIP2 stabilizes the GPCR:G protein complex. Such stiffening effects are absent in non-ionic detergent micelles, and therefore, more transitions have been observed in detergents. The inter-residue distances that change significantly upon GPCR activation are known as activation microswitches. The activation microswitches show different levels of activation in the cell membrane, in the pure POPC bilayer, and in detergents. Thus, the temporal heat map of different activation microswitches calculated from the MD simulations suggests a rheostat model of GPCR activation microswitches rather than the binary switch model. These simulation results connect the chemistry of cell membrane lipids to receptor activity, which is useful for the design of detergents mimicking the cell membrane.

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.

Allosteric interactions in the parathyroid hormone GPCR-arrestin complex formation

Peptide ligands of class B G-protein-coupled receptors act via a two-step binding process, but the essential mechanisms that link their extracellular binding to intracellular receptor-arrestin interactions are not fully understood. Using NMR, crosslinking coupled to mass spectrometry, signaling experiments and computational approaches on the parathyroid hormone (PTH) type 1 receptor (PTHR), we show that initial binding of the PTH C-terminal part constrains the conformation of the flexible PTH N-terminal signaling epitope before a second binding event occurs. A 'hot-spot' PTH residue, His9, that inserts into the PTHR transmembrane domain at this second step allosterically engages receptor-arrestin coupling. A conformational change in PTHR intracellular loop 3 permits favorable interactions with β-arrestin's finger loop. These results unveil structural determinants for PTHR-arrestin complex formation and reveal that the two-step binding mechanism proceeds via cooperative fluctuations between ligand and receptor, which extend to other class B G-protein-coupled receptors.

Structural insights into emergent signaling modes of G protein-coupled receptors

G protein-coupled receptors (GPCRs) represent the largest family of cell membrane proteins, with >800 GPCRs in humans alone, and recognize highly diverse ligands, ranging from photons to large protein molecules. Very important to human medicine, GPCRs are targeted by about 35% of prescription drugs. GPCRs are characterized by a seven-transmembrane α-helical structure, transmitting extracellular signals into cells to regulate major physiological processes via heterotrimeric G proteins and β-arrestins. Initially viewed as receptors whose signaling via G proteins is delimited to the plasma membrane, it is now recognized that GPCRs signal also at various intracellular locations, and the mechanisms and (patho)physiological relevance of such signaling modes are actively investigated. The propensity of GPCRs to adopt different signaling modes is largely encoded in the structural plasticity of the receptors themselves and of their signaling complexes. Here, we review emerging modes of GPCR signaling via endosomal membranes and the physiological implications of such signaling modes. We further summarize recent structural insights into mechanisms of GPCR activation and signaling. We particularly emphasize the structural mechanisms governing the continued GPCR signaling from endosomes and the structural aspects of the GPCR resensitization mechanism and discuss the recently uncovered and important roles of lipids in these processes.

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.

In Silico Drug Design for Purinergic GPCRs: Overview on Molecular Dynamics Applied to Adenosine and P2Y Receptors

Molecular modeling has contributed to drug discovery for purinergic GPCRs, including adenosine receptors (ARs) and P2Y receptors (P2YRs). Experimental structures and homology modeling have proven to be useful in understanding and predicting structure activity relationships (SAR) of agonists and antagonists. This review provides an excursus on molecular dynamics (MD) simulations applied to ARs and P2YRs. The binding modes of newly synthesized A1AR- and A3AR-selective nucleoside derivatives, potentially of use against depression and inflammation, respectively, have been predicted to recapitulate their SAR and the species dependence of A3AR affinity. P2Y12R and P2Y1R crystallographic structures, respectively, have provided a detailed understanding of the recognition of anti-inflammatory P2Y14R antagonists and a large group of allosteric and orthosteric antagonists of P2Y1R, an antithrombotic and neuroprotective target. MD of A2AAR (an anticancer and neuroprotective target), A3AR, and P2Y1R has identified microswitches that are putatively involved in receptor activation. The approach pathways of different ligands toward A2AAR and P2Y1R binding sites have also been explored. A1AR, A2AAR, and A3AR were utilizes to study allosteric phenomena, but locating the binding site of structurally diverse allosteric modulators, such as an A3AR enhancer LUF6000, is challenging. Ligand residence time, a predictor of in vivo efficacy, and the structural role of water were investigated through A2AAR MD simulations. Thus, new MD and other modeling algorithms have contributed to purinergic GPCR drug discovery.

Insights into adenosine A2A receptor activation through cooperative modulation of agonist and allosteric lipid interactions

The activation process of G protein-coupled receptors (GPCRs) has been extensively studied, both experimentally and computationally. In particular, Molecular Dynamics (MD) simulations have proven useful in exploring GPCR conformational space. The typical behaviour of class A GPCRs, when subjected to unbiased MD simulations from their crystallized inactive state, is to fluctuate between inactive and intermediate(s) conformations, even with bound agonist. Fully active conformation(s) are rarely stabilized unless a G protein is also bound. Despite several crystal structures of the adenosine A2a receptor (A2aR) having been resolved in complex with co-crystallized agonists and Gs protein, its agonist-mediated activation process is still not completely understood. In order to thoroughly examine the conformational landscape of A2aR activation, we performed unbiased microsecond-length MD simulations in quadruplicate, starting from the inactive conformation either in apo or with bound agonists: endogenous adenosine or synthetic NECA, embedded in two homogeneous phospholipid membranes: 1,2-dioleoyl-sn-glycerol-3-phosphoglycerol (DOPG) or 1,2-dioleoyl-sn-glycerol-3-phosphocholine (DOPC). In DOPC with bound adenosine or NECA, we observe transition to an intermediate receptor conformation consistent with the known adenosine-bound crystal state. In apo state in DOPG, two different intermediate conformations are obtained. One is similar to that observed with bound adenosine in DOPC, while the other is closer to the active state but not yet fully active. Exclusively, in DOPG with bound adenosine or NECA, we reproducibly identify receptor conformations with fully active features, which are able to dock Gs protein. These different receptor conformations can be attributed to the action/absence of agonist and phospholipid-mediated allosteric effects on the intracellular side of the receptor.

Novel roles of phosphoinositides in signaling, lipid transport, and disease

Phosphoinositides (PPIns) are lipid signaling molecules that act as master regulators of cellular signaling. Recent studies have revealed novel roles of PPIns in myriad cellular processes and multiple human diseases mediated by misregulation of PPIn signaling. This review will present a timely summary of recent discoveries in PPIn biology, specifically their role in regulating unexpected signaling pathways, modification of signaling outcomes downstream of integral membrane proteins, and novel roles in lipid transport. This has revealed new roles of PPIns in regulating membrane trafficking, immunity, cell polarity, and response to extracellular signals. A specific focus will be on novel opportunities to target PPIn metabolism for treatment of human diseases, including cancer, pathogen infection, developmental disorders, and immune disorders.

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.

G-Protein-Coupled Receptor-Membrane Interactions Depend on the Receptor Activation State

G-protein-coupled receptors (GPCRs) are the largest family of human membrane proteins and serve as primary targets of approximately one-third of currently marketed drugs. In particular, adenosine A1 receptor (A1 AR) is an important therapeutic target for treating cardiac ischemia-reperfusion injuries, neuropathic pain, and renal diseases. As a prototypical GPCR, the A1 AR is located within a phospholipid membrane bilayer and transmits cellular signals by changing between different conformational states. It is important to elucidate the lipid-protein interactions in order to understand the functional mechanism of GPCRs. Here, all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method were performed on both the inactive (antagonist bound) and active (agonist and G-protein bound) A1 AR, which was embedded in a 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) lipid bilayer. In the GaMD simulations, the membrane lipids played a key role in stabilizing different conformational states of the A1 AR. Our simulations further identified important regions of the receptor that interacted distinctly with the lipids in highly correlated manner. Activation of the A1 AR led to differential dynamics in the upper and lower leaflets of the lipid bilayer. In summary, GaMD enhanced simulations have revealed strongly coupled dynamics of the GPCR and lipids that depend on the receptor activation state. © 2019 Wiley Periodicals, Inc.

Energy Landscapes Reveal Agonist Control of G Protein-Coupled Receptor Activation via Microswitches

Agonist binding to G protein-coupled receptors (GPCRs) leads to conformational changes in the transmembrane region that activate cytosolic signaling pathways. Although high-resolution structures of different receptor states are available, atomistic details of allosteric signaling across the membrane remain elusive. We calculated free energy landscapes of β₂ adrenergic receptor activation using atomistic molecular dynamics simulations in an optimized string of swarms framework, which shed new light on how microswitches govern the equilibrium between conformational states. Contraction of the extracellular binding site in the presence of the agonist BI-167107 is obligatorily coupled to conformational changes in a connector motif located in the core of the transmembrane region. The connector is probabilistically coupled to the conformation of the intracellular region. An active connector promotes desolvation of a buried cavity, a twist of the conserved NPxxY motif, and an interaction between two conserved tyrosines in transmembrane helices 5 and 7 (Y–Y motif), which lead to a larger population of active-like states at the G protein binding site. This coupling is augmented by protonation of the strongly conserved Asp79^2.50. The agonist binding site hence communicates with the intracellular region via a cascade of locally connected microswitches. Characterization of these can be used to understand how ligands stabilize distinct receptor states and contribute to development drugs with specific signaling properties. The developed simulation protocol can likely be transferred to other class A GPCRs.

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

Reconstruction of apo A2A receptor activation pathways reveal ligand-competent intermediates and state-dependent cholesterol hotspots

G-protein coupled receptors (GPCRs) play a pivotal role in transmitting signals at the cellular level. Structural insights can be exploited to support GPCR structure-based drug discovery endeavours. Despite advances in GPCR crystallography, active state structures are scarce. Molecular dynamics (MD) simulations have been used to explore the conformational landscape of GPCRs. Efforts have been made to retrieve active state conformations starting from inactive structures, however to date this has not been possible without using an energy bias. Here, we reconstruct the activation pathways of the apo adenosine receptor (A2A), starting from an inactive conformation, by applying adaptive sampling MD combined with a goal-oriented scoring function. The reconstructed pathways reconcile well with experiments and help deepen our understanding of A2A regulatory mechanisms. Exploration of the apo conformational landscape of A2A reveals the existence of ligand-competent states, active intermediates and state-dependent cholesterol hotspots of relevance for drug discovery. To the best of our knowledge this is the first time an activation process has been elucidated for a GPCR starting from an inactive structure only, using a non-biased MD approach, opening avenues for the study of ligand binding to elusive yet pharmacologically relevant GPCR states.

Local Enrichment of Unsaturated Chains around the A2A Adenosine Receptor

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

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.

Allostery in G protein-coupled receptors investigated by molecular dynamics simulations

G-protein-coupled receptors (GPCRs) are allosteric signaling machines that trigger distinct functional responses depending on the particular conformational state they adopt upon binding. This so-called GPCR functional selectivity is prompted by ligands of different efficacy binding at orthosteric or allosteric sites on the receptor, as well as by interactions with intracellular protein partners or other receptor types. Molecular dynamics (MD) simulations can provide important mechanistic, thermodynamic, and kinetic insights into these interactions at a level of molecular detail that is necessary to rightly inform modern drug discovery. Here, we review the most recent MD contributions to understanding GPCR allostery, with an emphasis on their strengths and limitations.

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.