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

Probing the energy landscape of the lipid interactions of the serotonin1A receptor

G protein-coupled receptors (GPCRs) are lipid-regulated transmembrane proteins that play a central role in cell signaling and pharmacology. Although the role of membrane lipids in GPCR function is well established, the underlying GPCR-lipid interactions have not been thermodynamically characterized due to the complexity of these interactions. In this work, we estimate the energetics and dynamics of lipid association from coarse-grain simulations of the serotonin1A receptor embedded in a complex membrane. We show that lipids bind to the receptor with varying energetics of 1-4 kT, and timescales of 1-10 μs. The most favorable energetics and longest residence times are observed for cholesterol, glycosphingolipid GM1, phosphatidylethanolamine (PE) and phosphatidylserine (PS) lipids. Multi-exponential fitting of the contact probability suggests distinct dynamic regimes, corresponding to ps, ns and μs timescales, that we correlate with the annular, intermediate and non-annular lipid sites. The timescales of lipid binding correspond to high barrier heights, despite their relatively weaker energetics. Our results highlight that GPCR-lipid interactions are driven by both thermodynamic interactions and the dynamical features of lipid binding.

Bilayer lipids modulate ligand binding to atypical chemokine receptor 3

Chemokine receptors belong to the large class of G protein-coupled receptors (GPCRs) and are involved in a number of (patho)physiological processes. Previous studies highlighted the importance of membrane lipids for modulating GPCR structure and function. However, the underlying mechanisms of how lipids regulate GPCRs are often poorly understood. Here, we report that anionic lipid bilayers increase the binding affinity of the chemokine CXCL12 for the atypical chemokine receptor 3 (ACKR3) by modulating the CXCL12 binding kinetics. Notably, the anionic bilayer favors CXCL12 over the more positively charged chemokine CXCL11, which we explained by bilayer interactions orienting CXCL12 but not CXCL11 for productive ACKR3 binding. Furthermore, our data suggest a stabilization of active ACKR3 conformations in anionic bilayers. Taken together, the described regulation of chemokine selectivity of ACKR3 by the lipid bilayer proposes an extended version of the classical model of chemokine binding including the lipid environment of the receptor.

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.

Curvature Footprints of Transmembrane Proteins in Simulations with the Martini Force Field

Membranes play essential roles in biological systems and are tremendously diverse in the topologies and chemical and elastic properties that define their functions. In many cases, a given membrane may display considerable heterogeneity, with localized clusters of lipids and proteins exhibiting distinct characteristics compared to adjoining regions. These lipid-protein assemblies can span nanometers to micrometers and are associated with cellular processes such as transport and signaling. While lipid-protein assemblages are dynamic, they can be stabilized by coupling between local membrane composition and shape. Due to the inherent difficulty in resolving atomistic details of membrane proteins in their native lipid environments, these complexes are notoriously challenging to study experimentally; however, molecular dynamics (MD) simulations might be a viable alternative. Here, we aim to assess the utility of coarse-grained (CG) MD simulations with the Martini force field for studying membrane curvature induced by transmembrane (TM) proteins that are reported to generate local curvature. The direction and magnitude of curvature induced by five different TM proteins, as well as certain lipid-protein and protein-protein interactions, were found to be in good agreement with available reference data.

Cellular lipids regulate the conformational ensembles of the disordered intracellular loop 3 in β2-adrenergic receptor

The intracellular loops of G protein-coupled receptors (GPCRs) have been shown to play a key role in G protein coupling and selectivity. We recently showed that the intrinsically disordered third intracellular loop (ICL3) of β2-adrenergic receptor is dynamic and equilibrates between open and closed conformations to regulate the G protein coupling. In this study, using the extensive molecular dynamics simulations in multi-lipid bilayer models, we show that the lipid phosphatidylinositol 4,5-bisphosphate (PIP2) stabilizes the active state of β2-adrenergic receptor by keeping ICL3 in an open conformation. This stabilization results in a tilt of the receptor within the membrane. Additionally, the ganglioside lipid, GM3 interacts with extracellular loops, impacting the ligand binding site allosterically. This demonstrates the active role of the chemistry of lipids in stabilizing specific GPCR conformations.

Molecular dynamics simulations of lipid-protein interactions in SLC4 proteins

