Protein dynamics and the allosteric transitions of pentameric receptor channels

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.

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.

Mechanisms Underlying the Strong Inhibition of Muscle-Type Nicotinic Receptors by Tetracaine

Nicotinic acetylcholine (ACh) receptors (nAChRs) are included among the targets of a variety of local anesthetics, although the molecular mechanisms of blockade are still poorly understood. Some local anesthetics, such as lidocaine, act on nAChRs by different means through their ability to present as both charged and uncharged molecules. Thus, we explored the mechanisms of nAChR blockade by tetracaine, which at physiological pH is almost exclusively present as a positively charged local anesthetic. The nAChRs from Torpedo electroplaques were transplanted to Xenopus oocytes and the currents elicited by ACh (I_ACh s), either alone or co-applied with tetracaine, were recorded. Tetracaine reversibly blocked I_ACh , with an IC₅₀ (i.e., the concentration required to inhibit half the maximum I_ACh ) in the submicromolar range. Notably, at very low concentrations (0.1 μM), tetracaine reduced I_ACh in a voltage-dependent manner, the more negative potentials produced greater inhibition, indicating open-channel blockade. When the tetracaine concentration was increased to 0.7 μM or above, voltage-independent inhibition was also observed, indicating closed-channel blockade. The I_ACh inhibition by pre-application of just 0.7 μM tetracaine before superfusion of ACh also corroborated the notion of tetracaine blockade of resting nAChRs. Furthermore, tetracaine markedly increased nAChR desensitization, mainly at concentrations equal or higher than 0.5 μM. Interestingly, tetracaine did not modify desensitization when its binding within the channel pore was prevented by holding the membrane at positive potentials. Tetracaine-nAChR interactions were assessed by virtual docking assays, using nAChR models in the closed and open states. These assays revealed that tetracaine binds at different sites of the nAChR located at the extracellular and transmembrane domains, in both open and closed conformations. Extracellular binding sites seem to be associated with closed-channel blockade; whereas two sites within the pore, with different affinities for tetracaine, contribute to open-channel blockade and the enhancement of desensitization, respectively. These results demonstrate a concentration-dependent heterogeneity of tetracaine actions on nAChRs, and contribute to a better understanding of the complex modulation of muscle-type nAChRs by local anesthetics. Furthermore, the combination of functional and virtual assays to decipher nAChR-tetracaine interactions has allowed us to tentatively assign the main nAChR residues involved in these modulating actions.

Systematic Dissociation Pathway Searches Guided by Principal Component Modes

We introduce a novel method, Pathway Search guided by Internal Motions (PSIM), that efficiently finds molecular dissociation pathways of a ligand-receptor system with guidance from principal component (PC) modes obtained from molecular dynamics (MD) simulations. Modeling ligand-receptor dissociation pathways can provide insights into molecular recognition and has practical applications, including understanding kinetic mechanisms and barriers to binding/unbinding as well as design of drugs with desired kinetic properties. PSIM uses PC modes in multilayer internal coordinates to identify natural molecular motions that guide the search for conformational switches and unbinding pathways. The new multilayer internal coordinates overcome problems with Cartesian and classical internal coordinates that fail to smoothly present dihedral rotation or generate nonphysical distortions. We used HIV-1 protease, which has large-scale flap motions, as an example protein to demonstrate use of the multilayer internal coordinates. We provide examples of algorithms and implementation of PSIM with alanine dipeptide and chemical host-guest systems, 2-naphthyl ethanol-β-cyclodextrin and tetramethylammonium-cryptophane complexes. Tetramethylammonium-cryptophane has slow binding/unbinding kinetics. Its residence time, the length to dissociate tetramethylammonium from the host, is ∼14 s from experiments, and PSIM revealed 4 dissociation pathways in approximately 150 CPU h. We also searched the releasing pathways for the product glyceraldehyde-3-phosphate from tryptophan synthase, and one complete dissociation pathway was constructed after running multiple search iterations in approximately 300 CPU h. With guidance by internal PC modes from MD simulations, the PSIM method has advantages over simulation-based methods to search for dissociation pathways of molecular systems with slow noncovalent kinetic behavior.

Monod-Wyman-Changeux Analysis of Ligand-Gated Ion Channel Mutants

We present a framework for computing the gating properties of ligand-gated ion channel mutants using the Monod-Wyman-Changeux (MWC) model of allostery. We derive simple analytic formulas for key functional properties such as the leakiness, dynamic range, half-maximal effective concentration ([EC₅₀]), and effective Hill coefficient, and explore the full spectrum of phenotypes that are accessible through mutations. Specifically, we consider mutations in the channel pore of nicotinic acetylcholine receptor (nAChR) and the ligand binding domain of a cyclic nucleotide-gated (CNG) ion channel, demonstrating how each mutation can be characterized as only affecting a subset of the biophysical parameters. In addition, we show how the unifying perspective offered by the MWC model allows us, perhaps surprisingly, to collapse the plethora of dose-response data from different classes of ion channels into a universal family of curves.

