Mutational studies using a cation-conducting GABAA receptor reveal the selectivity determinants of the Cys-loop family of ligand-gated ion channels

The present study evaluates how four key amino acid residue positions (- 4' to - 1') within the M1-M2 linker of the GABA(A) receptor beta subunit influences ion selectivity of a cation-conducting GABA receptor. Cation selectivity was found to be highly dependent on the side-chains of the amino acid residues present. The critical factor for cation selectivity was the presence of a negatively charged Glu or Asp residue in the -1' position. Receptors containing the neutral amino acids Gln or Asn or a positively charged Arg residue were anion selective. In the presence of a -1' Glu residue, the amino acids in adjacent positions were also found to be important determinants of cation selectivity. Moreover, the length of the M1-M2 linker as well as the presence of a Pro residue within this segment also affected ion selectivity, suggesting that the local environment and three-dimensional position of the -1' Glu are essential determinants of cation permeation. Conversely, no specific amino acid residues were found to be essential for anion selectivity, suggesting that the basic architecture of the selectivity segment of this class of receptor channels is optimally suited for anion conduction.     

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

Charge selectivity of the Cys-loop family of ligand-gated ion channels

The determinants of charge selectivity of the Cys-loop family of ligand-gated ion channels have been studied for more than a decade. The investigations have mainly covered homomeric receptors e.g. the nicotinic acetylcholine receptor alpha7, the glycine receptor alpha1 and the serotonin receptor 5-HT(3A). Only recently, the determinants of charge selectivity of heteromeric receptors have been addressed for the GABA(A) receptor alpha2beta3gamma2. For all receptor subtypes, the selectivity determinants have been located to an intracellular linker between transmembrane domains M1 and M2. Two features of the M1-M2 linker appear to control ion selectivity. A central role for charged amino acid residues in selectivity has been almost universally observed. Furthermore, recent studies point to an important role of the size of the narrowest constriction in the pore. In the present review, these determinants of charge selectivity of the Cys-loop family of ligand-gated ion channels will be discussed in detail.     

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

Desensitization mechanism in prokaryotic ligand-gated ion channel

Crystal structures of Gloeobacter violaceus ligand-gated ion channel (GLIC), a proton-gated prokaryotic homologue of pentameric ligand-gated ion channel (LGIC) from G. violaceus, have provided high-resolution models of the channel architecture and its role in selective ion conduction and drug binding. However, it is still unclear which functional states of the LGIC gating scheme these crystal structures represent. Much of this uncertainty arises from a lack of thorough understanding of the functional properties of these prokaryotic channels. To elucidate the molecular events that constitute gating, we have carried out an extensive characterization of GLIC function and dynamics in reconstituted proteoliposomes by patch clamp measurements and EPR spectroscopy. We find that GLIC channels show rapid activation upon jumps to acidic pH followed by a time-dependent loss of conductance because of desensitization. GLIC desensitization is strongly coupled to activation and is modulated by voltage, permeant ions, pore-blocking drugs, and membrane cholesterol. Many of these properties are parallel to functions observed in members of eukaryotic LGIC. Conformational changes in loop C, measured by site-directed spin labeling and EPR spectroscopy, reveal immobilization during desensitization analogous to changes in LGIC and acetylcholine binding protein. Together, our studies suggest conservation of mechanistic aspects of desensitization among LGICs of prokaryotic and eukaryotic origin.     

Ginsenosides regulate ligand-gated ion channels from the outside

Treatment with ginsenosides, the major active ingredients of Panax ginseng, produces a variety of physiological effects on the central and peripheral nervous systems. Ginsenosides inhibit various types of ligand-gated ion channel but it is not clear whether they act from within or outside the cell since they are somewhat membrane-permeable. In the present study, we used the Xenopus oocyte gene expression system to determine from which side of the cell membrane the ginsenoside Rg3 (Rg3), and M4, a ginsenoside metabolite, act to regulate ligand-gated ion channel activity. Ligand-gated ion currents were measured using the two-electrode voltage clamp technique. Rg3 and M4 inhibited 5-HT3A and a3b4 nACh receptor-mediated ion currents when present outside of the cell but not when injected intracellularly. We also examined the effect of these agents on oocytes expressing the gustatory cGMP-gated ion channel, which is known to have a cGMP binding site on the intracellular side of the plasma membrane and is only activated by cytosolic cGMP. Rg3 inhibited cGMP-gated ion currents when applied extracellularly or to an outside-out patch clamp, but not when injected into the cytosol or when using an excised inside-out patch clamp. These results indicate that Rg3 and M4 regulate ligand-gated ion channel activity from the extracellular side.     

Heterogeneous behavior of metalloproteins toward metal ion binding and selectivity: insights from molecular dynamics studies

