A predicted binding site for cholesterol on the GABAA receptor

Forty Years Searching for Neurosteroid Binding Sites on GABAA Receptors

Brain inhibition is a vital process for controlling and sculpting the excitability of the central nervous system in healthy individuals. This level of control is provided over several timescales and involves the neurotransmitter GABA acting at inhibitory synapses to: rapidly inhibit neurons by activating the GABAA receptor; over a slower timescale, to tonically activate extrasynaptic GABAA receptors to provide a low level of background inhibition; and finally, to activate G-protein coupled GABAB receptors to control transmitter release by inhibiting presynaptic Ca2+ channels whilst providing postsynaptic inhibition via K+ channel activation. From this plethora of roles for GABA and its receptors, the GABAA receptor isoform is of major interest due to its dynamic functional plasticity, which in part, is due to being targeted by modulatory brain neurosteroids derived from sex and stress hormones. This family of neurosteroids can, depending on their structure, potentiate, activate and also inhibit the activity of GABAA receptors to affect brain inhibition. This review tracks the methods that have been deployed in probing GABAA receptors, and charts the sterling efforts made by several groups to locate the key neurosteroid binding sites that affect these important receptors. Increasing our knowledge of these binding sites will greatly facilitate our understanding of the physiological roles of neurosteroids and will help to advance their use as novel therapeutics to combat debilitating brain diseases.

GABA and astrocytic cholesterol determine the lipid environment of GABAAR in cultured cortical neurons

The γ-aminobutyric acid (GABA) type A receptor (GABAAR), a GABA activated pentameric chloride channel, mediates fast inhibitory neurotransmission in the brain. The lipid environment is critical for GABAAR function. How lipids regulate the channel in the cell membrane is not fully understood. Here we employed super resolution imaging of lipids to demonstrate that the agonist GABA induces a rapid and reversible membrane translocation of GABAAR to phosphatidylinositol 4,5-bisphosphate (PIP2) clusters in mouse primary cortical neurons. This translocation relies on nanoscopic separation of PIP2 clusters and lipid rafts (cholesterol-dependent ganglioside clusters). In a resting state, the GABAAR associates with lipid rafts and this colocalization is enhanced by uptake of astrocytic secretions. These astrocytic secretions enhance endocytosis and delay desensitization. Our findings suggest intercellular signaling from astrocytes regulates GABAAR location based on lipid uptake in neurons. The findings have implications for treating mood disorders associated with altered neural excitability.

Transcriptional Dysregulation of Cholesterol Synthesis Underlies Hyposensitivity to GABA in the Ventral Tegmental Area During Acute Alcohol Withdrawal

BACKGROUND

The ventral tegmental area (VTA) is a dopaminergic brain area that is critical in the development and maintenance of addiction. During withdrawal from chronic ethanol exposure, the response of VTA neurons to GABA (gamma-aminobutyric acid) is reduced through an epigenetically regulated mechanism. In the current study, a whole-genome transcriptomic approach was used to investigate the underlying molecular mechanism of GABA hyposensitivity in the VTA during withdrawal after chronic ethanol exposure.

METHODS

We performed RNA sequencing of the VTA of Sprague Dawley male rats withdrawn for 24 hours from a chronic ethanol diet as well as sequencing of the VTA of control rats fed the Lieber-DeCarli diet. RNA sequencing data were analyzed using weighted gene coexpression network analysis to identify modules that contained coexpressed genes. Validation was performed with quantitative polymerase chain reaction, gas chromatography-mass spectrometry, and electrophysiological assays.

RESULTS

Pathway and network analysis of weighted gene coexpression network analysis module 1 revealed a significant downregulation of genes associated with the cholesterol synthesis pathway. Consistent with this association, VTA cholesterol levels were significantly decreased during withdrawal. Chromatin immunoprecipitation indicated a decrease in levels of acetylated H3K27 at the transcriptional control regions of these genes. Electrophysiological studies in VTA slices demonstrated that GABA hyposensitivity during withdrawal was normalized by addition of exogenous cholesterol. In addition, inhibition of cholesterol synthesis produced GABA hyposensitivity, which was reversed by adding exogenous cholesterol to VTA slices.

CONCLUSIONS

These results suggest that decreased expression of cholesterol synthesis genes may regulate GABA hyposensitivity of VTA neurons during alcohol withdrawal. Increasing cholesterol levels in the brain may be a novel avenue for therapeutic intervention to reverse detrimental effects of chronic alcohol exposure.

Knockdown of the glutamate-gated chloride channel gene decreases emamectin benzoate susceptibility in the fall armyworm, Spodoptera frugiperda

Emamectin benzoate (EB), a derivative of avermectin, is the primary insecticide used to control the fall armyworm (FAW) in China. However, the specific molecular targets of EB against FAW remain unclear. In this study, we cloned the glutamate-gated chloride channel (GluCl) gene, which is known to be a primary molecular target for avermectin. We first investigated the transcript levels of SfGluCl in FAW and found that the expression level of SfGluCl in the head and nerve cord was significantly higher than that in other tissues. Furthermore, we found that the expression level of SfGluCl was significantly higher in eggs than that in other developmental stages, including larvae, pupae, and adults. Additionally, we identified three variable splice forms of SfGluCl in exons 3 and 9 and found that their splice frequencies remained unaffected by treatment with the LC50 of EB. RNAi mediated knockdown of SfGluCl showed a significant reduction of 42% and 65% after 48 and 72 h of dsRNA feeding, respectively. Importantly, knockdown of SfGluCl sifgnificantly reduced LC50 and LC90 EB treatment induced mortality of FAW larvae by 15% and 44%, respectively, compared to the control group feeding by dsEGFP. In contrast, there were no significant changes in the mortality of FAW larvae treated with the control insecticides chlorantraniliprole and spinetoram. Finally, molecular docking simulations revealed that EB bound to the large amino-terminal extracellular domain of SfGluCl by forming five hydrogen bonds, two alkyl hydrophobic interactions and one salt bridge. These findings strongly suggest that GluCl may serve as one of the molecular targets of EB in FAW, shedding light on the mode of action of this important insecticide.

