From hopanoids to cholesterol: Molecular clocks of pentameric ligand-gated ion channels

Pentameric ligand-gated ion channels (pLGICs) and their lipid microenvironments appear to have acquired mutually adaptive traits along evolution: 1) the three-ring architecture of their transmembrane (TM) region; 2) the ability of the outermost TM ring to convey lipid signals to the middle ring, which passes them on to the central pore ring, and 3) consensus motifs for sterol recognition in all pLGICs. Hopanoids are triterpenoid fossil lipids that constitute invaluable biomarkers for tracing evolution at the molecular scale. The cyanobacterium Gloeobacter violaceus is the oldest known living organism in which the X-ray structure of its pLGIC, GLIC, reveals the presence of the above attributes and, as discussed in this review, the ability to bind hopanoids. ELIC, the pLGIC from the bacillus Erwinia chrysanthemi is the only other known case to date. Both prokaryotes lack cholesterol but their pLGICs exhibit the same sterol motifs as mammalian pLGIC. This remarkable conservation suggests that the association of sterols and hopanoid surrogate molecules arose from the early need in prokaryotes to stabilize pLGIC TM regions by means of relatively rigid lipid molecules. The conservation of these phenotypic traits along such a long phylogenetic span leads us to suggest the possible co-evolution of these sterols with pLGICs.     

Water and ion permeability of a claudin model: A computational study

At present, the three‐dimensional structure of the multimeric paracellular claudin pore is unknown. Using extant biophysical data concerning the size of the pore and permeation of water and cations through it, two three‐dimensional models of the pore are constructed in silico. Molecular Dynamics (MD) calculations are then performed to compute water and sodium ion permeation fluxes under the influence of applied hydrostatic pressure. Comparison to experiment is made, under the assumption that the hydrostatic pressure applied in the simulations is equivalent to osmotic pressure induced in experimental measurements of water/ion permeability. One model, in which pore‐lining charged is distributed evenly over a selectivity filter section 10–16 Å in length, is found to be generally consistent with experimental data concerning the dependence of water and ion permeation on channel pore diameter, pore length, and the sign and magnitude of pore lining charge. The molecular coupling mechanism between water and ion flow under conditions where hydrostatic pressure is applied is computationally elucidated. Proteins 2016; 84:305–315. © 2016 Wiley Periodicals, Inc.     

Spontaneous cleavage of proteins at serine and threonine is facilitated by zinc

Old proteins are widely distributed in the body. Over time, they deteriorate and many spontaneous reactions, for example isomerisation of Asp and Asn, can be replicated by incubation of peptides under physiological conditions. One of the signatures of long‐lived proteins that has proven to be difficult to replicate in vitro is cleavage on the N‐terminal side of Ser residues, and this is important since cleavage at Ser, and also Thr, has been observed in a number of human proteins. In this study, the autolysis of Ser‐ and Thr‐containing peptides was investigated with particular reference to discovering factors that promote cleavage adjacent to Ser/Thr at neutral pH. It was found that zinc catalyses cleavage of the peptide bond on the N‐terminal side of Ser residues and further that this process is markedly accelerated if a His residue is adjacent to the Ser. NMR analysis indicated that the imidazole group co‐ordinates zinc and that once zinc is co‐ordinated, it can polarize the carbonyl group of the peptide bond in a manner analogous to that observed in the active site of the metalloexopeptidase, carboxypeptidase A. The hydroxyl side chain of Ser/Thr is then able to cleave the adjacent peptide bond. These observations enable an understanding of the origin of common truncations observed in long‐lived proteins, for example truncation on the N‐terminal side of Ser 8 in Abeta, Ser 19 in alpha B crystallin and Ser 66 in alpha A crystallin. The presence of zinc may therefore significantly affect the long‐term stability of cellular proteins.     

Molecular dynamics studies on structure and dynamics of phospholamban monomer and pentamer in membranes

