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

Assessing the Martini 3 protein model: A review of its path and potential

Coarse-grained (CG) protein models have become indispensable tools for studying many biological protein details, from conformational dynamics to the organization of protein macro-complexes, and even the interaction of proteins with other molecules. The Martini force field is one of the most widely used CG models for bio-molecular simulations, partly because of the enormous success of its protein model. With the recent release of a new and improved version of the Martini force field - Martini 3 - a new iteration of its protein model was also made available. The Martini 3 protein force field is an evolution of its Martini 2 counterpart, aimed at improving many of the shortcomings that had been previously identified. In this mini-review, we first provide a general overview of the model and then focus on the successful advances made in the short time since its release, many of which would not have been possible before. Furthermore, we discuss reported limitations, potential directions for model improvement and comment on what the likely future development and application avenues are.

The impact of volatile anesthetics and propofol on phosphatidylinositol 4,5-bisphosphate signaling

Phosphatidylinositol 4,5-bisphosphate (PIP2), as well as other anionic phospholipids, play a pivotal role in various cellular processes, including ion channel regulation, receptor trafficking, and intracellular signaling pathways. The binding of volatile anesthetics and propofol to PIP2 leads to alterations in PIP2-mediated signaling causing modulation of ion channels such as ɣ-aminobutyric acid type A (GABAA) receptors, voltage-gated calcium channels, and potassium channels through various mechanisms. Additionally, the interaction between anionic phospholipids and G protein-coupled receptors plays a critical role in various anesthetic pathways, with these anesthetic-induced changes impacting PIP2 levels which cause cascading effects on receptor trafficking, including GABAA receptor internalization. This comprehensive review of various mechanisms of interaction provides insights into the intricate interplay between PIP2 signaling and anesthetic-induced changes, shedding light on the molecular mechanisms underlying anesthesia.

Translational use of homing peptides: Tumor and placental targeting

HYPOTHESIS

Tissue-specific homing peptides have been shown to improve chemotherapeutic efficacy due to their trophism for tumor cells. Other sequences that selectively home to the placenta are providing new and safer therapeutics to treat complications in pregnancy. Our hypothesis is that the placental homing peptide RSGVAKS (RSG) may have binding affinity to cancer cells, and that insight can be gained into the binding mechanisms of RSG and the tumor homing peptide CGKRK to model membranes that mimic the primary lipid compositions of the respective cells.

EXPERIMENTS

Following cell culture studies on the binding efficacy of the peptides on a breast cancer cell line, a systematic translational characterization is delivered using ellipsometry, Brewster angle microscopy and neutron reflectometry of the extents, structures, and dynamics of the interactions of the peptides with the model membranes on a Langmuir trough.

FINDINGS

We start by revealing that RSG does indeed have binding affinity to breast cancer cells. The peptide is then shown to exhibit stronger interactions and greater penetration than CGKRK into both model membranes, combined with greater disruption to the lipid component. RSG also forms aggregates bound to the model membranes, yet both peptides bind to a greater extent to the placental than cancer model membranes. The results demonstrate the potential for varying local reservoirs of peptide within cell membranes that may influence receptor binding. The innovative nature of our findings motivates the urgent need for more studies involving multifaceted experimental platforms to explore the use of specific peptide sequences to home to different cellular targets.

GPCRs in the round: SMA-like copolymers and SMALPs as a platform for investigating GPCRs

G-protein-coupled receptors (GPCRs) are the largest family of membrane proteins, regulate a plethora of physiological responses and are the therapeutic target for 30-40% of clinically-prescribed drugs. They are integral membrane proteins deeply embedded in the plasma membrane where they activate intracellular signalling via coupling to G-proteins and β-arrestin. GPCRs are in intimate association with the bilayer lipids and that lipid environment regulates the signalling functions of GPCRs. This complex lipid 'landscape' is both heterogeneous and dynamic. GPCR function is modulated by bulk membrane properties including membrane fluidity, microdomains, curvature, thickness and asymmetry but GPCRs are also regulated by specific lipid:GPCR binding, including cholesterol and anionic lipids. Understanding the molecular mechanisms whereby GPCR signalling is regulated by lipids is a very active area of research currently. A major advance in membrane protein research in recent years was the application of poly(styrene-co-maleic acid) (SMA) copolymers. These spontaneously generate SMA lipid particles (SMALPs) encapsulating membrane protein in a nano-scale disc of cell membrane, thereby removing the historical need for detergent and preserving lipid:GPCR interaction. The focus of this review is how GPCR-SMALPs are increasing our understanding of GPCR structure and function at the molecular level. Furthermore, an increasing number of 'second generation' SMA-like copolymers have been reported recently. These are reviewed from the context of increasing our understanding of GPCR molecular mechanisms. Moreover, their potential as a novel platform for downstream biophysical and structural analyses is assessed and looking ahead, the translational application of SMA-like copolymers to GPCR drug discovery programmes in the future is considered.

