A gating mechanism of pentameric ligand-gated ion channels

Hidden complexity of α7 nicotinic acetylcholine receptor desensitization revealed by MD simulations and Markov state modeling

The α7 nicotinic acetylcholine receptor is a pentameric ligand-gated ion channel that plays an important role in neuronal signaling throughout the nervous system. Its implication in neurological disorders and inflammation has spurred the development of numerous compounds that enhance channel activation. However, the therapeutic potential of these compounds has been limited by the characteristically fast desensitization of the α7 receptor. Using recent high-resolution structures from cryo-EM, and all-atom molecular dynamic simulations augmented by Markov state modeling, here we explore the mechanism of α7 receptor desensitization and its implication on allosteric modulation. The results provide a precise characterization of the desensitization gate and illuminate the mechanism of ion-pore opening/closing with an agonist bound. In addition, the simulations reveal the existence of a short-lived, open-channel intermediate between the activated and desensitized states that rationalizes the paradoxical pharmacology of the L247T mutant and may be relevant to type-II allosteric modulation. This analysis provides an interpretation of the signal transduction mechanism and its regulation in α7 receptors.

Electrically silent mutants unravel the mechanism of binding-gating coupling in Cys-loop receptors

The transduction of extracellular chemical signals into intracellular events relies on the communication between neighboring domains of membrane receptors. In the particular case of Cys-loop receptor channels, five short stretches of amino acids, one per subunit, link the extracellular and transmembrane domains in such a way that the ion permeability of the latter and the affinity for neurotransmitters of the former become tied to each other. Here, using direct functional approaches, we set out to understand the molecular bases of this crucial interdependence through the characterization of total loss-of-current mutations at the interface between domains. Our results indicate that domain-domain proximity plays a previously unnoticed critical role inasmuch as inserting a single residue in each linker rendered the two domains independent of each other. In marked contrast, loss-of-current mutations that leave the linkers' length unaltered did not compromise the interdomain coupling, but rather, seemed to cause agonist-bound closed receptors to desensitize without appreciably opening.

A topological mechanism for robust and efficient global oscillations in biological networks

Long and stable timescales are often observed in complex biochemical networks, such as in emergent oscillations. How these robust dynamics persist remains unclear, given the many stochastic reactions and shorter time scales demonstrated by underlying components. We propose a topological model that produces long oscillations around the network boundary, reducing the system dynamics to a lower-dimensional current in a robust manner. Using this to model KaiC, which regulates the circadian rhythm in cyanobacteria, we compare the coherence of oscillations to that in other KaiC models. Our topological model localizes currents on the system edge, with an efficient regime of simultaneously increased precision and decreased cost. Further, we introduce a new predictor of coherence from the analysis of spectral gaps, and show that our model saturates a global thermodynamic bound. Our work presents a new mechanism and parsimonious description for robust emergent oscillations in complex biological networks.

Ligand-gated ion channels as potential biomarkers for ADT-mediated cognitive decline in prostate cancer patients

Men with prostate cancer are at increased risk of developing cognitive decline by the use of second-generation androgen signaling inhibitors. To date, reliable and sensitive biomarkers that could distinguish men at high risk of cognitive dysfunction under androgen deprivation therapy (ADT) have not been characterized. We used high-throughput transcriptional profiling utilizing human prostate cancer cell culture models mimicking ADT, biomarker selection using minimal common oncology data elements-cytoscape, and bioinformatic analyses employing Advaita® iPathwayGuide and DisGeNET for identification of disease-related gene associations. Validation analysis of genes was performed on brain neuronal and glial cells by quantitative real-time polymerase chain reaction assay. Our systematic analysis of androgen deprivation-associated genes involved multiple biological processes, including neuroactive ligand-receptor interaction, axon guidance, cytokine-cytokine receptor interaction, and metabolic and cancer signaling pathways. Genes associated with neuroreceptor ligand interaction, including gamma-aminobutyric acid (GABA) A and B receptors and nuclear core proteins, were identified as top upstream regulators. Functional enrichment and protein-protein interaction network analysis highlighted the role of ligand-gated ion channels (LGICs) and their receptors in cognitive dysfunction. Gene-disease association assigned forgetfulness, intellectual disability, visuospatial deficit, bipolar disorder, and other neurocognitive impairment with upregulation of type-1 angiotensin II receptor, brain-derived neurotrophic factor, GABA type B receptor subunit 2 (GABBR2), GABRA3, GABRA5, GABRB1, glycine receptor beta, glutamate ionotropic receptor N-methyl-D-aspartate receptor (NMDA) type subunit 1, glutamate ionotropic receptor NMDA type subunit 2D, 5-hydroxytryptamine receptor 1D, interferon beta 1, and nuclear receptor subfamily 3 group C member 1 as top differentially expressed genes. Validation studies of brain glial cells, neurons, and patients on ADT demonstrated the association of these genes with cognitive decline. Our findings highlight LGICs as potential biomarkers for ADT-mediated cognitive decline. Further validation of these biomarkers may lead to future practical clinical use.

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.

A release of local subunit conformational heterogeneity underlies gating in a muscle nicotinic acetylcholine receptor

Synaptic receptors respond to neurotransmitters by opening an ion channel across the post-synaptic membrane to elicit a cellular response. Here we use recent Torpedo acetylcholine receptor structures and functional measurements to delineate a key feature underlying allosteric communication between the agonist-binding extracellular and channel-gating transmembrane domains. Extensive mutagenesis at this inter-domain interface re-affirms a critical energetically coupled role for the principal α subunit β1-β2 and M2-M3 loops, with agonist binding re-positioning a key β1-β2 glutamate/valine to facilitate the outward motions of a conserved M2-M3 proline to open the channel gate. Notably, the analogous structures in non-α subunits adopt a locally active-like conformation in the apo state even though each L9' hydrophobic gate residue in each pore-lining M2 α-helix is closed. Agonist binding releases local conformational heterogeneity transitioning all five subunits into a conformationally symmetric open state. A release of conformational heterogeneity provides a framework for understanding allosteric communication in pentameric ligand-gated ion channels.

Phenols and GABAA receptors: from structure and molecular mechanisms action to neuropsychiatric sequelae

