A Structural Model of the Human α7 Nicotinic Receptor in an Open Conformation

Molecular Dynamics Refinement of Open State Serotonin 5-HT3A Receptor Structures

Pentameric ligand-gated ion channels play an important role in mediating fast neurotransmissions. As a member of this receptor family, cation-selective 5-HT₃ receptors are a clinical target for treating nausea and vomiting associated with chemotherapy and radiation therapy (Thompson and Lummis, 2006). Multiple cryo-electron microscopy (cryo-EM) structures of 5-HT₃ receptors have been determined in distinct functional states (e.g., open, closed, etc.) (Basak et al., 2018; Basak et al., 2018; Polovinkin et al., 2018; Zhang et al., 2015). However, recent work has shown that the transmembrane pores of the open 5-HT₃ receptor structures rapidly collapse and become artificially asymmetric in molecular dynamics (MD) simulations. To avoid this hydrophobic collapse, Dämgen and Biggin developed an equilibration protocol that led to a stable open state structure of the glycine receptor in MD simulations (Dämgen and Biggin, 2020). However, the protocol failed to yield open-like structures of the 5-HT₃ receptor in our simulations. Here, we present a refined equilibration protocol that involves the rearrangement of the transmembrane helices to achieve stable open state structures of the 5-HT₃ receptor that allow both water and ion permeation through the channel. Notably, channel gating is mediated through collective movement of the transmembrane helices, involving not only pore lining M2 helices but also their cross-talk with the adjacent M1 and M3 helices. Thus, the successful application of our refined equilibration protocol underscores the importance of the conformational coupling between the transmembrane helices in stabilizing open-like structures of the 5-HT₃ receptor.

Computational methods and theory for ion channel research

Ion channels are fundamental biological devices that act as gates in order to ensure selective ion transport across cellular membranes; their operation constitutes the molecular mechanism through which basic biological functions, such as nerve signal transmission and muscle contraction, are carried out. Here, we review recent results in the field of computational research on ion channels, covering theoretical advances, state-of-the-art simulation approaches, and frontline modeling techniques. We also report on few selected applications of continuum and atomistic methods to characterize the mechanisms of permeation, selectivity, and gating in biological and model channels.

A Computational Analysis of the Factors Governing the Dynamics of α7 nAChR and Its Homologs

The α7 nicotinic acetylcholine receptor is a homopentameric ion channel from the Cys-loop receptor superfamily targeted for psychiatric indications and inflammatory pain. Molecular dynamics studies of the receptor have focused on residue mobility and global conformational changes to address receptor function. However, a comparative analysis of α7 with its homologs that cannot trigger channel opening has not been made so far. To identify the residues involved in α7 activation, we ran triplicate 500-ns molecular dynamics simulations with an α7 extracellular domain homology model and two acetylcholine-binding protein homologs. We tested the effect of ligand binding and amino acid sequence on the structure and dynamics of the three proteins. We found that mobile regions identified based on root mean-square deviation and root mean-square fluctuation values are not always consistent among the individual α7 extracellular domain simulations. Comparison of the replica-average properties of the three proteins based on dynamic cross-correlation maps showed that ligand binding affects the coupling between the C-loop and the Cys-loop, vestibular loop, and β1-β2 loops. In addition, the main-immunogenic-region-like domain of α7 went through correlated motions with multiple domains of the receptor. These correlated motions were absent or diminished in α7 homologs, suggesting a unique role in α7 activation.

Thermodynamics and Kinetics of Ion Permeation in Wild-Type and Mutated Open Active Conformation of the Human α7 Nicotinic Receptor

