Exploring GPCR–Lipid Interactions by Molecular Dynamics Simulations: Excitements, Challenges, and the Way Forward

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.

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

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

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.

Assessing Protein-Protein Docking Protocols: Case Studies of G-Protein-Coupled Receptor Interactions

The interactions of G-protein-coupled receptors (GPCRs) with other proteins are critical in several cellular processes but resolving their structural dynamics remains challenging. An increasing number of GPCR complexes have been experimentally resolved but others including receptor variants are yet to be characterized, necessitating computational predictions of their interactions. Although integrative approaches with multi-scale simulations would provide rigorous estimates of their conformational dynamics, protein-protein docking remains a first tool of choice of many researchers due to the availability of open-source programs and easy to use web servers with reasonable predictive power. Protein-protein docking algorithms have limited ability to consider protein flexibility, environment effects, and entropy contributions and are usually a first step towards more integrative approaches. The two critical steps of docking: the sampling and scoring algorithms have improved considerably and their performance has been validated against experimental data. In this chapter, we provide an overview and generalized protocol of a few docking protocols using GPCRs as test cases. In particular, we demonstrate the interactions of GPCRs with extracellular protein ligands and an intracellular protein effectors (G-protein) predicted from docking approaches and test their limitations. The current chapter will help researchers critically assess docking protocols and predict experimentally testable structures of GPCR complexes.

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.

Dynamic Insights into the Self-Activation Pathway and Allosteric Regulation of the Orphan G-Protein-Coupled Receptor GPR52

Within over 800 members of G-protein-coupled receptors, there are numerous orphan receptors whose endogenous ligands are largely unknown, providing many opportunities for novel drug discovery. However, the lack of an in-depth understanding of the intrinsic working mechanism for orphan receptors severely limits the related rational drug design. The G-protein-coupled receptor 52 (GPR52) is a unique orphan receptor that constitutively increases cellular 5'-cyclic adenosine monophosphate (cAMP) levels without binding any exogenous agonists and has been identified as a promising therapeutic target for central nervous system disorders. Although recent structural biology studies have provided snapshots of both active and inactive states of GPR52, the mechanism of the conformational transition between these states remains unclear. Here, an acceptable self-activation pathway for GPR52 was proposed through 6 μs Gaussian accelerated molecular dynamics (GaMD) simulations, in which the receptor spontaneously transitions from the active state to that matching the inactive crystal structure. According to the three intermediate states of the receptor obtained by constructing a reweighted potential of mean force, how the allosteric regulation occurs between the extracellular orthosteric binding pocket and the intracellular G-protein-binding site is revealed. Combined with the independent gradient model, several important microswitch residues and the allosteric communication pathway that directly links the two regions are both identified. Transfer entropy calculations not only reveal the complex allosteric signaling within GPR52 but also confirm the unique role of ECL2 in allosteric regulation, which is mutually validated with the results of GaMD simulations. Overall, this work elucidates the allosteric mechanism of GPR52 at the atomic level, providing the most detailed information to date on the self-activation of the orphan receptor.

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.

Cofilin-Membrane Interactions: Electrostatic Effects in Phosphoinositide Lipid Binding

The actin cytoskeleton interacts with the cell membrane primarily through the indirect interactions of actin-binding proteins such as cofilin-1. The molecular mechanisms underlying the specific interactions of cofilin-1 with membrane lipids are still unclear. Here, we performed coarse-grain molecular dynamics simulations of cofilin-1 with complex lipid bilayers to analyze the specificity of protein-lipid interactions. We observed the maximal interactions with phosphoinositide (PIP) lipids, especially PIP₂ and PIP₃ lipids. A good match was observed between the residues predicted to interact and previous experimental studies. The clustering of PIP lipids around the membrane bound protein leads to an overall lipid demixing and gives rise to persistent membrane curvature. Further, through a series of control simulations, we observe that both electrostatics and geometry are critical for specificity of lipid binding. Our current study is a step towards understanding the physico-chemical basis of cofilin-PIP lipid interactions.

Non-ionic cholesterol-based additives for the stabilization of membrane proteins

We report herein the synthesis of two non-ionic amphiphiles with a cholesterol hydrophobic moiety that can be used as chemical additives for biochemical studies of membrane proteins. They were designed to show a high similarity with the planar steroid core of cholesterol and small-to-medium polar head groups attached at the C3 position of ring-A on the sterol skeleton. The two Chol-Tris and Chol-DG have a Tris-hydroxymethyl and a branched diglucose polar head group, respectively, which provide them sufficient water solubility when mixed with the "gold standard" detergent n-Dodecyl-β-D-Maltoside (DDM). The colloidal properties of these mixed micelles were investigated by means of surface tension (SFT) measurements and dynamic light scattering (DLS) experiments and showed the formation of globular micelles of about 8 nm in diameter with a critical micellar concentration of 0.20 mM for DDM:Chol-DG and 0.22 mM for DDM:Chol-Tris. We showed that mixed micelles do not alter the extraction potency of a G-protein coupled receptor (GPCR): the human adenosine A2A receptor (A_2AR). The thermostabilizing effect of the mixed micelles was confirmed on two GPCRs, A_2AR and the growth hormone secretagogue receptor (GHSR). Finally, these two mixed micelles were found suitable for the purification of an active form of A_2AR which remained able to bind two ligands of different class i.e. the specific agonist CGS-21680 and the specific inverse agonist ZM-241385. This suggests that Chol-Tris and Chol-DG may be used as a non-ionic alternative to the cholesteryl hemisuccinate (CHS) stabilizing agent.

