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

A millisecond coarse-grained simulation approach to decipher allosteric cannabinoid binding at the glycine receptor α1

Glycine receptors (GlyR) are regulated by small-molecule binding at several allosteric sites. Cannabinoids like tetrahydrocannabinol (THC) and N-arachidonyl-ethanol-amide (AEA) potentiate the GlyR response but their mechanism of action is not fully established. By combining millisecond coarse-grained (CG) MD simulations powered by Martini 3 with backmapping to all-atom representations, we have characterized the cannabinoid-binding site(s) at the zebrafish GlyR-α1 active state with atomic resolution. Based on hundreds of thousand ligand-binding events, we find that cannabinoids bind to the transmembrane domain of the receptor at both intrasubunit and intersubunit sites. For THC, the intrasubunit binding mode predicted in simulation is in excellent agreement with recent cryo-EM structures, while intersubunit binding recapitulates in full previous mutagenesis experiments. Intriguingly, AEA is predicted to bind at the same intersubunit site despite the strikingly different chemistry. Statistical analyses of the ligand-receptor interactions highlight potentially relevant residues for GlyR potentiation, offering experimentally testable predictions. The predictions for AEA have been validated by electrophysiology recordings of rationally designed mutants. The results highlight the existence of multiple cannabinoid-binding sites for the allosteric regulation of GlyR and put forward an effective strategy for the identification and structural characterization of allosteric binding sites.

Alchemical Transformations and Beyond: Recent Advances and Real-World Applications of Free Energy Calculations in Drug Discovery

Computational methods constitute efficient strategies for screening and optimizing potential drug molecules. A critical factor in this process is the binding affinity between candidate molecules and targets, quantified as binding free energy. Among various estimation methods, alchemical transformation methods stand out for their theoretical rigor. Despite challenges in force field accuracy and sampling efficiency, advancements in algorithms, software, and hardware have increased the application of free energy perturbation (FEP) calculations in the pharmaceutical industry. Here, we review the practical applications of FEP in drug discovery projects since 2018, covering both ligand-centric and residue-centric transformations. We show that relative binding free energy calculations have steadily achieved chemical accuracy in real-world applications. In addition, we discuss alternative physics-based simulation methods and the incorporation of deep learning into free energy calculations.

Allosteric Cholesterol Site in Glycine Receptors Characterized through Molecular Simulations

Glycine receptors are pentameric ligand-gated ion channels that conduct chloride ions across postsynaptic membranes to facilitate fast inhibitory neurotransmission. In addition to gating by the glycine agonist, interactions with lipids and other compounds in the surrounding membrane environment modulate their function, but molecular details of these interactions remain unclear, in particular, for cholesterol. Here, we report coarse-grained simulations in a model neuronal membrane for three zebrafish glycine receptor structures representing apparent resting, open, and desensitized states. We then converted the systems to all-atom models to examine detailed lipid interactions. Cholesterol bound to the receptor at an outer-leaflet intersubunit site, with a preference for the open and desensitized versus resting states, indicating that it can bias receptor function. Finally, we used short atomistic simulations and iterative amino acid perturbations to identify residues that may mediate allosteric gating transitions. Frequent cholesterol contacts in atomistic simulations clustered with residues identified by perturbation analysis and overlapped with mutations influencing channel function and pathology. Cholesterol binding at this site was also observed in a recently reported pig heteromeric glycine receptor. These results indicate state-dependent lipid interactions relevant to allosteric transitions of glycine receptors, including specific amino acid contacts applicable to biophysical modeling and pharmaceutical design.

State-dependent binding of cholesterol and an anionic lipid to the muscle-type Torpedo nicotinic acetylcholine receptor

