Channel opening motion of alpha7 nicotinic acetylcholine receptor as suggested by normal mode analysis

Comparative Analysis of Structural and Dynamical Features of Ribosome Upon Association With mRNA Reveals Potential Role of Ribosomal Proteins

Ribosomes play a critical role in maintaining cellular proteostasis. The binding of messenger RNA (mRNA) to the ribosome regulates kinetics of protein synthesis. To generate an understanding of the structural, mechanistic, and dynamical features of mRNA recognition in the ribosome, we have analysed mRNA-protein interactions through a structural comparison of the ribosomal complex in the presence and absence of mRNA. To do so, we compared the 3-Dimensional (3D) structures of components of the two assembly structures and analysed their structural differences because of mRNA binding, using elastic network models and structural network-based analysis. We observe that the head region of 30S ribosomal subunit undergoes structural displacement and subunit rearrangement to accommodate incoming mRNA. We find that these changes are observed in proteins that lie far from the mRNA-protein interface, implying allostery. Further, through perturbation response scanning, we show that the proteins S13, S19, and S20 act as universal sensors that are sensitive to changes in the inter protein network, upon binding of 30S complex with mRNA and other initiation factors. Our study highlights the significance of mRNA binding in the ribosome complex and identifies putative allosteric sites corresponding to alterations in structure and/or dynamics, in regions away from mRNA binding sites in the complex. Overall, our work provides fresh insights into mRNA association with the ribosome, highlighting changes in the interactions and dynamics of the ribosome assembly because of the binding.

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.

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'.

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.

Toward a Model for Activation of Orai Channel

Store-operated calcium release-activated calcium (CRAC) channels mediate a variety of cellular signaling functions. The CRAC channel pore-forming protein, Orai1, is a hexamer arranged with 3-fold symmetry. Despite its importance in moving Ca²⁺ ions into cells, a detailed mechanistic understanding of Orai1 activation is lacking. Herein, a working model is proposed for the putative open state of Orai from Drosophila melanogaster (dOrai), which involves a “twist-to-open” gating mechanism. The proposed model is supported by energetic, structural, and experimental evidence. Fluorescent imaging demonstrates that each subunit on the intracellular side of the pore is inherently strongly cross-linked, which is important for coupling to STIM1, the pore activator, and graded activation of the Orai1 channel. The proposed model thus paves the way for understanding key aspects of calcium signaling at a molecular level.

•

Mechanical coupling within the calcium channel pore is critical for its activation

•

Molecular modeling could disclose gating mechanism of ion channels at atomic level

•

The predicted open-state structure of the pore was further confirmed by experiments

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'.

Modeling ion permeation in wild-type and mutant human α7 nachr ion channels

Molecular dynamics simulations of wild type and two mutant (T248F and L251T) human [Formula: see text]7 nicotinic acetylcholine receptors (nAChR) have been performed. The channel transmembrane domains were modeled from the closed channel structure from torpedo ray (PDB ID 2BG9) and embedded in DPPC lipid bilayers, surrounded by physiological saline solution. An external electric field was used to obtain stable open channel structures. The adaptive biasing force (ABF) method was used to obtain potential of mean force (PMF) profiles for Na[Formula: see text] ion translocation through the wild type and mutant receptors. Based on the geometry and PMF profiles, the channel gate was found to be at one of the two hydrophobic conserved regions (V249-L251) near the lower end of the channel. The L251T mutation reduced the energetic barrier by 1.9[Formula: see text]kcal/mol, consistent with a slight increase in the channel radius in the bottleneck region. On the other hand, the T248F mutation caused a significant decrease in the channel radius (0.4 Å) and a substantial increase of 3.9[Formula: see text]kcal/mol in the energetic barrier. Ion permeation in all three structures was compared and found to be consistent with barrier height values. Using an external field in an incrementally increasing manner was found to be an effective way to obtain stable open, conducting channel structures.

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.

Molecular simulation study of the unbinding of α-conotoxin [ϒ4E]GID at the α7 and α4β2 neuronal nicotinic acetylcholine receptors

The α7 and α4β2 neuronal nicotinic receptors belonging to the family of ligand-gated ion channels are most prevalent in the brain, and are implicated in various neurodegenerative disorders. α-conotoxin GID (and its analogue [ϒ4E]GID) specifically inhibits these subtypes, with more affinity towards the human α7 (hα7) subtype, and is valuable in understanding the physiological roles of these receptors. In this study, we use umbrella-sampling molecular dynamics simulations to understand the mechanism of interaction between [ϒ4E]GID and the agonist binding pockets of the α4β2 and the hα7 receptors, and to estimate their relative binding affinities (ΔG_bind). The obtained ΔG_bind values indicate stronger interaction with the hα7 receptor, in agreement with previous experimental studies. Simulations also revealed different unbinding pathways between the two receptor subtypes, enabling identification of a number of interactions at locations far from the orthosteric binding site which may explain the difference in [ϒ4E]GID potency. The pathways identified will help in the design of novel conotoxins with increased potency at α4β2, for which there is currently no known highly potent conotoxin inhibitor. Computational mutational free energy analyses also revealed a number of possible single-site mutations to GID which might enhance its selective binding to α4β2 over α7.

RNA-DNA sequence differences in Saccharomyces cerevisiae