Smith-Lemli-Opitz syndrome: A pathophysiological manifestation of the Bloch hypothesis

The biosynthesis of cholesterol, an essential component of higher eukaryotic membranes, was worked out by Konrad Bloch (and Feodor Lynen) in the 1960s and they received the Nobel Prize around that time in recognition of their pioneering contributions. An elegant consequence of this was a hypothesis proposed by Konrad Bloch (the Bloch hypothesis) which suggests that each subsequent intermediate in the cholesterol biosynthesis pathway is superior in supporting membrane function in higher eukaryotes relative to its precursor. In this review, we discuss an autosomal recessive metabolic disorder, known as Smith-Lemli-Opitz syndrome (SLOS), associated with a defect in the Kandutsch-Russell pathway of cholesterol biosynthesis that results in accumulation of the immediate precursor of cholesterol in its biosynthetic pathway (7-dehydrocholesterol) and an altered cholesterol to total sterol ratio. Patients suffering from SLOS have several developmental, behavioral and cognitive abnormalities for which no drug is available yet. We characterize SLOS as a manifestation of the Bloch hypothesis and review its molecular etiology and current treatment. We further discuss defective Hedgehog signaling in SLOS and focus on the role of the serotonin_1A receptor, a representative neurotransmitter receptor belonging to the GPCR family, in SLOS. Notably, ligand binding activity and cellular signaling of serotonin_1A receptors are impaired in SLOS-like condition. Importantly, cellular localization and intracellular trafficking of the serotonin_1A receptor (which constitute an important determinant of a GPCR cellular function) are compromised in SLOS. We highlight some of the recent developments and emerging concepts in SLOS pathobiology and suggest that novel therapies based on trafficking defects of target receptors could provide new insight into treatment of SLOS.

Thermodynamic architecture and conformational plasticity of GPCRs

G-protein-coupled receptors (GPCRs) are ubiquitous integral membrane proteins involved in diverse cellular signaling processes. Here, we carry out a large-scale ensemble thermodynamic study of 45 ligand-free GPCRs employing a structure-based statistical mechanical framework. We find that multiple partially structured states co-exist in the GPCR native ensemble, with the TM helices 1, 6 and 7 displaying varied folding status, and shaping the conformational landscape. Strongly coupled residues are anisotropically distributed, accounting for only 13% of the residues, illustrating that a large number of residues are inherently dynamic. Active-state GPCRs are characterized by reduced conformational heterogeneity with altered coupling-patterns distributed throughout the structural scaffold. In silico alanine-scanning mutagenesis reveals that extra- and intra-cellular faces of GPCRs are coupled thermodynamically, highlighting an exquisite structural specialization and the fluid nature of the intramolecular interaction network. The ensemble-based perturbation methodology presented here lays the foundation for understanding allosteric mechanisms and the effects of disease-causing mutations in GCPRs.

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.

Structural dynamics of Smoothened (SMO) in the ciliary membrane and its interaction with membrane lipids

The Smoothened receptor (SMO, a 7 pass transmembrane domain, Class F GPCR family protein) plays a crucial role in the Hedgehog (HH) signaling pathway, which is involved in embryonic development and is implicated in various types of cancer throughout the animal kingdom. In the absence of HH signaling, SMO is inhibited by Patched 1 (PTC1; a 12 pass transmembrane domain protein), which is localized in the primary cilia. HH binding leads to the dislocation of PTC1 from the cilia, thus making way for SMO to localize in the primary cilia, as an essential prerequisite for its activation. We have carried out MARTINI coarse-grained molecular dynamics simulations of SMO in POPC and in ciliary membrane models, respectively, to study the interactions of SMO with cholesterol and other lipid molecules in the ciliary membrane, and to gain molecular-level insights into the role of the primary cilia in shaping the functional dynamics of SMO. We are able to identify the interaction of membrane cholesterols with definite sites and domains within SMO and relate them with known cholesterol-binding sequence and structure motifs. We show that cholesterol interactions with the transmembrane domain TMD, unlike those with the cysteine-rich domain (CRD) and the intracellular domain (ICD), are through residues belonging to known cholesterol-binding motifs. Notably, a few persistent interactions of cholesterol with lower TM cholesterol-binding domains are governed by the presence of multiple cholesterol-binding motifs. These analyses have also helped to identify and define a strict cholesterol consensus motif (CCM), which may well steer cholesterol into the hitherto identified binding sites within the TMD of SMO. We have also reported the interaction of phosphatidylinositol 4-phosphate with the intracellular region of transmembrane (TM) helices (TM1, TM3, TM4, and TM5), intracellular loop1, helix8, and Arg/Lys clusters of the ICD. Structural analysis of SMO domains shows significant changes in the CRD and ICD, during the course of the simulation. Further detailed analysis of the dynamics of the TMD reveals the movements of TM5, TM6, and TM7, linked with the helix8, which are possibly involved in shaping the conformational disposition of the ICD. The movement of these TM helices could possibly be a consequence of interactions involving the extracellular domain and extracellular loops. In addition, our analysis also shows that phosphatidylinositol-4-phosphate (PI4P), along with some ICD cholesterols, are implicated in anchoring SMO in the membrane.