γ-Aminobutyric acid type A receptors (GABAARs) are members of the pentameric ligand-gated ion channel (pLGIC) family, which are widespread throughout the invertebrate and vertebrate central nervous system. GABAARs are engaged in short-term changes of the neuronal concentrations of chloride (Cl-) and bicarbonate (HCO3 -) ions by their passive permeability through the ion channel pore. GABAARs are regulated by various structurally diverse phenolic substances ranging from simple phenols to complex polyphenols. The wide chemical and structural variability of phenols suggest similar and different binding sites on GABAARs, allowing them to manifest themselves as activators, inhibitors, or allosteric ligands of GABAAR function. Interest in phenols is associated with their great potential for GABAAR modulation, but also with their subsequent negative or positive role in neurological and psychiatric disorders. This review focuses on the GABAergic deficit hypotheses during neurological and psychiatric disorders induced by various phenols. We summarize the structure-activity relationship of general phenol groups concerning their differential roles in the manifestation of neuropsychiatric symptoms. We describe and analyze the role of GABAAR subunits in manifesting various neuropathologies and the molecular mechanisms underlying their modulation by phenols. Finally, we discuss how phenol drugs can modulate GABAAR activity via desensitization and resensitization. We also demonstrate a novel pharmacological approach to treat neuropsychiatric disorders via regulation of receptor phosphorylation/dephosphorylation.

Developing an Implicit Solvation Machine Learning Model for Molecular Simulations of Ionic Media

Molecular dynamics (MD) simulations of biophysical systems require accurate modeling of their native environment, i.e., aqueous ionic solution, as it critically impacts the structure and function of biomolecules. On the other hand, the models should be computationally efficient to enable simulations of large spatiotemporal scales. Here, we present the deep implicit solvation model for sodium chloride solutions that satisfies both requirements. Owing to the use of the neural network potential, the model can capture the many-body potential of mean force, while the implicit water treatment renders the model inexpensive. We demonstrate our approach first for pure ionic solutions with concentrations ranging from physiological to 2 M. We then extend the model to capture the effective ion interactions in the vicinity and far away from a DNA molecule. In both cases, the structural properties are in good agreement with all-atom MD, showcasing a general methodology for the efficient and accurate modeling of ionic media.

Modular Gating of Ion Transport by Postsynthetic Charge Transfer Complexation in a Metal-Organic Framework

Nature's design of biological ion channels that demonstrates efficient gating and selectivity brings to light a very promising model to mimic and design for achieving selective and tunable ion transport. Functionalized nanopores that permit modulation of the pore wall charges are a compelling approach to gain control over the ion transport mechanism through the pores. This makes way for employing a noncovalent supramolecular approach for attaining charge reversal of the MOF pore walls using donor-acceptor pairs that can demonstrate strong charge transfer interactions. Herein, robust Zr4+-based mesoporous MOF-808 was postsynthetically modified into an anion-selective nanochannel (MOF-808-MV) by modification with dicationic viologen-based motifs. Charge modulation and even reversal of the MOF-808-MV pore walls were then explored taking advantage of strong charge transfer interactions between the grafted dicationic viologen acceptor moieties and anionic, π-electron-rich donor guest molecules such as pyranine (PYR) and tetrathiafulvalene tetrabenzoic acid (TTF-TA). Tunability of the MOF pore charge from positive to neutral to negative was achieved via simple methodologies such as diffusion control in case of guest molecule like PYR and by pH modulation for pH-responsive guest like TTF-TA. This results in a concomitant modulation in the selectivity of the nanochannel, rendering it from anion-selective to ambipolar to cation-selective. Furthermore, as a real-time application of this ion channel, Na+ ion conductivity (σ = 3.5 × 10-5 S cm-1) was studied at ambient temperature.

Structural insights into opposing actions of neurosteroids on GABAA receptors

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

Conformational transitions and allosteric modulation in a heteromeric glycine receptor

Glycine Receptors (GlyRs) provide inhibitory neuronal input in the spinal cord and brainstem, which is critical for muscle coordination and sensory perception. Synaptic GlyRs are a heteromeric assembly of α and β subunits. Here we present cryo-EM structures of full-length zebrafish α1βBGlyR in the presence of an antagonist (strychnine), agonist (glycine), or agonist with a positive allosteric modulator (glycine/ivermectin). Each structure shows a distinct pore conformation with varying degrees of asymmetry. Molecular dynamic simulations found the structures were in a closed (strychnine) and desensitized states (glycine and glycine/ivermectin). Ivermectin binds at all five interfaces, but in a distinct binding pose at the β-α interface. Subunit-specific features were sufficient to solve structures without a fiduciary marker and to confirm the 4α:1β stoichiometry recently observed. We also report features of the extracellular and intracellular domains. Together, our results show distinct compositional and conformational properties of α1βGlyR and provide a framework for further study of this physiologically important channel.

Allosteric modulation of α1β3γ2 GABAA receptors by farnesol through the neurosteroid sites

In mammalian cells, all-trans farnesol, a 15-carbon isoprenol, is a product of the mevalonate pathway. It is the natural substrate of alcohol dehydrogenase and a substrate for CYP2E1, two enzymes implicated in ethanol metabolism. Studies have shown that farnesol is present in the human brain and inhibits voltage-gated Ca2+ channels at much lower concentrations than ethanol. Here we show that farnesol modulates the activity of γ-aminobutyric acid type A receptors (GABAARs), some of which also mediate the sedative activity of ethanol. Electrophysiology experiments performed in HEK cells expressing human α1β3γ2 or α6β3γ2 GABAARs revealed that farnesol increased chloride currents through positive allosteric modulation of these receptors and showed dependence on both the alcoholic functional group of farnesol and the length of the alkyl chain for activity. In silico studies using long-timescale unbiased all-atom molecular dynamics (MD) simulations of the human α1β3γ2 GABAA receptors revealed that farnesol modulates the channel by directly binding to the transmembrane neurosteroid-binding site, after partitioning into the surrounding membrane and reaching the receptor by lateral diffusion. Channel activation by farnesol was further characterized by several structural and dynamic variables, such as global twisting of the receptor's extracellular domain, tilting of the transmembrane M2 helices, radius, cross-sectional area, hydration status, and electrostatic potential of the channel pore. Our results expand the pharmacological activities of farnesol to yet another class of ion channels implicated in neurotransmission, thus providing a novel path for understanding and treatment of diseases involving GABAA receptor dysfunction.