The SLC4 family of secondary bicarbonate transporters is responsible for the transport of HCO3-, CO32-, Cl-, Na+, K+, NH3, and H+, which are necessary for regulation of pH and ion homeostasis. They are widely expressed in numerous tissues throughout the body and function in different cell types with different membrane properties. Potential lipid roles in SLC4 function have been reported in experimental studies, focusing mostly on two members of the family: AE1 (Cl-/HCO3- exchanger) and NBCe1 (Na+-CO32-cotransporter). Previous computational studies of the outward-facing state of AE1 with model lipid membranes revealed enhanced protein-lipid interactions between cholesterol (CHOL) and phosphatidylinositol bisphosphate (PIP2). However, the protein-lipid interactions in other members of the family and other conformation states are still poorly understood and this precludes the detailed studies of a potential regulatory role for lipids in the SLC4 family. In this work, we performed coarse-grained and atomistic molecular dynamics simulations on three members of the SLC4 family with different transport modes: AE1, NBCe1, and NDCBE (an Na+-CO32-/Cl- exchanger), in model HEK293 membranes consisting of CHOL, PIP2, phosphatidylcholine, phosphatidylethanolamine, phosphatidylserine, and sphingomyelin. The recently resolved inward-facing state of AE1 was also included in the simulations. Lipid-protein contact analysis of the simulated trajectories was performed with the ProLint server, which provides a multitude of visualization tools for illustration of areas of enhanced lipid-protein contact and identification of putative lipid binding sites within the protein matrix. We observed enrichment of CHOL and PIP2 around all proteins with subtle differences in their distribution depending on the protein type and conformation state. Putative binding sites were identified for CHOL, PIP2, phosphatidylcholine, and sphingomyelin in the three studied proteins, and their potential roles in the SLC4 transport function, conformational transition, and protein dimerization are discussed.

Cholesterol Dependence of the Conformational Changes in Metabotropic Glutamate Receptor 1

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

Electron videography of a lipid-protein tango

Biological phenomena, from enzymatic catalysis to synaptic transmission, originate in the structural transformations of biomolecules and biomolecular assemblies in liquid water. However, directly imaging these nanoscopic dynamics without probes or labels has been a fundamental methodological challenge. Here, we developed an approach for "electron videography"-combining liquid phase electron microscopy with molecular modeling-with which we filmed the nanoscale structural fluctuations of individual, suspended, and unlabeled membrane protein nanodiscs in liquid. Systematic comparisons with biochemical data and simulation indicate the graphene encapsulation involved can afford sufficiently reduced effects of the illuminating electron beam for these observations to yield quantitative fingerprints of nanoscale lipid-protein interactions. Our results suggest that lipid-protein interactions delineate dynamically modified membrane domains across unexpectedly long ranges. Moreover, they contribute to the molecular mechanics of the nanodisc as a whole in a manner specific to the protein within. Overall, this work illustrates an experimental approach to film, quantify, and understand biomolecular dynamics at the nanometer scale.

Analyzing lipid distributions and curvature in molecular dynamics simulations of complex membranes

We describe methods to analyze lipid distributions and curvature in membranes with complex lipid mixtures and embedded membrane proteins. We discuss issues involved in these analyses, available tools to calculate curvature preferences of lipids and proteins, and focus on tools developed in our group for visual analysis of lipid-protein interactions and the analysis of membrane curvature.

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

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

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.

High-throughput virtual screening of potential inhibitors of GPR52 using docking and biased sampling method for Huntington's disease therapy