Preferential electrostatic interactions of phosphatidic acid with arginines

Phosphatidic acid (PA) is an anionic lipid that preferentially interacts with proteins in a diverse set of cellular processes such as transport, apoptosis, and neurotransmission. One such interaction is that of the PA lipids with the proteins of voltage-sensitive ion channels. In comparison to several other similarly charged anionic lipids, PA lipids exhibit much stronger interactions. Intrigued and motivated by this finding, we sought out to gain deeper understanding into the electrostatic interactions of anionic lipids with charged proteins. Using the voltage sensor domain (VSD) of the KvAP channel as a model system, we performed long-timescale atomistic simulations to analyze the interactions of POPA, POPG, and POPI lipids with arginines (ARGs). Our simulations reveal two mechanisms. First, POPA is able to interact not only with surface ARGs but is able to snorkel and interact with a buried arginine. POPG and POPI lipids on the other hand show weak interactions even with both the surface and buried ARGs. Second, deprotonated POPA with -2 charge is able to break the salt-bridge connection between VSD protein segments and establish its own electrostatic bond with the ARG. Based on these findings, we propose a headgroup size hypothesis for preferential solvation of proteins by charged lipids. These findings may be valuable in understanding how PA lipids could be modulating kinematics of transmembrane proteins in cellular membranes.

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.

Palmitoylation of the Glucagon-like Peptide-1 Receptor Modulates Cholesterol Interactions at the Receptor-Lipid Microenvironment

The G protein-coupled receptor (GPCR) superfamily of cell surface receptors has been shown to be functionally modulated by post-translational modifications. The glucagon-like peptide receptor-1 (GLP-1R), which is a drug target in diabetes and obesity, undergoes agonist-dependent palmitoyl tail conjugation. The palmitoylation in the C-terminal domain of GLP-1R has been suggested to modulate the receptor-lipid microenvironment. In this work, we have performed coarse-grain molecular dynamics simulations of palmitoylated and nonpalmitoylated GLP-1R to analyze the differential receptor-lipid interactions. Interestingly, the placement and dynamics of the C-terminal domain of GLP-1R are found to be directly dependent on the palmitoyl tail. We observe that both cholesterol and phospholipids interact with the receptor but display differential interactions in the presence and absence of the palmitoyl tail. We characterize important cholesterol-binding sites and validate sites that have been previously reported in experimentally resolved structures of the receptor. We show that the receptor acts like a conduit for cholesterol flip-flop by stabilizing cholesterol in the membrane core. Taken together, our work represents an important step in understanding the molecular effects of lipid modifications in GPCRs.

Cellular Lipids Regulate the Conformational Ensembles of the Disordered Intracellular Loop 3 in β2 Adrenergic Receptor

The structurally disordered intracellular loops (ICLs) of G protein-coupled receptors (GPCRs) play a critical role in G protein coupling. In our previous work, we used a combination of FRET-based and computational methodologies to show that the third intracellular loop (ICL3) modulates the activity and G protein coupling selectivity in GPCRs. In the current study, we have uncovered the role of several lipid components in modulating the conformational ensemble of ICL3 of the β2-adrenergic receptor (β2AR). Our findings indicate that phosphatidylinositol 4,5-bisphosphate (PIP2) in the inner leaflet of the membrane bilayer acts as a stabilizing anchor for ICL3, opening the intracellular cavity to facilitate G protein coupling. This interaction between PIP2 and ICL3 causes tilting of β2AR within the cellular membrane. Notably, this tilting of the receptor is supported by ganglioside GM3 stabilizing the extracellular loops on the outer leaflet of the bilayer, thereby exerting an allosteric effect on the orthosteric ligand binding pocket. Our results underscore the significance of lipids in modulating GPCR activity, proposing an allosteric mechanism that occurs through the receptor's orientation within the membrane.

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

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