A cholesterol analog stabilizes the human β2-adrenergic receptor nonlinearly with temperature

In cell membranes, G protein-coupled receptors (GPCRs) interact with cholesterol, which modulates their assembly, stability, and conformation. Previous studies have shown how cholesterol modulates the structural properties of GPCRs at ambient temperature. Here, we characterized the mechanical, kinetic, and energetic properties of the human β2-adrenergic receptor (β2AR) in the presence and absence of the cholesterol analog cholesteryl hemisuccinate (CHS) at room temperature (25°C), at physiological temperature (37°C), and at high temperature (42°C). We found that CHS stabilized various structural regions of β2AR differentially, which changed nonlinearly with temperature. Thereby, the strongest effects were observed for structural regions that are important for receptor signaling. Moreover, at 37°C, but not at 25° or 42°C, CHS caused β2AR to increase and stabilize conformational substates to adopt to basal activity. These findings indicate that the nonlinear, temperature-dependent action of CHS in modulating the structural and functional properties of this GPCR is optimized for 37°C.

Effect of Hypoxia on the Function of the Human Serotonin1A Receptor

Cellular hypoxia causes numerous pathophysiological conditions associated with the disruption of oxygen homeostasis. Under oxygen-deficient conditions, cells adapt by controlling the cellular functions to facilitate the judicious use of available oxygen, such as cessation of cell growth and proliferation. In higher eukaryotes, the process of cholesterol biosynthesis is intimately coupled to the availability of oxygen, where the synthesis of one molecule of cholesterol requires 11 molecules of O₂. Cholesterol is an essential component of higher eukaryotic membranes and is crucial for the physiological functions of several membrane proteins and receptors. The serotonin_1A receptor, an important neurotransmitter G protein-coupled receptor associated with cognition and memory, has previously been shown to depend on cholesterol for its signaling and function. In this work, in order to explore the interdependence of oxygen levels, cholesterol biosynthesis, and the function of the serotonin_1A receptor, we developed a cellular hypoxia model to explore the function of the human serotonin_1A receptor heterologously expressed in Chinese hamster ovary cells. We observed cell cycle arrest at G1/S phase and the accumulation of lanosterol in cell membranes under hypoxic conditions, thereby validating our cellular model. Interestingly, we observed a significant reduction in ligand binding and disruption of downstream cAMP signaling of the serotonin_1A receptor under hypoxic conditions. To the best of our knowledge, our results represent the first report linking the function of the serotonin_1A receptor with hypoxia. From a broader perspective, these results contribute to our overall understanding of the molecular basis underlying neurological conditions often associated with hypoxia-induced brain dysfunction.

Molecular Mechanisms Underlying Caveolin-1 Mediated Membrane Curvature

Caveolin-1 is one of the main protein components of caveolae that acts as a mechanosensor at the cell membrane. The interactions of caveolin-1 with membranes have been shown to lead to complex effects such as curvature and the clustering of specific lipids. Here, we review the emerging concepts on the molecular interactions of caveolin-1, with a focus on insights from coarse-grain molecular dynamics simulations. Consensus structural models of caveolin-1 report a helix-turn-helix core motif with flanking domains of higher disorder that could be membrane composition dependent. Caveolin-1 appears to be mainly surface-bound and does not embed very deep in the membrane to which it is bound. The most interesting aspect of caveolin-1 membrane binding is the interplay of cholesterol clustering and membrane curvature. Although cholesterol has been reported to cluster in the vicinity of caveolin-1 by several approaches, simulations show that the clustering is maximal in membrane leaflet opposing the surface-bound caveolin-1. The intrinsic negative curvature of cholesterol appears to stabilize the negative curvature in the opposing leaflet. In fact, the simulations show that blocking cholesterol clustering (through artificial position restraints) blocks membrane curvature, and vice versa. Concomitant with cholesterol clustering is sphingomyelin clustering, again in the opposing leaflet, but in a concentration-dependent manner. The differential stress due to caveolin-1 binding and the inherent asymmetry of the membrane leaflets could be the determinant for membrane curvature and needs to be further probed. The review is an important step to reconcile the molecular level details emerging from simulations with the mesoscopic details provided by state of the art experimental approaches.

CG2AT2: an Enhanced Fragment-Based Approach for Serial Multi-scale Molecular Dynamics Simulations