Huntington's disease (HD) is a rare and progressive neurodegenerative disorder caused by polyglutamine (poly-Q) mutations of the huntingtin (HTT) gene resulting in chorea, cognitive, and psychiatric dysfunctions. Being a monogenic condition, reducing the levels of the mutated huntingtin protein (mHTT) holds promise as an effective therapeutic approach. GPR52, an orphan G-protein coupled receptor (GPCR), enriched in the striatum, is a novel target for slowing down the progression of HD by lowering the mHTT levels. Therefore, the study focuses on identifying potent small-molecule inhibitors for GPR52 using a combination of robust high-throughput virtual screening (HTVS) and pharmacokinetics profiling followed by fast pulling of ligand (FPL) and umbrella sampling (US) simulations. Initially, screening a library of 2,36,545 compounds was done against the binding pocket of GPR52. Based on binding affinity, stereochemical and non-bonded interactions, and pharmacokinetic profiling, 50 compounds were shortlisted. Selected hit compounds 1, 2, and 3 were subjected to FPL simulations with applied external bias potential to investigate their unique dissociation pathways and intermolecular interactions over time. Subsequently, the US simulations were performed on the selected hit compounds to estimate their binding free energy (ΔG). The analysis of the trajectories obtained from simulations revealed that the residues TYR34, TYR185, GLY187, ASP188, ILE189, SER299, PHE300, and THR303 within the active site of GPR52 were significant for efficient ligand binding through the formation of various hydrogen bond interactions and hydrophobic contacts. Out of the three hit compounds, compound 3 had the lowest ΔG of - 20.82 ± 0.44 kcal/mol. The study identified compounds 1, 2, and 3 as potential molecules that can be developed as GPR52 inhibitors holding promise for lowering mHTT levels.

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.

Pragmatic Coarse-Graining of Proteins: Models and Applications

The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.

Coarse-grained molecular dynamics simulations of lipid-protein interactions in SLC4 proteins

The SLC4 family of secondary bicarbonate transporters is responsible for the transport of HCO 3 -, CO 3 2- , Cl - , Na + , K + , NH 3 and H + necessary for regulation of pH and ion homeostasis. They are widely expressed in numerous tissues throughout the body and function in different cell types with different membrane properties. Potential lipid roles in SLC4 function have been reported in experimental studies, focusing mostly on two members of the family: AE1 (Cl - /HCO 3 - exchanger) and NBCe1 (Na + -CO 3 2- cotransporter). Previous computational studies of the outward facing (OF) state of AE1 with model lipid membranes revealed enhanced protein-lipid interactions between cholesterol (CHOL) and phosphatidylinositol bisphosphate (PIP2). However, the protein-lipid interactions in other members of the family and other conformation states are still poorly understood and this precludes the detailed studies of a potential regulatory role for lipids in the SLC4 family. In this work, we performed multiple 50 µs coarse-grained molecular dynamics simulations on three members of the SLC4 family with different transport modes: AE1, NBCe1 and NDCBE (a Na + -CO 3 2- /Cl - exchanger), in model HEK293 membranes consisting of CHOL, PIP2, phosphatidylcholine (POPC), phosphatidylethanolamine (POPE), phosphatidylserine (POPS), and sphingomyelin (POSM). The recently resolved inward-facing (IF) state of AE1 was also included in the simulations. Lipid-protein contact analysis of the simulated trajectories was performed with the ProLint server, which provides a multitude of visualization tools for illustration of areas of enhanced lipid-protein contact and identification of putative lipid binding sites within the protein matrix. We observed enrichment of CHOL and PIP2 around all proteins with subtle differences in their distribution depending on the protein type and conformation state. Putative binding sites were identified for CHOL, PIP2, POPC, and POSM in the three studied proteins and their potential roles in the SLC4 transport function, conformational transition and protein dimerization were discussed.

Statement of significance

The SLC4 protein family is involved in critical physiological processes like pH and blood pressure regulation and maintenance of ion homeostasis. Its members can be found in various tissues. A number of studies suggest possible lipid regulation of the SLC4 function. However, the protein-lipid interactions in the SLC4 family are still poorly understood. Here we make use of long coarse-grained molecular dynamics simulations to assess the protein-lipid interactions in three SLC4 proteins with different transport modes, AE1, NBCe1, and NDCBE. We identify putative lipid binding sites for several lipid types of potential mechanistic importance, discuss them in the framework of the known experimental data and provide a necessary basis for further studies on lipid regulation of SLC4 function.

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

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

Cholesterol in Class C GPCRs: Role, Relevance, and Localization

G-protein coupled receptors (GPCRs), one of the largest superfamilies of cell-surface receptors, are heptahelical integral membrane proteins that play critical roles in virtually every organ system. G-protein-coupled receptors operate in membranes rich in cholesterol, with an imbalance in cholesterol level within the vicinity of GPCR transmembrane domains affecting the structure and/or function of many GPCRs, a phenomenon that has been linked to several diseases. These effects of cholesterol could result in indirect changes by altering the mechanical properties of the lipid environment or direct changes by binding to specific sites on the protein. There are a number of studies and reviews on how cholesterol modulates class A GPCRs; however, this area of study is yet to be explored for class C GPCRs, which are characterized by a large extracellular region and often form constitutive dimers. This review highlights specific sites of interaction, functions, and structural dynamics involved in the cholesterol recognition of the class C GPCRs. We summarize recent data from some typical family members to explain the effects of membrane cholesterol on the structural features and functions of class C GPCRs and speculate on their corresponding therapeutic potential.