Structural basis of dimerization of chemokine receptors CCR5 and CXCR4

G protein-coupled receptors (GPCRs) are prominent drug targets responsible for extracellular-to-intracellular signal transduction. GPCRs can form functional dimers that have been poorly characterized so far. Here, we show the dimerization mechanism of the chemokine receptors CCR5 and CXCR4 by means of an advanced free-energy technique named coarse-grained metadynamics. Our results reproduce binding events between the GPCRs occurring in the minute timescale, revealing a symmetric and an asymmetric dimeric structure for each of the three investigated systems, CCR5/CCR5, CXCR4/CXCR4, and CCR5/CXCR4. The transmembrane helices TM4-TM5 and TM6-TM7 are the preferred binding interfaces for CCR5 and CXCR4, respectively. The identified dimeric states differ in the access to the binding sites of the ligand and G protein, indicating that dimerization may represent a fine allosteric mechanism to regulate receptor activity. Our study offers structural basis for the design of ligands able to modulate the formation of CCR5 and CXCR4 dimers and in turn their activity, with therapeutic potential against HIV, cancer, and immune-inflammatory diseases.

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

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

Molecular Biophysics of Class A G Protein Coupled Receptors-Lipids Interactome at a Glance-Highlights from the A2A Adenosine Receptor

G protein-coupled receptors (GPCRs) are embedded in phospholipid membrane bilayers with cholesterol representing 34% of the total lipid content in mammalian plasma membranes. Membrane lipids interact with GPCRs structures and modulate their function and drug-stimulated signaling through conformational selection. It has been shown that anionic phospholipids form strong interactions between positively charged residues in the G protein and the TM5-TM6-TM 7 cytoplasmic interface of class A GPCRs stabilizing the signaling GPCR-G complex. Cholesterol with a high content in plasma membranes can be identified in more specific sites in the transmembrane region of GPCRs, such as the Cholesterol Consensus Motif (CCM) and Cholesterol Recognition Amino Acid Consensus (CRAC) motifs and other receptor dependent and receptor state dependent sites. Experimental biophysical methods, atomistic (AA) MD simulations and coarse-grained (CG) molecular dynamics simulations have been applied to investigate these interactions. We emphasized here the impact of phosphatidyl inositol-4,5-bisphosphate (PtdIns(4,5)P₂ or PIP₂), a minor phospholipid component and of cholesterol on the function-related conformational equilibria of the human A_2A adenosine receptor (A_2AR), a representative receptor in class A GPCR. Several GPCRs of class A interacted with PIP₂ and cholesterol and in many cases the mechanism of the modulation of their function remains unknown. This review provides a helpful comprehensive overview for biophysics that enter the field of GPCRs-lipid systems.

Stabilization of pre-existing neurotensin receptor conformational states by β-arrestin-1 and the biased allosteric modulator ML314

The neurotensin receptor 1 (NTS1) is a G protein-coupled receptor (GPCR) with promise as a drug target for the treatment of pain, schizophrenia, obesity, addiction, and various cancers. A detailed picture of the NTS1 structural landscape has been established by X-ray crystallography and cryo-EM and yet, the molecular determinants for why a receptor couples to G protein versus arrestin transducers remain poorly defined. We used 13CεH3-methionine NMR spectroscopy to show that binding of phosphatidylinositol-4,5-bisphosphate (PIP2) to the receptor's intracellular surface allosterically tunes the timescale of motions at the orthosteric pocket and conserved activation motifs - without dramatically altering the structural ensemble. β-arrestin-1 further remodels the receptor ensemble by reducing conformational exchange kinetics for a subset of resonances, whereas G protein coupling has little to no effect on exchange rates. A β-arrestin biased allosteric modulator transforms the NTS1:G protein complex into a concatenation of substates, without triggering transducer dissociation, suggesting that it may function by stabilizing signaling incompetent G protein conformations such as the non-canonical state. Together, our work demonstrates the importance of kinetic information to a complete picture of the GPCR activation landscape.

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.

Setting up multicolour TIRF microscopy down to the single molecule level