Structure-function Studies of GABA (A) Receptors and Related computer-aided Studies

The γ-aminobutyric acid type A receptor (GABA (A) receptor) is a membrane protein activated by the neurotransmitter GABA. Structurally, this major inhibitory neurotransmitter receptor in the human central nervous system is a pentamer that can be built from a selection of 19 subunits consisting of α(1,2,3,4,5 or 6), β (1,2 or 3), γ (1,2 or 3), ρ (1,2 or 3), and δ, π, θ, and ε. This creates several possible pentameric arrangements, which also influence the pharmacological and physiological properties of the receptor. The complexity and heterogeneity of the receptors are further increased by the addition of short and long splice variants in several subunits and the existence of multiple allosteric binding sites and expansive ligands that can bind to the receptors. Therefore, a comprehensive understanding of the structure and function of the receptors is required to gain novel insights into the consequences of receptor dysfunction and subsequent drug development studies. Notably, advancements in computational-aided studies have facilitated the elucidation of residual interactions and exploring energy binding, which may otherwise be challenging to investigate. In this review, we aim to summarize the current understanding of the structure and function of GABA (A) receptors obtained from advancements in computational-aided applications.

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

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

Fluorescence microscopy imaging of a neurotransmitter receptor and its cell membrane lipid milieu

Hampered by the diffraction phenomenon, as expressed in 1873 by Abbe, applications of optical microscopy to image biological structures were for a long time limited to resolutions above the ∼200 nm barrier and restricted to the observation of stained specimens. The introduction of fluorescence was a game changer, and since its inception it became the gold standard technique in biological microscopy. The plasma membrane is a tenuous envelope of 4 nm-10 nm in thickness surrounding the cell. Because of its highly versatile spectroscopic properties and availability of suitable instrumentation, fluorescence techniques epitomize the current approach to study this delicate structure and its molecular constituents. The wide spectral range covered by fluorescence, intimately linked to the availability of appropriate intrinsic and extrinsic probes, provides the ability to dissect membrane constituents at the molecular scale in the spatial domain. In addition, the time resolution capabilities of fluorescence methods provide complementary high precision for studying the behavior of membrane molecules in the time domain. This review illustrates the value of various fluorescence techniques to extract information on the topography and motion of plasma membrane receptors. To this end I resort to a paradigmatic membrane-bound neurotransmitter receptor, the nicotinic acetylcholine receptor (nAChR). The structural and dynamic picture emerging from studies of this prototypic pentameric ligand-gated ion channel can be extrapolated not only to other members of this superfamily of ion channels but to other membrane-bound proteins. I also briefly discuss the various emerging techniques in the field of biomembrane labeling with new organic chemistry strategies oriented to applications in fluorescence nanoscopy, the form of fluorescence microscopy that is expanding the depth and scope of interrogation of membrane-associated phenomena.

Positive Allosteric Modulators of Glycine Receptors and Their Potential Use in Pain Therapies

Glycine receptors are ligand-gated ion channels that mediate synaptic inhibition throughout the mammalian spinal cord, brainstem, and higher brain regions. They have recently emerged as promising targets for novel pain therapies due to their ability to produce antinociception by inhibiting nociceptive signals within the dorsal horn of the spinal cord. This has greatly enhanced the interest in developing positive allosteric modulators of glycine receptors. Several pharmaceutical companies and research facilities have attempted to identify new therapeutic leads by conducting large-scale screens of compound libraries, screening new derivatives from natural sources, or synthesizing novel compounds that mimic endogenous compounds with antinociceptive activity. Advances in structural techniques have also led to the publication of multiple high-resolution structures of the receptor, highlighting novel allosteric binding sites and providing additional information for previously identified binding sites. This has greatly enhanced our understanding of the functional properties of glycine receptors and expanded the structure activity relationships of novel pharmacophores. Despite this, glycine receptors are yet to be used as drug targets due to the difficulties in obtaining potent, selective modulators with favorable pharmacokinetic profiles that are devoid of side effects. This review presents a summary of the structural basis for how current compounds cause positive allosteric modulation of glycine receptors and discusses their therapeutic potential as analgesics. SIGNIFICANCE STATEMENT: Chronic pain is a major cause of disability, and in Western societies, this will only increase as the population ages. Despite the high level of prevalence and enormous socioeconomic burden incurred, treatment of chronic pain remains limited as it is often refractory to current analgesics, such as opioids. The National Institute for Drug Abuse has set finding effective, safe, nonaddictive strategies to manage chronic pain as their top priority. Positive allosteric modulators of glycine receptors may provide a therapeutic option.

Cholesterol-induced suppression of Kir2 channels is mediated by decoupling at the inter-subunit interfaces

Cholesterol is a major regulator of multiple types of ion channels. Although there is increasing information about cholesterol binding sites, the molecular mechanisms through which cholesterol binding alters channel function are virtually unknown. In this study, we used a combination of Martini coarse-grained simulations, a network theory-based analysis, and electrophysiology to determine the effect of cholesterol on the dynamic structure of the Kir2.2 channel. We found that increasing membrane cholesterol reduced the likelihood of contact between specific regions of the cytoplasmic and transmembrane domains of the channel, most prominently at the subunit-subunit interfaces of the cytosolic domains. This decrease in contact was mediated by pairwise interactions of specific residues and correlated to the stoichiometry of cholesterol binding events. The predictions of the model were tested by site-directed mutagenesis of two identified residues-V265 and H222-and high throughput electrophysiology.