Alterations of RNA sequences and structures, such as those from editing and alternative splicing, result in two or more RNA transcripts from a DNA template. It was thought that in yeast, RNA editing only occurs in tRNAs. Here, we found that Saccharomyces cerevisiae have all 12 types of RNA–DNA sequence differences (RDDs) in the mRNA. We showed these sequence differences are propagated to proteins, as we identified peptides encoded by the RNA sequences in addition to those by the DNA sequences at RDD sites. RDDs are significantly enriched at regions with R-loops. A screen of yeast mutants showed that RDD formation is affected by mutations in genes regulating R-loops. Loss-of-function mutations in ribonuclease H, senataxin, and topoisomerase I that resolve RNA–DNA hybrids lead to increases in RDD frequency. Our results demonstrate that RDD is a conserved process that diversifies transcriptomes and proteomes and provide a mechanistic link between R-loops and RDDs.

A combined coarse-grained and all-atom simulation of TRPV1 channel gating and heat activation

Coarse-grained modeling and all-atom molecular dynamics simulation provide insight into the mechanism for heat activation of TRPV1 gating.

The transient receptor potential (TRP) channels act as key sensors of various chemical and physical stimuli in eukaryotic cells. Despite years of study, the molecular mechanisms of TRP channel activation remain unclear. To elucidate the structural, dynamic, and energetic basis of gating in TRPV1 (a founding member of the TRPV subfamily), we performed coarse-grained modeling and all-atom molecular dynamics (MD) simulation based on the recently solved high resolution structures of the open and closed form of TRPV1. Our coarse-grained normal mode analysis captures two key modes of collective motions involved in the TRPV1 gating transition, featuring a quaternary twist motion of the transmembrane domains (TMDs) relative to the intracellular domains (ICDs). Our transition pathway modeling predicts a sequence of structural movements that propagate from the ICDs to the TMDs via key interface domains (including the membrane proximal domain and the C-terminal domain), leading to sequential opening of the selectivity filter followed by the lower gate in the channel pore (confirmed by modeling conformational changes induced by the activation of ICDs). The above findings of coarse-grained modeling are robust to perturbation by lipids. Finally, our MD simulation of the ICD identifies key residues that contribute differently to the nonpolar energy of the open and closed state, and these residues are predicted to control the temperature sensitivity of TRPV1 gating. These computational predictions offer new insights to the mechanism for heat activation of TRPV1 gating, and will guide our future electrophysiology and mutagenesis studies.

Allosteric Modulation of Alpha7 Nicotinic Receptors: Mechanistic Insight through Metadynamics and Essential Dynamics

Increasing attention has recently been devoted to allosteric modulators, as they can provide inherent advantages over classic receptor agonists. In the field of nicotinic receptors (nAChRs), the main advantage is that allosteric modulators can trigger pharmacological responses, limiting receptor desensitization. Most of the known allosteric ligands are "positive allosteric modulators" (PAMs), which increase both sensitivity to receptor agonists and current amplitude. Intriguingly, some allosteric modulators are also able to activate the α7 receptor (α7-nAChR) even in the absence of orthosteric agonists. These compounds have been named "ago-allosteric modulators" and GAT107 has been studied in depth because of its unique mechanism of action. We here investigate by molecular dynamics simulations, metadynamics, and essential dynamics the activation mechanism of α7-nAChR, in the presence of different nicotinic modulators. We determine the free energy profiles associated with the closed-to-open motion of the loop C, and we highlight mechanistic differences observed in the presence of different modulators. In particular, we demonstrate that GAT107 triggers conformational motions and cross-talk similar to those observed when the α7-nACh receptor is in complex with both an agonist and an allosteric modulator.

Pathways and Barriers for Ion Translocation through the 5-HT3A Receptor Channel

Pentameric ligand gated ion channels (pLGICs) are ionotropic receptors that mediate fast intercellular communications at synaptic level and include either cation selective (e.g., nAChR and 5-HT₃) or anion selective (e.g., GlyR, GABA_A and GluCl) membrane channels. Among others, 5-HT₃ is one of the most studied members, since its first cloning back in 1991, and a large number of studies have successfully pinpointed protein residues critical for its activation and channel gating. In addition, 5-HT₃ is also the target of a few pharmacological treatments due to the demonstrated benefits of its modulation in clinical trials. Nonetheless, a detailed molecular analysis of important protein features, such as the origin of its ion selectivity and the rather low conductance as compared to other channel homologues, has been unfeasible until the recent crystallization of the mouse 5-HT₃A receptor. Here, we present extended molecular dynamics simulations and free energy calculations of the whole 5-HT₃A protein with the aim of better understanding its ion transport properties, such as the pathways for ion permeation into the receptor body and the complex nature of the selectivity filter. Our investigation unravels previously unpredicted structural features of the 5-HT₃A receptor, such as the existence of alternative intersubunit pathways for ion translocation at the interface between the extracellular and the transmembrane domains, in addition to the one along the channel main axis. Moreover, our study offers a molecular interpretation of the role played by an arginine triplet located in the intracellular domain on determining the characteristic low conductance of the 5-HT₃A receptor, as evidenced in previous experiments. In view of these results, possible implications on other members of the superfamily are suggested.

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

