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

A bupropion modulatory site in the Gloeobacter violaceus ligand-gated ion channel

Bupropion is an atypical antidepressant and smoking cessation drug that causes adverse effects such as insomnia, irritability, and anxiety. Bupropion inhibits dopamine and norepinephrine reuptake transporters and eukaryotic cation-conducting pentameric ligand-gated ion channels, such as nicotinic acetylcholine and serotonin type 3A receptors, at clinically relevant concentrations. Here, we demonstrate that bupropion also inhibits a prokaryotic homolog of pentameric ligand-gated ion channels, the Gloeobacter violaceus ligand-gated ion channel (GLIC). Using the GLIC as a model, we used molecular docking to predict binding sites for bupropion. Bupropion was found to bind to several sites within the transmembrane domain, with the predominant site being localized to the interface between transmembrane segments M1 and M3 of two adjacent subunits. Residues W213, T214, and W217 in the first transmembrane segment, M1, and F267 and I271 in the third transmembrane segment, M3, most frequently reside within a 4 Å distance from bupropion. We then used single amino acid substitutions at these positions and two-electrode voltage-clamp recordings to determine their impact on bupropion inhibitory effects. The substitution T214F alters bupropion potency by shifting the half-maximal inhibitory concentration to a 13-fold higher value compared to wild-type GLIC. Residue T214 is found within a previously identified binding pocket for neurosteroids and lipids in the GLIC. This intersubunit binding pocket is structurally conserved and almost identical to a binding pocket described for neurosteroids in γ-aminobutyric acid type A receptors. Our data thus suggest that the T214 that lines a previously identified lipophilic binding pocket in GLIC and γ-aminobutyric acid type A receptors is also a modulatory site for bupropion interaction with the GLIC.

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

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

Targeting TRP channels: recent advances in structure, ligand binding, and molecular mechanisms

Transient receptor potential (TRP) channels are a large and diverse family of transmembrane ion channels that are widely expressed, have important physiological roles, and are associated with many human diseases. These proteins are actively pursued as promising drug targets, benefitting greatly from advances in structural and mechanistic studies of TRP channels. At the same time, the complex, polymodal activation and regulation of TRP channels have presented formidable challenges. In this short review, we summarize recent progresses toward understanding the structural basis of TRP channel function, as well as potential ligand binding sites that could be targeted for therapeutics. A particular focus is on the current understanding of the molecular mechanisms of TRP channel activation and regulation, where many fundamental questions remain unanswered. We believe that a deeper understanding of the functional mechanisms of TRP channels will be critical and likely transformative toward developing successful therapeutic strategies targeting these exciting proteins. This endeavor will require concerted efforts from computation, structural biology, medicinal chemistry, electrophysiology, pharmacology, drug safety and clinical studies.

Lipid nanodisc scaffold and size alter the structure of a pentameric ligand-gated ion channel

Lipid nanodiscs have become a standard tool for studying membrane proteins, including using single particle cryo-electron microscopy (cryo-EM). We find that reconstituting the pentameric ligand-gated ion channel (pLGIC), Erwinia ligand-gated ion channel (ELIC), in different nanodiscs produces distinct structures by cryo-EM. The effect of the nanodisc on ELIC structure extends to the extracellular domain and agonist binding site. Additionally, molecular dynamic simulations indicate that nanodiscs of different size impact ELIC structure and that the nanodisc scaffold directly interacts with ELIC. These findings suggest that the nanodisc plays a crucial role in determining the structure of pLGICs, and that reconstitution of ion channels in larger nanodiscs may better approximate a lipid membrane environment.

Identifying Membrane Protein-Lipid Interactions with Lipidomic Lipid Exchange-Mass Spectrometry