Druggable Lipid Binding Sites in Pentameric Ligand-Gated Ion Channels and Transient Receptor Potential Channels

Lipids modulate the function of many ion channels, possibly through direct lipid-protein interactions. The recent outpouring of ion channel structures by cryo-EM has revealed many lipid binding sites. Whether these sites mediate lipid modulation of ion channel function is not firmly established in most cases. However, it is intriguing that many of these lipid binding sites are also known sites for other allosteric modulators or drugs, supporting the notion that lipids act as endogenous allosteric modulators through these sites. Here, we review such lipid-drug binding sites, focusing on pentameric ligand-gated ion channels and transient receptor potential channels. Notable examples include sites for phospholipids and sterols that are shared by anesthetics and vanilloids. We discuss some implications of lipid binding at these sites including the possibility that lipids can alter drug potency or that understanding protein-lipid interactions can guide drug design. Structures are only the first step toward understanding the mechanism of lipid modulation at these sites. Looking forward, we identify knowledge gaps in the field and approaches to address them. These include defining the effects of lipids on channel function in reconstituted systems using asymmetric membranes and measuring lipid binding affinities at specific sites using native mass spectrometry, fluorescence binding assays, and computational approaches.

Dysregulation of Neuronal Nicotinic Acetylcholine Receptor-Cholesterol Crosstalk in Autism Spectrum Disorder

Autism spectrum disorder (ASD) is a set of complex neurodevelopmental diseases that include impaired social interaction, delayed and disordered language, repetitive or stereotypic behavior, restricted range of interests, and altered sensory processing. The underlying causes of the core symptoms remain unclear, as are the factors that trigger their onset. Given the complexity and heterogeneity of the clinical phenotypes, a constellation of genetic, epigenetic, environmental, and immunological factors may be involved. The lack of appropriate biomarkers for the evaluation of neurodevelopmental disorders makes it difficult to assess the contribution of early alterations in neurochemical processes and neuroanatomical and neurodevelopmental factors to ASD. Abnormalities in the cholinergic system in various regions of the brain and cerebellum are observed in ASD, and recently altered cholesterol metabolism has been implicated at the initial stages of the disease. Given the multiple effects of the neutral lipid cholesterol on the paradigm rapid ligand-gated ion channel, the nicotinic acetylcholine receptor, we explore in this review the possibility that the dysregulation of nicotinic receptor-cholesterol crosstalk plays a role in some of the neurological alterations observed in ASD.

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.

Spontaneous lipid binding to the nicotinic acetylcholine receptor in a native membrane

The nicotinic acetylcholine receptor (nAChR) and other pentameric ligand-gated ion channels are native to neuronal membranes with an unusual lipid composition. While it is well-established that these receptors can be significantly modulated by lipids, the underlying mechanisms have been primarily studied in model membranes with few lipid species. Here, we use coarse-grained molecular dynamics simulation to probe specific binding of lipids in a complex quasi-neuronal membrane. We ran a total of 50 μs of simulations of a single nAChR in a membrane composed of 36 species of lipids. Competition between multiple lipid species produces a complex distribution. We find that overall, cholesterol selects for concave inter-subunit sites and polyunsaturated fatty acids select for convex M4 sites, while monounsaturated and saturated lipids are unenriched in the nAChR boundary. We propose the "density-threshold affinity" as a metric calculated from continuous density distributions, which reduces to a standard affinity in two-state binding. We find that the density-threshold affinity for M4 weakens with chain rigidity, which suggests that flexible chains may help relax packing defects caused by the conical protein shape. For any site, PE headgroups have the strongest affinity of all phospholipid headgroups, but anionic lipids still yield moderately high affinities for the M4 sites as expected. We observe cooperative effects between anionic headgroups and saturated chains at the M4 site in the inner leaflet. We also analyze affinities for individual anionic headgroups. When combined, these insights may reconcile several apparently contradictory experiments on the role of anionic phospholipids in modulating nAChR.

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

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

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

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

Tracking Membrane Protein Dynamics in Real Time

Abstract

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

Graphic Abstract

Cholesterol content in the membrane promotes key lipid-protein interactions in a pentameric serotonin-gated ion channel

Pentameric ligand-gated ion channels (pLGICs), embedded in the lipid membranes of nerve cells, mediate fast synaptic transmission and are major pharmaceutical targets. Because of their complexity and the limited knowledge of their structure, their working mechanisms have still to be fully unraveled at the molecular level. Over the past few years, evidence that the lipid membrane may modulate the function of membrane proteins, including pLGICs, has emerged. Here, we investigate, by means of molecular dynamics simulations, the behavior of the lipid membrane at the interface with the 5-HT_3A receptor (5-HT_3AR), a representative pLGIC which is the target of nausea-suppressant drugs, in a nonconductive state. Three lipid compositions are studied, spanning different concentrations of the phospholipids, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine, and of cholesterol, hence a range of viscosities. A variety of lipid interactions and persistent binding events to different parts of the receptor are revealed in the investigated models, providing snapshots of the dynamical environment at the membrane-receptor interface. Some of these events result in lipid intercalation within the transmembrane domain, and others reach out to protein key sections for signal transmission and receptor activation, such as the Cys-loop and the M2-M3 loop. In particular, phospholipids, with their long hydrophobic tails, play an important role in these interactions, potentially providing a bridge between these two structures. A higher cholesterol content appears to promote lipid persistent binding to the receptor.

Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development

Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.

Thermophoretic analysis of ligand-specific conformational states of the inhibitory glycine receptor embedded in copolymer nanodiscs