Investigating biological mechanisms in ever greater detail requires continuous advances in microscopy techniques and setups. Total internal reflection fluorescence (TIRF) microscopy is a well-established technique for visualizing processes on the cell membrane. TIRF allows studies down to the single molecule level, mainly in single-colour applications. Instead, multicolour setups are still limited. Here, we describe our strategies for implementing a multi-channel TIRF microscopy system capable of simultaneous two-channel excitation and detection, starting from a single-colour commercial setup. First, we report some applications at high molecule density and then focus on the challenges we faced for achieving the single molecule level simultaneously in different channels, showing that rigorous optimizations on the setup are needed to increase its sensitivity up to this point, from camera setting to background minimization. We also discuss our strategies regarding crucial points of fluorescent labelling for this type of experiment: labelling strategy, kind of probe, efficiency, and orthogonality of the reaction, all of which are aspects that can influence the achievable results. This work may provide useful guidelines for setting up advanced single-molecule multi-channel TIRF experiments to obtain insights into interaction mechanisms on the cell membrane of living cells.

Lipid Modulation of a Class B GPCR: Elucidating the Modulatory Role of PI(4,5)P2 Lipids

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂) lipids have been shown to stabilize an active conformation of class A G-protein coupled receptors (GPCRs) through a conserved binding site, not present in class B GPCRs. For class B GPCRs, previous molecular dynamics (MD) simulation studies have shown PI(4,5)P₂ interacting with the Glucagon receptor (GCGR), which constitutes an important target for diabetes and obesity therapeutics. In this work, we applied MD simulations supported by native mass spectrometry (nMS) to study lipid interactions with GCGR. We demonstrate how tail composition plays a role in modulating the binding of PI(4,5)P₂ lipids to GCGR. Specifically, we find the PI(4,5)P₂ lipids to have a higher affinity toward the inactive conformation of GCGR. Interestingly, we find that in contrast to class A GPCRs, PI(4,5)P₂ appear to stabilize the inactive conformation of GCGR through a binding site conserved across class B GPCRs but absent in class A GPCRs. This suggests differences in the regulatory function of PI(4,5)P₂ between class A and class B GPCRs.

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.

Membrane phosphoinositides regulate GPCR-β-arrestin complex assembly and dynamics

Binding of arrestin to phosphorylated G protein-coupled receptors (GPCRs) is crucial for modulating signaling. Once internalized, some GPCRs remain complexed with β-arrestins, while others interact only transiently; this difference affects GPCR signaling and recycling. Cell-based and in vitro biophysical assays reveal the role of membrane phosphoinositides (PIPs) in β-arrestin recruitment and GPCR-β-arrestin complex dynamics. We find that GPCRs broadly stratify into two groups, one that requires PIP binding for β-arrestin recruitment and one that does not. Plasma membrane PIPs potentiate an active conformation of β-arrestin and stabilize GPCR-β-arrestin complexes by promoting a fully engaged state of the complex. As allosteric modulators of GPCR-β-arrestin complex dynamics, membrane PIPs allow for additional conformational diversity beyond that imposed by GPCR phosphorylation alone. For GPCRs that require membrane PIP binding for β-arrestin recruitment, this provides a mechanism for β-arrestin release upon translocation of the GPCR to endosomes, allowing for its rapid recycling.

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.

Crystal structure of a constitutive active mutant of adenosine A2A receptor

The adenosine A_2A receptor (A_2AAR) is a prototypical member of the class A subfamily of G-protein-coupled receptors (GPCRs) that is widely distributed in various tissues and organs of the human body, and participates in many important signal-regulation processes. We have previously summarized a common activation pathway of class A GPCRs in which a series of conserved residues/motifs undergo conformational change during extracellular agonist binding and finally induce the coupling of intracellular G protein. Through this mechanism we have successfully predicted several novel constitutive active or inactive mutations for A_2AAR. To reveal the molecular mechanism of mutation-induced constitutive activity, we determined the structure of a typical mutant I92N complexed with the agonist UK-432097. The mutated I92N forms a hydrophilic interaction network with nearby residues including Trp^6.48 of the CWxP motif, which is absent in wild-type A_2AAR. Although the mutant structure is similar overall to the previously determined intermediate-state A_2AAR structure (PDB ID 3qak) [Xu, Wu, Katritch, Han, Jacobson, Gao, Cherezov & Stevens (2011). Science, 332, 322-327 ▸], molecular dynamics simulations suggest that the I92N mutant stabilizes the metastable intermediate state through the hydrophilic interaction network and favors the conformational transition of the receptor towards the active state. This research provides a structural template towards the special pharmacological outcome triggered by conformational mutation and sheds light on future structural or pharmaco-logical studies among class A GPCRs.