Phospholamban (PLB) is an integral membrane protein of 52 residues that regulates the activity of the sarcoplasmic reticulum calcium pump in cardiac muscle cells through reversible phosphorylation of Ser16. To explore its possible conformations and dynamics in a monomeric state, we have performed comparative molecular dynamics simulations of unphosphorylated and phosphorylated PLB (pPLB) with various orientations in POPC membranes. The simulations indicate that dynamics of the cytoplasmic domain is highly dependent on its interactions with membranes, that is, large conformational changes in the absence of membrane interactions, but very restricted dynamics in their presence. pPLB shows more structural flexibility in its cytoplasmic domain, which is consistent with experimental observations. We have also performed a simulation of a PLB pentameric structure (the so‐called bellflower model), recently determined in micelles, to investigate its behaviors in a POPC membrane. The cytoplasmic domain in each monomer shows uncorrelated dynamics and undergoes large conformational changes toward the membrane surface during the simulation, which supports the so‐called pinwheel model of the PLB pentamer structure. The hydrophobic nature of the pentameric pore excludes water molecules in the pore region, which illustrates that the pore appears to be an energetic barrier for ion and water translocation. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Analysis of Hyperekplexia Mutations Identifies Transmembrane Domain Rearrangements That Mediate Glycine Receptor Activation

Pentameric ligand-gated ion channels (pLGICs) mediate numerous physiological processes and are therapeutic targets for a wide range of clinical indications. Elucidating the structural differences between their closed and open states may help in designing improved drugs that bias receptors toward the desired conformational state. We recently showed that two new hyperekplexia mutations, Q226E and V280M, induced spontaneous activity in α1 glycine receptors. Gln-226, located near the top of transmembrane (TM) 1, is closely apposed to Arg-271 at the top of TM2 in the neighboring subunit. Using mutant cycle analysis, we inferred that Q226E induces activation via an enhanced electrostatic attraction to Arg-271. This would tilt the top of TM2 toward TM1 and hence away from the pore axis to open the channel. We also concluded that the increased side chain volume of V280M, in the TM2-TM3 loop, exerts a steric repulsion against Ile-225 at the top of TM1 in the neighboring subunit. We infer that this steric repulsion would tilt the top of TM3 radially outwards against the stationary TM1 and thus provide space for TM2 to relax away from the pore axis to create an open channel. Because the transmembrane domain movements inferred from this functional analysis are consistent with the structural differences evident in the x-ray atomic structures of closed and open state bacterial pLGICs, we propyose that the model of pLGIC activation as outlined here may be broadly applicable across the eukaryotic pLGIC receptor family.     

Brownian dynamics simulation of ion flow through porin channels

Bacterial porins, which allow the passage of solutes across the outer bacterial membrane, are structurally well characterized. They therefore lend themselves to detailed studies of the determinants of ion flow through transmembraneous channels. In a comparative study, we have performed Brownian dynamics simulations to obtain statistically significant transfer efficiencies for cations and anions through matrix porin OmpF, osmoporin OmpK36, phosphoporin PhoE and two OmpF charge mutants.The simulations show that the electrostatic potential at the highly charged channel constriction serves to enhance ion permeability of either cations or anions, dependent on the type of porin. At the same time translocation of counterions is not severely impeded. At the constriction, cations and anions follow distinct trajectories, due to the segregation of basic and acidic protein residues.Simulated ion selectivity and relative conductance agree well with experimental values, and are dependent crucially on the charge constellation at the pore constriction. The experimentally observed decrease in ion selectivity and single channel conductance with increasing ionic strength is well reproduced and can be attributed to electrostatic shielding of the pore lining.     

Common binding sites for cholesterol and neurosteroids on a pentameric ligand-gated ion channel