Phylogenetic analyses of 5-hydroxytryptamine 3 (5-HT3) receptors in Metazoa

The 5-hydroxytrptamine 3 (5-HT3) receptor is a member of the 'Cys-loop' family and the only pentameric ligand gated ion channel among the serotonin receptors. 5-HT3 receptors play an important role in controlling growth, development, and behaviour in animals. Several 5-HT3 receptor antagonists are used to treat diseases (e.g., irritable bowel syndrome, nausea and emesis). Humans express five different subunits (A-E) enabling a variety of heteromeric receptors to form but all contain 5HT3A subunits. However, the information available about the 5-HT3 receptor subunit occurrence among the metazoan lineages is minimal. In the present article we searched for 5-HT3 receptor subunit homologs from different phyla in Metazoa. We identified more than 1000 5-HT3 receptor subunits in Metazoa in different phyla and undertook simultaneous phylogenetic analysis of 526 5HT3A, 358 5HT3B, 239 5HT3C, 70 5HT3D, and 173 5HT3E sequences. 5-HT3 receptor subunits were present in species belonging to 11 phyla: Annelida, Arthropoda, Chordata, Cnidaria, Echinodermata, Mollusca, Nematoda, Orthonectida, Platyhelminthes, Rotifera and Tardigrada. All subunits were most often identified in Chordata phylum which was strongly represented in searches. Using multiple sequence alignment, we investigated variations in the ligand binding region of the 5HT3A subunit protein sequences in the metazoan lineage. Several critical amino acid residues important for ligand binding (common structural features) are commonly present in species from Nematoda and Platyhelminth gut parasites through to Chordata. Collectively, this better understanding of the 5-HT3 receptor evolutionary patterns raises possibilities of future pharmacological challenges facing Metazoa including effects on parasitic and other species in ecosystems that contain 5-HT3 receptor ligands.

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.

Illumination of a progressive allosteric mechanism mediating the glycine receptor activation

Pentameric ligand-gated ion channel mediate signal transduction at chemical synapses by transiting between resting and open states upon neurotransmitter binding. Here, we investigate the gating mechanism of the glycine receptor fluorescently labeled at the extracellular-transmembrane interface by voltage-clamp fluorometry (VCF). Fluorescence reports a glycine-elicited conformational change that precedes pore opening. Low concentrations of glycine, partial agonists or specific mixtures of glycine and strychnine trigger the full fluorescence signal while weakly activating the channel. Molecular dynamic simulations of a partial agonist bound-closed Cryo-EM structure show a highly dynamic nature: a marked structural flexibility at both the extracellular-transmembrane interface and the orthosteric site, generating docking properties that recapitulate VCF data. This work illuminates a progressive propagating transition towards channel opening, highlighting structural plasticity within the mechanism of action of allosteric effectors.

Wordom update 2: A user-friendly program for the analysis of molecular structures and conformational ensembles

We present the second update of Wordom, a user-friendly and efficient program for manipulation and analysis of conformational ensembles from molecular simulations. The actual update expands some of the existing modules and adds 21 new modules to the update 1 published in 2011. The new adds can be divided into three sets that: 1) analyze atomic fluctuations and structural communication; 2) explore ion-channel conformational dynamics and ionic translocation; and 3) compute geometrical indices of structural deformation. Set 1 serves to compute correlations of motions, find geometrically stable domains, identify a dynamically invariant core, find changes in domain-domain separation and mutual orientation, perform wavelet analysis of large-scale simulations, process the output of principal component analysis of atomic fluctuations, perform functional mode analysis, infer regions of mechanical rigidity, analyze overall fluctuations, and perform the perturbation response scanning. Set 2 includes modules specific for ion channels, which serve to monitor the pore radius as well as water or ion fluxes, and measure functional collective motions like receptor twisting or tilting angles. Finally, set 3 includes tools to monitor structural deformations by computing angles, perimeter, area, volume, β-sheet curvature, radial distribution function, and center of mass. The ring perception module is also included, helpful to monitor supramolecular self-assemblies. This update places Wordom among the most suitable, complete, user-friendly, and efficient software for the analysis of biomolecular simulations. The source code of Wordom and the relative documentation are available under the GNU general public license at http://wordom.sf.net.

α1 Proline 277 Residues Regulate GABAAR Gating through M2-M3 Loop Interaction in the Interface Region

Cys-loop receptors are a superfamily of transmembrane, pentameric receptors that play a crucial role in mammalian CNS signaling. Physiological activation of these receptors is typically initiated by neurotransmitter binding to the orthosteric binding site, located at the extracellular domain (ECD), which leads to the opening of the channel pore (gate) at the transmembrane domain (TMD). Whereas considerable knowledge on molecular mechanisms of Cys-loop receptor activation was gathered for the acetylcholine receptor, little is known with this respect about the GABA_A receptor (GABA_AR), which mediates cellular inhibition. Importantly, several static structures of GABA_AR were recently described, paving the way to more in-depth molecular functional studies. Moreover, it has been pointed out that the TMD-ECD interface region plays a crucial role in transduction of conformational changes from the ligand binding site to the channel gate. One of the interface structures implicated in this transduction process is the M2-M3 loop with a highly conserved proline (P277) residue. To address this issue specifically for α₁β₂γ_2L GABA_AR, we choose to substitute proline α₁P277 with amino acids with different physicochemical features such as electrostatic charge or their ability to change the loop flexibility. To address the functional impact of these mutations, we performed macroscopic and single-channel patch-clamp analyses together with modeling. Our findings revealed that mutation of α₁P277 weakly affected agonist binding but was critical for all transitions of GABA_AR gating: opening/closing, preactivation, and desensitization. In conclusion, we provide evidence that conservative α₁P277 at the interface is strongly involved in regulating the receptor gating.

Exploring the Conformational Impact of Glycine Receptor TM1-2 Mutations Through Coarse-Grained Analysis and Atomistic Simulations