Lipids can play important roles in modulating membrane protein structure and function. However, it is challenging to identify natural lipids bound to membrane proteins in complex bilayers. Here, we developed lipidomic lipid exchange-mass spectrometry (LX-MS) to study the lipid affinity for membrane proteins on a lipidomic scale. We first mix membrane protein nanodiscs with empty nanodiscs that have no embedded membrane proteins. After allowing lipids to passively exchange between the two populations, we separate the two types of nanodiscs and perform lipidomic analysis on each with liquid chromatography and MS. Enrichment of lipids in the membrane protein nanodiscs reveals the affinity of individual lipids for binding the target membrane protein. We apply this approach to study three membrane proteins. With the Escherichia coli ammonium transporter AmtB and aquaporin AqpZ in nanodiscs with E. coli polar lipid extracts, we detected binding of cardiolipin and phosphatidyl-glycerol lipids to the proteins. With the acetylcholine receptor in nanodiscs with brain polar lipid extracts, we discovered a complex set of lipid interactions that depended on the head group and tail composition. Overall, lipidomic LX-MS provides a detailed understanding of the lipid-binding affinity and thermodynamics for membrane proteins in complex bilayers and provides a unique perspective on the chemical environment surrounding membrane proteins.

Structural insights into opposing actions of neurosteroids on GABAA receptors

γ-Aminobutyric acid type A (GABAA) receptors mediate fast inhibitory signaling in the brain and are targets of numerous drugs and endogenous neurosteroids. A subset of neurosteroids are GABAA receptor positive allosteric modulators; one of these, allopregnanolone, is the only drug approved specifically for treating postpartum depression. There is a consensus emerging from structural, physiological and photolabeling studies as to where positive modulators bind, but how they potentiate GABA activation remains unclear. Other neurosteroids are negative modulators of GABAA receptors, but their binding sites remain debated. Here we present structures of a synaptic GABAA receptor bound to allopregnanolone and two inhibitory sulfated neurosteroids. Allopregnanolone binds at the receptor-bilayer interface, in the consensus potentiator site. In contrast, inhibitory neurosteroids bind in the pore. MD simulations and electrophysiology support a mechanism by which allopregnanolone potentiates channel activity and suggest the dominant mechanism for sulfated neurosteroid inhibition is through pore block.

Screening for bilayer-active and likely cytotoxic molecules reveals bilayer-mediated regulation of cell function

A perennial problem encountered when using small molecules (drugs) to manipulate cell or protein function is to assess whether observed changes in function result from specific interactions with a desired target or from less specific off-target mechanisms. This is important in laboratory research as well as in drug development, where the goal is to identify molecules that are unlikely to be successful therapeutics early in the process, thereby avoiding costly mistakes. We pursued this challenge from the perspective that many bioactive molecules (drugs) are amphiphiles that alter lipid bilayer elastic properties, which may cause indiscriminate changes in membrane protein (and cell) function and, in turn, cytotoxicity. Such drug-induced changes in bilayer properties can be quantified as changes in the monomer↔dimer equilibrium for bilayer-spanning gramicidin channels. Using this approach, we tested whether molecules in the Pathogen Box (a library of 400 drugs and drug-like molecules with confirmed activity against tropical diseases released by Medicines for Malaria Venture to encourage the development of therapies for neglected tropical diseases) are bilayer modifiers. 32% of the molecules in the Pathogen Box were bilayer modifiers, defined as molecules that at 10 µM shifted the monomer↔dimer equilibrium toward the conducting dimers by at least 50%. Correlation analysis of the molecules' reported HepG2 cell cytotoxicity to bilayer-modifying potency, quantified as the shift in the gramicidin monomer↔dimer equilibrium, revealed that molecules producing <25% change in the equilibrium had significantly lower probability of being cytotoxic than molecules producing >50% change. Neither cytotoxicity nor bilayer-modifying potency (quantified as the shift in the gramicidin monomer↔dimer equilibrium) was well predicted by conventional physico-chemical descriptors (hydrophobicity, polar surface area, etc.). We conclude that drug-induced changes in lipid bilayer properties are robust predictors of the likelihood of membrane-mediated off-target effects, including cytotoxicity.

The second PI(3,5)P2 binding site in the S0 helix of KCNQ1 stabilizes PIP2-at the primary PI1 site with potential consequences on intermediate-to-open state transition