Cholesterol is an essential component of cell membranes, and is required for mammalian pentameric ligand-gated ion channel (pLGIC) function. Computational studies suggest direct interactions between cholesterol and pLGICs but experimental evidence identifying specific binding sites is limited. In this study, we mapped cholesterol binding to Gloeobacter ligand-gated ion channel (GLIC), a model pLGIC chosen for its high level of expression, existing crystal structure, and previous use as a prototypic pLGIC. Using two cholesterol analogue photolabeling reagents with the photoreactive moiety on opposite ends of the sterol, we identified two cholesterol binding sites: an intersubunit site between TM3 and TM1 of adjacent subunits and an intrasubunit site between TM1 and TM4. In both the inter- and intrasubunit sites, cholesterol is oriented such that the 3‑OH group points toward the center of the transmembrane domains rather than toward either the cytosolic or extracellular surfaces. We then compared this binding to that of the cholesterol metabolite, allopregnanolone, a neurosteroid that allosterically modulates pLGICs. The same binding pockets were identified for allopregnanolone and cholesterol, but the binding orientation of the two ligands was markedly different, with the 3‑OH group of allopregnanolone pointing to the intra- and extracellular termini of the transmembrane domains rather than to their centers. We also found that cholesterol increases, whereas allopregnanolone decreases the thermal stability of GLIC. These data indicate that cholesterol and neurosteroids bind to common hydrophobic pockets in the model pLGIC, GLIC, but that their effects depend on the orientation and specific molecular interactions unique to each sterol.     

Pentagon packing models for "all-pentamer" virus structures.

A connection is made between 1) the observed structures of virus capsids whose capsomers are all pentamers and 2) the mathematical problem of determination of the largest size of a given number of equal regular spherical pentagons that can be packed on the surface of the unit sphere without overlapping. It is found that papillomaviruses provide the conjectured solution to the spherical pentagon packing problem for 72 pentagons. Thus, a study of some virus structures has given additional insight into a mathematical problem. At the same time this mathematical problem enables prediction of an octahedral form of papillomavirus particles consisting of 24 pentamers. It is also found that the various tubular and spherical "all-pentamer" virus structures identified so far can be represented by closet-packing arrangements of equal morphological units composed of equal regular pentagons on a cylinder and on a sphere.     

Exploration of the structural features defining the conduction properties of a synthetic ion channel.

The finite-difference Poisson-Boltzmann methodology was applied to a series of parallel, alpha-helical bundle models of the designed ion channel peptide Ac-(LSSLLSL)3-CONH2. This method is able to fully describe the current-voltage curves for this channel and quantitatively explains their cation selectivity and rectification. We examined a series of energy-minimized models representing different aggregation states, side-chain rotamers, and helical rotations, as well as an ensemble of structures from a molecular dynamics trajectory. Potential energies were computed for single, permeating K+ and Cl- ions at a series of positions along a central pathway through the models. A variable-electric-field Nernst-Planck electrodiffusion model was used, with two adjustable parameters representing the diffusion coefficients of K+ and Cl- to scale the individual ion current magnitudes. The ability of a given DelPhi potential profile to fit the experimental data depended strongly on the magnitude of the desolvation of the permeating ion. Below a pore radius of 3.8 A, the predicted profiles showed large energy barriers, and the experimental data could be fit only with unrealistically high values for the K+ and Cl- diffusion coefficients. For pore radii above 3.8 A, the desolvation energies were 2kT or less. The electrostatic calculations were sensitive to positioning of the Ser side chains, with the best fits associated with maximum exposure of the Ser side-chain hydroxyls to the pore. The backbone component was shown to be the major source of asymmetry in the DelPhi potential profiles. Only two of the energy-minimized structures were able to explain the experimental data, whereas an average of the dynamics structures gave excellent agreement with experimental results. Thus this method provides a promising approach to prediction of current-voltage curves from three-dimensional structures of ion channel proteins.     

Ab initio molecular orbital study of the interaction of Li, Na and K with the pore components of ion channels: Consideration of the size, structure and selectivity of the pore of the channels