Nicotinic acetylcholine receptors (nAchRs) are ligand-gated ion channels that regulate chemical transmission at the neuromuscular junction. Structural information is available at low resolution from open and closed forms of an eukaryotic receptor, and at high resolution from other members of the same structural family, two prokaryotic orthologs and an eukaryotic GluCl channel. Structures of human channels however are still lacking. Homology modeling and Molecular Dynamics simulations are valuable tools to predict structures of unknown proteins, however, for the case of human nAchRs, they have been unsuccessful in providing a stable open structure so far. This is due to different problems with the template structures: on one side the homology with prokaryotic species is too low, while on the other the open eukaryotic GluCl proved itself unstable in several MD studies and collapsed to a dehydrated, non-conductive conformation, even when bound to an agonist. Aim of this work is to obtain, by a mixing of state-of-the-art homology and simulation techniques, a plausible prediction of the structure (still unknown) of the open state of human α7 nAChR complexed with epibatidine, from which it is possible to start structural and functional test studies. To prevent channel closure we employ a restraint that keeps the transmembrane pore open, and obtain in this way a stable, hydrated conformation. To further validate this conformation, we run four long, unbiased simulations starting from configurations chosen at random along the restrained trajectory. The channel remains stable and hydrated over the whole runs. This allows to assess the stability of the putative open conformation over a cumulative time of 1 μs, 800 ns of which are of unbiased simulation. Mostly based on the analysis of pore hydration and size, we suggest that the obtained structure has reasonable chances to be (at least one of the possible) structures of the channel in the open conformation.

The nicotinic acetylcholine receptor and its prokaryotic homologues: Structure, conformational transitions & allosteric modulation

Pentameric ligand-gated ion channels (pLGICs) play a central role in intercellular communications in the nervous system by converting the binding of a chemical messenger - a neurotransmitter - into an ion flux through the postsynaptic membrane. Here, we present an overview of the most recent advances on the signal transduction mechanism boosted by X-ray crystallography of both prokaryotic and eukaryotic homologues of the nicotinic acetylcholine receptor (nAChR) in conjunction with time-resolved analyses based on single-channel electrophysiology and Molecular Dynamics simulations. The available data consistently point to a global mechanism of gating that involves a large reorganization of the receptor mediated by two distinct quaternary transitions: a global twisting and a radial expansion/contraction of the extracellular domain. These transitions profoundly modify the organization of the interface between subunits, which host several sites for orthosteric and allosteric modulatory ligands. The same mechanism may thus mediate both positive and negative allosteric modulations of pLGICs ligand binding at topographically distinct sites. The emerging picture of signal transduction is expected to pave the way to new pharmacological strategies for the development of allosteric modulators of nAChR and pLGICs in general. This article is part of the Special Issue entitled 'The Nicotinic Acetylcholine Receptor: From Molecular Biology to Cognition'.

Protein dynamics and the allosteric transitions of pentameric receptor channels

The recent application of molecular dynamics (MD) methodology to investigate the allosteric transitions of the acetylcholine receptor and its prokaryotic and eukaryotic pentameric homologs has yielded new insights into the mechanisms of signal transduction by these receptors. Combined with available data on X-ray structures, MD techniques enable description of the dynamics of the conformational change at the atomic level, intra-molecular propagation of this signal transduction mechanism as a concerted stepwise process at physiological timescales and the control of this process by allosteric modulators, thereby offering new perspectives for drug design.

Constrained geometric simulation of the nicotinic acetylcholine receptor

Constrained geometric simulations have been performed for the recently published closed-channel state of the nicotinic acetylcholine receptor. These simulations support the theory that correlated motion in the flexible β-sheet structure of the extracellular domain helps to communicate a "conformational wave", spreading from the acetylcholine binding pocket. Furthermore, we have identified key residues that act at the interface between subunits and between domains that could potentially facilitate rapid communication between the binding site and the transmembrane gate.

Pentameric Ligand-gated Ion Channels : Insights from Computation

Pentameric ligand-gated ion channels (pLGICs) conduct upon the binding of an agonist and are fundamental to neurotransmission. New insights into the complex mechanisms underlying pLGIC gating, ion selectivity, and modulation have recently been gained via a series of crystal structures in prokaryotes and C .elegans, as well as computational studies relying on these structures. Here we review contributions from a variety of computational approaches, including normal mode analysis, automated docking, and fully atomistic molecular dynamics simulation. Examples from our own research, particularly concerning interactions with general anesthetics and lipids, are used to illustrate predictive results complementary to crystallographic studies.

Conformational changes in acetylcholine binding protein investigated by temperature accelerated molecular dynamics

Despite the large number of studies available on nicotinic acetylcholine receptors, a complete account of the mechanistic aspects of their gating transition in response to ligand binding still remains elusive. As a first step toward dissecting the transition mechanism by accelerated sampling techniques, we study the ligand-induced conformational changes of the acetylcholine binding protein (AChBP), a widely accepted model for the full receptor extracellular domain. Using unbiased Molecular Dynamics (MD) and Temperature Accelerated Molecular Dynamics (TAMD) simulations we investigate the AChBP transition between the apo and the agonist-bound state. In long standard MD simulations, both conformations of the native protein are stable, while the agonist-bound structure evolves toward the apo one if the orientation of few key sidechains in the orthosteric cavity is modified. Conversely, TAMD simulations initiated from the native conformations are able to produce the spontaneous transition. With respect to the modified conformations, TAMD accelerates the transition by at least a factor 10. The analysis of some specific residue-residue interactions points out that the transition mechanism is based on the disruption/formation of few key hydrogen bonds. Finally, while early events of ligand dissociation are observed already in standard MD, TAMD accelerates the ligand detachment and, at the highest TAMD effective temperature, it is able to produce a complete dissociation path in one AChBP subunit.

Acetylcholine promotes binding of α-conotoxin MII at α3 β2 nicotinic acetylcholine receptors