A mechanism for acetylcholine receptor gating based on structure, coupling, phi, and flip

Gupta et al. use single-channel electrophysiology to investigate the gating mechanism of acetylcholine receptor ion channels. They propose that channel opening starts at the M2–M3 linker and ligand-binding sites and proceeds through four brief intermediate conformations before ending with the collapse of a gate bubble.

Nicotinic acetylcholine receptors are allosteric proteins that generate membrane currents by isomerizing (“gating”) between resting and active conformations under the influence of neurotransmitters. Here, to explore the mechanisms that link the transmitter-binding sites (TBSs) with the distant gate, we use mutant cycle analyses to measure coupling between residue pairs, phi value analyses to sequence domain rearrangements, and current simulations to reproduce a microsecond shut component (“flip”) apparent in single-channel recordings. Significant interactions between amino acids separated by >15 Å are rare; an exception is between the αM2–M3 linkers and the TBSs that are ∼30 Å apart. Linker residues also make significant, local interactions within and between subunits. Phi value analyses indicate that without agonists, the linker is the first region in the protein to reach the gating transition state. Together, the phi pattern and flip component suggest that a complete, resting↔active allosteric transition involves passage through four brief intermediate states, with brief shut events arising from sojourns in all or a subset. We derive energy landscapes for gating with and without agonists, and propose a structure-based model in which resting→active starts with spontaneous rearrangements of the M2–M3 linkers and TBSs. These conformational changes stabilize a twisted extracellular domain to promote transmembrane helix tilting, gate dilation, and the formation of a “bubble” that collapses to initiate ion conduction. The energy landscapes suggest that twisting is the most energetically unfavorable step in the resting→active conformational change and that the rate-limiting step in the reverse process is bubble formation.

A Gating Mechanism of the Serotonin 5-HT3 Receptor

Our recently solved high-resolution structure of the serotonin 5-HT3 receptor (5-HT3R) delivered the first detailed structural insights for a mammalian pentameric ligand-gated ion channel. Based on this structure, we here performed a total of 2.8-μs all-atom molecular dynamics simulations to unravel at atomic detail how neurotransmitter binding on the extracellular domain induces sequential conformational transitions in the receptor, opening an ion channel and translating a chemical signal into electrical impulses across the membrane. We found that serotonin binding first induces distinct conformational fluctuations at the side chain of W156 in the highly conserved ligand-binding cage, followed by tilting-twisting movements of the extracellular domain which couple to the transmembrane TM2 helices, opening the hydrophobic gate at L260 and forming a continuous transmembrane water pathway. The structural transitions in the receptor's transmembrane part finally couple to the intracellular MA helix bundle, opening lateral ports for ion passage.

Function of the M1 π-helix in endplate receptor activation and desensitization

KEY POINTS

A conserved proline in M1 causes a kink between α and π helical segments. The kink is under greater tension in the resting versus active conformation. The kink and the agonist do not interact directly. The π-helix separates the gating functions of the extracellular and transmembrane domains. Mutations of the conserved proline and propofol increase desensitization.

ABSTRACT

Nicotinic acetylcholine receptors (AChRs) switch on/off to generate transient membrane currents (C↔O; closed-open 'gating') and enter/recover from long-lived, refractory states (O↔D; 'desensitization'). The M1 transmembrane helix of the muscle endplate AChR is linked to a β-strand of the extracellular domain that extends to a neurotransmitter binding site. We used electrophysiology to measure the effects of mutations of amino acids that are located at a proline kink in M1 that separates π and α helices, in both α (N217, V218 and P221) and non-α subunits. In related receptors, the kink is straighter and more stable in O vs. C structures (gating is 'spring-loaded'). None of the AChR kink mutations had a measureable effect on agonist affinity but many influenced the allosteric gating constant substantially. Side chains in the M1 α-helix experience extraordinarily large energy differences between C and O structures, probably because of a ∼2 Å displacement and tilt of M2 relative to M1. There is a discrete break in the character of the gating transition state between αN217 and αV218, indicating that the π-helix is a border between extracellular- and transmembrane-domain function. Mutations of the conserved M1 proline, and the anaesthetic propofol, increase a rate constant for desensitization. The results suggest that straightening of the M1 proline kink triggers AChR desensitization.