Pentameric Ligand-gated Ion Channels : Insights from Computation

Neurosteroids and their potential as a safer class of general anesthetics

Neurosteroids (NS) are a class of steroids that are synthesized within the central nervous system (CNS). Various NS can either enhance or inhibit CNS excitability and they play important biological roles in brain development, brain function and as mediators of mood. One class of NS, 3α-hydroxy-pregnane steroids such as allopregnanolone (AlloP) or pregnanolone (Preg), inhibits neuronal excitability; these endogenous NS and their analogues have been therapeutically applied as anti-depressants, anti-epileptics and general anesthetics. While NS have many favorable properties as anesthetics (e.g. rapid onset, rapid recovery, minimal cardiorespiratory depression, neuroprotection), they are not currently in clinical use, largely due to problems with formulation. Recent advances in understanding NS mechanisms of action and improved formulations have rekindled interest in development of NS as sedatives and anesthetics. In this review, the synthesis of NS, and their mechanism of action will be reviewed with specific emphasis on their binding sites and actions on γ-aminobutyric acid type A (GABAA) receptors. The potential advantages of NS analogues as sedative and anesthetic agents will be discussed.

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.

Carboxylate binding prefers two cations to one

Almost all studies of specific ion binding by carboxylates (-COO^-) have considered only a single cation, but clustering of ions and ligands is a common phenomenon. We apply density functional theory to investigate how variations in the number of acetate ligands in binding to two monovalent cations affects ion binding preferences. We study a series of monovalent (Li⁺, Na⁺, K⁺, Cs⁺) ions relevant to experimental work on many topics, including ion channels, battery storage, water purification and solar cells. We find that the preferred optimal structure has 3 acetates except for Cs⁺, which has 2 acetates. The optimal coordination of the cation by the carboxylate O atoms is 4 for both Na⁺ and K⁺, and 3 for Li⁺ and Cs⁺. There is a 4-fold coordination minimum just a few kcal mol^-1 higher than the optimal 3-fold structure for Li⁺. For two cations, multiple minima occur in the vicinity of the lowest free energy state. We find that, for Li, Na and K, the preferred optimal structure with two cations is favored over a mixture of single cation complexes, providing a basis for understanding ionic cluster formation that is relevant for engineering proteins and other materials for rapid, selective ion transport.

Molecular Mingling: Multimodal Predictions of Ligand Promiscuity in Pentameric Ligand-Gated Ion Channels

Background: Human pentameric ligand-gated ion channels (pLGICs) comprise nicotinic acetylcholine receptors (nAChRs), 5-hydroxytryptamine type 3 receptors (5-HT3Rs), zinc-activated channels (ZAC), γ-aminobutyric acid type A receptors (GABAARs) and glycine receptors (GlyRs). They are recognized therapeutic targets of some of the most prescribed drugs like general anesthetics, anxiolytics, smoking cessation aids, antiemetics and many more. Currently, approximately 100 experimental structures of pLGICs with ligands bound exist in the protein data bank (PDB). These atomic-level 3D structures enable the generation of a comprehensive binding site inventory for the superfamily and the in silico prediction of binding site properties. Methods: A panel of high throughput in silico methods including pharmacophore screening, conformation analysis and descriptor calculation was applied to a selection of allosteric binding sites for which in vitro screens are lacking. Variant abundance near binding site forming regions and computational docking complement the approach. Results: The structural data reflects known and novel binding sites, some of which may be unique to individual receptors, while others are broadly conserved. The membrane spanning domain, comprising four highly conserved segments, contains ligand interaction sites for which in vitro assays suitable for high throughput screenings are critically lacking. This is also the case for structurally more variable novel sites in the extracellular domain. Our computational results suggest that the phytocannabinoid Δ9-tetrahydrocannabinol (Δ9-THC) can utilize multiple pockets which are likely to exist on most superfamily members. Conclusion: With this study, we explore the potential for polypharmacology among pLGICs. Our data suggest that ligands can display two forms of promiscuity to an extent greater than what has been realized: 1) Ligands can interact with homologous sites in many members of the superfamily, which bears toxicological relevance. 2) Multiple pockets in distinct localizations of individual receptor subtypes share common ligands, which counteracts efforts to develop selective agents. Moreover, conformational states need to be considered for in silico drug screening, as certain binding sites display considerable flexibility. In total, this work contributes to a better understanding of polypharmacology across pLGICs and provides a basis for improved structure guided in silico drug development and drug derisking.

Photopharmacology of Ion Channels through the Light of the Computational Microscope

The optical control and investigation of neuronal activity can be achieved and carried out with photoswitchable ligands. Such compounds are designed in a modular fashion, combining a known ligand of the target protein and a photochromic group, as well as an additional electrophilic group for tethered ligands. Such a design strategy can be optimized by including structural data. In addition to experimental structures, computational methods (such as homology modeling, molecular docking, molecular dynamics and enhanced sampling techniques) can provide structural insights to guide photoswitch design and to understand the observed light-regulated effects. This review discusses the application of such structure-based computational methods to photoswitchable ligands targeting voltage- and ligand-gated ion channels. Structural mapping may help identify residues near the ligand binding pocket amenable for mutagenesis and covalent attachment. Modeling of the target protein in a complex with the photoswitchable ligand can shed light on the different activities of the two photoswitch isomers and the effect of site-directed mutations on photoswitch binding, as well as ion channel subtype selectivity. The examples presented here show how the integration of computational modeling with experimental data can greatly facilitate photoswitchable ligand design and optimization. Recent advances in structural biology, both experimental and computational, are expected to further strengthen this rational photopharmacology approach.