Coarse-grained molecular dynamics provides a means for simulating the assembly and interactions of macromolecular complexes at a reduced level of representation, thereby allowing both longer timescale and larger sized simulations. Here, we describe an enhanced fragment-based protocol for converting macromolecular complexes from coarse-grained to atomistic resolution, for further refinement and analysis. While the focus is upon systems that comprise an integral membrane protein embedded in a phospholipid bilayer, the technique is also suitable for membrane-anchored and soluble protein/nucleotide complexes. Overall, this provides a method for generating an accurate and well-equilibrated atomic-level description of a macromolecular complex. The approach is evaluated using a diverse test set of 11 system configurations of varying size and complexity. Simulations are assessed in terms of protein stereochemistry, conformational drift, lipid/protein interactions, and lipid dynamics.

Caveolin induced membrane curvature and lipid clustering: two sides of the same coin?

Caveolin-1 (cav-1) is a multi-domain membrane protein that is a key player in cell signaling, endocytosis and mechanoprotection. It is the principle component of cholesterol-rich caveolar domains and has been reported to induce membrane curvature. The molecular mechanisms underlying the interactions of cav-1 with complex membranes, leading to modulation of membrane topology and the formation of cholesterol-rich domains, remain elusive. In this study, we aim to understand the effect of lipid composition by analyzing the interactions of cav-1 with complex membrane bilayers comprised of about sixty lipid types. We have performed a series of coarse-grain molecular dynamics simulations using the Martini force-field with a cav-1 protein construct (residue 82-136) that includes the membrane binding domains and a palmitoyl tail. We observe that cav-1 induces curvature in this complex membrane, though it is restricted to a nanometer length scale. Concurrently, we observe a clustering of cholesterol, sphingolipids and other lipid molecules leading to the formation of nanodomains. Direct microsecond timescale interactions are observed for specific lipids such as cholesterol, phosphatidylserine and phosphatidylethanolamine lipid types. The results indicate that there is an interplay between membrane topology and lipid species. Our work is a step toward understanding how lipid composition and organization regulate the formation of caveolae, in the context of endocytosis and cell signaling.

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.

Cholesterol in GPCR Structures: Prevalence and Relevance

Bound cholesterol molecules are emerging as important hallmarks of GPCR structures. In this commentary, we analyze their statistical prevalence and biological relevance.

Conformational plasticity and dynamic interactions of the N-terminal domain of the chemokine receptor CXCR1

The dynamic interactions between G protein-coupled receptors (GPCRs) and their cognate protein partners are central to several cell signaling pathways. For example, the association of CXC chemokine receptor 1 (CXCR1) with its cognate chemokine, interleukin-8 (IL8 or CXCL8) initiates pathways leading to neutrophil-mediated immune responses. The N-terminal domain of chemokine receptors confers ligand selectivity, but unfortunately the conformational dynamics of this intrinsically disordered region remains unresolved. In this work, we have explored the interaction of CXCR1 with IL8 by microsecond time scale coarse-grain simulations, complemented by atomistic models and NMR chemical shift predictions. We show that the conformational plasticity of the apo-receptor N-terminal domain is restricted upon ligand binding, driving it to an open C-shaped conformation. Importantly, we corroborated the dynamic complex sampled in our simulations against chemical shift perturbations reported by previous NMR studies and show that the trends are similar. Our results indicate that chemical shift perturbation is often not a reporter of residue contacts in such dynamic associations. We believe our results represent a step forward in devising a strategy to understand intrinsically disordered regions in GPCRs and how they acquire functionally important conformational ensembles in dynamic protein-protein interfaces.

Lipid metabolism and lipid signals in aging and longevity

Lipids play crucial roles in regulating aging and longevity. In the past few decades, a series of genetic pathways have been discovered to regulate lifespan in model organisms. Interestingly, many of these regulatory pathways are linked to lipid metabolism and lipid signaling. Lipid metabolic enzymes undergo significant changes during aging and are regulated by different longevity pathways. Lipids also actively modulate lifespan and health span as signaling molecules. In this review, we summarize recent insights into the roles of lipid metabolism and lipid signaling in aging and discuss lipid-related interventions in promoting longevity.

Cholesterol as a modulator of cannabinoid receptor CB2 signaling