The Phosphatidylinositol 3-phosphate 5-kinase Type III PIKfyve is the main source for selectively generated phosphatidylinositol 3,5-bisphosphate (PI(3,5)P₂), a known regulator of membrane protein trafficking. PI(3,5)P₂ facilitates the cardiac KCNQ1/KCNE1 channel plasma membrane abundance and therewith increases the macroscopic current amplitude. Functional-physical interaction of PI(3,5)P₂ with membrane proteins and its structural impact is not sufficiently understood. This study aimed to identify molecular interaction sites and stimulatory mechanisms of the KCNQ1/KCNE1 channel via the PIKfyve-PI(3,5)P₂ axis. Mutational scanning at the intracellular membrane leaflet and nuclear magnetic resonance (NMR) spectroscopy identified two PI(3,5)P₂ binding sites, the known PIP₂ site PS1 and the newly identified N-terminal α-helix S0 as relevant for functional PIKfyve effects. Cd²⁺ coordination to engineered cysteines and molecular modeling suggest that repositioning of S0 stabilizes the channel s open state, an effect strictly dependent on parallel binding of PI(3,5)P₂ to both sites.

Structure and function meet at the nicotinic acetylcholine receptor-lipid interface

The nicotinic acetylcholine receptor (nAChR) is a transmembrane protein that mediates fast intercellular communication in response to the endogenous neurotransmitter acetylcholine. It is the best characterized and archetypal molecule in the superfamily of pentameric ligand-gated ion channels (pLGICs). As a typical transmembrane macromolecule, it interacts extensively with its vicinal lipid microenvironment. Experimental evidence provides a wealth of information on receptor-lipid crosstalk: the nAChR exerts influence on its immediate membrane environment and conversely, the lipid moiety modulates ligand binding, affinity state transitions and gating of ion translocation functions of the receptor protein. Recent cryogenic electron microscopy (cryo-EM) studies have unveiled the occurrence of sites for phospholipids and cholesterol on the lipid-exposed regions of neuronal and electroplax nAChRs, confirming early spectroscopic and affinity labeling studies demonstrating the close contact of lipid molecules with the receptor transmembrane segments. This new data provides structural support to the postulated "lipid sensor" ability displayed by the outer ring of M4 transmembrane domains and their modulatory role on nAChR function, as we postulated a decade ago. Borrowing from the best characterized nAChR, the electroplax (muscle-type) receptor, and exploiting new structural information on the neuronal nAChR, it is now possible to achieve an improved depiction of these sites. In combination with site-directed mutagenesis, single-channel electrophysiology, and molecular dynamics studies, the new structural information delivers a more comprehensive portrayal of these lipid-sensitive loci, providing mechanistic explanations for their ability to modulate nAChR properties and raising the possibility of targetting them in disease.

Membrane phospholipids control gating of the mechanosensitive potassium leak channel TREK1

Tandem pore domain (K2P) potassium channels modulate resting membrane potentials and shape cellular excitability. For the mechanosensitive subfamily of K2Ps, the composition of phospholipids within the bilayer strongly influences channel activity. To examine the molecular details of K2P lipid modulation, we solved cryo-EM structures of the TREK1 K2P channel bound to either the anionic lipid phosphatidic acid (PA) or the zwitterionic lipid phosphatidylethanolamine (PE). At the extracellular face of TREK1, a PA lipid inserts its hydrocarbon tail into a pocket behind the selectivity filter, causing a structural rearrangement that recapitulates mutations and pharmacology known to activate TREK1. At the cytoplasmic face, PA and PE lipids compete to modulate the conformation of the TREK1 TM4 gating helix. Our findings demonstrate two distinct pathways by which anionic lipids enhance TREK1 activity and provide a framework for a model that integrates lipid gating with the effects of other mechanosensitive K2P modulators.

Open-channel structure of a pentameric ligand-gated ion channel reveals a mechanism of leaflet-specific phospholipid modulation

Pentameric ligand-gated ion channels (pLGICs) mediate synaptic transmission and are sensitive to their lipid environment. The mechanism of phospholipid modulation of any pLGIC is not well understood. We demonstrate that the model pLGIC, ELIC (Erwinia ligand-gated ion channel), is positively modulated by the anionic phospholipid, phosphatidylglycerol, from the outer leaflet of the membrane. To explore the mechanism of phosphatidylglycerol modulation, we determine a structure of ELIC in an open-channel conformation. The structure shows a bound phospholipid in an outer leaflet site, and structural changes in the phospholipid binding site unique to the open-channel. In combination with streamlined alchemical free energy perturbation calculations and functional measurements in asymmetric liposomes, the data support a mechanism by which an anionic phospholipid stabilizes the activated, open-channel state of a pLGIC by specific, state-dependent binding to this site.