Sphingosine-1-phosphate receptor 2 plays a dual role depending on the stage of cell differentiation in renal epithelial cells

Epithelial renal cells have the ability to adopt different cellular phenotypes through epithelial-mesenchymal transition (EMT) and mesenchymal-epithelial transition (MET). These processes are increasingly recognized as important repair factors following acute renal tubular injury. Sphingosine-1-phosphate (S1P) is a bioactive sphingolipid with impact on proliferation, growth, migration, and differentiation which has significant implication in various diseases including cancer and kidney fibrosis. Here we demonstrated that S1P can exert by activating S1P receptor 2 (S1PR2) different functions depending on the stage of cell differentiation. We observed that the differences in the migratory profile of Madin-Darby canine kidney (MDCK) cells depend both on their stage of cell differentiation and the activity of S1PR2, a receptor that can either promote or inhibit the migratory process. Meanwhile in non-differentiated cells S1PR2 activation avoids migration, it is essential on fully differentiated cells. This is the first time that an antagonist effect of S1PR2 was reported for the same cell type. Moreover, in fully differentiated cells, S1PR2 activation is crucial for the progression of EMT - characterized by adherent junctions disassembly, β-catenin and SNAI2 nuclear translocation and vimentin expression- and depends on ERK 1/2 activation and nuclear translocation. These findings provide a new perspective about the different S1PR2 functions depending on the stage of cell differentiation that can be critical to the modulation of renal epithelial cell plasticity, potentially paving the way for innovative research with pathophysiologic relevance.

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.

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.

Lipid-mediated prestin organization in outer hair cell membranes and its implications in sound amplification

Prestin is a high-density motor protein in the outer hair cells (OHCs), whose conformational response to acoustic signals alters the shape of the cell, thereby playing a major role in sound amplification by the cochlea. Despite recent structures, prestin's intimate interactions with the membrane, which are central to its function remained unresolved. Here, employing a large set (collectively, more than 0.5 ms) of coarse-grained molecular dynamics simulations, we demonstrate the impact of prestin's lipid-protein interactions on its organization at densities relevant to the OHCs and its effectiveness in reshaping OHCs. Prestin causes anisotropic membrane deformation, which mediates a preferential membrane organization of prestin where deformation patterns by neighboring copies are aligned constructively. The resulting reduced membrane rigidity is hypothesized to maximize the impact of prestin on OHC reshaping. These results demonstrate a clear case of protein-protein cooperative communication in membrane, purely mediated by interactions with lipids.

Intermolecular Interactions in G Protein-Coupled Receptor Allosteric Sites at the Membrane Interface from Molecular Dynamics Simulations and Quantum Chemical Calculations

Allosteric modulators are called promising candidates in G protein-coupled receptor (GPCR) drug development by displaying subtype selectivity and more specific receptor modulation. Among the allosteric sites known to date, cavities at the receptor-lipid interface represent an uncharacteristic binding location that raises many questions about the ligand interactions and stability, the binding site structure, and how all of these are affected by lipid molecules. In this work, we analyze interactions in the allosteric sites of the PAR2, C5aR1, and GCGR receptors in three lipid compositions using molecular dynamics simulations. In addition, we performed quantum chemical calculations involving the symmetry-adapted perturbation theory (SAPT) and the natural population analysis to quantify the strength of intermolecular interactions. We show that besides classical hydrogen bonds, weak polar interactions such as O-HC, O-Br, and long-range electrostatics with the backbone amides contribute to the stability of allosteric modulators at the receptor-lipid interface. The allosteric cavities are detectable in various membrane compositions. The availability of polar atoms for interactions in such cavities can be assessed by water molecules from simulations. Although ligand-lipid interactions are weak, lipid tails play a role in ligand binding pose stability and the size of allosteric cavities. We discuss physicochemical aspects of ligand binding at the receptor-lipid interface and suggest a compound library enriched by weak donor groups for ligand search in such sites.