The glycine receptor (GlyR), a member of the pentameric ligand-gated ion channel family (pLGIC), displays remarkable variations in the affinity and efficacy of the full agonist glycine and the partial agonist taurine depending on the cell system used. Despite detailed insights in the GlyR three-dimensional structure and activation mechanism, little is known about conformational rearrangements induced by these agonists. Here, we characterized the conformational states of the α1 GlyR upon binding of glycine and taurine by microscale thermophoresis expressed in HEK293 cells and Xenopus oocytes after solubilization in amphipathic styrene-maleic acid copolymer nanodiscs. Our results show that glycine and taurine induce different conformational transitions of the GlyR upon ligand binding. In contrast, the variability of agonist affinity is not mediated by an altered conformational change. Thus, our data shed light on specific agonist induced conformational features and mechanisms of pLGIC upon ligand binding determining receptor activation in native environments.

Identification of N-acyl amino acids that are positive allosteric modulators of glycine receptors

Glycine receptors (GlyRs) mediate inhibitory neurotransmission within the spinal cord and play a crucial role in nociceptive signalling. This makes them primary targets for the development of novel chronic pain therapies. Endogenous lipids have previously been shown to modulate glycine receptors and produce analgesia in pain models, however little is known about what chemical features mediate these effects. In this study, we characterised lipid modulation of GlyRs by screening a library of N-acyl amino acids across all receptor subtypes and determined chemical features crucial for their activity. Acyl-glycine's with a C18 carbon tail were found to produce the greatest potentiation, and require a cis double bond within the central region of the carbon tail (ω6 - ω9) to be active. At 1 µM, C18 ω6,9 glycine potentiated glycine induced currents in α₃ and α₃β receptors by over 50%, and α₁, α₂, α₁β and α₂β receptors by over 100%. C18 ω9 glycine (N-oleoyl glycine) significantly enhance glycine induced peak currents and cause a dose-dependent shift in the glycine concentration response. In the presence of 3 µM C18 ω9 glycine, the EC_5o of glycine at the α₁ receptor was reduced from 17 µM to 10 µM. This study has identified several acyl-amino acids which are positive allosteric modulators of GlyRs and make promising lead compounds for the development of novel chronic pain therapies.

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

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

Functional marriage in plasma membrane: Critical cholesterol level-optimal protein activity

In physiology, homeostasis refers to the condition where a system exhibits an optimum functional level. In contrast, any variation from this optimum is considered as a dysfunctional or pathological state. In this review, we address the proposal that a critical cholesterol level in the plasma membrane is required for the proper functioning of transmembrane proteins. Thus, membrane cholesterol depletion or enrichment produces a loss or gain of direct cholesterol-protein interaction and/or changes in the physical properties of the plasma membrane, which affect the basal or optimum activity of transmembrane proteins. Whether or not this functional switching is a generalized mechanism exhibited for all transmembrane proteins, or if it works just for an exclusive group of them, is an open question and an attractive subject to explore at a basic, pharmacological and clinical level.

Microsecond-timescale simulations suggest 5-HT-mediated preactivation of the 5-HT3A serotonin receptor

Aided by efforts to improve their speed and efficiency, molecular dynamics (MD) simulations provide an increasingly powerful tool to study the structure-function relationship of pentameric ligand-gated ion channels (pLGICs). However, accurate reporting of the channel state and observation of allosteric regulation by agonist binding with MD remains difficult due to the timescales necessary to equilibrate pLGICs from their artificial and crystalized conformation to a more native, membrane-bound conformation in silico. Here, we perform multiple all-atom MD simulations of the homomeric 5-hydroxytryptamine 3A (5-HT3A) serotonin receptor for 15 to 20 μs to demonstrate that such timescales are critical to observe the equilibration of a pLGIC from its crystalized conformation to a membrane-bound conformation. These timescales, which are an order of magnitude longer than any previous simulation of 5-HT3A, allow us to observe the dynamic binding and unbinding of 5-hydroxytryptamine (5-HT) (i.e., serotonin) to the binding pocket located on the extracellular domain (ECD) and allosteric regulation of the transmembrane domain (TMD) from synergistic 5-HT binding. While these timescales are not long enough to observe complete activation of 5-HT3A, the allosteric regulation of ion gating elements by 5-HT binding is indicative of a preactive state, which provides insight into molecular mechanisms that regulate channel activation from a resting state. This mechanistic insight, enabled by microsecond-timescale MD simulations, will allow a careful examination of the regulation of pLGICs at a molecular level, expanding our understanding of their function and elucidating key structural motifs that can be targeted for therapeutic regulation.

Modulation of Transient Receptor Potential C Channel Activity by Cholesterol

Changes of cholesterol level in the plasma membrane of cells have been shown to modulate ion channel function. The proposed mechanisms underlying these modulations include association of cholesterol to a single binding site at a single channel conformation, association to a highly flexible cholesterol binding site adopting multiple poses, and perturbation of lipid rafts. These perturbations have been shown to induce reversible targeting of mammalian transient receptor potential C (TRPC) channels to the cholesterol-rich membrane environment of lipid rafts. Thus, the observed inhibition of TRPC channels by methyl-β-cyclodextrin (MβCD), which induces cholesterol efflux from the plasma membrane, may result from disruption of lipid rafts. This perturbation was also shown to disrupt multimolecular signaling complexes containing TRPC channels. The Drosophila TRP and TRP-like (TRPL) channels belong to the TRPC channel subfamily. When the Drosophila TRPL channel was expressed in S2 or HEK293 cells and perfused with MβCD, the TRPL current was abolished in less than 100 s, fitting well the fast kinetic phase of cholesterol sequestration experiments in cells. It was thus suggested that the fast kinetics of TRPL channel suppression by MβCD arise from disruption of lipid rafts. Accordingly, lipid raft perturbation by cholesterol sequestration could give clues to the function of lipid environment in TRPC channel activity and its mechanism.