α-Conotoxin MII (α-CTxMII) is a 16-residue peptide with the sequence GCCSNPVCHLEHSNLC, containing Cys2-Cys8 and Cys3-Cys16 disulfide bonds. This peptide, isolated from the venom of the marine cone snail Conus magus, is a potent and selective antagonist of neuronal nicotinic acetylcholine receptors (nAChRs). To evaluate the impact of channel-ligand interactions on ligand-binding affinity, homology models of the heteropentameric α3β2-nAChR were constructed. The models were created in MODELLER with the aid of experimentally characterized structures of the Torpedo marmorata-nAChR (Tm-nAChR, PDB ID: 2BG9) and the Aplysia californica-acetylcholine binding protein (Ac-AChBP, PDB ID: 2BR8) as templates for the α3- and β2-subunit isoforms derived from rat neuronal nAChR primary amino acid sequences. Molecular docking calculations were performed with AutoDock to evaluate interactions of the heteropentameric nAChR homology models with the ligands acetylcholine (ACh) and α-CTxMII. The nAChR homology models described here bind ACh with binding energies commensurate with those of previously reported systems, and identify critical interactions that facilitate both ACh and α-CTxMII ligand binding. The docking calculations revealed an increased binding affinity of the α3β2-nAChR for α-CTxMII with ACh bound to the receptor, and this was confirmed through two-electrode voltage clamp experiments on oocytes from Xenopus laevis. These findings provide insights into the inhibition and mechanism of electrostatically driven antagonist properties of the α-CTxMIIs on nAChRs.

Allosteric regulation of pentameric ligand-gated ion channels: an emerging mechanistic perspective

Pentameric ligand-gated ion channels (pLGICs) play a central role in intercellular communications in the nervous system by converting the binding of a chemical messenger—a neurotransmitter—into an ion flux through the postsynaptic membrane. They are oligomeric assemblies that provide prototypical examples of allosterically regulated integral membrane proteins. Here, we present an overview of the most recent advances on the signal transduction mechanism based on the X-ray structures of both prokaryotic and invertebrate eukaryotic pLGICs and atomistic Molecular Dynamics simulations. The present results suggest that ion gating involves a large structural reorganization of the molecule mediated by two distinct quaternary transitions, a global twisting and the blooming of the extracellular domain, which can be modulated by ligand binding at the topographically distinct orthosteric and allosteric sites. The emerging model of gating is consistent with a wealth of functional studies and will boost the development of novel pharmacological strategies.

NMR structures of the human α7 nAChR transmembrane domain and associated anesthetic binding sites

The α7 nicotinic acetylcholine receptor (nAChR), assembled as homomeric pentameric ligand-gated ion channels, is one of the most abundant nAChR subtypes in the brain. Despite its importance in memory, learning and cognition, no structure has been determined for the α7 nAChR TM domain, a target for allosteric modulators. Using solution state NMR, we determined the structure of the human α7 nAChR TM domain (PDB ID: 2MAW) and demonstrated that the α7 TM domain formed functional channels in Xenopus oocytes. We identified the associated binding sites for the anesthetics halothane and ketamine; the former cannot sensitively inhibit α7 function, but the latter can. The α7 TM domain folds into the expected four-helical bundle motif, but the intra-subunit cavity at the extracellular end of the α7 TM domain is smaller than the equivalent cavity in the α4β2 nAChRs (PDB IDs: 2LLY; 2LM2). Neither drug binds to the extracellular end of the α7 TM domain, but two halothane molecules or one ketamine molecule binds to the intracellular end of the α7 TM domain. Halothane and ketamine binding sites are partially overlapped. Ketamine, but not halothane, perturbed the α7 channel-gate residue L9'. Furthermore, halothane did not induce profound dynamics changes in the α7 channel as observed in α4β2. The study offers a novel high-resolution structure for the human α7 nAChR TM domain that is invaluable for developing α7-specific therapeutics. It also provides evidence to support the hypothesis: only when anesthetic binding perturbs the channel pore or alters the channel motion, can binding generate functional consequences.

A gating mechanism of pentameric ligand-gated ion channels

Pentameric ligand-gated ion channels (pLGICs) play a central role in intercellular communication in the nervous system and are involved in fundamental processes such as attention, learning, and memory. They are oligomeric protein assemblies that convert a chemical signal into an ion flux through the postsynaptic membrane, but the molecular mechanism of gating ions has remained elusive. Here, we present atomistic molecular dynamics simulations of the prokaryotic channels from Gloeobacter violaceus (GLIC) and Erwinia chrysanthemi (ELIC), whose crystal structures are thought to represent the active and the resting states of pLGICs, respectively, and of the eukaryotic glutamate-gated chloride channel from Caenorhabditis elegans (GluCl), whose open-channel structure was determined complexed with the positive allosteric modulator ivermectin. Structural observables extracted from the trajectories of GLIC and ELIC are used as progress variables to analyze the time evolution of GluCl, which was simulated in the absence of ivermectin starting from the structure with bound ivermectin. The trajectory of GluCl with ivermectin removed shows a sequence of structural events that couple agonist unbinding from the extracellular domain to ion-pore closing in the transmembrane domain. Based on these results, we propose a structural mechanism for the allosteric communication leading to deactivation/activation of the GluCl channel. This model of gating emphasizes the coupling between the quaternary twisting and the opening/closing of the ion pore and is likely to apply to other members of the pLGIC family.

Structural analysis and ion translocation mechanisms of the muscle-type acetylcholine receptor channel

PURPOSE

The aim of this work is to analyze the conformational changes in the acetylcholine receptor caused by channel opening and to investigate the electrostatic profile during ion translocation through the channel.