Molecular Simulations and Drug Discovery of Adenosine Receptors

G protein-coupled receptors (GPCRs) represent the largest family of human membrane proteins. Four subtypes of adenosine receptors (ARs), the A1AR, A2AAR, A2BAR and A3AR, each with a unique pharmacological profile and distribution within the tissues in the human body, mediate many physiological functions and serve as critical drug targets for treating numerous human diseases including cancer, neuropathic pain, cardiac ischemia, stroke and diabetes. The A1AR and A3AR preferentially couple to the Gi/o proteins, while the A2AAR and A2BAR prefer coupling to the Gs proteins. Adenosine receptors were the first subclass of GPCRs that had experimental structures determined in complex with distinct G proteins. Here, we will review recent studies in molecular simulations and computer-aided drug discovery of the adenosine receptors and also highlight their future research opportunities.

Interaction With the Lipid Membrane Influences Fentanyl Pharmacology

Overdose deaths from fentanyl have reached epidemic proportions in the USA and are increasing worldwide. Fentanyl is a potent opioid agonist that is less well reversed by naloxone than morphine. Due to fentanyl’s high lipophilicity and elongated structure we hypothesised that its unusual pharmacology may be explained by its interactions with the lipid membrane on route to binding to the μ-opioid receptor (MOPr). Through coarse-grained molecular dynamics simulations, electrophysiological recordings and cell signalling assays, we determined how fentanyl and morphine access the orthosteric pocket of MOPr. Morphine accesses MOPr via the aqueous pathway; first binding to an extracellular vestibule, then diffusing into the orthosteric pocket. In contrast, fentanyl may take a novel route; first partitioning into the membrane, before accessing the orthosteric site by diffusing through a ligand-induced gap between the transmembrane helices. In electrophysiological recordings fentanyl-induced currents returned after washout, suggesting fentanyl deposits in the lipid membrane. However, mutation of residues forming the potential MOPr transmembrane access site did not alter fentanyl’s pharmacological profile in vitro. A high local concentration of fentanyl in the lipid membrane, possibly in combination with a novel lipophilic binding route, may explain the high potency and lower susceptibility of fentanyl to reversal by naloxone.

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.

Allosteric modulation of the adenosine A2A receptor by cholesterol

Cholesterol is a major component of the cell membrane and commonly regulates membrane protein function. Here, we investigate how cholesterol modulates the conformational equilibria and signaling of the adenosine A2A receptor (A2AR) in reconstituted phospholipid nanodiscs. This model system conveniently excludes possible effects arising from cholesterol-induced phase separation or receptor oligomerization and focuses on the question of allostery. GTP hydrolysis assays show that cholesterol weakly enhances the basal signaling of A2AR while decreasing the agonist EC50. Fluorine nuclear magnetic resonance (19F NMR) spectroscopy shows that this enhancement arises from an increase in the receptor's active state population and a G-protein-bound precoupled state. 19F NMR of fluorinated cholesterol analogs reveals transient interactions with A2AR, indicating a lack of high-affinity binding or direct allosteric modulation. The combined results suggest that the observed allosteric effects are largely indirect and originate from cholesterol-mediated changes in membrane properties, as shown by membrane fluidity measurements and high-pressure NMR.

Structures of neurokinin 1 receptor in complex with Gq and Gs proteins reveal substance P binding mode and unique activation features