Pentameric ligand-gated ion channels (PLGICs) are a family of proteins that convert chemical signals into ion fluxes through cellular membranes. Their structures are highly conserved across all kingdoms from bacteria to eukaryotes. Beyond their classical roles in neurotransmission and neurological disorders, PLGICs have been recently related to cell proliferation and cancer. Here, we focus on the best characterized eukaryotic channel, the glycine receptor (GlyR), to investigate its mutational patterns in genomic-wide tumor screens and compare them with mutations linked to hyperekplexia (HPX), a Mendelian neuromotor disease that disrupts glycinergic currents. Our analysis highlights that cancer mutations significantly accumulate across TM1 and TM2, partially overlapping with HPX changes. Based on 3D-clustering, conservation, and phenotypic data, we select three mutations near the pore, expected to impact GlyR conformation, for further study by molecular dynamics (MD). Using principal components from experimental GlyR ensembles as framework, we explore the motions involved in transitions from the human closed and desensitized structures and how they are perturbed by mutations. Our MD simulations show that WT GlyR spontaneously explores opening and re-sensitization transitions that are significantly impaired by mutations, resulting in receptors with altered permeability and desensitization properties in agreement with HPX functional data.

Supramolecular Mechanosensitive Potassium Channel Formed by Fluorinated Amphiphilic Cyclophane

Inspired by mechanosensitive potassium channels found in nature, we developed a fluorinated amphiphilic cyclophane composed of fluorinated rigid aromatic units connected via flexible hydrophilic octa(ethylene glycol) chains. Microscopic and emission spectroscopic studies revealed that the cyclophane could be incorporated into the hydrophobic layer of the lipid bilayer membranes and self-assembled to form a supramolecular transmembrane ion channel. Current recording measurements using cyclophane-containing planer lipid bilayer membranes successfully demonstrated an efficient transmembrane ion transport. We also demonstrated that the ion transport property was sensitive to the mechanical forces applied to the membranes. In addition, ion transport assays using pH-sensitive fluorescence dye revealed that the supramolecular channel possesses potassium ion selectivity. We also performed all-atom hybrid quantum-mechanical/molecular mechanical simulations to assess the channel structures at atomic resolution and the mechanism of selective potassium ion transport. This research demonstrated the first example of a synthetic mechanosensitive potassium channel, which would open a new door to sensing and manipulating biologically important processes and purification of key materials in industries.

Database of glutamate-gated chloride (GluCl) subunits across 125 nematode species: patterns of gene accretion and sequence diversification

Glutamate-gated chloride channels belong to the Cys-loop receptor superfamily. Glutamate-gated chloride channels are activated by glutamate and form substrates for the antiparasitic drugs from the avermectin family. Glutamate-gated chloride channels are pentameric, and each subunit contains an N-terminal extracellular domain that binds glutamate and 4 helical transmembrane domains, which contain binding sites for avermectin drugs. In order to provide more insight into phylum-wide patterns of glutamate-gated chloride subunit gene expansion and sequence diversity across nematodes, we have developed a database of predicted glutamate-gated chloride subunit genes from 125 nematode species. Our analysis into this dataset described assorted patterns of species-specific glutamate-gated chloride gene counts across different nematodes as well as sequence diversity in key residues thought to be involved in avermectin binding.

Antagonistic Mechanism of α-Conotoxin BuIA toward the Human α3β2 Nicotinic Acetylcholine Receptor

Nicotinic acetylcholine receptors (nAChRs) are pentameric ligand-gated ion channels that are abundantly expressed in the central and peripheral nervous systems, playing an important role in mediating neurotransmitter release and inter-synaptic signaling. Dysfunctional nAChRs are associated with neurological disorders, and studying the structure and function of nAChRs is essential for development of drugs or strategies for treatment of related diseases. α-Conotoxins are selective antagonists of the nAChR and are an important class of drug leads. So far, the antagonistic mechanism of α-conotoxins toward the nAChRs is still unclear. In this study, we built an α3β2 nAChR homology model and investigated its conformational transition mechanism upon binding with a highly potent inhibitor, α-conotoxin BuIA, through μs molecular dynamic simulations and site-directed mutagenesis studies. The results suggested that the α3β2 nAChR underwent global conformational transitions and was stabilized into a closed state with three hydrophobic gates present in the transmembrane domain by BuIA. Finally, the probable antagonistic mechanism of BuIA was proposed. Overall, the closed-state model of the α3β2 nAChR bound with BuIA is not only essential for understanding the antagonistic mechanism of α-conotoxins but also particularly valuable for development of therapeutic inhibitors in future.

Towards the Idea of Molecular Brains

How can single cells without nervous systems perform complex behaviours such as habituation, associative learning and decision making, which are considered the hallmark of animals with a brain? Are there molecular systems that underlie cognitive properties equivalent to those of the brain? This review follows the development of the idea of molecular brains from Darwin’s “root brain hypothesis”, through bacterial chemotaxis, to the recent discovery of neuron-like r-protein networks in the ribosome. By combining a structural biology view with a Bayesian brain approach, this review explores the evolutionary labyrinth of information processing systems across scales. Ribosomal protein networks open a window into what were probably the earliest signalling systems to emerge before the radiation of the three kingdoms. While ribosomal networks are characterised by long-lasting interactions between their protein nodes, cell signalling networks are essentially based on transient interactions. As a corollary, while signals propagated in persistent networks may be ephemeral, networks whose interactions are transient constrain signals diffusing into the cytoplasm to be durable in time, such as post-translational modifications of proteins or second messenger synthesis. The duration and nature of the signals, in turn, implies different mechanisms for the integration of multiple signals and decision making. Evolution then reinvented networks with persistent interactions with the development of nervous systems in metazoans. Ribosomal protein networks and simple nervous systems display architectural and functional analogies whose comparison could suggest scale invariance in information processing. At the molecular level, the significant complexification of eukaryotic ribosomal protein networks is associated with a burst in the acquisition of new conserved aromatic amino acids. Knowing that aromatic residues play a critical role in allosteric receptors and channels, this observation suggests a general role of π systems and their interactions with charged amino acids in multiple signal integration and information processing. We think that these findings may provide the molecular basis for designing future computers with organic processors.

Markov state models of proton- and pore-dependent activation in a pentameric ligand-gated ion channel