METHODS

A computational model of the human muscle-type acetylcholine receptor (AChR) was built and used to analyze channel structure and its interactions with different ions. Using the Torpedo AChR crystal structure as a homologous template, the 3D structure of the human muscle-type AChR was reconstructed.

RESULTS

This first model is optimized and an open structure of the channel is generated using Normal Mode Analysis in order to assess morphologic and energetic differences between open and closed structures. In addition, the issue of ion translocation is investigated in further detail. Results elucidate different aspects of the channel: channel gate structure, channel interactions with translocating ions, differences between muscle-type AChR and previous neuronal-type AChR models.

CONCLUSIONS

The model constructed here is ideal for further computational studies on muscle-type AChR and its pathologic mutations.

Emerging approaches to probing ion channel structure and function

Ion channels, as membrane proteins, are the sensors of the cell. They act as the first line of communication with the world beyond the plasma membrane and transduce changes in the external and internal environments into unique electrical signals to shape the responses of excitable cells. Because of their importance in cellular communication, ion channels have been intensively studied at the structural and functional levels. Here, we summarize the diverse approaches, including molecular and cellular, chemical, optical, biophysical, and computational, used to probe the structural and functional rearrangements that occur during channel activation (or sensitization), inactivation (or desensitization), and various forms of modulation. The emerging insights into the structure and function of ion channels by multidisciplinary approaches allow the development of new pharmacotherapies as well as new tools useful in controlling cellular activity.

GABA binding to an insect GABA receptor: a molecular dynamics and mutagenesis study

RDL receptors are GABA-activated inhibitory Cys-loop receptors found throughout the insect CNS. They are a key target for insecticides. Here, we characterize the GABA binding site in RDL receptors using computational and electrophysiological techniques. A homology model of the extracellular domain of RDL was generated and GABA docked into the binding site. Molecular dynamics simulations predicted critical GABA binding interactions with aromatic residues F206, Y254, and Y109 and hydrophilic residues E204, S176, R111, R166, S176, and T251. These residues were mutated, expressed in Xenopus oocytes, and their functions assessed using electrophysiology. The data support the binding mechanism provided by the simulations, which predict that GABA forms many interactions with binding site residues, the most significant of which are cation-π interactions with F206 and Y254, H-bonds with E204, S205, R111, S176, T251, and ionic interactions with R111 and E204. These findings clarify the roles of a range of residues in binding GABA in the RDL receptor, and also show that molecular dynamics simulations are a useful tool to identify specific interactions in Cys-loop receptors.

End-plate acetylcholine receptor: structure, mechanism, pharmacology, and disease

The synapse is a localized neurohumoral contact between a neuron and an effector cell and may be considered the quantum of fast intercellular communication. Analogously, the postsynaptic neurotransmitter receptor may be considered the quantum of fast chemical to electrical transduction. Our understanding of postsynaptic receptors began to develop about a hundred years ago with the demonstration that electrical stimulation of the vagus nerve released acetylcholine and slowed the heart beat. During the past 50 years, advances in understanding postsynaptic receptors increased at a rapid pace, owing largely to studies of the acetylcholine receptor (AChR) at the motor endplate. The endplate AChR belongs to a large superfamily of neurotransmitter receptors, called Cys-loop receptors, and has served as an exemplar receptor for probing fundamental structures and mechanisms that underlie fast synaptic transmission in the central and peripheral nervous systems. Recent studies provide an increasingly detailed picture of the structure of the AChR and the symphony of molecular motions that underpin its remarkably fast and efficient chemoelectrical transduction.

Tightening of the ATP-binding sites induces the opening of P2X receptor channels

The opening of ligand-gated ion channels in response to agonist binding is a fundamental process in biology. In ATP-gated P2X receptors, little is known about the molecular events that couple ATP binding to channel opening. In this paper, we identify structural changes of the ATP site accompanying the P2X2 receptor activation by engineering extracellular zinc bridges at putative mobile regions as revealed by normal mode analysis. We provide evidence that tightening of the ATP sites shaped like open 'jaws' induces opening of the P2X ion channel. We show that ATP binding favours jaw tightening, whereas binding of a competitive antagonist prevents gating induced by this movement. Our data reveal the inherent dynamic of the binding jaw, and provide new structural insights into the mechanism of P2X receptor activation.

Gating mechanism of a P2X4 receptor developed from normal mode analysis and molecular dynamics simulations

P2X receptors are trimeric ATP-gated cation channels participating in diverse physiological processes. How ATP binding triggers channel opening remains unclear. Here the gating mechanism of a P2X receptor was studied by normal mode analysis and molecular dynamics (MD) simulations. Based on the resting-state crystal structure, a normal mode involving coupled motions of three β-strands (β1, β13, and β14) at the trimeric interface of the ligand-binding ectodomain and the pore-lining helix (TM2) in the transmembrane domain (TMD) was identified. The resulting widening of the fenestrations above the TMD and opening of the transmembrane pore produce known signatures of channel activation. In MD simulations, ATP was initially placed in the putative binding pocket (defined by four charged residues located in β1, β13 and β14) in two opposite orientations, with the adenine either proximal or distal to the TMD. In the proximal orientation, the triphosphate group extends outward to draw in the four charged residues, leading to closure of β13/β14 toward β1. The adenine ring, wedged between β1 and β13, acts as a fulcrum for the β14 lever, turning a modest closure around the triphosphate group into significant opening of the pre-TM2 loop. The motions of these β-strands are similar to those in the putative channel-activation normal mode. In the distal orientation, the ATP stabilizes the trimeric interface and the closure of the pre-TM2 loop, possibly representing desensitization. Our computational studies produced the first complete model, supported by experimental data, for how ATP binding triggers activation of a P2X receptor.