Untangling Direct and Domain-Mediated Interactions Between Nicotinic Acetylcholine Receptors in DHA-Rich Membranes

At the neuromuscular junction (NMJ), the nicotinic acetylcholine receptor (nAChR) self-associates to give rise to rapid muscle movement. While lipid domains have maintained nAChR aggregates in vitro, their specific roles in nAChR clustering are currently unknown. In the present study, we carried out coarse-grained molecular dynamics simulations (CG-MD) of 1-4 nAChR molecules in two membrane environments: one mixture containing domain-forming, homoacidic lipids, and a second mixture consisting of heteroacidic lipids. Spontaneous dimerization of nAChRs was up to ten times more likely in domain-forming membranes; however, the effect was not significant in four-protein systems, suggesting that lipid domains are less critical to nAChR oligomerization when protein concentration is higher. With regard to lipid preferences, nAChRs consistently partitioned into liquid-disordered domains occupied by the omega-3 ([Formula: see text]-3) fatty acid, docosahexaenoic acid (DHA); enrichment of DHA boundary lipids increased with protein concentration, particularly in homoacidic membranes. This result suggests dimer formation blocks access of saturated chains and cholesterol, but not polyunsaturated chains, to boundary lipid sites.

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.

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

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

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

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

Involvement of the GABAA receptor α subunit in the mode of action of etifoxine

Etifoxine (EFX) is a non-benzodiazepine psychoactive drug which exhibits anxiolytic effects through a dual mechanism, by directly binding to GABA_A receptors (GABA_ARs) and to the mitochondrial 18-kDa translocator protein, resulting in the potentiation of the GABAergic function. The β subunit subtype plays a key role in the EFX-GABA_AR interaction, however this does not explain the anxiolytic effects of this drug. Here, we combined behavioral and electrophysiological experiments to challenge the role of the GABA_AR α subunit in the EFX mode of action. After single administrations of anxiolytic doses (25-50 mg/kg, intraperitoneal), EFX did not induce any neurological nor locomotor impairments, unlike the benzodiazepine bromazepam (0.5-1 mg/kg, intraperitoneal). We established the EFX pharmacological profile on heteropentameric GABA_ARs constructed with α1 to α6 subunit expressed in Xenopus oocyte. Unlike what is known for benzodiazepines, neither the γ nor δ subunits influenced EFX-mediated potentiation of GABA-evoked currents. EFX acted first as a partial agonist on α2β3γ2S, α3β3γ2S, α6β3γ2S and α6β3δ GABA_ARs, but not on α1β3γ2S, α4β3γ2S, α4β3δ nor α5β3γ2S GABA_ARs. Moreover, EFX exhibited much higher positive allosteric modulation towards α2β3γ2S, α3β3γ2S and α6β3γ2S than for α1β3γ2S, α4β3γ2S and α5β3γ2S GABA_ARs. At 20 μM, corresponding to brain concentration at anxiolytic doses, EFX increased GABA potency to the highest extent for α3β3γ2S GABA_ARs. We built a docking model of EFX on α3β3γ2S GABA_ARs, which is consistent with a binding site located between α and β subunits in the extracellular domain. In conclusion, EFX preferentially potentiates α2β3γ2S and α3β3γ2S GABA_ARs, which might support its advantageous anxiolytic/sedative balance.

The nicotinic acetylcholine receptor: a typical 'allosteric machine'

The concept of allosteric interaction was initially proposed to account for the inhibitory feedback mechanism mediated by bacterial regulatory enzymes. In contrast with the classical mechanism of competitive, steric, interaction between ligands for a common site, allosteric interactions take place between topographically distinct sites and are mediated by a discrete and reversible conformational change of the protein. The concept was soon extended to membrane receptors for neurotransmitters and shown to apply to the signal transduction process which, in the case of the acetylcholine nicotinic receptor (nAChR), links the ACh binding site to the ion channel. Pharmacological effectors, referred to as allosteric modulators, such as Ca²⁺ ions and ivermectin, were discovered that enhance the transduction process when they bind to sites distinct from the orthosteric ACh site and the ion channel. The recent X-ray and electron microscopy structures, at atomic resolution, of the resting and active conformations of several homologues of the nAChR, in combination with atomistic molecular dynamics simulations reveal a stepwise quaternary transition in the transduction process with tertiary changes modifying the boundaries between subunits. These interfaces host orthosteric and allosteric modulatory sites which structural organization changes in the course of the transition. The nAChR appears as a typical allosteric machine. The model emerging from these studies has led to the conception and development of several new pharmacological agents.This article is part of a discussion meeting issue 'Allostery and molecular machines'.

GABAA receptor family: overview on structural characterization

The pentameric γ-aminobutyric acid type A receptors are ion channels activated by ligands, which intervene in the rapid inhibitory transmission in the mammalian CNS. Due to their rich pharmacology and therapeutic potential, it is essential to understand their structure and function thoroughly. This deep characterization was hampered by the lack of experimental structural information for many years. Thus, computational techniques have been extensively combined with experimental data, in order to undertake the study of γ-aminobutyric acid type A receptors and their interaction with drugs. Here, we review the exciting journey made to assess the structures of these receptors and outline major outcomes. Finally, we discuss the brand new structure of the α1β2γ2 subtype and the amazing advances it brings to the field.

Boundary lipids of the nicotinic acetylcholine receptor: Spontaneous partitioning via coarse-grained molecular dynamics simulation