Ligand-gated ion channels conduct currents in response to chemical stimuli, mediating electrochemical signaling in neurons and other excitable cells. For many channels, the details of gating remain unclear, partly due to limited structural data and simulation timescales. Here, we used enhanced sampling to simulate the pH-gated channel GLIC, and construct Markov state models (MSMs) of gating. Consistent with new functional recordings, we report in oocytes, our analysis revealed differential effects of protonation and mutation on free-energy wells. Clustering of closed- versus open-like states enabled estimation of open probabilities and transition rates, while higher-order clustering affirmed conformational trends in gating. Furthermore, our models uncovered state- and protonation-dependent symmetrization. This demonstrates the applicability of MSMs to map energetic and conformational transitions between ion-channel functional states, and how they reproduce shifts upon activation or mutation, with implications for modeling neuronal function and developing state-selective drugs.

Nicotinic receptors: From protein allostery to computational neuropharmacology

We propose an extension and further development of the Monod-Wyman-Changeux model for allosteric transitions of regulatory proteins to brain communications and specifically to neurotransmitters receptors, with the nicotinic acetylcholine receptor (nAChR) as a model of ligand-gated ion channels. The present development offers an expression of the change of the gating isomerization constant caused by pharmacological ligand binding in terms of its value in the absence of ligands and several "modulation factors", which vary with orthosteric ligand binding (agonists/antagonists), allosteric ligand binding (positive allosteric modulators/negative allosteric modulators) and receptor desensitization. The new - explicit - formulation of such "modulation factors", provides expressions for the pharmacological attributes of potency, efficacy, and selectivity for the modulatory ligands (including endogenous neurotransmitters) in terms of their binding affinity for the active, resting, and desensitized states of the receptor. The current formulation provides ways to design neuroactive compounds with a controlled pharmacological profile, opening the field of computational neuro-pharmacology.

Mutational analysis to explore long-range allosteric couplings involved in a pentameric channel receptor pre-activation and activation

Pentameric ligand-gated ion channels (pLGICs) mediate chemical signaling through a succession of allosteric transitions that are yet not completely understood as intermediate states remain poorly characterized by structural approaches. In a previous study on the prototypic bacterial proton-gated channel GLIC, we generated several fluorescent sensors of the protein conformation that report a fast transition to a pre-active state, which precedes the slower process of activation with pore opening. Here, we explored the phenotype of a series of allosteric mutations, using simultaneous steady-state fluorescence and electrophysiological measurements over a broad pH range. Our data, fitted to a three-state Monod-Wyman-Changeux model, show that mutations at the subunit interface in the extracellular domain (ECD) principally alter pre-activation, while mutations in the lower ECD and in the transmembrane domain principally alter activation. We also show that propofol alters both transitions. Data are discussed in the framework of transition pathways generated by normal mode analysis (iModFit). It further supports that pre-activation involves major quaternary compaction of the ECD, and suggests that activation involves principally a reorganization of a ‘central gating region’ involving a contraction of the ECD β-sandwich and the tilt of the channel lining M2 helix.

Advances in hyperekplexia and other startle syndromes

Startle, a basic alerting reaction common to all mammals, is described as a sudden involuntary movement of the body evoked by all kinds of sudden and unexpected stimulus. Startle syndromes are heterogeneous groups of disorders with abnormal and exaggerated responses to startling events, including hyperekplexia, stimulus-induced disorders, and neuropsychiatric startle syndromes. Hyperekplexia can be attributed to a genetic, idiopathic, or symptomatic cause. Excluding secondary factors, hereditary hyperekplexia, a rare neurogenetic disorder with highly genetic heterogeneity, is characterized by neonatal hypertonia, exaggerated startle response provoked by the sudden external stimuli, and followed by a short period of general stiffness. It mainly arises from defects of inhibitory glycinergic neurotransmission. GLRA1 is the major pathogenic gene of hereditary hyperekplexia, along with many other genes involved in the function of glycinergic inhibitory synapses. While about 40% of patients remain negative genetic findings. Clonazepam, which can specifically upgrade the GABARA1 chloride channels, is the main and most effective administration for hereditary hyperekplexia patients. In this review, with the aim at enhancing the recognition and prompting potential treatment for hyperekplexia, we focused on discussing the advances in hereditary hyperekplexia genetics and the expound progress in pathogenic mechanisms of the glycinergic-synapse-related pathway and then followed by a brief overview of other common startle syndromes.

Biomimetic calcium-inactivated ion/molecular channel

A phosphopeptide-modified nanochannel was prepared based on a conical polymeric nanopore. It shows a reversible Ca2+-induced inactivation effect toward the ion flow and molecular transport, resulting from Ca2+ binding-caused surface charge neutralization and hydrophilicity reduction, and Ca2+ removal by the competitive binding.

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.

Is the Antidepressant Activity of Selective Serotonin Reuptake Inhibitors Mediated by Nicotinic Acetylcholine Receptors?

It is generally assumed that selective serotonin reuptake inhibitors (SSRIs) induce antidepressant activity by inhibiting serotonin (5-HT) reuptake transporters, thus elevating synaptic 5-HT levels and, finally, ameliorates depression symptoms. New evidence indicates that SSRIs may also modulate other neurotransmitter systems by inhibiting neuronal nicotinic acetylcholine receptors (nAChRs), which are recognized as important in mood regulation. There is a clear and strong association between major depression and smoking, where depressed patients smoke twice as much as the normal population. However, SSRIs are not efficient for smoking cessation therapy. In patients with major depressive disorder, there is a lower availability of functional nAChRs, although their amount is not altered, which is possibly caused by higher endogenous ACh levels, which consequently induce nAChR desensitization. Other neurotransmitter systems have also emerged as possible targets for SSRIs. Studies on dorsal raphe nucleus serotoninergic neurons support the concept that SSRI-induced nAChR inhibition decreases the glutamatergic hyperstimulation observed in stress conditions, which compensates the excessive 5-HT overflow in these neurons and, consequently, ameliorates depression symptoms. At the molecular level, SSRIs inhibit different nAChR subtypes by noncompetitive mechanisms, including ion channel blockade and induction of receptor desensitization, whereas α9α10 nAChRs, which are peripherally expressed and not directly involved in depression, are inhibited by competitive mechanisms. According to the functional and structural results, SSRIs bind within the nAChR ion channel at high-affinity sites that are spread out between serine and valine rings. In conclusion, SSRI-induced inhibition of a variety of nAChRs expressed in different neurotransmitter systems widens the complexity by which these antidepressants may act clinically.