Ab initio molecular orbital calculations were made for the various types of structures of the pore of the ion channel and the results were applied to the permeability model by Hille, an extension of the Eyring rate theory. In Hille's model, ion passage through the channel is regarded as a kinetic process. Accordingly, it is thought that the interaction energy between cation and ligand, the easier the passage is, due to the lower activation energy. The calculated interaction energy was in the order Li+ greater than Na+ greater than K+ for all models. The optimum size of the pore determined from the interaction energy depends on the structure of the filter. The size for the pentagon was largest, followed by the hexagon and tetragon. On the other hand, the size depends hardly at all on the kind of ligand molecules. In the case of the tetragon, the sizes for the Na and K channels were nearly the same as those estimated from the model building and inhibitor-blocking experiment. The interaction energy between the ionized carboxyl group and the cation was extremely large, clearly reflecting the experimental fact that the carboxyl group in the pore has an important role in making the passage of the cation through the channel easier by dehydrating the water molecules. By analysis of the interaction energy, it was revealed that the contribution of the electrostatic energy was predominant, although the contributions of the other effects might not be negligible. Among these effects, the value of the charge transfer energy is largest, and this is noteworthy in connection with the selective transmission of cations through the overlap of orbitals. It is concluded that the quantum-chemical indices such as interaction energy and the electronic charge calculated by the sophisticated ab initio method help to shed light on the nature of the pore of the ion channel.     

Reconstitution of a chloroplast protein import channel.

The chloroplastic outer envelope protein OEP75 with a molecular weight of 75 kDa probably forms the central pore of the protein import machinery of the outer chloroplastic membrane. Patch-clamp analysis shows that heterologously expressed, purified and reconstituted OEP75 constitutes a voltage-gated ion channel with a unit conductance of Lambda = 145pS. Activation of the OEP75 channel in vitro is completely dependent on the magnitude and direction of the voltage gradient. Therefore, movements of protein charges of parts of OEP75 in the membrane electric field are required either for pore formation or its opening. In the presence of precursor protein from only one side of the bilayer, strong flickering and partial closing of the channel was observed, indicating a specific interaction of the precursor with OEP75. The comparatively low ionic conductance of OEP75 is compatible with a rather narrow aqueous pore (dporeapproximately equal to 8-9 A). Provided that protein and ion translocation occur through the same pore, this implies that the environment of the polypeptide during the transit is mainly hydrophilic and that protein translocation requires almost complete unfolding of the precursor.     

Molecular dissection of Cl--selective Cys-loop receptor points to components that are dispensable or essential for channel activity

Cys-loop receptors are pentameric ligand-gated ion channels (pLGICs) that bind neurotransmitters to open an intrinsic transmembrane ion channel pore. The recent crystal structure of a prokaryotic pLGIC from the cyanobacterium Gloeobacter violaceus (GLIC) revealed that it naturally lacks an N-terminal extracellular α helix and an intracellular domain that are typical of eukaryotic pLGICs. GLIC does not respond to neurotransmitters acting at eukaryotic pLGICs but is activated by protons. To determine whether the structural differences account for functional differences, we used a eukaryotic chimeric acetylcholine-glutamate pLGIC that was modified to carry deletions corresponding to the sequences missing in the prokaryotic homolog GLIC. Deletions made in the N-terminal extracellular α helix did not prevent the expression of receptor subunits and the appearance of receptor assemblies on the cell surface but abolished the capability of the receptor to bind α-bungarotoxin (a competitive antagonist) and to respond to the neurotransmitter. Other truncated chimeric receptors that lacked the intracellular domain did bind ligands; displayed robust acetylcholine-elicited responses; and shared with the full-length chimeric receptor similar anionic selectivity, effective open pore diameter, and unitary conductance. We suggest that the integrity of the N-terminal α helix is crucial for ligand accommodation because it stabilizes the intersubunit interfaces adjacent to the neurotransmitter-binding pocket(s). We also conclude that the intracellular domain of the chimeric acetylcholine-glutamate receptor does not modulate the ion channel conductance and is not involved in positioning of the pore-lining helices in the conformation necessary for coordinating a Cl- ion within the intracellular vestibule of the ion channel pore.     

The SecY complex forms a channel capable of ionic discrimination

Protein translocation across the bacterial membrane occurs at the SecY complex or channel. The resting SecY channel is impermeable to small molecules owing to a plug domain that creates a seal. Here, we report that a channel loosely sealed, or with a plug locked open, does not, however, lead to general membrane permeability. Instead, strong selectivity towards small monovalent anions, especially chloride, is observed. Mutations in the pore ring‐structure increase both the translocation activity of the channel and its ionic conductance, however the selectivity is maintained. The same ionic specificity also occurs at the onset of protein translocation and across the archaeal SecY complex. Thus, the ion‐conducting characteristic of the channel seems to be conserved as a normal consequence of protein translocation. We propose that the pore ring‐structure forms a selectivity filter, allowing cells to tolerate channels with imperfect plugs.     