Reconstituted nicotinic acetylcholine receptors (nAChRs) exhibit significant gain-of-function upon addition of cholesterol to reconstitution mixtures, and cholesterol affects the organization of nAChRs within domain-forming membranes, but whether nAChR partitions to cholesterol-rich liquid-ordered ("raft" or l_o) domains or cholesterol-poor liquid-disordered (l_do) domains is unknown. We use coarse-grained molecular dynamics simulations to observe spontaneous interactions of cholesterol, saturated lipids, and polyunsaturated (PUFA) lipids with nAChRs. In binary Dipalmitoylphosphatidylcholine:Cholesterol (DPPC:CHOL) mixtures, both CHOL and DPPC acyl chains were observed spontaneously entering deep "non-annular" cavities in the nAChR TMD, particularly at the subunit interface and the β subunit center, facilitated by the low amino acid density in the cryo-EM structure of nAChR in a native membrane. Cholesterol was highly enriched in the annulus around the TMD, but this effect extended over (at most) 5-10 Å. In domain-forming ternary mixtures containing PUFAs, the presence of a single receptor did not significantly affect the likelihood of domain formation. nAChR partitioned to any cholesterol-poor l_do domain that was present, regardless of whether the l_do or l_o domain lipids had PC or PE headgroups. Enrichment of PUFAs among boundary lipids was positively correlated with their propensity for demixing from cholesterol-rich phases. Long n-3 chains (tested here with Docosahexaenoic Acid, DHA) were highly enriched in annular and non-annular embedded sites, partially displacing cholesterol and completely displacing DPPC, and occupying sites even deeper within the bundle. Shorter n-6 chains were far less effective at displacing cholesterol from non-annular sites.

Common binding sites for cholesterol and neurosteroids on a pentameric ligand-gated ion channel

Cholesterol is an essential component of cell membranes, and is required for mammalian pentameric ligand-gated ion channel (pLGIC) function. Computational studies suggest direct interactions between cholesterol and pLGICs but experimental evidence identifying specific binding sites is limited. In this study, we mapped cholesterol binding to Gloeobacter ligand-gated ion channel (GLIC), a model pLGIC chosen for its high level of expression, existing crystal structure, and previous use as a prototypic pLGIC. Using two cholesterol analogue photolabeling reagents with the photoreactive moiety on opposite ends of the sterol, we identified two cholesterol binding sites: an intersubunit site between TM3 and TM1 of adjacent subunits and an intrasubunit site between TM1 and TM4. In both the inter- and intrasubunit sites, cholesterol is oriented such that the 3‑OH group points toward the center of the transmembrane domains rather than toward either the cytosolic or extracellular surfaces. We then compared this binding to that of the cholesterol metabolite, allopregnanolone, a neurosteroid that allosterically modulates pLGICs. The same binding pockets were identified for allopregnanolone and cholesterol, but the binding orientation of the two ligands was markedly different, with the 3‑OH group of allopregnanolone pointing to the intra- and extracellular termini of the transmembrane domains rather than to their centers. We also found that cholesterol increases, whereas allopregnanolone decreases the thermal stability of GLIC. These data indicate that cholesterol and neurosteroids bind to common hydrophobic pockets in the model pLGIC, GLIC, but that their effects depend on the orientation and specific molecular interactions unique to each sterol.

Molecular Dynamics Simulations of Kir2.2 Interactions with an Ensemble of Cholesterol Molecules

Cholesterol is a major regulator of multiple types of ion channels, but the specific mechanisms and the dynamics of its interactions with the channels are not well understood. Kir2 channels were shown to be sensitive to cholesterol through direct interactions with "cholesterol-sensitive" regions on the channel protein. In this work, we used Martini coarse-grained simulations to analyze the long (μs) timescale dynamics of cholesterol with Kir2.2 channels embedded into a model membrane containing POPC phospholipid with 30 mol% cholesterol. This approach allows us to simulate the dynamic, unbiased migration of cholesterol molecules from the lipid membrane environment to the protein surface of Kir2.2 and explore the favorability of cholesterol interactions at both surface sites and recessed pockets of the channel. We found that the cholesterol environment surrounding Kir channels forms a complex milieu of different short- and long-term interactions, with multiple cholesterol molecules concurrently interacting with the channel. Furthermore, utilizing principles from network theory, we identified four discrete cholesterol-binding sites within the previously identified cholesterol-sensitive region that exist depending on the conformational state of the channel-open or closed. We also discovered that a twofold decrease in the cholesterol level of the membrane, which we found earlier to increase Kir2 activity, results in a site-specific decrease of cholesterol occupancy at these sites in both the open and closed states: cholesterol molecules at the deepest of these discrete sites shows no change in occupancy at different cholesterol levels, whereas the remaining sites showed a marked decrease in occupancy.

Challenges and approaches to understand cholesterol-binding impact on membrane protein function: an NMR view

Experimental evidence for a direct role of lipids in determining the structure, dynamics, and function of membrane proteins leads to the term 'functional lipids'. In particular, the sterol molecule cholesterol modulates the activity of many membrane proteins. The precise nature of cholesterol-binding sites and the consequences of modulation of local membrane micro-viscosity by cholesterol, however, is often unknown. Here, we review the current knowledge of the interaction of cholesterol with transmembrane proteins, with a special focus on structural aspects of the interaction derived from nuclear magnetic resonance approaches. We highlight examples of the importance of cholesterol modulation of membrane protein function, discuss the specificity of cholesterol binding, and review the proposed binding motifs from a molecular perspective. We conclude with a short perspective on what could be future trends in research efforts targeted towards a better understanding of cholesterol/membrane protein interactions.

Orthosteric and benzodiazepine cavities of the α1β2γ2 GABAA receptor: insights from experimentally validated in silico methods