About one-third of the existing proteins require metal ions as cofactors for their catalytic activities and structural complexities. While many of them bind only to a specific metal, others bind to multiple (different) metal ions. However, the exact mechanism of their metal preference has not been deduced to clarity. In this study, we used molecular dynamics (MD) simulations to investigate whether a cognate metal (bound to the structure) can be replaced with other similar metal ions. We have chosen seven different proteins (phospholipase A₂, sucrose phosphatase, pyrazinamidase, cysteine dioxygenase (CDO), plastocyanin, monoclonal anti-CD4 antibody Q425, and synaptotagmin 1 C₂B domain) bound to seven different divalent metal ions (Ca²⁺, Mg²⁺, Zn²⁺, Fe²⁺, Cu²⁺, Ba²⁺, and Sr²⁺, respectively). In total, 49 MD simulations each of 50 ns were performed and each trajectory was analyzed independently. Results demonstrate that in some cases, cognate metal ions can be exchanged with similar metal ions. On the contrary, some proteins show binding affinity specifically to their cognate metal ions. Surprisingly, two proteins CDO and plastocyanin which are known to bind Fe²⁺ and Cu²⁺, respectively, do not exhibit binding affinity to any metal ion. Furthermore, the study reveals that in some cases, the active site topology remains rigid even without cognate metals, whereas, some require them for their active site stability. Thus, it will be interesting to experimentally verify the accuracy of these observations obtained computationally. Moreover, the study can help in designing novel active sites for proteins to sequester metal ions particularly of toxic nature.     

Sequence dependence of amyloid fibril formation: insights from molecular dynamics simulations

The clarification of the physico-chemical determinants underlying amyloid deposition is critical for our understanding of misfolding diseases. With this purpose we have performed a systematic all-atom molecular dynamics (MD) study of a series of single point mutants of the de novo designed amyloidogenic peptide STVIIE. Sixteen different 50ns long simulations using explicit solvent have been carried out starting from four different conformations of a polymeric six-stranded beta-sheet. The simulations have provided evidence for the influence of a small number of site-specific hydrophobic interactions on the packing and stabilization of nascent aggregates, as well as the interplay between side-chain interactions and the net charge of the molecule on the strand arrangement of polymeric beta-sheets. This MD analysis has also shed light into the origin of the position dependence on mutation of beta-sheet polymerization that was found experimentally for this model system. Our results suggest that MD can be applied to detect critical positions for beta-sheet aggregation within a given amyloidogenic stretch. Studies similar to the one presented here can guide site-directed mutations or the design of drugs that specifically disrupt the key stabilizing interactions of beta-sheet aggregates.     

Phosphorylation of ligand-gated ion channels: a possible mode of synaptic plasticity.

Most neurotransmitter receptors examined to date have been shown either to be regulated by protein phosphorylation or to contain consensus sequences for phosphorylation by protein kinases. Neurotransmitter receptors that mediate rapid synaptic transmission in the nervous system are the ligand-gated ion channels and include the nicotinic acetylcholine receptors of muscle and nerve and the excitatory and inhibitory amino acid receptors: the glutamate, GABAA, and glycine receptors. These receptors are multimeric proteins composed of homologous subunits which each span the membrane several times and contain a large intracellular loop that is a mosaic of consensus sites for protein phosphorylation. Recent evidence has suggested that extracellular signals released from the presynaptic neuron, such as neurotransmitters and neuropeptides as well as an extracellular matrix protein, regulate the phosphorylation of ligand-gated ion channels. The functional effects of phosphorylation are varied and include the regulation of receptor desensitization rate, subunit assembly, and receptor aggregation at the synapse. These results suggest that phosphorylation of neurotransmitter receptors represents a major mechanism in the regulation of their function and may play an important role in synaptic plasticity.     

"Adaptive" behavior of ligand-gated ion channels: constraints by thermodynamics.

The calcium-induced calcium release channel of the cardiac sarcoplasmic reticulum has been reported to inactivate in a novel manner (termed "adaptation"), which permits reactivation by exposure to successively higher concentrations of calcium. I examined the limitations placed by thermodynamics on the possible kinetic mechanisms for such behavior. The mechanism suggested by Gyorke and Fill, in which the affinity of a calcium-binding site decreases during adaptation, is not thermodynamically feasible for a passive system, but requires an external input of free energy. Possible sources of such energy are 1) metabolic energy, which is excluded by the fact that adaptation was observed in isolated channels in the absence of ATP, or 2) coupling of ion permeation to gating, for which there is currently no evidence. I derived a general limit on the thermodynamic feasibility of a sequence of channel activations and adaptations, irrespective of channel kinetics, from the requirement that the free energy must decrease during the spontaneous evolution of the system from the state existing immediately after a step increase in [Ca2+] to the state of maximum open probability that follows. The opening of the channel must involve an increase in free energy, which must be compensated by the free energy released by the incremental binding of calcium. This requirement leads to a complicated system of inequalities, which was simplified and manipulated algebraically into the form of a linear programming problem. Numerical solution of this problem showed that the sequence of adaptations of the SR channel observed by Gyorke and Fill requires the presence of at least 10 calcium-binding sites on the channel if it is to occur in the absence of exogenous sources of free energy. This indicates either that a large number of calcium-binding sites participate in the regulation of the SR calcium release channel, or that the existing data are significantly flawed with respect to the low open probability in the resting state, the importance of "calcium spike" artifacts from flash photolysis, or both.     