Rapid simulation of protein motion: merging flexibility, rigidity and normal mode analyses

Protein function frequently involves conformational changes with large amplitude on timescales which are difficult and computationally expensive to access using molecular dynamics. In this paper, we report on the combination of three computationally inexpensive simulation methods--normal mode analysis using the elastic network model, rigidity analysis using the pebble game algorithm, and geometric simulation of protein motion--to explore conformational change along normal mode eigenvectors. Using a combination of ElNemo and First/Froda software, large-amplitude motions in proteins with hundreds or thousands of residues can be rapidly explored within minutes using desktop computing resources. We apply the method to a representative set of six proteins covering a range of sizes and structural characteristics and show that the method identifies specific types of motion in each case and determines their amplitude limits.

Molecular dynamics

Molecular dynamics (MD) simulation holds the promise of revealing the mechanisms of biological processes in their ultimate detail. It is carried out by computing the interaction forces acting on each atom and then propagating the velocities and positions of the atoms by numerical integration of Newton's equations of motion. In this review, we present an overview of how the MD simulation can be conducted to address computational toxicity problems. The study cases will cover a standard MD simulation performed to investigate the overall flexibility of a cytochrome P450 (CYP) enzyme and a set of more advanced MD simulations to examine the barrier to ion conduction in a human α7 nicotinic acetylcholine receptor (nAChR).

The harmonic analysis of cylindrically symmetric proteins: a comparison of Dronpa and a DNA sliding clamp

The harmonic analysis of two types of proteins with cylindrical symmetry is performed by the Standard Force Field Normal Mode Analysis and by the elastic network model. For both proteins the global elastic modes are assigned to their characteristic topologies. Dronpa is a rigid β-barrel structure, presenting the twisting, bending and breathing motion of a cylindrical rod. The β sliding clamp of Escherichia coli is a hexagonal β-wheel, consisting of rigid segments. In its spectrum four classes of vibrations are identified which are characteristic of an elastic torus. Correlation diagrams and RMSF analysis are compared. The results provide not only a comprehensive validation of the use of both methods to describe the elastic behavior according to the low-frequency normal modes, but also depict the correlated motions of β-barrel and β-wheel proteins. The harmonic flexibility of the Dronpa protein is compared to the principal components of molecular dynamics (MD) simulation. A functionally important localized cleft opening mode is found, which is not detected by harmonic analysis.

Intramembrane proton binding site linked to activation of bacterial pentameric ion channel

Prokaryotic orthologs of eukaryotic Cys-loop receptor channels recently emerged as structural and mechanistic surrogates to investigate this superfamily of intercellular signaling proteins. Here, we examine proton activation of the prokaryotic ortholog GLIC using patch clamp electrophysiology, mutagenesis, and molecular dynamics (MD) simulations. Whole-cell current recordings from human embryonic kidney (HEK) 293 cells expressing GLIC show half-maximal activation at pH 6, close to the pK(a) of histidine, implicating the three native His residues in proton sensing linked to activation. The mutation H235F abolishes proton activation, H277Y is without effect, and all nine mutations of His-127 prevent expression on the cell surface. In the GLIC crystal structure, His-235 on transmembrane (TM) α-helix 2, hydrogen bonds to the main chain carbonyl oxygen of Ile-259 on TM α-helix 3. MD simulations show that when His-235 is protonated, the hydrogen bond persists, and the channel remains in the open conformation, whereas when His-235 is deprotonated, the hydrogen bond dissociates, and the channel closes. Mutations of the proximal Tyr-263, which also links TM α-helices 2 and 3 via a hydrogen bond, alter proton sensitivity over a 1.5 pH unit range. MD simulations show that mutations of Tyr-263 alter the hydrogen bonding capacity of His-235. The overall findings show that His-235 in the TM region of GLIC is a novel proton binding site linked to channel activation.

Competitive binding at a nicotinic receptor transmembrane site of two α7-selective positive allosteric modulators with differing effects on agonist-evoked desensitization

Positive allosteric modulators (PAMs) of nicotinic acetylcholine receptors (nAChRs) have attracted considerable interest as a novel area of therapeutic drug discovery. Two types of α7-selective PAMs have been identified (type I and type II). Whilst both potentiate peak agonist-induced responses, they have different effects on the rate of agonist-induced receptor desensitization. Type I PAMs have little or no effect on the rapid rate of desensitization that is characteristic of α7 nAChRs, whereas type II PAMs cause dramatic slowing of receptor desensitization. Previously, we have obtained evidence indicating that PNU-120596, a type II PAM, causes potentiation by interacting with an allosteric transmembrane site. In contrast, other studies have demonstrated the importance of the ‘M2–M3 segment’ in modulating the effects of the type I PAM NS1738 and have led to the proposal that NS1738 may interact with the extracellular N-terminal domain. Here, our aim has been to compare the mechanism of allosteric potentiation of α7 nAChRs by NS1738 and PNU-120596. Functional characterization of a series of mutated α7 nAChRs indicates that mutation of amino acids within a proposed intrasubunit transmembrane cavity have a broadly similar effect on these two PAMs. In addition, we have employed a functional assay designed to examine the ability of ligands to act competitively at either the orthosteric or allosteric binding site of α7 nAChRs. These data, together with computer docking simulations, lead us to conclude that both the type I PAM NS1738 and the type II PAM PNU-120596 bind competitively at a mutually exclusive intrasubunit transmembrane site.