Regulation of PD-L1 through direct binding of cholesterol to CRAC motifs

Cholesterol, an essential molecule for cell structure, function, and viability, plays crucial roles in the development, progression, and survival of cancer cells. Earlier studies have shown that cholesterol-lowering drugs can inhibit the high expression of programmed-death ligand 1 (PD-L1) that contributes to immunoevasion in cancer cells. However, the regulatory mechanism of cell surface PD-L1 abundance by cholesterol is still controversial. Here, using nuclear magnetic resonance and biochemical techniques, we demonstrated that cholesterol can directly bind to the transmembrane domain of PD-L1 through two cholesterol-recognition amino acid consensus (CRAC) motifs, forming a sandwich-like architecture and stabilizing PD-L1 to prevent downstream degradation. Mutations at key binding residues prohibit PD-L1-cholesterol interactions, decreasing the cellular abundance of PD-L1. Our results reveal a unique regulatory mechanism that controls the stability of PD-L1 in cancer cells, providing an alternative method to overcome PD-L1-mediated immunoevasion in cancers.

Ptchd1 mediates opioid tolerance via cholesterol-dependent effects on μ-opioid receptor trafficking

Repeated exposure to opioids causes tolerance, which limits their analgesic utility and contributes to overdose and abuse liability. However, the molecular mechanisms underpinning tolerance are not well understood. Here, we used a forward genetic screen in Caenorhabditis elegans for unbiased identification of genes regulating opioid tolerance which revealed a role for PTR-25/Ptchd1. We found that PTR-25/Ptchd1 controls μ-opioid receptor trafficking and that these effects were mediated by the ability of PTR-25/Ptchd1 to control membrane cholesterol content. Electrophysiological studies showed that loss of Ptchd1 in mice reduced opioid-induced desensitization of neurons in several brain regions and the peripheral nervous system. Mice and C. elegans lacking Ptchd1/PTR-25 display similarly augmented responses to opioids. Ptchd1 knockout mice fail to develop analgesic tolerance and have greatly diminished somatic withdrawal. Thus, we propose that Ptchd1 plays an evolutionarily conserved role in protecting the μ-opioid receptor against overstimulation.

In Silico Assessment of the Lipid Fingerprint Signature of ATP2, the Essential P4-ATPase of Malaria Parasites

ATP2, a putative type 4 P-type ATPase, is a phosphatidylinositol-4-phosphate (PI4P)-regulated phospholipid transporter with an interesting potential as an antimalarial drug target due to its conservation across Plasmodium species and its essential role in the life cycle of Plasmodium falciparum. Despite its importance, the exact mechanism of its action and regulation is still not fully understood. In this study we used coarse-grained molecular dynamics (CG-MD) to elucidate the lipid–protein interactions between a heterogeneous lipid membrane containing phosphatidylinositol and Plasmodium chabaudi ATP2 (PcATP2), an ortholog of P. falciparum ATP2. Our study reveals structural information of the lipid fingerprint of ATP2, and provides structural information on the potential phosphatidylinositol allosteric binding site. Moreover, we identified a set of evolutionary conserved residues that may play a key role in the binding and stabilization of lipids in the binding pocket.

Systematic simulation of the interactions of pleckstrin homology domains with membranes

Pleckstrin homology (PH) domains can recruit proteins to membranes by recognition of phosphatidylinositol phosphate (PIP) lipids. Several family members are linked to diseases including cancer. We report the systematic simulation of the interactions of 100 mammalian PH domains with PIP-containing membranes. The observed PIP interaction hotspots recapitulate crystallographic binding sites and reveal a number of insights: (i) The β1 and β2 strands and their connecting loop constitute the primary PIP interaction site but are typically supplemented by interactions at the β3-β4 and β5-β6 loops; (ii) we reveal exceptional cases such as the Exoc8 PH domain; (iii) PH domains adopt different membrane-bound orientations and induce clustering of anionic lipids; and (iv) beyond family-level insights, our dataset sheds new light on individual PH domains, e.g., by providing molecular detail of secondary PIP binding sites. This work provides a global view of PH domain/membrane association involving multivalent association with anionic lipids.

The role of cholesterol binding in the control of cholesterol by the Scap–Insig system