Neurotoxins from Snake Venoms and α-Conotoxin ImI Inhibit Functionally Active Ionotropic γ-Aminobutyric Acid (GABA) Receptors

Ionotropic receptors of γ-aminobutyric acid (GABA_AR) regulate neuronal inhibition and are targeted by benzodiazepines and general anesthetics. We show that a fluorescent derivative of α-cobratoxin (α-Ctx), belonging to the family of three-finger toxins from snake venoms, specifically stained the α1β3γ2 receptor; and at 10 μm α-Ctx completely blocked GABA-induced currents in this receptor expressed in Xenopus oocytes (IC₅₀ = 236 nm) and less potently inhibited α1β2γ2 ≈ α2β2γ2 > α5β2γ2 > α2β3γ2 and α1β3δ GABA_ARs. The α1β3γ2 receptor was also inhibited by some other three-finger toxins, long α-neurotoxin Ls III and nonconventional toxin WTX. α-Conotoxin ImI displayed inhibitory activity as well. Electrophysiology experiments showed mixed competitive and noncompetitive α-Ctx action. Fluorescent α-Ctx, however, could be displaced by muscimol indicating that most of the α-Ctx-binding sites overlap with the orthosteric sites at the β/α subunit interface. Modeling and molecular dynamic studies indicated that α-Ctx or α-bungarotoxin seem to interact with GABA_AR in a way similar to their interaction with the acetylcholine-binding protein or the ligand-binding domain of nicotinic receptors. This was supported by mutagenesis studies and experiments with α-conotoxin ImI and a chimeric Naja oxiana α-neurotoxin indicating that the major role in α-Ctx binding to GABA_AR is played by the tip of its central loop II accommodating under loop C of the receptors.     

P2X1 and P2X3 receptors form stable trimers: a novel structural motif of ligand-gated ion channels.

P2X receptors are cation channels gated by extracellular ATP. The seven known P2X isoforms possess no sequence homology with other proteins. Here we studied the quaternary structure of P2X receptors by chemical cross-linking and blue native PAGE. P2X1 and P2X3 were N-terminally tagged with six histidine residues to allow for non-denaturing receptor isolation from cRNA-injected, [35S]methionine-labeled oocytes. The His-tag did not change the electrophysiological properties of the P2X1 receptor. His-P2X1 was found to carry four N-glycans per polypeptide chain, only one of which acquired Endo H resistance en route to the plasma membrane. 3, 3'-Dithiobis(sulfosuccinimidylpropionate) (DTSSP) and two of three bifunctional analogues of the P2X receptor antagonist pyridoxalphosphate-6-azophenyl-2',4'-disulfonic acid (PPADS) cross-linked digitonin-solubilized His-P2X1 and His-P2X3 quantitatively to homo-trimers. Likewise, when analyzed by blue native PAGE, P2X receptors purified in digitonin or dodecyl-beta-D-maltoside migrated entirely as non-covalently linked homo-trimers, whereas the alpha2 beta gamma delta nicotinic acetylcholine receptor (used as a positive control) migrated as the expected pentamer. P2X monomers remained undetected soon after synthesis, indicating that trimerization occurred in the endoplasmic reticulum. The plasma membrane form of His-P2X1 was also identified as a homo-trimer. If n-octylglucoside was used for P2X receptor solubilization, homo-hexamers were observed, suggesting that trimers can aggregate to form larger complexes. We conclude that trimers represent an essential element of P2X receptor structure. Keywords: blue native PAGE/cross-linking/P2X receptor/quaternary structure.     

Dominant-negative effects in prion diseases: insights from molecular dynamics simulations on mouse prion protein chimeras