γ-aminobutyric acid-type A (GABA_A) receptors mediate fast synaptic inhibition in the central nervous system of mammals. They are modulated via several sites by numerous compounds, which include GABA, benzodiazepines, ethanol, neurosteroids and anaesthetics among others. Due to their potential as targets of novel drugs, a detailed knowledge of their structure-function relationships is needed. Here, we present the model of the α₁β₂γ₂ subtype GABA_A receptor in the APO state and in complex with selected ligands, including agonists, antagonists and allosteric modulators. The model is based on the crystallographic structure of the human β₃ homopentamer GABA_A receptor. The complexes were refined using atomistic molecular dynamics simulations. This allowed a broad description of the binding modes and the detection of important interactions in agreement with experimental information. From the best of our knowledge, this is the only model of the α₁β₂γ₂ GABA_A receptor that represents altogether the desensitized state of the channel and comprehensively describes the interactions of ligands of the orthosteric and benzodiazepines binding sites in agreement with the available experimental data. Furthermore, it is able to explain small differences regarding the binding of a variety of chemically divergent ligands. Finally, this new model may pave the way for the design of focused experimental studies that will allow a deeper description of the receptor.

Identification of binding sites contributing to volatile anesthetic effects on GABA type A receptors

Most general anesthetics enhance GABA type A (GABA_A) receptor activity at clinically relevant concentrations. Sites of action of volatile anesthetics on the GABA_A receptor remain unknown, whereas sites of action of many intravenous anesthetics have been identified in GABA_A receptors by using photolabeling. Here, we used photoactivatable analogs of isoflurane (AziISO) and sevoflurane (AziSEVO) to locate their sites on α₁β₃γ_2L and α₁β₃ GABA_A receptors. As with isoflurane and sevoflurane, AziISO and AziSEVO enhanced the currents elicited by GABA. AziISO and AziSEVO each labeled 10 residues in α₁β₃ receptors and 9 and 8 residues, respectively, in α₁β₃γ_2L receptors. Photolabeled residues were concentrated in transmembrane domains and located in either subunit interfaces or in the interface between the extracellular domain and the transmembrane domain. The majority of these transmembrane residues were protected from photolabeling with the addition of excess parent anesthetic, which indicated specificity. Binding sites were primarily located within α+/β- and β+/α- subunit interfaces, but residues in the α+/γ- interface were also identified, which provided a basis for differential receptor subtype sensitivity. Isoflurane and sevoflurane did not always share binding sites, which suggests an unexpected degree of selectivity.-Woll, K. A., Zhou, X., Bhanu, N. V., Garcia, B. A., Covarrubias, M., Miller, K. W., Eckenhoff, R. G. Identification of binding sites contributing to volatile anesthetic effects on GABA type A receptors.

Crystal structures of a GABAA-receptor chimera reveal new endogenous neurosteroid-binding sites

γ-Aminobutyric acid receptors (GABA_ARs) are vital for controlling excitability in the brain. This is emphasized by the numerous neuropsychiatric disorders that result from receptor dysfunction. A critical component of most native GABA_ARs is the α subunit. Its transmembrane domain is the target for many modulators, including endogenous brain neurosteroids that impact anxiety, stress and depression, and for therapeutic drugs, such as general anesthetics. Understanding the basis for the modulation of GABA_AR function requires high-resolution structures. Here we present the first atomic structures of a GABA_AR chimera at 2.8-Å resolution, including those bound with potentiating and inhibitory neurosteroids. These structures define new allosteric binding sites for these modulators that are associated with the α-subunit transmembrane domain. Our findings will enable the exploitation of neurosteroids for therapeutic drug design to regulate GABA_ARs in neurological disorders.

Allosteric modulation as a unifying mechanism for receptor function and regulation

Four major receptor families enable cells to respond to chemical and physical signals from their proximal environment. The ligand- and voltage-gated ion channels, G-protein-coupled receptors, nuclear hormone receptors and receptor tyrosine kinases are all allosteric proteins that carry multiple, spatially distinct, yet conformationally linked ligand-binding sites. Recent studies point to common mechanisms governing the allosteric transitions of these receptors, including the impact of oligomerization, pre-existing and functionally distinct conformational ensembles, intrinsically disordered regions, and the occurrence of allosteric modulatory sites. Importantly, synthetic allosteric modulators are being discovered for these receptors, providing an enriched, yet challenging, landscape for novel therapeutics.

Regulation of Ca2+-Sensitive K+ Channels by Cholesterol and Bile Acids via Distinct Channel Subunits and Sites

Cholesterol (CLR) conversion into bile acids (BAs) in the liver constitutes the major pathway for CLR elimination from the body. Moreover, these steroids regulate each other's metabolism. While the roles of CLR and BAs in regulating metabolism and tissue function are well known, research of the last two decades revealed the existence of specific protein receptors for CLR or BAs in tissues with minor contribution to lipid metabolism, raising the possibility that these lipids serve as signaling molecules throughout the body. Among other lipids, CLR and BAs regulate ionic current mediated by the activity of voltage- and Ca²⁺-gated, K⁺ channels of large conductance (BK channels) and, thus, modulate cell physiology and participate in tissue pathophysiology. Initial work attributed modification of BK channel function by CLR or BAs to the capability of these steroids to directly interact with bilayer lipids and thus alter the physicochemical properties of the bilayer with eventual modification of BK channel function. Based on our own work and that of others, we now review evidence that supports direct interactions between CLR or BA and specific BK protein subunits, and the consequence of such interactions on channel activity and organ function, with a particular emphasis on arterial smooth muscle. For each steroid type, we will also briefly discuss several mechanisms that may underlie modification of channel steady-state activity. Finally, we will present novel computational data that provide a chemical basis for differential recognition of CLR vs lithocholic acid by distinct BK channel subunits and recognition sites.

Tryptophan and Cysteine Mutations in M1 Helices of α1β3γ2L γ-Aminobutyric Acid Type A Receptors Indicate Distinct Intersubunit Sites for Four Intravenous Anesthetics and One Orphan Site