Filling the Gap in Understanding the Mechanism of GABAAR and Propofol Using Computational Approaches

γ-Aminobutyric acid type-A receptors (GABA_ARs) play a critical role in neural transmission by mediating the inhibitory neural firing and are the target of many psychiatric drugs. Among them, propofol is one of the most widely used and important general anesthetics in clinics. Recent advances in structural biology revealed the structure of a human GABA_AR in both open and closed states. Yet, the detailed mechanism of the receptor and propofol remains to be fully understood. Therefore, in this study, based on the previous successes in structural biology, a variety of computational techniques were applied to fill the gap between previous experimental studies. This study investigated the ion-conducting mechanism of GABA_AR, predicted the possible binding mechanism of propofol, and revealed a new motion mechanism of transmembrane domain (TMD) helices. We hope that this study may contribute to future studies on ion-channel receptors, general anesthetics, and drug development.

Signal transduction through Cys-loop receptors is mediated by the nonspecific bumping of closely apposed domains

One of the most fundamental questions in the field of Cys-loop receptors (pentameric ligand-gated ion channels, pLGICs) is how the affinity for neurotransmitters and the conductive/nonconductive state of the transmembrane pore are correlated despite the ∼60-Å distance between the corresponding domains. Proposed mechanisms differ, but they all converge into the idea that interactions between wild-type side chains across the extracellular-transmembrane-domain (ECD-TMD) interface are crucial for this phenomenon. Indeed, the successful design of fully functional chimeras that combine intact ECD and TMD modules from different wild-type pLGICs has commonly been ascribed to the residual conservation of sequence that exists at the level of the interfacial loops even between evolutionarily distant parent channels. Here, using mutagenesis, patch-clamp electrophysiology, and radiolabeled-ligand binding experiments, we studied the effect of eliminating this residual conservation of sequence on ion-channel function and cell-surface expression. From our results, we conclude that proper state interconversion ("gating") does not require conservation of sequence-or even physicochemical properties-across the ECD-TMD interface. Wild-type ECD and TMD side chains undoubtedly interact with their surroundings, but the interactions between them-straddling the interface-do not seem to be more important for gating than those occurring elsewhere in the protein. We propose that gating of pLGICs requires, instead, that the overall structure of the interfacial loops be conserved, and that their relative orientation and distance be the appropriate ones for changes in one side to result in changes in the other, in a phenomenon akin to the nonspecific "bumping" of closely apposed domains.

Molecular Simulations of Hydrophobic Gating of Pentameric Ligand Gated Ion Channels: Insights into Water and Ions

Ion channels are proteins which form gated nanopores in biological membranes. Many channels exhibit hydrophobic gating, whereby functional closure of a pore occurs by local dewetting. The pentameric ligand gated ion channels (pLGICs) provide a biologically important example of hydrophobic gating. Molecular simulation studies comparing additive vs polarizable models indicate predictions of hydrophobic gating are robust to the model employed. However, polarizable models suggest favorable interactions of hydrophobic pore-lining regions with chloride ions, of relevance to both synthetic carriers and channel proteins. Electrowetting of a closed pLGIC hydrophobic gate requires too high a voltage to occur physiologically but may inform designs for switchable nanopores. Global analysis of ∼200 channels yields a simple heuristic for structure-based prediction of (closed) hydrophobic gates. Simulation-based analysis is shown to provide an aid to interpretation of functional states of new channel structures. These studies indicate the importance of understanding the behavior of water and ions within the nanoconfined environment presented by ion channels.

Cholesterol content in the membrane promotes key lipid-protein interactions in a pentameric serotonin-gated ion channel

Pentameric ligand-gated ion channels (pLGICs), embedded in the lipid membranes of nerve cells, mediate fast synaptic transmission and are major pharmaceutical targets. Because of their complexity and the limited knowledge of their structure, their working mechanisms have still to be fully unraveled at the molecular level. Over the past few years, evidence that the lipid membrane may modulate the function of membrane proteins, including pLGICs, has emerged. Here, we investigate, by means of molecular dynamics simulations, the behavior of the lipid membrane at the interface with the 5-HT_3A receptor (5-HT_3AR), a representative pLGIC which is the target of nausea-suppressant drugs, in a nonconductive state. Three lipid compositions are studied, spanning different concentrations of the phospholipids, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine, and of cholesterol, hence a range of viscosities. A variety of lipid interactions and persistent binding events to different parts of the receptor are revealed in the investigated models, providing snapshots of the dynamical environment at the membrane-receptor interface. Some of these events result in lipid intercalation within the transmembrane domain, and others reach out to protein key sections for signal transmission and receptor activation, such as the Cys-loop and the M2-M3 loop. In particular, phospholipids, with their long hydrophobic tails, play an important role in these interactions, potentially providing a bridge between these two structures. A higher cholesterol content appears to promote lipid persistent binding to the receptor.

Pair of Residue Substitutions at the Outer Mouth of the Channel Pore Act as Inputs for a Boolean Logic "OR" Gate Based on the Glycine Receptor

The glycine receptor (GlyR) is a ligand-activated chloride channel, whose mutations are the major cause of hereditary hyperekplexia. The hyperekplexia-causing R271Q mutation, which is located at the extracellular outer mouth of the channel pore, dramatically impairs the GlyR function manifesting a reduced sensitivity toward glycine. This study reports that a second mutation, S273D, rescues the function of the R271Q GlyR to that of the wild-type (WT) GlyR. Surprisingly, the S273D mutation, when introduced to the WT GlyR, does not further increase the receptor function. In other words, the compromised function of the 271Q 273S GlyR (i.e., the R271Q GlyR) can be rescued to WT levels by the introduction of either, or both, of the Q271R and S273D substitutions. From the perspective of Boolean logic gates, the Q271R and S273D substitutions act as inputs for an OR gate based on the GlyR. Further experiments revealed that the negative-charge carried by the 273 residue is essential for the expression of the OR gate and that the expression of the OR gate is residue-position-specific. In addition, mechanistic investigation implied that the 273 residue influences the 271 residue, which might underpin the unique nonadditive OR gate relationship between these two residues. Such an ion-channel-based OR gate, expressing output in the form of electrical current, could potentially be developed to digitally manipulate neuronal activity.