Mutations in the prion protein (PrP) can cause spontaneous prion diseases in humans (Hu) and animals. In transgenic mice, mutations can determine the susceptibility to the infection of different prion strains. Some of these mutations also show a dominant-negative effect, thus halting the replication process by which wild type mouse (Mo) PrP is converted into Mo scrapie. Using all-atom molecular dynamics (MD) simulations, here we studied the structure of HuPrP, MoPrP, 10?Hu/MoPrP chimeras, and 1 Mo/sheepPrP chimera in explicit solvent. Overall, ～2?μs of MD were collected. Our findings suggest that the interactions between α1 helix and N-terminal of α3 helix are critical in prion propagation, whereas the β2–α2 loop conformation plays a role in the dominant-negative effect.

An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:4.     

Recognition dynamics of dopamine to human Monoamine oxidase B: role of Leu171/Gln206 and conserved water molecules in the active site cavity

The human Monoamine oxidase (hMAO) metabolizes several biogenic amine neurotransmitters and is involved in different neurological disorders. Extensive MD simulation studies of dopamine-docked hMAO B structures have revealed the stabilization of amino-terminal of the substrate by a direct and water-mediated interaction of catalytic tyrosines, Gln206, and Leu171 residues. The catechol ring of the substrate is stabilized by Leu171(C–H)?π(Dop)?(H–C) Ile199 interaction. Several conserved water molecules are observed to play a role in the recognition of substrate to the enzyme, where W1 and W2 associate in dopamine– FAD interaction, reversible dynamics of W3 and W4 influenced the coupling of Tyr435 to Trp432 and FAD, and W5 and W8 stabilized the catalytic Tyr188/398 residues. The W6, W7, and W8 water centers are involved in the recognition of catalytic residues and FAD with the N⁺- site of dopamine through hydrogen bonding interaction. The recognition of substrate to gating residues is made through W9, W10, and W11 water centers. Beside the interplay of water molecules, the catalytic aromatic cage has also been stabilized by π?water, π?C–H, and π?π interactions. The topology of conserved water molecular sites along with the hydration dynamics of catalytic residues, FAD, and dopamine has added a new feature on the substrate binding chemistry in hMAO B which may be useful for substrate analog inhibitor design.     

Ligand unbinding pathways from the vitamin D receptor studied by molecular dynamics simulations

Molecular dynamics simulation techniques have been used to study the unbinding pathways of 1α,25-dihydroxyvitamin D₃ from the ligand-binding pocket of the vitamin D receptor (VDR). The pathways observed in a large number of relatively short (<200 ps) random acceleration molecular dynamics (RAMD) trajectories were found to be in fair agreement, both in terms of pathway locations and deduced relative preferences, compared to targeted molecular dynamics (TMD) and streered molecular dynamics simulations (SMD). However, the high-velocity ligand expulsions of RAMD tend to favor straight expulsion trajectories and the observed relative frequencies of different pathways were biased towards the probability of entering a particular exit channel. Simulations indicated that for VDR the unbinding pathway between the H1–H2 loop and the β-sheet between H5 and H6 is more favorable than the pathway located between the H1–H2 loop and H3. The latter pathway has been suggested to be the most likely unbinding path for thyroid hormone receptors (TRs) and a likely path for retinoic acid receptor. Ligand entry/exit through these two pathways would not require displacement of H12 from its agonistic position. Differences in the packing of the H1, H2, H3 and β-sheet region explain the changed relative preference of the two unbinding pathways in VDR and TRs. Based on the crystal structures of the ligand binding domains of class 2 nuclear receptors, whose members are VDR and TRs, this receptor class can be divided in two groups according to the packing of the H1, H2, H3 and β-sheet region. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Role of 6-Mercaptopurine in the potential therapeutic targets DNA base pairs and G-quadruplex DNA: insights from quantum chemical and molecular dynamics simulations

The theoretical studies on DNA with the anticancer drug 6-Mercaptopurine (6-MP) are investigated using theoretical methods to shed light on drug designing. Among the DNA base pairs considered, 6-MP is stacked with GC with the highest interaction energy of –46.19 kcal/mol. Structural parameters revealed that structure of the DNA base pairs is deviated from the planarity of the equilibrium position due to the formation of hydrogen bonds and stacking interactions with 6-MP. These deviations are verified through the systematic comparison between X–H bond contraction and elongation and the associated blue shift and red shift values by both NBO analysis and vibrational analysis. Bent’s rule is verified for the C–H bond contraction in the 6-MP interacted base pairs. The AIM results disclose that the higher values of electron density (ρ) and Laplacian of electron density (?²ρ) indicate the increased overlap between the orbitals that represent the strong interaction and positive values of the total electron density show the closed-shell interaction. The relative sensitivity of the chemical shift values for the DNA base pairs with 6-MP is investigated to confirm the hydrogen bond strength. Molecular dynamics simulation studies of G-quadruplex DNA d(TGGGGT)₄ with 6-MP revealed that the incorporation of 6-MP appears to cause local distortions and destabilize the G-quadruplex DNA.