The neurokinin 1 receptor (NK₁R) is involved in inflammation and pain transmission. This pathophysiologically important G protein–coupled receptor is predominantly activated by its cognate agonist substance P (SP) but also by the closely related neurokinins A and B. Here, we report cryo–electron microscopy structures of SP-bound NK₁R in complex with its primary downstream signal mediators, G_q and G_s. Our structures reveal how a polar network at the extracellular, solvent-exposed receptor surface shapes the orthosteric pocket and that NK₁R adopts a noncanonical active-state conformation with an interface for G protein binding, which is distinct from previously reported structures. Detailed comparisons with antagonist-bound NK₁R crystal structures reveal that insurmountable antagonists induce a distinct and long-lasting receptor conformation that sterically blocks SP binding. Together, our structures provide important structural insights into ligand and G protein promiscuity, the lack of basal signaling, and agonist- and antagonist-induced conformations in the neurokinin receptor family.

Modulation of adenosine A2a receptor oligomerization by receptor activation and PIP2 interactions

GPCRs have been shown to form oligomers, which generate distinctive signaling outcomes. However, the structural nature of the oligomerization process remains uncertain. We have characterized oligomeric configurations of the adenosine A2a receptor (A2aR) by combining large-scale molecular dynamics simulations with Markov state models. These oligomeric structures may also serve as templates for studying oligomerization of other class A GPCRs. Our simulation data revealed that receptor activation results in enhanced oligomerization, more diverse oligomer populations, and a more connected oligomerization network. The active state conformation of the A2aR shifts protein-protein association interfaces to those involving intracellular loop ICL3 and transmembrane helix TM6. Binding of PIP₂ to A2aR stabilizes protein-protein interactions via PIP₂-mediated association interfaces. These results indicate that A2aR oligomerization is responsive to the local membrane lipid environment. This, in turn, suggests a modulatory effect on A2aR whereby a given oligomerization profile favors the dynamic formation of specific supramolecular signaling complexes.

Insights into the Role of Membrane Lipids in the Structure, Function and Regulation of Integral Membrane Proteins

Membrane proteins exist within the highly hydrophobic membranes surrounding cells and organelles, playing key roles in cellular function. It is becoming increasingly clear that the membrane does not just act as an appropriate environment for these proteins, but that the lipids that make up these membranes are essential for membrane protein structure and function. Recent technological advances in cryogenic electron microscopy and in advanced mass spectrometry methods, as well as the development of alternative membrane mimetic systems, have allowed experimental study of membrane protein–lipid complexes. These have been complemented by computational approaches, exploiting the ability of Molecular Dynamics simulations to allow exploration of membrane protein conformational changes in membranes with a defined lipid content. These studies have revealed the importance of lipids in stabilising the oligomeric forms of membrane proteins, mediating protein–protein interactions, maintaining a specific conformational state of a membrane protein and activity. Here we review some of the key recent advances in the field of membrane protein–lipid studies, with major emphasis on respiratory complexes, transporters, channels and G-protein coupled receptors.

A novel high-throughput screen for identifying lipids that stabilise membrane proteins in detergent based solution

Membrane proteins have a range of crucial biological functions and are the target of about 60% of all prescribed drugs. For most studies, they need to be extracted out of the lipid-bilayer, e.g. by detergent solubilisation, leading to the loss of native lipids, which may disturb important protein-lipid/bilayer interactions and thus functional and structural integrity. Relipidation of membrane proteins has proven extremely successful for studying challenging targets, but the identification of suitable lipids can be expensive and laborious. Therefore, we developed a screen to aid the high-throughput identification of beneficial lipids. The screen covers a large lipid space and was designed to be suitable for a range of stability assessment methods. Here, we demonstrate its use as a tool for identifying stabilising lipids for three membrane proteins: a bacterial pyrophosphatase (Tm-PPase), a fungal purine transporter (UapA) and a human GPCR (A_2AR). A_2AR is stabilised by cholesteryl hemisuccinate, a lipid well known to stabilise GPCRs, validating the approach. Additionally, our screen also identified a range of new lipids which stabilised our test proteins, providing a starting point for further investigation and demonstrating its value as a novel tool for membrane protein research. The pre-dispensed screen will be made commercially available to the scientific community in future and has a number of potential applications in the field.

Structures of rhodopsin in complex with G-protein-coupled receptor kinase 1