Scap and Insig, two proteins embedded in the membrane of the endoplasmic reticulum (ER), regulate the synthesis of cholesterol in animal cells by forming a dimer in the presence of high concentrations of cholesterol. Cryo-electron microscopic structures for the Scap–Insig dimer show a sterol-binding site at the dimer interface, but none of the structures include cholesterol itself. Here, a molecular docking approach developed to characterise cholesterol binding to the transmembrane (TM) regions of membrane proteins is used to characterise cholesterol binding to sites on the TM surface of the dimer and to the interfacial binding site. Binding of cholesterol is also observed at sites on the extra-membranous luminal domains of Scap, but the properties of these sites suggest that they will be unoccupied in vivo. Comparing the structure of Scap in the dimer with that predicted by AlphaFold for monomeric Scap suggests that dimer formation could result in relocation of TM helix 7 of Scap and of the loop between TM6 and 7, and that this could be the key change on Scap that signals that there is a high concentration of cholesterol in the ER.

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.

Regulation of Cholesterol Binding to the Receptor Patched1 by its interactions With the Ligand Sonic Hedgehog (Shh)

The Hedgehog (Hh) signaling pathway is essential in cell development and regeneration, which is activated by the ligand Sonic hedgehog (Shh). The binding of Shh to its receptor Patched1 (PTCH1) releases the inhibitory effect on the downstream protein Smoothened (SMO), a G-protein-coupled-receptor (GPCR) protein. Cholesterol was supposed to function as a secondary messenger between PTCH1 and SMO. However, the molecular mechanism of this regulation process is still unclear. Therefore, microsecond coarse-grained molecular dynamics simulations were performed to investigate the protein-lipid interactions of the PTCH1 monomer and dimer-Shh complex. It was observed that the binding of cholesterols to the monomer is more stable than that to the dimer-Shh complex. It is regulated by the enrichment of Ganglioside lipids around proteins and the conformation of Y446, a residue in the sterol-sensing domain (SSD). The regulation of Shh on the dynamics of PTCH1 was further analyzed to explore the allosteric communication pathways between the Shh and the SSD. Our study provides structural and dynamic details of an additional perspective on the regulation of Hh signaling pathway through the lipid micro-environments of PTCH1.

The Bak core dimer focuses triacylglycerides in the membrane

Apoptosis, the intrinsic programmed cell death process, is mediated by the Bcl-2 family members Bak and Bax. Activation via formation of symmetric core dimers and oligomerization on the mitochondrial outer membrane (MOM) leads to permeabilization and cell death. Although this process is linked to the MOM, the role of the membrane in facilitating such pores is poorly understood. We recently described Bak core domain dimers, revealing lipid binding sites and an initial role of lipids in oligomerization. Here we describe simulations that identified localized clustering and interaction of triacylglycerides (TAGs) with a minimized Bak dimer construct. Coalescence of TAGs occurred beneath this Bak dimer, mitigating dimer-induced local membrane thinning and curvature in representative coarse-grain MOM and model membrane systems. Furthermore, the effects observed as a result of coarse-grain TAG cluster formation was concentration dependent, scaling from low physiological MOM concentrations to those found in other organelles. We find that increasing the TAG concentration in liposomes mimicking the MOM decreased the ability of activated Bak to permeabilize these liposomes. These results suggest that the presence of TAGs within a Bak-lipid membrane preserves membrane integrity and is associated with reduced membrane stress, suggesting a possible role of TAGs in Bak-mediated apoptosis.

Insights into lipid-protein interactions from computer simulations

Lipid-protein interactions play an important direct role in the function of many membrane proteins. We argue they are key players in membrane structure, modulate membrane proteins in more subtle ways than direct binding, and are important for understanding the mechanism of classes of hydrophobic drugs. By directly comparing membrane proteins from different families in the same, complex lipid mixture, we found a unique lipid environment for every protein. Extending this work, we identified both differences and similarities in the lipid environment of GPCRs, dependent on which family they belong to and in some cases their conformational state, with particular emphasis on the distribution of cholesterol. More recently, we have been studying modes of coupling between protein conformation and local membrane properties using model proteins. In more applied approaches, we have used similar methods to investigate specific hypotheses on interactions of lipid and lipid-like molecules with ion channels. We conclude this perspective with some considerations for future work, including a new more sophisticated coarse-grained force field (Martini 3), an interactive visual exploration framework, and opportunities to improve sampling.

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

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

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.