The ability of the Torpedo nicotinic acetylcholine receptor (nAChR) to undergo agonist-induced conformational transitions requires the presence of cholesterol and/or anionic lipids. Here we use recently solved structures along with multiscale molecular dynamics simulations to examine lipid binding to the nAChR in bilayers that have defined effects on nAChR function. We examine how phosphatidic acid and cholesterol, lipids that support conformational transitions, individually compete for binding with phosphatidylcholine, a lipid that does not. We also examine how the two lipids work synergistically to stabilize an agonist-responsive nAChR. We identify rapidly exchanging lipid binding sites, including both phospholipid sites with a high affinity for phosphatidic acid and promiscuous cholesterol binding sites in the grooves between adjacent transmembrane α-helices. A high affinity cholesterol site is confirmed in the inner leaflet framed by a key tryptophan residue on the MX α-helix. Our data provide insight into the dynamic nature of lipid-nAChR interactions and set the stage for a detailed understanding of the mechanisms by which lipids facilitate nAChR function at the neuromuscular junction.

Cryo-EM structures of prokaryotic ligand-gated ion channel GLIC provide insights into gating in a lipid environment

GLIC, a proton-activated prokaryotic ligand-gated ion channel, served as a model system for understanding the eukaryotic counterparts due to their structural and functional similarities. Despite extensive studies conducted on GLIC, the molecular mechanism of channel gating in the lipid environment requires further investigation. Here, we present the cryo-EM structures of nanodisc-reconstituted GLIC at neutral and acidic pH in the resolution range of 2.6 - 3.4 Å. In our apo state at pH 7.5, the extracellular domain (ECD) displays conformational variations compared to the existing apo structures. At pH 4.0, three distinct conformational states (C1, C2 and O states) are identified. The protonated structures exhibit a compacted and counter-clockwise rotated ECD compared with our apo state. A gradual widening of the pore in the TMD is observed upon reducing the pH, with the widest pore in O state, accompanied by several layers of water pentagons. The pore radius and molecular dynamics (MD) simulations suggest that the O state represents an open conductive state. We also observe state-dependent interactions between several lipids and proteins that may be involved in the regulation of channel gating. Our results provide comprehensive insights into the importance of lipids impact on gating.

Discovery of lipid binding sites in a ligand-gated ion channel by integrating simulations and cryo-EM

Ligand-gated ion channels transduce electrochemical signals in neurons and other excitable cells. Aside from canonical ligands, phospholipids are thought to bind specifically to the transmembrane domain of several ion channels. However, structural details of such lipid contacts remain elusive, partly due to limited resolution of these regions in experimental structures. Here, we discovered multiple lipid interactions in the channel GLIC by integrating cryo-electron microscopy and large-scale molecular simulations. We identified 25 bound lipids in the GLIC closed state, a conformation where none, to our knowledge, were previously known. Three lipids were associated with each subunit in the inner leaflet, including a buried interaction disrupted in mutant simulations. In the outer leaflet, two intrasubunit sites were evident in both closed and open states, while a putative intersubunit site was preferred in open-state simulations. This work offers molecular details of GLIC-lipid contacts particularly in the ill-characterized closed state, testable hypotheses for state-dependent binding, and a multidisciplinary strategy for modeling protein-lipid interactions.

Startle Disease: New Molecular Insights into an Old Neurological Disorder

Startle disease (SD) is characterized by enhanced startle responses, generalized muscle stiffness, unexpected falling, and fatal apnea episodes due to disturbed feedback inhibition in the spinal cord and brainstem of affected individuals. Mutations within the glycine receptor (GlyR) subunit and glycine transporter 2 (GlyT2) genes have been identified in individuals with SD. Impaired inhibitory neurotransmission in SD is due to pre- and/or postsynaptic GlyR or presynaptic GlyT2 dysfunctions. Previous research has focused on mutated GlyRs and GlyT2 that impair ion channel/transporter function or trafficking. With insights provided by recently solved cryo-electron microscopy and X-ray structures of GlyRs, a detailed picture of structural transitions important for receptor gating has emerged, allowing a deeper understanding of SD at the molecular level. Moreover, studies on novel SD mutations have demonstrated a higher complexity of SD, with identification of additional clinical signs and symptoms and interaction partners representing key players for fine-tuning synaptic processes. Although our knowledge has steadily improved during the last years, changes in synaptic localization and GlyR or GlyT2 homeostasis under disease conditions are not yet completely understood. Combined proteomics, interactomics, and high-resolution microscopy techniques are required to reveal alterations in receptor dynamics at the synaptic level under disease conditions.

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.

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.

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.