Molecular studies of human pentameric ligand-gated ion channels (LGICs) expressed in neurons and at neuromuscular junctions are of utmost importance in the development of therapeutic strategies for neurological disorders. We focus here on the nicotinic acetylcholine receptor nAChR-α7, a homopentameric channel widely expressed in the human brain, with a proven role in a wide spectrum of disorders including schizophrenia and Alzheimer’s disease. By exploiting an all-atom structural model of the full (transmembrane and extracellular) protein in the open, agonist-bound conformation we recently developed, we evaluate the free energy and the mean first passage time of single-ion permeation using molecular dynamics simulations and the milestoning method with Voronoi tessellation. The results for the wild-type channel provide the first available mapping of the potential of mean force in the full-length α7 nAChR, reveal its expected cationic nature, and are in good agreement with simulation data for other channels of the LGIC family and with experimental data on nAChRs. We then investigate the role of a specific mutation directly related to ion selectivity in LGICs, the E-1′ → A-1′ substitution at the cytoplasmatic selectivity filter. We find that the mutation strongly affects sodium and chloride permeation in opposite directions, leading to a complete inversion of selectivity, at variance with the limited experimental results available that classify this mutant as cationic. We thus provide structural determinants for the observed cationic-to-anionic inversion, revealing a key role of the protonation state of residue rings far from the mutation, in the proximity of the hydrophobic channel gate.

Nicotinic receptor pharmacology in silico: Insights and challenges

Nicotinic acetylcholine receptors (nAChR) are homo- or hetero-pentameric ligand-gated ion channels of the Cys-loop superfamily and play important roles in the nervous system and muscles. Studies on nAChR benefit from in silico modeling due to the lack of high-resolution structures for most receptor subtypes and challenges in experiments addressing the complex mechanism of activation involving allosteric sites. Although there is myriad of computational modeling studies on nAChR, the multitude of the methods and parameters used in these studies makes modeling nAChR a daunting task, particularly for the non-experts in the field. To address this problem, the modeling literature on Torpedo nAChR and α7 nAChR were focused on as examples of heteromeric and homomeric nAChR, and the key in silico modeling studies between the years 1995-2019 were concisely reviewed. This was followed by a critical analysis of these studies by comparing the findings with each other and with the emerging experimental and computational data on nAChR. Based on these critical analyses, suggestions were made to guide the future researchers in the field of in silico modeling of nAChR. This article is part of the special issue on 'Contemporary Advances in Nicotine Neuropharmacology'.

A Refined Open State of the Glycine Receptor Obtained via Molecular Dynamics Simulations

Pentameric ligand-gated ion channels are key players in mediating fast neurotransmission. Glycine receptors are chloride-selective members of this receptor family that mediate inhibitory synaptic transmission and are implicated in neurological disorders including autism and hyperekplexia. They have been structurally characterized by both X-ray crystallography and cryoelectron microscopy (cryo-EM) studies, with the latter giving rise to what was proposed as a possible open state. However, recent work has questioned the physiological relevance of this open state structure, since it rapidly collapses in molecular dynamics simulations. Here, we show that the collapse can be avoided by a careful equilibration protocol that reconciles the more problematic regions of the original density map and gives a stable open state that shows frequent selective chloride permeation. The protocol developed in this work provides a means to refine open-like structures of the whole pentameric ligand-gated ion channel superfamily and reconciles the previous issues with the cryo-EM structure.

In Silico Modeling of the α7 Nicotinic Acetylcholine Receptor: New Pharmacological Challenges Associated with Multiple Modes of Signaling

The α7 nicotinic acetylcholine receptor is a homopentameric ion-channel of the Cys-loop superfamily characterized by its low probability of opening, high calcium permeability, and rapid desensitization. The α7 receptor has been targeted for the treatment of the cognitive symptoms of schizophrenia, depression, and Alzheimer's disease, but it is also involved in inflammatory modulation as a part of the cholinergic anti-inflammatory pathway. Despite its functional importance, in silico studies of the α7 receptor cannot produce a general model explaining the structural features of receptor activation, nor predict the mode of action for various ligand classes. Two particular problems in modeling the α7 nAChR are the absence of a high-resolution structure and the presence of five potentially nonequivalent orthosteric ligand binding sites. There is wide variability regarding the templates used for homology modeling, types of ligands investigated, simulation methods, and simulation times. However, a systematic survey focusing on the methodological similarities and differences in modeling α7 has not been done. In this work, we make a critical analysis of the modeling literature of α7 nAChR by comparing the findings of computational studies with each other and with experimental studies under the main topics of structural studies, ligand binding studies, and comparisons with other nAChR. In light of our findings, we also summarize current problems in the field and make suggestions for future studies concerning modeling of the α7 receptor.