Simulating ion channel activation mechanisms using swarms of trajectories

Atomic-level studies of protein activity represent a significant challenge as a result of the complexity of conformational changes occurring on wide-ranging timescales, often greatly exceeding that of even the longest simulations. A prime example is the elucidation of protein allosteric mechanisms, where localized perturbations transmit throughout a large macromolecule to generate a response signal. For example, the conversion of chemical to electrical signals during synaptic neurotransmission in the brain is achieved by specialized membrane proteins called pentameric ligand-gated ion channels. Here, the binding of a neurotransmitter results in a global conformational change to open an ion-conducting pore across the nerve cell membrane. X-ray crystallography has produced static structures of the open and closed states of the proton-gated GLIC pentameric ligand-gated ion channel protein, allowing for atomistic simulations that can uncover changes related to activation. We discuss a range of enhanced sampling approaches that could be used to explore activation mechanisms. In particular, we describe recent application of an atomistic string method, based on Roux's "swarms of trajectories" approach, to elucidate the sequence and interdependence of conformational changes during activation. We illustrate how this can be combined with transition analysis and Brownian dynamics to extract thermodynamic and kinetic information, leading to understanding of what controls ion channel function. © 2019 Wiley Periodicals, Inc.

Microsecond-timescale simulations suggest 5-HT-mediated preactivation of the 5-HT3A serotonin receptor

Aided by efforts to improve their speed and efficiency, molecular dynamics (MD) simulations provide an increasingly powerful tool to study the structure-function relationship of pentameric ligand-gated ion channels (pLGICs). However, accurate reporting of the channel state and observation of allosteric regulation by agonist binding with MD remains difficult due to the timescales necessary to equilibrate pLGICs from their artificial and crystalized conformation to a more native, membrane-bound conformation in silico. Here, we perform multiple all-atom MD simulations of the homomeric 5-hydroxytryptamine 3A (5-HT3A) serotonin receptor for 15 to 20 μs to demonstrate that such timescales are critical to observe the equilibration of a pLGIC from its crystalized conformation to a membrane-bound conformation. These timescales, which are an order of magnitude longer than any previous simulation of 5-HT3A, allow us to observe the dynamic binding and unbinding of 5-hydroxytryptamine (5-HT) (i.e., serotonin) to the binding pocket located on the extracellular domain (ECD) and allosteric regulation of the transmembrane domain (TMD) from synergistic 5-HT binding. While these timescales are not long enough to observe complete activation of 5-HT3A, the allosteric regulation of ion gating elements by 5-HT binding is indicative of a preactive state, which provides insight into molecular mechanisms that regulate channel activation from a resting state. This mechanistic insight, enabled by microsecond-timescale MD simulations, will allow a careful examination of the regulation of pLGICs at a molecular level, expanding our understanding of their function and elucidating key structural motifs that can be targeted for therapeutic regulation.

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.

Functional movements of the GABA type A receptor

We have performed a Monte Carlo simulation of the GABA type A receptor. We have analysed the configurations and developed a correlation tensor method to predict receptor gating.

Large-Scale Conformational Changes and Protein Function: Breaking the in silico Barrier

Large-scale conformational changes are essential to link protein structures with their function at the cell and organism scale, but have been elusive both experimentally and computationally. Over the past few years developments in cryo-electron microscopy and crystallography techniques have started to reveal multiple snapshots of increasingly large and flexible systems, deemed impossible only short time ago. As structural information accumulates, theoretical methods become central to understand how different conformers interconvert to mediate biological function. Here we briefly survey current in silico methods to tackle large conformational changes, reviewing recent examples of cross-validation of experiments and computational predictions, which show how the integration of different scale simulations with biological information is already starting to break the barriers between the in silico, in vitro, and in vivo worlds, shedding new light onto complex biological problems inaccessible so far.

Proline Residues in the Transmembrane/Extracellular Domain Interface Loops Have Different Behaviors in 5-HT3 and nACh Receptors

Cys-loop receptors are important drug targets that are involved in signaling in the nervous system. The binding of neurotransmitters in the extracellular region of these receptors triggers an allosteric activation mechanism, the full details of which remain elusive, although structurally flexible loops in the interface between the extracellular region of Cys-loop receptors and the pore-forming transmembrane domain are known to play an important role. Here we explore the roles of three largely conserved Pro residues in two of these loops, the Cys-loop and M2-M3 loop, in 5-HT₃A and α7 nACh receptors. The data from natural and noncanonical amino acid mutagenesis suggest that in both proteins a Pro is essential in the Cys-loop, probably because of its enhanced ability to form a cis peptide bond, although other factors are also involved. The important characteristics of Pros in the M2-M3 loop, however, differ in these two receptors: in the 5-HT₃ receptor, the Pros can be replaced by some charged amino acids resulting in EC₅₀s similar to those of wild-type receptors, while such substitutions in the nACh receptor ablate function. Ala substitution at one of these Pros also has different effects in the two receptors. Thus, our data show that even highly conserved residues can have distinct behaviors in related Cys-loop receptors.

Identification of the Initial Steps in Signal Transduction in the α4β2 Nicotinic Receptor: Insights from Equilibrium and Nonequilibrium Simulations

Nicotinic acetylcholine receptors (nAChRs) modulate synaptic transmission in the nervous system. These receptors have emerged as therapeutic targets in drug discovery for treating several conditions, including Alzheimer's disease, pain, and nicotine addiction. In this in silico study, we use a combination of equilibrium and nonequilibrium molecular dynamics simulations to map dynamic and structural changes induced by nicotine in the human α4β2 nAChR. They reveal a striking pattern of communication between the extracellular binding pockets and the transmembrane domains (TMDs) and show the sequence of conformational changes associated with the initial steps in this process. We propose a general mechanism for signal transduction for Cys-loop receptors: the mechanistic steps for communication proceed firstly through loop C in the principal subunit, and are subsequently transmitted, gradually and cumulatively, to loop F of the complementary subunit, and then to the TMDs through the M2-M3 linker.