G-protein-coupled receptor (GPCR) kinases (GRKs) selectively phosphorylate activated GPCRs, thereby priming them for desensitization¹. Although it is unclear how GRKs recognize these receptors^2-4, a conserved region at the GRK N terminus is essential for this process^5-8. Here we report a series of cryo-electron microscopy single-particle reconstructions of light-activated rhodopsin (Rho*) bound to rhodopsin kinase (GRK1), wherein the N terminus of GRK1 forms a helix that docks into the open cytoplasmic cleft of Rho*. The helix also packs against the GRK1 kinase domain and stabilizes it in an active configuration. The complex is further stabilized by electrostatic interactions between basic residues that are conserved in most GPCRs and acidic residues that are conserved in GRKs. We did not observe any density for the regulator of G-protein signalling homology domain of GRK1 or the C terminus of rhodopsin. Crosslinking with mass spectrometry analysis confirmed these results and revealed dynamic behaviour in receptor-bound GRK1 that would allow the phosphorylation of multiple sites in the receptor tail. We have identified GRK1 residues whose mutation augments kinase activity and crosslinking with Rho*, as well as residues that are involved in activation by acidic phospholipids. From these data, we present a general model for how a small family of protein kinases can recognize and be activated by hundreds of different GPCRs.

Allosteric modulation of ghrelin receptor signaling by lipids

The membrane is an integral component of the G protein-coupled receptor signaling machinery. Here we demonstrate that lipids regulate the signaling efficacy and selectivity of the ghrelin receptor GHSR through specific interactions and bulk effects. We find that PIP2 shifts the conformational equilibrium of GHSR away from its inactive state, favoring basal and agonist-induced G protein activation. This occurs because of a preferential binding of PIP2 to specific intracellular sites in the receptor active state. Another lipid, GM3, also binds GHSR and favors G protein activation, but mostly in a ghrelin-dependent manner. Finally, we find that not only selective interactions but also the thickness of the bilayer reshapes the conformational repertoire of GHSR, with direct consequences on G protein selectivity. Taken together, this data illuminates the multifaceted role of the membrane components as allosteric modulators of how ghrelin signal could be propagated.

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).

Bacterial Analogs to Cholesterol Affect Dimerization of Proteorhodopsin and Modulates Preferred Dimer Interface

Hopanoids, the bacterial analogues of sterols, are ubiquitous in bacteria and play a significant role in organismal survival under stressful environments. Unlike sterols, hopanoids have a high degree of variation in the size and chemical nature of the substituent attached to the ring moiety, leading to different effects on the structure and dynamics of biological membranes. While it is understood that hopanoids can indirectly tune membrane physical properties, little is known on the role that hopanoids may play in affecting the organization and behavior of bacterial membrane proteins. In this work we used coarse-grained molecular dynamics simulations to characterize the effects of two hopanoids, diploptene (DPT) and bacteriohopanetetrol (BHT), on the oligomerization of proteorhodopsin (PR) in a model membrane composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phophoethanolamine (POPE) and 1-palmitoyl-2-oleoyl-sn-3-phosphoglycerol (POPG). PR is a bacterial membrane protein that functions as a light-activated proton pump. We chose PR based on its ability to adopt a distribution of oligomeric states in different membrane environments. Furthermore, the efficiency of proton pumping in PR is intimately linked to its organization into oligomers. Our results reveal that both BHT and DPT indirectly affect dimerization by tuning membrane properties in a fashion that is concentration-dependent. Variation in their interaction with PR in the membrane-embedded and the cytoplasmic regions leads to distinctly different effects on the plasticity of the dimer interface. BHT has the ability to intercalate between monomers in the dimeric interface, whereas DPT shifts dimerization interactions via packing of the interleaflet region of the membrane. Our results show a direct relationship between hopanoid structure and lateral organization of PR, providing a first glimpse at how these bacterial analogues to eukaryotic sterols produce very similar biophysical effects within the cell membrane.

Activation of G-protein-coupled receptors is thermodynamically linked to lipid solvation

Preferential lipid solvation of the G-protein-coupled A_2A adenosine receptor (A_2AR) is evaluated from 35 μs of all-atom molecular dynamics simulation. A coarse-grained transition matrix algorithm is developed to overcome slow equilibration of the first solvation shell, obtaining estimates of the free energy of solvation by different lipids for the receptor in different activation states. Results indicate preference for solvation by unsaturated chains, which favors the active receptor. A model for lipid-dependent G-protein-coupled receptor activity is proposed in which the chemical potential of lipids in the bulk membrane modulates receptor activity. The entropies associated with moving saturated and unsaturated lipids from bulk to A_2AR's first solvation shell are evaluated. Overall, the acyl chains are more disordered (i.e., obtain a favorable entropic contribution) when partitioning to the receptor surface, and this effect is augmented for the saturated chains, which are relatively more ordered in bulk.