Molecular dynamics simulations of water within models of ion channels.

The transbilayer pores formed by ion channel proteins contain extended columns of water molecules. The dynamic properties of such waters have been suggested to differ from those of water in its bulk state. Molecular dynamics simulations of ion channel models solvated within and at the mouths of their pores are used to investigate the dynamics and structure of intra-pore water. Three classes of channel model are investigated: a) parallel bundles of hydrophobic (Ala20) alpha-helices; b) eight-stranded hydrophobic (Ala10) antiparallel beta-barrels; and c) parallel bundles of amphipathic alpha-helices (namely, delta-toxin, alamethicin, and nicotinic acetylcholine receptor M2 helix). The self-diffusion coefficients of water molecules within the pores are reduced significantly relative to bulk water in all of the models. Water rotational reorientation rates are also reduced within the pores, particularly in those pores formed by alpha-helix bundles. In the narrowest pore (that of the Ala20 pentameric helix bundle) self-diffusion coefficients and reorientation rates of intra-pore waters are reduced by approximately an order of magnitude relative to bulk solvent. In Ala20 helix bundles the water dipoles orient antiparallel to the helix dipoles. Such dipole/dipole interaction between water and pore may explain how water-filled ion channels may be formed by hydrophobic helices. In the bundles of amphipathic helices the orientation of water dipoles is modulated by the presence of charged side chains. No preferential orientation of water dipoles relative to the pore axis is observed in the hydrophobic beta-barrel models.     

Phosphorylation mediated structural and functional changes in pentameric ligand-gated ion channels: Implications for drug discovery

Pentameric ligand-gated ion channels (pLGICs) mediate numerous physiological processes, including fast neurotransmission in the brain. They are targeted by a large number of clinically-important drugs and disruptions to their function are associated with many neurological disorders. The phosphorylation of pLGICs can result in a wide range of functional consequences. Indeed, many neurological disorders result from pLGIC phosphorylation. For example, chronic pain is caused by the protein kinase A-mediated phosphorylation of α3 glycine receptors and nicotine addiction is mediated by the phosphorylation of α4- or α7-containing nicotinic receptors. A recent study demonstrated that phosphorylation can induce a global conformational change in a pLGIC that propagates to the neurotransmitter-binding site. Here we present evidence that phosphorylation-induced global conformational changes may be a universal phenomenon in pLGICs. This raises the possibility of designing drugs to specifically treat disease-modified pLGICs. This review summarizes some of the opportunities available in this area.     

The proton ladder,a static mechanism for ion/proton coports and counterports

An ion/proton counterport is formed simply by locating a chain of ionizable residues connected by a proton conducting path near a passive ion pore which spans the membrane. The electric coupling between the ion in transit through the pore and the residues can ensure that for each ion passing through the pore in one direction a proton is driven along the chain of ionizable residues (the proton ladder) in the same or in the opposite direction. The mechanism is symmetrical in that a trans-membrane ion gradient may drive protons against their electrochemical potential gradient or a proton gradient may drive ions against theirs. The mechanism is applicable to cation or anion channels and to coports or counterports. No mechanical motion is required other than the motion of the ions and the protons. Monte Carlo computer simulations are performed on the model and its predicted properties are listed. The new type of counterport model is compared with currently used models. Offprint requests to: D. T Edmonds     

Arabidopsis thaliana glutamate receptor ion channel function demonstrated by ion pore transplantation