Modeling lipid-protein interactions for coarse-grained lipid and Cα protein models

Biological membranes that play major roles in diverse functions are composed of numerous lipids and proteins, making them an important target for coarse-grained (CG) molecular dynamics (MD) simulations. Recently, we have developed the CG implicit solvent lipid force field (iSoLF) that has a resolution compatible with the widely used Cα protein representation [D. Ugarte La Torre and S. Takada, J. Chem. Phys. 153, 205101 (2020)]. In this study, we extended it and developed a lipid-protein interaction model that allows the combination of the iSoLF and the Cα protein force field, AICG2+. The hydrophobic-hydrophilic interaction is modeled as a modified Lennard-Jones potential in which parameters were tuned partly to reproduce the experimental transfer free energy and partly based on the free energy profile normal to the membrane surface from previous all-atom MD simulations. Then, the obtained lipid-protein interaction is tested for the configuration and placement of transmembrane proteins, water-soluble proteins, and peripheral proteins, showing good agreement with prior knowledge. The interaction is generally applicable and is implemented in the publicly available software, CafeMol.

Styrene-maleic acid copolymer effects on the function of the GPCR rhodopsin in lipid nanoparticles

Styrene-maleic acid (SMA) copolymers solubilize biological membranes to form lipid nanoparticles (SMALPs) that contain membrane proteins surrounded by native lipids, thus enabling the use of a variety of biophysical techniques for structural and functional studies. The question of whether SMALPs provide a truly natural environment or SMA solubilization affects the functional properties of membrane proteins, however, remains open. We address this question by comparing the photoactivation kinetics of rhodopsin, a G-protein-coupled receptor in the disk membranes of rod cells, in native membrane and SMALPs prepared at different molar ratios between SMA(3:1) and rhodopsin. Time-resolved absorption spectroscopy combined with complex kinetic analysis reveals kinetic and mechanistic differences between the native membrane and SMA-stabilized environment. The results suggest a range of molar ratios for nanoparticles suitable for kinetic studies.

Lipid distributions and transleaflet cholesterol migration near heterogeneous surfaces in asymmetric bilayers

Specific and nonspecific protein-lipid interactions in cell membranes have important roles in an abundance of biological functions. We have used coarse-grained (CG) molecular dynamics (MD) simulations to assess lipid distributions and cholesterol flipping dynamics around surfaces in a model asymmetric plasma membrane containing one of six structurally distinct entities: aquaporin-1 (AQP1), the bacterial β-barrel outer membrane proteins OmpF and OmpX, the KcsA potassium channel, the WALP23 peptide and a carbon nanotube (CNT). Our findings revealed varied lipid partitioning and cholesterol flipping times around the different solutes and putative cholesterol binding sites in AQP1 and KcsA. The results suggest that protein-lipid interactions can be highly variable, and that surface-dependent lipid profiles are effectively manifested in CG simulations with the Martini force field.

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.

The Role of the Membrane in Transporter Folding and Activity

The synthesis, folding, and function of membrane transport proteins are critical factors for defining cellular physiology. Since the stability of these proteins evolved amidst the lipid bilayer, it is no surprise that we are finding that many of these membrane proteins demonstrate coupling of their structure or activity in some way to the membrane. More and more transporter structures are being determined with some information about the surrounding membrane, and computational modeling is providing further molecular details about these solvation structures. Thus, the field is moving towards identifying which molecular mechanisms - lipid interactions, membrane perturbations, differential solvation, and bulk membrane effects - are involved in linking membrane energetics to transporter stability and function. In this review, we present an overview of these mechanisms and the growing evidence that the lipid bilayer is a major determinant of the fold, form, and function of membrane transport proteins in membranes.

In Silico Identification of Cholesterol Binding Motifs in the Chemokine Receptor CCR3