Molecular dynamics simulations of dihydro-β-erythroidine bound to the human α4β2 nicotinic acetylcholine receptor

BACKGROUND AND PURPOSE

The heteromeric α4β2 nicotinic acetylcholine receptor (nAChR) is abundant in the human brain and is associated with a range of CNS disorders. This nAChR subtype has been recently crystallised in a conformation that was proposed to represent a desensitised state. Here, we investigated the conformational transition mechanism of this nAChR from a desensitised to a closed/resting state.

EXPERIMENTAL APPROACH

The competitive antagonist dihydro-β-erythroidine (DHβE) was modelled by replacement of the agonist nicotine in the α4β2 nAChR experimental structure. DHβE is used both in vitro and in vivo for its ability to block α4β2 nAChRs. This system was studied by three molecular dynamics simulations with a combined simulation time of 2.6 μs. Electrophysiological studies of mutated receptors were performed to validate the simulation results.

KEY RESULTS

The relative positions of the extracellular and transmembrane domains in the models are distinct from those of the desensitised state structure and are compatible with experimental structures of Cys-loop receptors captured in a closed/resting state.

CONCLUSIONS AND IMPLICATIONS

Our model suggests that the side chains of α4 L257 (9') and α4 L264 (16') are the main constrictions in the transmembrane pore. The involvement of position 9' in channel gating is well established, but position 16' was only previously identified as a gate for the bacterial channels, ELIC and GLIC. L257 but not L264 was found to influence the slow component of desensitisation. The structure of the antagonist-bound state proposed here should be valuable for the development of therapeutic or insecticide compounds.

Mutagenesis computer experiments in pentameric ligand-gated ion channels: the role of simulation tools with different resolution

Pentameric ligand-gated ion channels (pLGICs) are an important class of widely expressed membrane neuroreceptors, which play a crucial role in fast synaptic communications and are involved in several neurological conditions. They are activated by the binding of neurotransmitters, which trigger the transmission of an electrical signal via facilitated ion flux. They can also be activated, inhibited or modulated by a number of drugs. Mutagenesis electrophysiology experiments, with natural or unnatural amino acids, have provided a large body of functional data that, together with emerging structural information from X-ray spectroscopy and cryo-electron microscopy, are helping unravel the complex working mechanisms of these neuroreceptors. Computer simulations are complementing these mutagenesis experiments, with insights at various levels of accuracy and resolution. Here, we review how a selection of computational tools, including first principles methods, classical molecular dynamics and enhanced sampling techniques, are contributing to construct a picture of how pLGICs function and can be pharmacologically targeted to treat the disorders they are responsible for.

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

Allostery in Its Many Disguises: From Theory to Applications

Allosteric regulation plays an important role in many biological processes, such as signal transduction, transcriptional regulation, and metabolism. Allostery is rooted in the fundamental physical properties of macromolecular systems, but its underlying mechanisms are still poorly understood. A collection of contributions to a recent interdisciplinary CECAM (Center Européen de Calcul Atomique et Moléculaire) workshop is used here to provide an overview of the progress and remaining limitations in the understanding of the mechanistic foundations of allostery gained from computational and experimental analyses of real protein systems and model systems. The main conceptual frameworks instrumental in driving the field are discussed. We illustrate the role of these frameworks in illuminating molecular mechanisms and explaining cellular processes, and describe some of their promising practical applications in engineering molecular sensors and informing drug design efforts.

"Are We There Yet?" Quest for Treatment of Refractory Epilepsy

Biochemical Autoregulatory Gene Therapy for Focal Epilepsy Lieb A, Qiu Y, Dixon CL, Heller JP, Walker MC, Schorge S, Kullmann DM. Nat Med. 2018;24:1324-1329. Despite the introduction of more than one dozen new antiepileptic drugs in the past 20 years, approximately one-third of people who develop epilepsy continue to have seizures on mono- or polytherapy. Viral-vector-mediated gene transfer offers the opportunity to design a rational treatment that builds on mechanistic understanding of seizure generation and that can be targeted to specific neuronal populations in epileptogenic foci. Several such strategies have shown encouraging results in different animal models, although clinical translation is limited by possible effects on circuits underlying cognitive, mnemonic, sensory, or motor function. Here, we describe an autoregulatory antiepileptic gene therapy, which relies on neuronal inhibition in response to elevations in extracellular glutamate. It is effective in a rodent model of focal epilepsy and is well tolerated, thus lowering the barrier to clinical translation.

A structural look at GABAA receptor mutations linked to epilepsy syndromes

Understanding the genetic variation in GABA_A receptor subunit genes (GABRs), GABRA1-6, GABRB1-3, GABRG1-3 and GABRD, in individuals affected by epilepsy may improve the diagnosis and treatment of epilepsy syndromes through identification of disease-associated variants. However, the lack of functional analysis and validation of many novel and previously reported familial and de novo mutations have made it challenging to address meaningful gene associations with epilepsy syndromes. GABA_A receptors belong to the Cys-loop receptor family. Even though GABA_A receptor mutant residues are widespread among different GABRs, their frequent occurrence in important structural domains that share common functional features suggests associations between structure and function.

Trans- Cis Proline Switches in a Pentameric Ligand-Gated Ion Channel: How They Are Affected by and How They Affect the Biomolecular Environment

Pentameric ligand-gated ion channels (pLGICs) are important neuroreceptors, embedded in neuronal membranes, that mediate fast synaptic transmission. The molecular details of their working mechanisms have still to be fully unravelled due to their complexity and limited structural information available. Here we focus on a potential molecular switch in a prototypical pLGIC, the serotonin-activated 5-HT₃ receptor, consisting of the trans- cis isomerization of a proline at the interface between the extracellular and transmembrane domain. Mutagenesis electrophysiology experiments previously showed that if such isomerization could not take place, the channel would not open, but the hypothetical role of this mechanism as key to channel gating is still debated. We investigate this switch within the receptor with molecular dynamics and enhanced sampling simulations. We analyze how the isomerization free energy landscape is affected by the receptor environment in comparison to simplified models. Moreover, we reveal how the isomerization, in turn, affects the structural and electrostatic properties of the receptor at the extracellular-transmembrane domain interface, e.g., by tuning the ion selectivity filter.