Mutations of the nACh Receptor M4 Helix Reveal Different Phenotypes in Different Expression Systems: Could Lipids be Responsible?

The role of the outermost helix (M4) in the pentameric ligand-gated ion channel (pLGIC) family is currently not fully understood. It is known that M4 is important for receptor assembly, possibly via interactions with neighboring M1 and M3 helices. M4 can also transmit information on the lipid content of the membrane to the gating mechanism, and it may form a link to the extracellular domain via the Cys-loop. Our previous study examining the α4β2 nACh receptor M4 helix using HEK cells indicated M4 here is more sensitive to change than those of other pLGIC. Many of these other studies, however, were performed in Xenopus oocytes. Here we examine the nine previously identified nonfunctional α4β2 nACh receptor M4 mutant receptors using this system. The data reveal that seven of these mutant receptors do function when expressed in oocytes, with only 2, the conserved Asp at the intracellular end of M4 and a Phe in the center, having a similar phenotype (nonfunctional) in both HEK cells and oocytes. The oocyte data are more consistent with studies in other pLGIC and demonstrate the importance of the expression system used. Of the many differences between these two expression systems, we suggest that the different lipid content of the plasma membrane is a possible candidate for explaining these discrepancies.

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.

Polyunsaturated fatty acids inhibit a pentameric ligand-gated ion channel through one of two binding sites

Polyunsaturated fatty acids (PUFAs) inhibit pentameric ligand-gated ion channels (pLGICs) but the mechanism of inhibition is not well understood. The PUFA, docosahexaenoic acid (DHA), inhibits agonist responses of the pLGIC, ELIC, more effectively than palmitic acid, similar to the effects observed in the GABA_A receptor and nicotinic acetylcholine receptor. Using photo-affinity labeling and coarse-grained molecular dynamics simulations, we identified two fatty acid binding sites in the outer transmembrane domain (TMD) of ELIC. Fatty acid binding to the photolabeled sites is selective for DHA over palmitic acid, and specific for an agonist-bound state. Hexadecyl-methanethiosulfonate modification of one of the two fatty acid binding sites in the outer TMD recapitulates the inhibitory effect of PUFAs in ELIC. The results demonstrate that DHA selectively binds to multiple sites in the outer TMD of ELIC, but that state-dependent binding to a single intrasubunit site mediates DHA inhibition of ELIC.

A Single Mutation in the Outer Lipid-Facing Helix of a Pentameric Ligand-Gated Ion Channel Affects Channel Function Through a Radially-Propagating Mechanism

Pentameric ligand-gated ion channels (pLGICs) mediate fast synaptic transmission and are crucial drug targets. Their gating mechanism is triggered by ligand binding in the extracellular domain that culminates in the opening of a hydrophobic gate in the transmembrane domain. This domain is made of four α-helices (M1 to M4). Recently the outer lipid-facing helix (M4) has been shown to be key to receptor function, however its role in channel opening is still poorly understood. It could act through its neighboring helices (M1/M3), or via the M4 tip interacting with the pivotal Cys-loop in the extracellular domain. Mutation of a single M4 tyrosine (Y441) to alanine renders one pLGIC-the 5-HT3A receptor-unable to function despite robust ligand binding. Using Y441A as a proxy for M4 function, we here predict likely paths of Y441 action using molecular dynamics, and test these predictions with functional assays of mutant receptors in HEK cells and Xenopus oocytes using fluorescent membrane potential sensitive dye and two-electrode voltage clamp respectively. We show that Y441 does not act via the M4 tip or Cys-loop, but instead connects radially through M1 to a residue near the ion channel hydrophobic gate on the pore-lining helix M2. This demonstrates the active role of the M4 helix in channel opening.