Ionotropic glutamate receptors (iGluRs) are ligand-gated cation channels that mediate fast excitatory neurotransmission in the mammalian central nervous system. In the model plant Arabidopsis thaliana, a large family of 20 genes encoding proteins that share similarities with animal iGluRs in sequence and predicted secondary structure has been discovered. Members of this family, termed AtGLRs (A. thaliana glutamate receptors), have been implicated in root development, ion transport, and several metabolic and signalling pathways. However, there is still no direct proof of ligand-gated ion channel function of any AtGLR subunit. We used a domain transplantation technique to directly test whether the putative ion pore domains of AtGLRs can conduct ions. To this end, we transplanted the ion pore domains of 17 AtGLR subunits into rat α-amino-3-hydroxy-5-methyl-4-isoxazole propionate (GluR1) and kainate (GluR6) receptor subunits and tested the resulting chimaeras for ion channel function in the Xenopus oocyte expression system. We show that AtGLR1.1 and AtGLR1.4 have functional Na⁺-, K⁺-, and Ca²⁺-permeable ion pore domains. The properties of currents through the AtGLR1.1 ion pore match those of glutamate-activated currents, depolarisations, and glutamate-triggered Ca²⁺ influxes observed in plant cells. We conclude that some AtGLRs have functional non-selective cation pores.     

Extended dipolar chain model for ion channels: electrostriction effects and the translocational energy barrier.

We reinvestigate the dipolar chain model for an ion channel. Our goal is to account for the influence that ion-induced electrostriction of channel water has on the translocational energy barriers experienced by different ions in the channel. For this purpose, we refine our former model by relaxing the positional constraint on the ion and the water dipoles and by including Lennard-Jones contributions in addition to the electrostatic interactions. The positions of the ion and the waters are established by minimization of the free energy. As before, interaction with the external medium is described via the image forces. Application to alkali cations show that the short range interactions modulate the free energy profiles leading to a selectivity sequence for translocation. We study the influence of some structural parameters on this sequence and compare our theoretical predictions with observed results for gramicidin.     

Unbiased simulations reveal the inward-facing conformation of the human serotonin transporter and Na(+) ion release

Monoamine transporters are responsible for termination of synaptic signaling and are involved in depression, control of appetite, and anxiety amongst other neurological processes. Despite extensive efforts, the structures of the monoamine transporters and the transport mechanism of ions and substrates are still largely unknown. Structural knowledge of the human serotonin transporter (hSERT) is much awaited for understanding the mechanistic details of substrate translocation and binding of antidepressants and drugs of abuse. The publication of the crystal structure of the homologous leucine transporter has resulted in homology models of the monoamine transporters. Here we present extended molecular dynamics simulations of an experimentally supported homology model of hSERT with and without the natural substrate yielding a total of more than 1.5 μs of simulation of the protein dimer. The simulations reveal a transition of hSERT from an outward-facing occluded conformation to an inward-facing conformation in a one-substrate-bound state. Simulations with a second substrate in the proposed symport effector site did not lead to conformational changes associated with translocation. The central substrate binding site becomes fully exposed to the cytoplasm leaving both the Na(+)-ion in the Na2-site and the substrate in direct contact with the cytoplasm through water interactions. The simulations reveal how sodium is released and show indications of early events of substrate transport. The notion that ion dissociation from the Na2-site drives translocation is supported by experimental studies of a Na2-site mutant. Transmembrane helices (TMs) 1 and 6 are identified as the helices involved in the largest movements during transport.     

Differences in hydration structure near hydrophobic and hydrophilic amino acids.

We use molecular dynamics to simulate recent neutron scattering experiments on aqueous solutions of N-acetyl-leucine-amide and N-acetyl-glutamine-amide, and break down the total scattering function into contributions from solute-solute, solute-water, water-water, and intramolecular correlations. We show that the shift of the main diffraction peak to smaller angle that is observed for leucine, but not for glutamine, is attributable primarily to alterations in water-water correlations relative to bulk. The perturbation of the water hydrogen-bonded network extends roughly two solvation layers from the hydrophobic side chain surface, and is characterized by a distribution of hydrogen bonded ring sizes that are more planar and are dominated by pentagons in particular than those near the hydrophilic side chain. The different structural organization of water near the hydrophobic solute that gives rise to the inward shift in the main neutron diffraction peak under ambient conditions may also provide insight into the same directional shift for pure liquid water as it is cooled and supercooled.