Activation of a Bacterial Mechanosensitive Channel, MscL, Underlies the Membrane Permeabilization of Dual-Targeting Antibacterial Compounds

Resistance to antibiotics is a serious and worsening threat to human health worldwide, and there is an urgent need to develop new antibiotics that can avert it. One possible solution is the development of compounds that possess multiple modes of action, requiring at least two mutations to acquire resistance. Compound SCH-79797 both avoids resistance and has two mechanisms of action: one inhibiting the folate pathway, and a second described as “membrane permeabilization”; however, the mechanism by which membranes from bacterial cells, but not the host, are disrupted has remained mysterious. The opening of the bacterial mechanosensitive channel of large conductance, MscL, which ordinarily serves the physiological role of osmotic emergency release valves gated by hypoosmotic shock, has been previously demonstrated to affect bacterial membrane permeabilization. MscL allows the rapid permeabilization of both ions and solutes through the opening of the largest known gated pore, which has a diameter of 30 Å. We found that SCH-79797 and IRS-16, a more potent derivative, directly bind to the MscL channel and produce membrane permeabilization as a result of its activation. These findings suggest that possessing or adding an MscL-activating component to an antibiotic compound could help to lower toxicity and evade antibiotic resistance.

Recent Insight into Lipid Binding and Lipid Modulation of Pentameric Ligand-Gated Ion Channels

Pentameric ligand-gated ion channels (pLGICs) play a leading role in synaptic communication, are implicated in a variety of neurological processes, and are important targets for the treatment of neurological and neuromuscular disorders. Endogenous lipids and lipophilic compounds are potent modulators of pLGIC function and may help shape synaptic communication. Increasing structural and biophysical data reveal sites for lipid binding to pLGICs. Here, we update our evolving understanding of pLGIC–lipid interactions highlighting newly identified modes of lipid binding along with the mechanistic understanding derived from the new structural data.

Regulation of Gramicidin Channel Function Solely by Changes in Lipid Intrinsic Curvature

Membrane protein function is regulated by the lipid bilayer composition. In many cases the changes in function correlate with changes in the lipid intrinsic curvature (c 0), and c 0 is considered a determinant of protein function. Yet, water-soluble amphiphiles that cause either negative or positive changes in curvature have similar effects on membrane protein function, showing that changes in lipid bilayer properties other than c 0 are important-and may be dominant. To further investigate the mechanisms underlying the bilayer regulation of protein function, we examined how maneuvers that alter phospholipid head groups effective "size"-and thereby c 0-alter gramicidin (gA) channel function. Using dioleoylphospholipids and planar bilayers, we varied the head groups' physical volume and the electrostatic repulsion among head groups (and thus their effective size). When 1,2-dioleyol-sn-glycero-3-phosphocholine (DOPC), was replaced by 1,2-dioleyol-sn-glycero-3-phosphoethanolamine (DOPE) with a smaller head group (causing a more negative c 0), the channel lifetime (τ) is decreased. When the pH of the solution bathing a 1,2-dioleyol-sn-glycero-3-phosphoserine (DOPS) bilayer is decreased from 7 to 3 (causing decreased head group repulsion and a more negative c 0), τ is decreased. When some DOPS head groups are replaced by zwitterionic head groups, τ is similarly decreased. These effects do not depend on the sign of the change in surface charge. In DOPE:DOPC (3:1) bilayers, pH changes from 5→9 to 5→0 (both increasing head group electrostatic repulsion, thereby causing a less negative c 0) both increase τ. Nor do the effects depend on the use of planar, hydrocarbon-containing bilayers, as similar changes were observed in hydrocarbon-free lipid vesicles. Altering the interactions among phospholipid head groups may alter also other bilayer properties such as thickness or elastic moduli. Such changes could be excluded using capacitance measurements and single channel measurements on gA channels of different lengths. We conclude that changes gA channel function caused by changes in head group effective size can be predicted from the expected changes in c 0.