State-dependent protein-lipid interactions of a pentameric ligand-gated ion channel in a neuronal membrane

Pentameric ligand-gated ion channels (pLGICs) are receptor proteins that are sensitive to their membrane environment, but the mechanism for how lipids modulate function under physiological conditions in a state dependent manner is not known. The glycine receptor is a pLGIC whose structure has been resolved in different functional states. Using a realistic model of a neuronal membrane coupled with coarse-grained molecular dynamics simulations, we demonstrate that some key lipid-protein interactions are dependent on the receptor state, suggesting that lipids may regulate the receptor's conformational dynamics. Comparison with existing structural data confirms known lipid binding sites, but we also predict further protein-lipid interactions including a site at the communication interface between the extracellular and transmembrane domain. Moreover, in the active state, cholesterol can bind to the binding site of the positive allosteric modulator ivermectin. These protein-lipid interaction sites could in future be exploited for the rational design of lipid-like allosteric drugs.

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

Endophilin recruitment drives membrane curvature generation through coincidence detection of GPCR loop interactions and negative lipid charge

Endophilin plays key roles during endocytosis of cellular receptors, including generating membrane curvature to drive internalization. Electrostatic interactions between endophilin's BIN/Amphiphysin/Rvs domain and anionic membrane lipids have been considered the major driving force in curvature generation. However, the SH3 domain of endophilin also interacts with the proline-rich third intracellular loop (TIL) of various G-protein-coupled receptors (GPCRs), and it is unclear whether this interaction has a direct role in generating membrane curvature during endocytosis. To examine this, we designed model membranes with a membrane density of 1400 receptors per μm2 represented by a covalently conjugated TIL region from the β1-adrenergic receptor. We observed that TIL recruits endophilin to membranes composed of 95 mol% of zwitterionic lipids via the SH3 domain. More importantly, endophilin recruited via TIL tubulates vesicles and gets sorted onto highly curved membrane tubules. These observations indicate that the cellular membrane bending and curvature sensing activities of endophilin can be facilitated through detection of the TIL of activated GPCRs in addition to binding to anionic lipids. Furthermore, we show that TIL electrostatically interacts with membranes composed of anionic lipids. Therefore, anionic lipids can modulate TIL/SH3 domain binding. Overall, our findings imply that an interplay between TIL, charged membrane lipids, BAR domain, and SH3 domain could exist in the biological system and that these components may act in coordination to regulate the internalization of cellular receptors.

Artificial Intelligence-Based Application to Explore Inhibitors of Neurodegenerative Diseases

Neuroinflammation is a common factor in neurodegenerative diseases, and it has been demonstrated that galectin-3 activates microglia and astrocytes, leading to inflammation. This means that inhibition of galectin-3 may become a new strategy for the treatment of neurodegenerative diseases. Based on this motivation, the objective of this study is to explore an integrated new approach for finding lead compounds that inhibit galectin-3, by combining universal artificial intelligence algorithms with traditional drug screening methods. Based on molecular docking method, potential compounds with high binding affinity were screened out from Chinese medicine database. Manifold artificial intelligence algorithms were performed to validate the docking results and further screen compounds. Among all involved predictive methods, the deep learning-based algorithm made 500 modeling attempts, and the square correlation coefficient of the best trained model on the test sets was 0.9. The XGBoost model reached a square correlation coefficient of 0.97 and a mean square error of only 0.01. We switched to the ZINC database and performed the same experiment, the results showed that the compounds in the former database showed stronger affinity. Finally, we further verified through molecular dynamics simulation that the complex composed of the candidate ligand and the target protein showed stable binding within 100 ns of simulation time. In summary, combined with the application based on artificial intelligence algorithms, we unearthed the active ingredients 1,2-Dimethylbenzene and Typhic acid contained in Crataegus pinnatifida and Typha angustata might be the effective inhibitors of neurodegenerative diseases. The high prediction accuracy of the models shows that it has practical application value on small sample data sets such as drug screening.

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

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

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.