Molecular dynamics: a powerful tool for studying the medicinal chemistry of ion channel modulators

Molecular dynamics (MD) simulations allow researchers to investigate the behavior of desired biological targets at ever-decreasing costs with ever-increasing precision. Among the biological macromolecules, ion channels are remarkable transmembrane proteins, capable of performing special biological processes and revealing a complex regulatory matrix, including modulation by small molecules, either endogenous or exogenous. Recently, given the developments in ion channel structure determination and accessibility of bio-computational techniques, MD and related tools are becoming increasingly popular in the intense research area regarding ligand-channel interactions. This review synthesizes and presents the most important fields of MD involvement in investigating channel-molecule interactions, including, but not limited to, deciphering the binding modes of ligands to their ion channel targets and the mechanisms through which chemical compounds exert their effect on channel function. Special attention is devoted to the importance of more elaborate methods, such as free energy calculations, while principles regarding drug design and discovery are highlighted. Several technical aspects involving the creation and simulation of channel-molecule MD systems (ligand parameterization, proper membrane setup, system building, etc.) are also presented.

Elephants in the Dark: Insights and Incongruities in Pentameric Ligand-gated Ion Channel Models

The superfamily of pentameric ligand-gated ion channels (pLGICs) comprises key players in electrochemical signal transduction across evolution, including historic model systems for receptor allostery and targets for drug development. Accordingly, structural studies of these channels have steadily increased, and now approach 250 depositions in the protein data bank. This review contextualizes currently available structures in the pLGIC family, focusing on morphology, ligand binding, and gating in three model subfamilies: the prokaryotic channel GLIC, the cation-selective nicotinic acetylcholine receptor, and the anion-selective glycine receptor. Common themes include the challenging process of capturing and annotating channels in distinct functional states; partially conserved gating mechanisms, including remodeling at the extracellular/transmembrane-domain interface; and diversity beyond the protein level, arising from posttranslational modifications, ligands, lipids, and signaling partners. Interpreting pLGIC structures can be compared to describing an elephant in the dark, relying on touch alone to comprehend the many parts of a monumental beast: each structure represents a snapshot in time under specific experimental conditions, which must be integrated with further structure, function, and simulations data to build a comprehensive model, and understand how one channel may fundamentally differ from another.

Functional characterization of neurotransmitter activation and modulation in a nematode model ligand-gated ion channel

The superfamily of pentameric ligand-gated ion channels includes neurotransmitter receptors that mediate fast synaptic transmission in vertebrates, and are targets for drugs including alcohols, anesthetics, benzodiazepines, and anticonvulsants. However, the mechanisms of ion channel opening, gating, and modulation in these receptors leave many open questions, despite their pharmacological importance. Subtle conformational changes in both the extracellular and transmembrane domains are likely to influence channel opening, but have been difficult to characterize given the limited structural data available for human membrane proteins. Recent crystal structures of a modified Caenorhabditis elegans glutamate-gated chloride channel (GluCl) in multiple states offer an appealing model system for structure-function studies. However, the pharmacology of the crystallographic GluCl construct is not well established. To establish the functional relevance of this system, we used two-electrode voltage-clamp electrophysiology in Xenopus oocytes to characterize activation of crystallographic and native-like GluCl constructs by L-glutamate and ivermectin. We also tested modulation by ethanol and other anesthetic agents, and used site-directed mutagenesis to explore the role of a region of Loop F which was implicated in ligand gating by molecular dynamics simulations. Our findings indicate that the crystallographic construct functionally models concentration-dependent agonism and allosteric modulation of pharmacologically relevant receptors. Specific substitutions at residue Leu174 in loop F altered direct L-glutamate activation, consistent with computational evidence for this region's role in ligand binding. These insights demonstrate conservation of activation and modulation properties in this receptor family, and establish a framework for GluCl as a model system, including new possibilities for drug discovery. In this study, we elucidate the validity of a modified glutamate-gated chloride channel (GluClcryst ) as a structurally accessible model for GABAA receptors. In contrast to native-like controls, GluClcryst exhibits classical activation by its neurotransmitter ligand L-glutamate. The modified channel is also sensitive to allosteric modulators associated with human GABAA receptors, and to site-directed mutations predicted to alter channel opening.

Beyond Membrane Protein Structure: Drug Discovery, Dynamics and Difficulties

Most of the previous content of this book has focused on obtaining the structures of membrane proteins. In this chapter we explore how those structures can be further used in two key ways. The first is their use in structure based drug design (SBDD) and the second is how they can be used to extend our understanding of their functional activity via the use of molecular dynamics. Both aspects now heavily rely on computations. This area is vast, and alas, too large to consider in depth in a single book chapter. Thus where appropriate we have referred the reader to recent reviews for deeper assessment of the field. We discuss progress via the use of examples from two main drug target areas; G-protein coupled receptors (GPCRs) and ion channels. We end with a discussion of some of the main challenges in the area.