Testing the applicability of Nernst-Planck theory in ion channels: comparisons with Brownian dynamics simulations

The macroscopic Nernst-Planck (NP) theory has often been used for predicting ion channel currents in recent years, but the validity of this theory at the microscopic scale has not been tested. In this study we systematically tested the ability of the NP theory to accurately predict channel currents by combining and comparing the results with those of Brownian dynamics (BD) simulations. To thoroughly test the theory in a range of situations, calculations were made in a series of simplified cylindrical channels with radii ranging from 3 to 15 Å, in a more complex 'catenary' channel, and in a realistic model of the mechanosensitive channel MscS. The extensive tests indicate that the NP equation is applicable in narrow ion channels provided that accurate concentrations and potentials can be input as the currents obtained from the combination of BD and NP match well with those obtained directly from BD simulations, although some discrepancies are seen when the ion concentrations are not radially uniform. This finding opens a door to utilising the results of microscopic simulations in continuum theory, something that is likely to be useful in the investigation of a range of biophysical and nano-scale applications and should stimulate further studies in this direction.

Agonist activation of alpha7 nicotinic acetylcholine receptors via an allosteric transmembrane site

Conventional nicotinic acetylcholine receptor (nAChR) agonists, such as acetylcholine, act at an extracellular "orthosteric" binding site located at the interface between two adjacent subunits. Here, we present evidence of potent activation of α7 nAChRs via an allosteric transmembrane site. Previous studies have identified a series of nAChR-positive allosteric modulators (PAMs) that lack agonist activity but are able to potentiate responses to orthosteric agonists, such as acetylcholine. It has been shown, for example, that TQS acts as a conventional α7 nAChR PAM. In contrast, we have found that a compound with close chemical similarity to TQS (4BP-TQS) is a potent allosteric agonist of α7 nAChRs. Whereas the α7 nAChR antagonist metyllycaconitine acts competitively with conventional nicotinic agonists, metyllycaconitine is a noncompetitive antagonist of 4BP-TQS. Mutation of an amino acid (M253L), located in a transmembrane cavity that has been proposed as being the binding site for PAMs, completely blocks agonist activation by 4BP-TQS. In contrast, this mutation had no significant effect on agonist activation by acetylcholine. Conversely, mutation of an amino acid located within the known orthosteric binding site (W148F) has a profound effect on agonist potency of acetylcholine (resulting in a shift of ∼200-fold in the acetylcholine dose-response curve), but had little effect on the agonist dose-response curve for 4BP-TQS. Computer docking studies with an α7 homology model provides evidence that both TQS and 4BP-TQS bind within an intrasubunit transmembrane cavity. Taken together, these findings provide evidence that agonist activation of nAChRs can occur via an allosteric transmembrane site.

Decrypting the sequence of structural events during the gating transition of pentameric ligand-gated ion channels based on an interpolated elastic network model

Despite many experimental and computational studies of the gating transition of pentameric ligand-gated ion channels (pLGICs), the structural basis of how ligand binding couples to channel gating remains unknown. By using a newly developed interpolated elastic network model (iENM), we have attempted to compute a likely transition pathway from the closed- to the open-channel conformation of pLGICs as captured by the crystal structures of two prokaryotic pLGICs. The iENM pathway predicts a sequence of structural events that begins at the ligand-binding loops and is followed by the displacements of two key loops (loop 2 and loop 7) at the interface between the extracellular and transmembrane domain, the tilting/bending of the pore-lining M2 helix, and subsequent movements of M4, M3 and M1 helices in the transmembrane domain. The predicted order of structural events is in broad agreement with the Φ-value analysis of α subunit of nicotinic acetylcholine receptor mutants, which supports a conserved core mechanism for ligand-gated channel opening in pLGICs. Further perturbation analysis has supported the critical role of certain intra-subunit and inter-subunit interactions in dictating the above sequence of events.

Relationships between structural dynamics and functional kinetics in oligomeric membrane receptors

Recent efforts to broaden understanding of the molecular mechanisms of membrane receptors in signal transduction make use of rate-equilibrium free-energy relationships (REFERs), previously applied to chemical reactions, enzyme kinetics, and protein folding. For oligomeric membrane receptors, we distinguish between a), the Leffler parameter α_L, to characterize the global transition state for the interconversion between conformations; and b), the Fersht parameter, ϕ_F, to assign the degree of progression of individual residue positions at the transition state. For both α_L and ϕ_F, insights are achieved by using harmonic energy profiles to reflect the dynamic nature of proteins, as illustrated with single-channel results reported for normal and mutant nicotinic receptors. We also describe new applications of α_L based on published results. For large-conductance calcium-activated potassium channels, data are satisfactorily fit with an α_L value of 0.65, in accord with REFERs. In contrast, results reported for the flip conformational state of glycine and nicotinic receptors are in disaccord with REFERs, since they yield α_L values outside the usual limits of 0–1. Concerning published ϕ_F values underlying the conformational wave hypothesis for nicotinic receptors, we note that interpretations may be complicated by variations in the width of harmonic energy profiles.

Unresponsive correlated motion in alpha7 nAChR to halothane binding explains its functional insensitivity to volatile anesthetics