CC motif chemokine receptor 3 (CCR3) is a Class A G protein-coupled receptor (GPCR) mainly responsible for the cellular trafficking of eosinophils. As such, it plays key roles in inflammatory conditions, such as asthma and arthritis, and the metastasis of many deadly forms of cancer. However, little is known about how CCR3 functionally interacts with its bilayer environment. Here, we investigate cholesterol binding sites in silico through Coarse-Grained Molecular Dynamics (MD) and Pylipid analysis using an extensively validated homology model based on the crystal structure of CCR5. These simulations identified several cholesterol binding sites containing Cholesterol Recognition/Interaction Amino Acid Consensus motif (CRAC) and its inversion CARC motifs in CCR3. One such site, a CARC site in TM1, in conjunction with aliphatic residues in TM7, emerged as a candidate for future investigation based on the cholesterol residency time within the binding pocket. This site forms the core of a cholesterol binding site previously observed in computational studies of CCR2 and CCR5. Most importantly, these cholesterol binding sites are conserved in other chemokine receptors and may provide clues to cholesterol regulation mechanisms in this subfamily of Class A GPCRs.

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.

ProLint: a web-based framework for the automated data analysis and visualization of lipid-protein interactions

The functional activity of membrane proteins is carried out in a complex lipid environment. Increasingly, it is becoming clear that lipids are an important player in regulating or generally modulating their activity. A routinely used method to gain insight into this interplay between lipids and proteins are Molecular Dynamics (MD) simulations, since they allow us to study interactions at atomic or near-atomic detail as a function of time. A major bottleneck, however, is analyzing and visualizing lipid-protein interactions, which, in practice, is a time-demanding task. Here, we present ProLint (www.prolint.ca), a webserver that completely automates analysis of MD generated files and visualization of lipid-protein interactions. Analysis is modular allowing users to select their preferred method, and visualization is entirely interactive through custom built applications that enable a detailed qualitative and quantitative exploration of lipid-protein interactions. ProLint also includes a database of published MD results that have been processed through the ProLint workflow and can be visualized by anyone regardless of their level of experience with MD. The automated analysis, feature-rich visualization, database integration, and open-source distribution with an easy to install process, will allow ProLint to become a routine workflow in lipid-protein interaction studies.

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

Interfacial binding sites for cholesterol on GABAA receptors and competition with neurosteroids

γ-Aminobutyric acid type A (GABA_A) receptors in the brain are located in the outer membranes of brain cells where the concentration of cholesterol is high. Of the 25 available high-resolution structures available for GABA_A receptors, none were determined in the presence of cholesterol, but four include resolved molecules of cholesterol hemisuccinate (CHS). Here, a molecular docking procedure is used to sweep the transmembrane (TM) surfaces of the receptors for cholesterol binding sites. Cholesterol docking poses determined in this way match 89% of the resolved CHS when CHS molecules deemed unlikely to represent typical bound cholesterols are excluded. The receptors are pentameric, and their TM surfaces consist of a set of five facets, each including pairs of TM helices from two adjacent subunits. Each facet contains hydrophobic hollows running from one side of the membrane to the other, within which are six potential binding sites for cholesterol, three on each side of the membrane. High-resolution structures of GABA_A receptors with bound neurosteroids show that neurosteroids bind in these cholesterol binding sites, so the binding of neurosteroids and cholesterol will be competitive.

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.

Dynamical Correlations Reveal Allosteric Sites in G Protein-Coupled Receptors

G protein-coupled Receptors (GPCRs) play a central role in many physiological processes and, consequently, constitute important drug targets. In particular, the search for allosteric drugs has recently drawn attention, since they could be more selective and lead to fewer side effects. Accordingly, computational tools have been used to estimate the druggability of allosteric sites in these receptors. In spite of many successful results, the problem is still challenging, particularly the prediction of hydrophobic sites in the interface between the protein and the membrane. In this work, we propose a complementary approach, based on dynamical correlations. Our basic hypothesis was that allosteric sites are strongly coupled to regions of the receptor that undergo important conformational changes upon activation. Therefore, using ensembles of experimental structures, normal mode analysis and molecular dynamics simulations we calculated correlations between internal fluctuations of different sites and a collective variable describing the activation state of the receptor. Then, we ranked the sites based on the strength of their coupling to the collective dynamics. In the β2 adrenergic (β2AR), glucagon (GCGR) and M2 muscarinic receptors, this procedure allowed us to correctly identify known allosteric sites, suggesting it has predictive value. Our results indicate that this dynamics-based approach can be a complementary tool to the existing toolbox to characterize allosteric sites in GPCRs.