Closed-Locked and Apo-Resting State Structures of the Human α7 Nicotinic Receptor: A Computational Study

Nicotinic acetylcholine receptors, belonging to the Cys-loop superfamily of ligand-gated ion channels (LGICs), are membrane proteins present in neurons and at neuromuscular junctions. They are responsible for signal transmission, and their function is regulated by neurotransmitters, agonists, and antagonists drugs. A detailed knowledge of their conformational transition in response to ligand binding is critical to understanding the basis of ligand-receptor interaction, in view of new pharmacological approaches to control receptor activity. However, the scarcity of experimentally derived structures of human channels makes this perspective extremely challenging. To contribute overcoming this issue, we have recently reported structural models for the open and the desensitized states of the human α7 nicotinic receptor. Here, we provide all-atom structural models of the same receptor in two different nonconductive states. The first structure, built via homology modeling and relaxed with extensive Molecular Dynamics simulations, represents the receptor bound to the natural antagonist α-conotoxin ImI. After comparison with available experimental data and computational models of other eukaryotic LGICs, we deem it consistent with the "closed-locked" state. The second model, obtained with simulations from the spontaneous relaxation of the open, agonist-bound α7 structure after ligand removal, recapitulates the characteristics of the apo-resting state of the receptor. These results add to our previous work on the active and desensitized state conformations, contributing to the structural characterization of the conformational landscape of the human α7 receptor and suggesting benchmarks to discriminate among conformations found in experiments or in simulations of LGICs. In particular key interactions at the interface between the extracellular domain and the transmembrane domain are identified, that could be critical to the α7 receptor function.

Discovery, synthesis, biological evaluation and structure-based optimization of novel piperidine derivatives as acetylcholine-binding protein ligands

The homomeric α7 nicotinic receptor (α7 nAChR) is widely expressed in the human brain that could be activated to suppress neuroinflammation, oxidative stress and neuropathic pain. Consequently, a number of α7 nAChR agonists have entered clinical trials as anti-Alzheimer's or anti-psychotic therapies. However, high-resolution crystal structure of the full-length α7 receptor is thus far unavailable. Since acetylcholine-binding protein (AChBP) from Lymnaea stagnalis is most closely related to the α-subunit of nAChRs, it has been used as a template for the N-terminal domain of α-subunit of nAChR to study the molecular recognition process of nAChR-ligand interactions, and to identify ligands with potential nAChR-like activities.Here we report the discovery and optimization of novel acetylcholine-binding protein ligands through screening, structure-activity relationships and structure-based design. We manually screened in-house CNS-biased compound library in vitro and identified compound 1, a piperidine derivative, as an initial hit with moderate binding affinity against AChBP (17.2% inhibition at 100 nmol/L). During the 1st round of optimization, with compound 2 (21.5% inhibition at 100 nmol/L) as the starting point, 13 piperidine derivatives with different aryl substitutions were synthesized and assayed in vitro. No apparent correlation was demonstrated between the binding affinities and the steric or electrostatic effects of aryl substitutions for most compounds, but compound 14 showed a higher affinity (K_i=105.6 nmol/L) than nicotine (K_i=777 nmol/L). During the 2nd round of optimization, we performed molecular modeling of the putative complex of compound 14 with AChBP, and compared it with the epibatidine-AChBP complex. The results suggested that a different piperidinyl substitution might confer a better fit for epibatidine as the reference compound. Thus, compound 15 was designed and identified as a highly affinitive acetylcholine-binding protein ligand. In this study, through two rounds of optimization, compound 15 (K_i=2.8 nmol/L) has been identified as a novel, piperidine-based acetylcholine-binding protein ligand with a high affinity.