Neuronal nicotinic acetylcholine receptors (nAChRs) have been implicated as targets for general anesthetics, but the functional responses to anesthetic modulation vary considerably among different subtypes of nAChRs. Inhaled general anesthetics, such as halothane, could effectively inhibit the channel activity of the alpha4beta2 nAChR but not the homologous alpha7 nAChR. To understand why alpha7 is insensitive to inhaled general anesthetics, we performed multiple sets of 20 ns molecular dynamics (MD) simulations on the closed- and open-channel alpha7 in the absence and presence of halothane and critically compared the results with those from our studies on the alpha4beta2 nAChR (Liu et al. J. Phys. Chem. B 2009, 113, 12581 and Liu et al. J. Phys. Chem. B 2010, 114, 626). Several halothane binding sites with fairly high binding affinities were identified in both closed- and open-channel alpha7, suggesting that a lack of sensitive functional responses of the alpha7 nAChR to halothane in the previous experiments was unlikely due to a lack of halothane interaction with alpha7. The binding affinities of halothane in alpha7 seemed to be protein conformation-dependent. Overall, halothane affinity was higher in the closed-channel alpha7. Halothane binding to alpha7 did not induce profound changes in alpha7 structure and dynamics that could be related to the channel function. In contrast, correlated motion of the open-channel alpha4beta2 was reduced substantially in the presence of halothane, primarily due to the more susceptible nature of beta2 to anesthetic modulation. The amphiphilic extracellular and transmembrane domain interface of the beta2 subunit is attractive to halothane and is susceptible to halothane perturbation, which may be why alpha4beta2 is functionally more sensitive to halothane than alpha7.

Normal mode gating motions of a ligand-gated ion channel persist in a fully hydrated lipid bilayer model

We have previously used molecular modeling and normal-mode analyses combined with experimental data to visualize a plausible model of a transmembrane ligand-gated ion channel. We also postulated how the gating motion of the channel may be affected by the presence of various ligands, especially anesthetics. As is typical for normal-mode analyses, those studies were performed in vacuo to reduce the computational complexity of the problem. While such calculations constitute an efficient way to model the large scale structural flexibility of transmembrane proteins, they can be criticized for neglecting the effects of an explicit phospholipid bilayer or hydrated environment. Here, we show the successful calculation of normal-mode motions for our model of a glycine α-1 receptor, now suspended in a fully hydrated lipid bilayer. Despite the almost uniform atomic density, the introduction of water and lipid does not grossly distort the overall gating motion. Normal-mode analysis revealed that even a fully immersed glycine α-1 receptor continues to demonstrate an iris-like channel gating motion as a low-frequency, high-amplitude natural harmonic vibration consistent with channel gating. Furthermore, the introduction of periodic boundary conditions allows the examination of simultaneous harmonic vibrations of lipid in synchrony with the protein gating motions that are compatible with reasonable lipid bilayer perturbations. While these perturbations tend to influence the overall protein motion, this work provides continued support for the iris-like motion model that characterizes gating within the family of ligand-gated ion channels.

Nicotinic acetylcholine receptor transmembrane mutations convert ivermectin from a positive to a negative allosteric modulator

Ivermectin is a macrocyclic lactone that acts as a positive allosteric modulator of alpha7 nicotinic acetylcholine receptors (nAChRs) but has no modulatory activity on 5-hydroxytryptamine (5-HT) type 3 (5-HT(3)) receptors. By examining the influence of ivermectin on subunit chimeras containing domains from the nAChR alpha7 subunit and the 5-HT3A subunit, we have concluded that the transmembrane domains play a critical role in influencing allosteric modulation by ivermectin. A series of mutations located within the alpha-helical transmembrane domains of the alpha7 subunit were examined, and seven were found to have significant effects on allosteric modulation by ivermectin. Four mutations (A225D, Q272V, T456Y, and C459Y) caused a significant reduction in the potency of ivermectin as an allosteric potentiator. Compared with wild-type alpha7 nAChRs, potentiation by ivermectin was reduced dramatically (by 89-97%) by these mutations. Somewhat unexpectedly, three mutations (S222M, M253L, and S276V located in TM1, TM2, and TM3) converted ivermectin from a positive allosteric modulator into an antagonist. Levels of inhibition of 56, 84, and 89% were observed on M253L, S276V, and S222M, respectively. Antagonism by ivermectin was insurmountable and had no effect on EC(50) of acetylcholine, indicating that it is acting noncompetitively. The seven mutations that influence allosteric modulation by ivermectin are located near a predicted intrasubunit transmembrane cavity. Computer docking simulations provide support for the hypothesis that ivermectin binds in close proximity to this cavity. We conclude that transmembrane mutations in alpha7 nAChRs are able to convert ivermectin from a positive to a negative allosteric modulator.

The role of Loop F in the activation of the GABA receptor

Functional studies of the ligand gated ion channel family (nicotinic acetylcholine, serotonin Type 3, glycine and GABA receptors) along with the crystal structure of the acetylcholine binding protein (AChBP) and molecular dynamics simulations of the nAChR structure have resulted in a structural model in which the agonist-binding pocket comprises six loops (A-F) contributed by adjacent subunits. It is presumed that the binding of agonist results in a local structural rearrangement that is then transduced to the gate, causing the pore to open. Efforts are underway to better define the specific roles of the six binding loops. Several studies have suggested Loop F may play a direct role in linking the structural rearrangement within the binding pocket to the gate, although other investigations have indicated Loop F may be crucial for locking the agonist molecule into the binding site. This review will focus on the controversy surrounding the role of Loop F during GABA receptor activation.