BACKGROUND

γ-Aminobutyric acid type A (GABAA) receptors mediate important effects of intravenous general anesthetics. Photolabel derivatives of etomidate, propofol, barbiturates, and a neurosteroid get incorporated in GABAA receptor transmembrane helices M1 and M3 adjacent to intersubunit pockets. However, photolabels have not been consistently targeted at heteromeric αβγ receptors and do not form adducts with all contact residues. Complementary approaches may further define anesthetic sites in typical GABAA receptors.

METHODS

Two mutation-based strategies, substituted tryptophan sensitivity and substituted cysteine modification-protection, combined with voltage-clamp electrophysiology in Xenopus oocytes, were used to evaluate interactions between four intravenous anesthetics and six amino acids in M1 helices of α1, β3, and γ2L GABAA receptor subunits: two photolabeled residues, α1M236 and β3M227, and their homologs.

RESULTS

Tryptophan substitutions at α1M236 and positional homologs β3L231 and γ2L246 all caused spontaneous channel gating and reduced γ-aminobutyric acid EC50. Substituted cysteine modification experiments indicated etomidate protection at α1L232C and α1M236C, R-5-allyl-1-methyl-5-(m-trifluoromethyl-diazirinylphenyl) barbituric acid protection at β3M227C and β3L231C, and propofol protection at α1M236C and β3M227C. No alphaxalone protection was evident at the residues the authors explored, and none of the tested anesthetics protected γ2I242C or γ2L246C.

CONCLUSIONS

All five intersubunit transmembrane pockets of GABAA receptors display similar allosteric linkage to ion channel gating. Substituted cysteine modification and protection results were fully concordant with anesthetic photolabeling at α1M236 and β3M227 and revealed overlapping noncongruent sites for etomidate and propofol in β-α interfaces and R-5-allyl-1-methyl-5-(m-trifluoromethyl-diazirinylphenyl) barbituric acid and propofol in α-β and γ-β interfaces. The authors' results identify the α-γ transmembrane interface as a potentially unique orphan modulator site.

Towards understanding the molecular basis of ion channel modulation by lipids: Mechanistic models and current paradigms

Research on ion channel modulation has become a hot topic because of the key roles these membrane proteins play in both prokaryotic and eukaryotic organisms. In this respect, lipid modulation adds to the overall modulatory mechanisms as a potential via to find new pharmacological targets for drug design based on interfering with lipid/channel interactions. However, our knowledge in this field is scarce and often circumscribed to the sites where lipids bind and/or its final functional consequences. To fully understand this process it is necessary to improve our knowledge on its molecular basis, from the binding sites to the signalling pathways that derive in structural and functional effects on the ion channel. In this review, we have compiled information about such mechanisms and established a classification into four different modes of action. Afterwards, we have revised in more detail the lipid modulation of Cys-loop receptors and of the potassium channel KcsA, which were chosen as model channels modulated by specific lipids. This article is part of a Special Issue entitled: Membrane Lipid Therapy: Drugs Targeting Biomembranes edited by Pablo V. Escribá.

Simulation of Ligand Binding to Membrane Proteins

Membrane proteins are involved in a large variety of functions. Most of these protein functions are regulated by ligand binding with diverse modes of action: agonists, partial agonists, antagonists, and allosteric modulators, potentiators and inhibitors. From the pharmacological point of view, membrane proteins are one if not the major target for drug development. However, experimental structure determination of membrane proteins in complex or in free form still represents a great challenge. Molecular dynamics (MD) simulations commonly reach the microsecond scale on membrane systems. This numerical tool is mature enough to predict and add molecular details on the different ligand-binding modes. In the present chapter, I will present the different steps to design, simulate, and analyze a MD simulation system containing a protein embedded in a membrane and surrounded by water and ligand. As an illustration, the simulation of the ligand-gated ion channel γ-aminobutyric acid type A receptor (GABA_AR) surrounded by one of its allosteric potentiators, bromoform, will be presented and discussed.

A Novel Bifunctional Alkylphenol Anesthetic Allows Characterization of γ-Aminobutyric Acid, Type A (GABAA), Receptor Subunit Binding Selectivity in Synaptosomes

Propofol, an intravenous anesthetic, is a positive modulator of the GABAA receptor, but the mechanistic details, including the relevant binding sites and alternative targets, remain disputed. Here we undertook an in-depth study of alkylphenol-based anesthetic binding to synaptic membranes. We designed, synthesized, and characterized a chemically active alkylphenol anesthetic (2-((prop-2-yn-1-yloxy)methyl)-5-(3-(trifluoromethyl)-3H-diazirin-3-yl)phenol, AziPm-click (1)), for affinity-based protein profiling (ABPP) of propofol-binding proteins in their native state within mouse synaptosomes. The ABPP strategy captured ∼4% of the synaptosomal proteome, including the unbiased capture of five α or β GABAA receptor subunits. Lack of γ2 subunit capture was not due to low abundance. Consistent with this, independent molecular dynamics simulations with alchemical free energy perturbation calculations predicted selective propofol binding to interfacial sites, with higher affinities for α/β than γ-containing interfaces. The simulations indicated hydrogen bonding is a key component leading to propofol-selective binding within GABAA receptor subunit interfaces, with stable hydrogen bonds observed between propofol and α/β cavity residues but not γ cavity residues. We confirmed this by introducing a hydrogen bond-null propofol analogue as a protecting ligand for targeted-ABPP and observed a lack of GABAA receptor subunit protection. This investigation demonstrates striking interfacial GABAA receptor subunit selectivity in the native milieu, suggesting that asymmetric occupancy of heteropentameric ion channels by alkylphenol-based anesthetics is sufficient to induce modulation of activity.