Signaling through integral membrane G protein-coupled receptors (GPCRs) is influenced by lipid composition of cell membranes. By using novel high affinity ligands of human cannabinoid receptor CB₂, we demonstrate that cholesterol increases basal activation levels of the receptor and alters the pharmacological categorization of these ligands. Our results revealed that (2-(6-chloro-2-((2,2,3,3-tetramethylcyclopropane-1-carbonyl)imino)benzo[d]thiazol-3(2H)-yl)ethyl acetate ligand (MRI-2646) acts as a partial agonist of CB₂ in membranes devoid of cholesterol and as a neutral antagonist or a partial inverse agonist in cholesterol-containing membranes. The differential effects of a specific ligand on activation of CB₂ in different types of membranes may have implications for screening of drug candidates in a search of modulators of GPCR activity. MD simulation suggests that cholesterol exerts an allosteric effect on the intracellular regions of the receptor that interact with the G-protein complex thereby altering the recruitment of G protein.

Role of Cholesterol in Transmembrane Dimerization of the ErbB2 Growth Factor Receptor

The association of ErbB2 growth factor receptors is critical for cell growth and potentiates tumor proliferation in several cancer types. An important aspect in ErbB2 association is the role of lipids such as cholesterol, especially since their metabolism is often reprogrammed in cancer cells. Here, we have coupled metadynamics with coarse-grain simulations to identify cholesterol effects in the transmembrane dimerization of ErbB2 receptors. Overall, cholesterol interactions are observed with the receptor that directly tunes the association energetics. Several dimer conformations are identified both in the presence and absence of cholesterol, although the dimer regime appears to be more favorable in the presence of cholesterol. We observe an overall modulation of the underlying energy profile and the symmetric active and inactive conformational states are not distinguished in the presence of cholesterol. We show that cholesterol binds to the receptor transmembrane domain at a site (CRAC motif) that overlaps with the dimer interface (SmXXXSm motif). The competition between the transmembrane interactions and cholesterol interactions decides the final conformational landscape. Our work is an important step toward characterizing cholesterol effects in ErbB2 membrane receptor function.

Lipid rafts as platforms for sphingosine 1-phosphate metabolism and signalling

Spontaneous segregation of cholesterol and sphingolipids as a liquid-ordered phase leads to their clustering in selected membrane areas, the lipid rafts. These specialized membrane domains enriched in gangliosides, sphingomyelin, cholesterol and selected proteins involved in signal transduction, organize and determine the function of multiprotein complexes involved in several aspects of signal transduction, thus regulating cell homeostasis. Sphingosine 1-phosphate, an important biologically active mediator, is involved in several signal transduction processes regulating a plethora of cell functions and, not only several of its downstream effectors tend to localize in lipid rafts, some of the enzymes involved in its pathway, of receptors involved in its signalling and its transporters have been often found in these membrane microdomains. Considering this, in this review we address what is currently known regarding the relationship between sphingosine 1-phosphate metabolism and signalling and plasma membrane lipid rafts.

Membrane cholesterol regulates endocytosis and trafficking of the serotonin1A receptor: Insights from acute cholesterol depletion

Endocytosis and intracellular trafficking constitute important regulatory features associated with G protein-coupled receptor (GPCR) function. GPCR endocytosis involves several remodeling events at the plasma membrane orchestrated by a concerted interplay of a large number of proteins and membrane lipids. Although considerable literature exists on the protein framework underlying GPCR endocytosis, the role of membrane lipids in this process remains largely unexplored. In order to explore the role of membrane cholesterol (an essential and important lipid in higher eukaryotes) in GPCR endocytosis, we monitored the effect of acute cholesterol depletion using methyl-β-cyclodextrin (MβCD) on endocytosis and intracellular trafficking of the serotonin_1A receptor, an important neurotransmitter GPCR. Our results show that the serotonin_1A receptor exhibits agonist-induced clathrin-mediated endocytosis with a concentration-dependent inhibition in internalization with increasing concentrations of MβCD, which was restored upon cholesterol replenishment. Interestingly, subsequent to internalization under these conditions, serotonin_1A receptors were re-routed toward lysosomal degradation, instead of endosomal recycling observed under normal conditions, thereby implicating membrane cholesterol in modulation of intracellular trafficking of the receptor. This raises the possibility of a novel cholesterol-dependent role of intracellular sorting proteins in GPCR trafficking. These results differ from our previous observations on the endocytosis of the serotonin_1A receptor upon statin-induced chronic cholesterol depletion, in terms of endocytic pathway. We conclude that analysis of complex cellular trafficking events such as GPCR endocytosis under acute and chronic cholesterol depletion conditions should be carried out with caution due to fundamental differences underlying these processes.

Interplay of cholesterol, membrane bilayers and the AT1R: A cholesterol consensus motif on AT1R is revealed

Hypertension, mediated by the Angiotensin II receptor type 1 (AT1R), is still the major cause of premature death despite the discovery of novel therapeutics, highlighting the importance of an in depth understanding of the drug-AT1R recognition mechanisms coupled with the impact of the membrane environment on the interaction of drugs with AT1R. Herein, we examine the interplay of cholesterol-lipid-candesartan and the AT1R using Molecular Dynamics simulations of a model membrane consisting of 60:40 mol%. DPPC:cholesterol, candesartan and the AT1R, mimicking the physiological cholesterol concentration in sarcolemma membranes. The simulations of the model membrane of 60:40 mol%. DPPC:cholesterol were further validated using DOSY NMR experiments. Interestingly, our results suggest a significant role of cholesterol in the AT1R function imposed through a Cholesterol Consensus Motif (CCM) in the receptor, which could be crucial in the drug binding process. Candesartan diffusion towards AT1R through incorporation into lipid bilayers, appears to be retarded by the presence of cholesterol. However, its direct approach towards AT1R may be facilitated through the mobility induced on the N-terminus by the cholesterol binding on the CCM these novel insights could pave the way towards the development of more potent pharmaceutical agents to combat hypertension more effectively.

Structure, dynamics and lipid interactions of serotonin receptors: excitements and challenges

Serotonin (5-hydroxytryptamine, 5-HT) is an intrinsically fluorescent neurotransmitter found in organisms spanning a wide evolutionary range. Serotonin exerts its diverse actions by binding to distinct cell membrane receptors which are classified into many groups. Serotonin receptors are involved in regulating a diverse array of physiological signaling pathways and belong to the family of either G protein-coupled receptors (GPCRs) or ligand-gated ion channels. Serotonergic signaling appears to play a key role in the generation and modulation of various cognitive and behavioral functions such as sleep, mood, pain, anxiety, depression, aggression, and learning. Serotonin receptors act as drug targets for a number of diseases, particularly neuropsychiatric disorders. The signaling mechanism and efficiency of serotonin receptors depend on their amazing ability to rapidly access multiple conformational states. This conformational plasticity, necessary for the wide variety of functions displayed by serotonin receptors, is regulated by binding to various ligands. In this review, we provide a succinct overview of recent developments in generating and analyzing high-resolution structures of serotonin receptors obtained using crystallography and cryo-electron microscopy. Capturing structures of distinct conformational states is crucial for understanding the mechanism of action of these receptors, which could provide important insight for rational drug design targeting serotonin receptors. We further provide emerging information and insight from studies on interactions of membrane lipids (such as cholesterol) with serotonin receptors. We envision that a judicious combination of analysis of high-resolution structures and receptor-lipid interaction would allow a comprehensive understanding of GPCR structure, function and dynamics, thereby leading to efficient drug discovery.

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.

Alfa2-adrenoblockers attenuate the elevated plasma cholesterol, anxiety levels and restore impaired spatial memory of rats under the chronic noise exposure

Noise is considered one of the environmental hazards that negatively affect health. It can cause damage to the auditory, neurological, hormonal and cardiovascular systems, in addition to impairing psychological and cognitive functions. Considering the significance of vascular disturbances and oxidative stress in the development of the aforementioned negative effects, the purpose of our investigation was to study the level of high density lipoprotein-cholesterol (HDL-Cl), low density lipoprotein-cholesterol (LDL-Cl), and total cholesterol (TCl) in plasma, in addition to the behavioral characteristics of white rats, and the effects of the α₂-adrenoblockers beditin and mesedin to reveal their antiatherogenic effect during noise exposure. The "Open field" and "Y-maze" tests were used in order to evaluate the behavioral states of the rats. Investigations were carried out on albino rats divided into 4 groups. The 1^st group of rats served as a control. The 2^nd, 3^rd and 4^th groups were exposed to 91 dBA of noise; the duration of exposure was 8 h per day for 60 days. The 3^rd group was injected with beditin and the 4^th group with mesedin, both intraperitoneally and repeatedly. According to our results, the chronic exposure to high-volume noise leads to the increase of plasma TCl and LDL-Cl concentrations and the decrease of HDL-Cl levels, resulting in increase of the atherogenic coefficient, which is estimated to be one of the main cardiovascular disease risk factors. The "Open field" and "Y-maze" tests revealed that chronic noise exposure caused disturbances in the behavioral activity, a noise duration-dependent delay in movement and orientation, increased anxiety and deficit in the animals' spatial memory. The administration of α₂-adrenoblockers to the noise-exposed animals had a regulatoryeffects of varying intensities, depending on the medication used and the studied parameters under the conditions of chronic acoustic stress.

Cell Cycle Dependent Modulation of Membrane Dipole Potential and Neurotransmitter Receptor Activity: Role of Membrane Cholesterol

The cell cycle is a sequential multistep process essential for growth and proliferation of cells that make up multicellular organisms. A number of nuclear and cytoplasmic proteins are known to modulate the cell cycle. Yet, the role of lipids, membrane organization, and physical properties in cell cycle progression remains largely elusive. Membrane dipole potential is an important physicochemical property and originates due to the electrostatic potential difference within the membrane because of nonrandom arrangement of amphiphile dipoles and water molecules at the membrane interface. In this work, we explored the modulation of membrane dipole potential in various stages of the cell cycle in CHO-K1 cells. Our results show that membrane dipole potential is highest in the G1 phase relative to S and G2/M phases. This was accompanied by regulation of membrane cholesterol content in the cell cycle. The highest cholesterol content was found in the G1 phase with a considerable reduction in cholesterol in S and G2/M phases. Interestingly, we noted a similarity in the dependence of membrane dipole potential and cholesterol with progress of the cell cycle. In addition, we observed an increase in neutral lipid (which contains esterified cholesterol) content as cells progressed from the G1 to G2/M phase via the S phase of the cell cycle. Importantly, we further observed a cell cycle dependent reduction in ligand binding activity of serotonin_1A receptors expressed in CHO-K1 cells. To the best of our knowledge, these results constitute the first report of cell cycle dependent modulation of membrane dipole potential and activity of a neurotransmitter receptor belonging to the G protein-coupled receptor family. We envision that understanding the basis of cell cycle events from a biophysical perspective would result in a deeper appreciation of the cell cycle and its regulation in relation to cellular function.

Molecular evolution of a collage of cholesterol interaction motifs in transmembrane helix V of the serotonin1A receptor

The human serotonin_1A receptor is a representative member of the superfamily of G protein-coupled receptors (GPCRs) and an important drug target for neurological disorders. Using a combination of biochemical, biophysical and molecular dynamics simulation approaches, we and others have shown that membrane cholesterol modulates the organization, dynamics and function of vertebrate serotonin_1A receptors. Previous studies have shown that the cytoplasmic portion of transmembrane helix V (TM V) and the extramembraneous intracellular loop 3 are critical for G-protein coupling, phosphorylation and desensitization of the receptor. We have recently resolved a collage of putative cholesterol interaction motifs from the amino acid sequence overlapping this region. In this paper, we explore the sequence plasticity of this fragment that may have adapted to altered membrane lipidome, after vertebrates evolved from primordial invertebrates. Since invertebrates have lower levels of membrane cholesterol relative to vertebrates, we compared TM V sequence fragments from invertebrate serotonin₁ receptors with vertebrate orthologs to infer the sequence plasticity in TM V. We report that the average number of cholesterol interaction motifs in TM V for diverse phyla represents an increasing trend that could mirror vertebrate evolution from primordial invertebrates. By statistical modeling, we propose that the collage of cholesterol interaction motifs in TM V of the human serotonin_1A receptor may have evolved from rudimentary collages, reminiscent of primordial invertebrate orthologs. Taken together, we propose that a repertoire of cholesterol-philic nonsynonymous substitutions may have enhanced collage complexity in TM V during vertebrate evolution.

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.

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

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

Cholesterol interaction motifs in G protein-coupled receptors: Slippery hot spots?

G protein-coupled receptors (GPCRs) are cell membrane associated signaling hubs that orchestrate a multitude of cellular functions upon binding to a diverse variety of extracellular ligands. Since GPCRs are integral membrane proteins with seven-transmembrane domain architecture, their function, organization and dynamics are intimately regulated by membrane lipids, such as cholesterol. Cholesterol is an extensively studied lipids in terms of its effects on GPCR structure and function. One of the possible mechanisms underlying modulation of GPCR function by cholesterol is via specific interaction of GPCRs with membrane cholesterol. These interactions of GPCRs with membrane cholesterol are often attributed to structural features of GPCRs that could facilitate their preferential association with cholesterol. In this backdrop, cholesterol interaction motifs represent putative interaction sites on GPCRs that could facilitate cholesterol-sensitive function of these receptors. In this review, we provide an overview of cholesterol interaction motifs found in GPCRs, which have been identified through a combination of crystallography, bioinformatics analysis, and functional studies. In addition, we will highlight, using specific examples, why mere presence of a cholesterol interaction motif at a given site may not directly implicate its role in interaction with membrane cholesterol. We therefore believe that experimental approaches, followed by functional analysis of cholesterol sensitivity of GPCRs, would provide a better understanding of the role played by these motifs in cholesterol-sensitive function. We envision that a comprehensive knowledge of cholesterol interaction sites in GPCRs would allow us to develop a better understanding of GPCR structure-function paradigm, and could be useful in future therapeutics. This article is categorized under: Models of Systems Properties and Processes > Mechanistic Models Analytical and Computational Methods > Computational Methods Laboratory Methods and Technologies > Macromolecular Interactions, Methods.

Statin-Induced Chronic Cholesterol Depletion Switches GPCR Endocytosis and Trafficking: Insights from the Serotonin1A Receptor

Endocytosis is a key regulatory mechanism adopted by G protein-coupled receptors (GPCRs) to modulate downstream signaling responses within a stringent spatiotemporal regime. Although the role of membrane lipids has been extensively studied in the context of the function, organization, and dynamics of GPCRs, their role in receptor endocytosis remains largely unexplored. Cholesterol, the predominant sterol in higher eukaryotes, plays a crucial role in maintaining the structure and organization of cell membranes and is involved in essential cellular processes in health and disease. The serotonin_1A receptor is a representative GPCR involved in neuronal development and in neuropsychiatric disorders such as anxiety and depression. We recently combined quantitative flow cytometric and confocal microscopic approaches to demonstrate that the serotonin_1A receptor undergoes clathrin-mediated endocytosis upon agonist stimulation and subsequently traffics along the endosomal recycling pathway. In this work, we show that statin-induced chronic cholesterol depletion switches the endocytic pathway of the serotonin_1A receptor from clathrin- to caveolin-mediated endocytosis. Interestingly, under these conditions, a significant proportion of endocytosed receptors is rerouted toward lysosomal degradation. To the best of our knowledge, these results constitute one of the first comprehensive reports on the role of membrane cholesterol in GPCR endocytosis and trafficking. These results are significant in our overall understanding of the modulatory effects of membrane lipids on GPCR endocytosis and trafficking and could provide novel insight in developing therapeutic interventions against neuropsychiatric disorders such as depression.

Role of cholesterol-mediated effects in GPCR heterodimers

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

Imaging of intermittent lipid-receptor interactions reflects changes in live cell membranes upon agonist-receptor binding

Protein-lipid interactions in cellular membranes modulate central cellular functions, are often transient in character, but occur too intermittently to be readily observable. We introduce transient state imaging (TRAST), combining sensitive fluorescence detection of fluorophore markers with monitoring of their dark triplet state transitions, allowing imaging of such protein-lipid interactions. We first determined the dark state kinetics of the biomembrane fluorophore 7-nitrobenz-2-oxa-1,3-diazole-4-yl (NBD) in lipid vesicles, and how its triplet state is quenched by spin-labels in the same membranes. We then monitored collisional quenching of NBD-lipid derivatives by spin-labelled stearic acids in live cell plasma membranes, and of NBD-lipid derivatives by spin-labelled G-Protein Coupled Receptors (GPCRs). We could then resolve transient interactions between the GPCRs and different lipids, how these interactions changed upon GPCR activation, thereby demonstrating a widely applicable means to image and characterize transient molecular interactions in live cell membranes in general, not within reach via traditional fluorescence readouts.

Computational Nanoscopy of Tight Junctions at the Blood-Brain Barrier Interface

The selectivity of the blood-brain barrier (BBB) is primarily maintained by tight junctions (TJs), which act as gatekeepers of the paracellular space by blocking blood-borne toxins, drugs, and pathogens from entering the brain. The BBB presents a significant challenge in designing neurotherapeutics, so a comprehensive understanding of the TJ architecture can aid in the design of novel therapeutics. Unraveling the intricacies of TJs with conventional experimental techniques alone is challenging, but recently developed computational tools can provide a valuable molecular-level understanding of TJ architecture. We employed the computational methods toolkit to investigate claudin-5, a highly expressed TJ protein at the BBB interface. Our approach started with the prediction of claudin-5 structure, evaluation of stable dimer conformations and nanoscale assemblies, followed by the impact of lipid environments, and posttranslational modifications on these claudin-5 assemblies. These led to the study of TJ pores and barriers and finally understanding of ion and small molecule transport through the TJs. Some of these in silico, molecular-level findings, will need to be corroborated by future experiments. The resulting understanding can be advantageous towards the eventual goal of drug delivery across the BBB. This review provides key insights gleaned from a series of state-of-the-art nanoscale simulations (or computational nanoscopy studies) performed on the TJ architecture.

Mechanism of μ-Opioid Receptor-Magnesium Interaction and Positive Allosteric Modulation

In the era of opioid abuse epidemics, there is an increased demand for understanding how opioid receptors can be allosterically modulated to guide the development of more effective and safer opioid therapies. Among the modulators of the μ-opioid (MOP) receptor, which is the pharmacological target for the majority of clinically used opioid drugs, are monovalent and divalent cations. Specifically, the monovalent sodium cation (Na⁺) has been known for decades to affect MOP receptor signaling by reducing agonist binding, whereas the divalent magnesium cation (Mg²⁺) has been shown to have the opposite effect, notwithstanding the presence of sodium chloride. Although ultra-high-resolution opioid receptor crystal structures have revealed a specific Na⁺ binding site and molecular dynamics (MD) simulation studies have supported the idea that this monovalent ion reduces agonist binding by stabilizing the receptor inactive state, the putative binding site of Mg²⁺ on the MOP receptor, as well as the molecular determinants responsible for its positive allosteric modulation of the receptor, are unknown. In this work, we carried out tens of microseconds of all-atom MD simulations to investigate the simultaneous binding of Mg²⁺ and Na⁺ cations to inactive and active crystal structures of the MOP receptor embedded in an explicit lipid-water environment and confirmed adequate sampling of Mg²⁺ ion binding with a grand canonical Monte Carlo MD method. Analyses of these simulations shed light on 1) the preferred binding sites of Mg²⁺ on the MOP receptor, 2) details of the competition between Mg²⁺ and Na⁺ cations for specific sites, 3) estimates of binding affinities, and 4) testable hypotheses of the molecular mechanism underlying the positive allosteric modulation of the MOP receptor by the Mg²⁺ cation.

Differential Dynamics Underlying the Gln27Glu Population Variant of the β2-Adrenergic Receptor

The β₂-adrenergic receptor (β₂AR) is a membrane-bound G-protein-coupled receptor and an important drug target for asthma. Clinical studies report that the population variant Gln27Glu is associated with a differential response to common asthma drugs, such as albuterol, isoproterenol and terbutaline. Interestingly, the 27th amino acid is positioned on the N-terminal region that is the most flexible and consequently the least studied part of the receptor. In this study, we probe the molecular origin of the differential drug binding by performing structural modeling and simulations of the wild-type (Gln) and variant (Glu) receptors followed by ensemble docking with the ligands, albuterol, isoproterenol and terbutaline. In line with clinical studies, the ligands were observed to interact preferentially with the Glu variant. Our results indicate that the Glu residue at the 27th position perturbs the network of electrostatic interactions that connects the N-terminal region to the binding site in the wild-type receptor. As a result, the Glu variant is observed to bind better to the three ligands tested in this study. Our study provides a structural basis to explain the variable drug response associated with the 27th position polymorphism in the β₂AR and is a starting step to identify genotype-specific therapeutics.