Probing the pore of ClC-0 by substituted cysteine accessibility method using methane thiosulfonate reagents

ClC channels are a family of protein molecules containing two ion-permeation pores. Although these transmembrane proteins are important for a variety of physiological functions, their molecular operations are only superficially understood. High-resolution X-ray crystallography techniques have recently revealed the structures of two bacterial ClC channels, but whether vertebrate ClC channel pores are similar to those of bacterial homologues is not clear. To study the pore architecture of the Torpedo ClC-0 channel, we employed the substituted-cysteine-accessibility method (SCAM) and used charged methane thiosulfonate (MTS) compounds to modify the introduced cysteine. Several conclusions were derived from this approach. First, the MTS modification pattern from Y512C to E526C in ClC-0, which corresponds to residues forming helix R in bacterial ClC channels, is indeed consistent with the suggested helical structure. Second, the ClC-0 pore is more accessible to the negatively charged than to the positively charged MTS compound, a pore property that is regulated by the intrinsic electrostatic potential in the pore. Finally, attempts to modify the introduced cysteine at positions intracellular to the selectivity filter did not result in larger MTS modification rates for the open-state channel, suggesting that the fast gate of ClC-0 cannot be located at a position intracellular to the Cl- selectivity filter. Thus, the proposal that the glutamate side chain is the fast gate of the channel is applicable to ClC-0, revealing a structural and functional conservation of ClC channels between bacterial and vertebrate species.     

A BEST example of channel structure annotation by molecular simulation

An increasing number of ion channel structures are being determined. This generates a need for computational tools to enable functional annotation of channel structures. However, several studies of ion channel and model pores have indicated that the physical dimensions of a pore are not always a reliable indicator of its conductive status. This is due to the unusual behavior of water within nano-confined spaces, resulting in a phenomenon referred to as “hydrophobic gating”. We have recently demonstrated how simulating the behavior of water within an ion channel pore can be used to predict its conductive status. In this addendum to our study, we apply this method to compare the recently solved structure of a mutant of the bestrophin chloride channel BEST1 with that of the wild-type channel. Our results support the hypothesis of a hydrophobic gate within the narrow neck of BEST1. This provides further validation that this simulation approach provides the basis for an accurate and computationally efficient tool for the functional annotation of ion channel structures.     

Analysis of the selectivity filter of the voltage-gated sodium channel NavRh

NaChBac is a bacterial voltage-gated sodium (Na_v) channel that shows sequence similarity to voltage-gated calcium channels. To understand the ion-permeation mechanism of Na_v channels, we combined molecular dynamics simulation, structural biology and electrophysiological approaches to investigate the recently determined structure of Na_vRh, a marine bacterial NaChBac ortholog. Two Na⁺ binding sites are identified in the selectivity filter (SF) in our simulations: The extracellular Na⁺ ion first approaches site 1 constituted by the side groups of Ser181 and Glu183, and then spontaneously arrives at the energetically more favorable site 2 formed by the carbonyl oxygens of Leu179 and Thr178. In contrast, Ca²⁺ ions are prone to being trapped by Glu183 at site 1, which then blocks the entrance of both Na⁺ and Ca²⁺ to the vestibule of the SF. In addition, Na⁺ permeates through the selective filter in an asymmetrical manner, a feature that resembles that of the mammalian Na_v orthologs. The study reported here provides insights into the mechanism of ion selectivity on Na⁺ over Ca²⁺ in mammalian Na_v channels.     

Pore Hydration States of KcsA Potassium Channels in Membranes

Water-filled hydrophobic cavities in channel proteins serve as gateways for transfer of ions across membranes, but their properties are largely unknown. We determined water distributions along the conduction pores in two tetrameric channels embedded in lipid bilayers using neutron diffraction: potassium channel KcsA and the transmembrane domain of M2 protein of influenza A virus. For the KcsA channel in the closed state, the distribution of water is peaked in the middle of the membrane, showing water in the central cavity adjacent to the selectivity filter. This water is displaced by the channel blocker tetrabutyl-ammonium. The amount of water associated with the channel was quantified, using neutron diffraction and solid state NMR. In contrast, the M2 proton channel shows a V-shaped water profile across the membrane, with a narrow constriction at the center, like the hourglass shape of its internal surface. These two types of water distribution are therefore very different in their connectivity to the bulk water. The water and protein profiles determined here provide important evidence concerning conformation and hydration of channels in membranes and the potential role of pore hydration in channel gating.     

Structural basis of pore formation by the bacterial toxin pneumolysin

The bacterial toxin pneumolysin is released as a soluble monomer that kills target cells by assembling into large oligomeric rings and forming pores in cholesterol-containing membranes. Using cryo-EM and image processing, we have determined the structures of membrane-surface bound (prepore) and inserted-pore oligomer forms, providing a direct observation of the conformational transition into the pore form of a cholesterol-dependent cytolysin. In the pore structure, the domains of the monomer separate and double over into an arch, forming a wall sealing the bilayer around the pore. This transformation is accomplished by substantial refolding of two of the four protein domains along with deformation of the membrane. Extension of protein density into the bilayer supports earlier predictions that the protein inserts beta hairpins into the membrane. With an oligomer size of up to 44 subunits in the pore, this assembly creates a transmembrane channel 260 A in diameter lined by 176 beta strands.     

Structure of the SthK Carboxy-Terminal Region Reveals a Gating Mechanism for Cyclic Nucleotide-Modulated Ion Channels

Cyclic nucleotide-sensitive ion channels are molecular pores that open in response to cAMP or cGMP, which are universal second messengers. Binding of a cyclic nucleotide to the carboxyterminal cyclic nucleotide binding domain (CNBD) of these channels is thought to cause a conformational change that promotes channel opening. The C-linker domain, which connects the channel pore to this CNBD, plays an important role in coupling ligand binding to channel opening. Current structural insight into this mechanism mainly derives from X-ray crystal structures of the C-linker/CNBD from hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels. However, these structures reveal little to no conformational changes upon comparison of the ligand-bound and unbound form. In this study, we take advantage of a recently identified prokaryote ion channel, SthK, which has functional properties that strongly resemble cyclic nucleotide-gated (CNG) channels and is activated by cAMP, but not by cGMP. We determined X-ray crystal structures of the C-linker/CNBD of SthK in the presence of cAMP or cGMP. We observe that the structure in complex with cGMP, which is an antagonist, is similar to previously determined HCN channel structures. In contrast, the structure in complex with cAMP, which is an agonist, is in a more open conformation. We observe that the CNBD makes an outward swinging movement, which is accompanied by an opening of the C-linker. This conformation mirrors the open gate structures of the K_v1.2 channel or MthK channel, which suggests that the cAMP-bound C-linker/CNBD from SthK represents an activated conformation. These results provide a structural framework for better understanding cyclic nucleotide modulation of ion channels, including HCN and CNG channels.     

Cumulative hydropathic topology of a voltage-gated sodium channel at atomic resolution

Voltage-gated sodium channels (NavChs) are biological pores that control the flow of sodium ions through the cell membrane. In humans, mutations in genes encoding NavChs can disrupt physiological cellular activity thus leading to a wide spectrum of diseases. Here, we present a topological connection between the functional architecture of a NavAb bacterial channel and accumulation of atomic hydropathicity around its pore. This connection is established via a scaling analysis methodology that elucidates how intrachannel hydropathic density variations translate into hydropathic dipole field configurations along the pore. Our findings suggest the existence of a nonrandom cumulative hydropathic topology that is organized parallel to the membrane surface so that pore's stability, as well as, gating behavior are guaranteed. Given the biophysical significance of the hydropathic effect, our study seeks to provide a computational framework for studying cumulative hydropathic topological properties of NavChs and pore-forming proteins in general.     

Potassium, sodium, calcium and glutamate-gated channels: pore architecture and ligand action

In the last decade, the idea of common organization of certain ion channel families exhibiting diverse physiological and pharmacological properties has received strong experimental support. Transmembrane topologies and patterns of the pore-facing residues are conserved in P-loop channels that include high-selective cation channels and certain ligand-gated channels. X-ray structures of bacterial K+ channels, KcsA, MthK and KvAP, help to understand structure-function relationships of other P-loop channels. Data on binding sites and mechanisms of action of ligands of K+, Na+, Ca2+ and glutamate gated ion channels are considered in view of their possible structural similarity to the bacterial K+ channels. Emphasized are structural determinants of ligand-receptor interactions within the channels and mechanisms of state-dependent action of the ligands.     

Intracellular regions of potassium channels: Kv2.1 and heag

Intracellular regions of voltage-gated potassium channels often comprise the largest part of the channel protein, and yet the functional role of these regions is not fully understood. For the Kv2.1 channel, although there are differences in activation kinetics between rat and human channels, there are, for instance, no differences in movement of the S4 region between the two channels, and indeed our mutagenesis studies have identified interacting residues in both the N- and C -terminal intracellular regions that are responsible for these functional effects. Furthermore, using FRET with fluorescent-tagged Kv2.1 channels, we have shown movement of the C-termini relative to the N-termini during activation. Such interactions and movements of the intracellular regions of the channel appear to form part of the channel gating machinery. Heag1 and heag2 channels also display differing activation properties, despite their considerable homology. By a chimeric approach, we have shown that these differences in activation kinetics are determined by multiple interacting regions in the N-terminus and membrane-spanning regions. Furthermore, alanine mutations of many residues in the C-terminal cyclic nucleotide binding domain affect activation kinetics. The data again suggest interacting regions between N- and C- termini that participate in the conformational changes during channel activation. Using a mass-spectrometry approach, we have identified α-tubulin and a heat shock protein as binding to the C-terminus of the heag2 channel, and α-tubulin itself has functional effects on channel activation kinetics. Clearly, the intracellular regions of these ion channels (and most likely many other ion channels too) are important regions in determining channel function. EBSA Satellite Meeting: Ion channels, Leeds, July 2007.     

Structural insights into the gating mechanisms of TRPV channels

Transient Receptor Potential channels from the vanilloid subfamily (TRPV) are a group of cation channels modulated by a variety of endogenous stimuli as well as a range of natural and synthetic compounds. Their roles in human health make them of keen interest, particularly from a pharmacological perspective. However, despite this interest, the complexity of these channels has made it difficult to obtain high resolution structures until recently. With the cryo-EM resolution revolution, TRPV channel structural biology has blossomed to produce dozens of structures, covering every TRPV family member and a variety of approaches to examining channel modulation. Here, we review all currently available TRPV structures and the mechanistic insights into gating that they reveal.     

Mechanisms of up-regulation of single calcium channels by serotonin in Helix pomatia neurons

Action of serotonin (5-HT) on single Ca(2+) channel activity was studied in identified neurons of snail Helix pomatia. Only one type of Ca(2+) channels of 5 pS unitary conductance was determined under patch-clamp cell-attached mode. Kinetic analysis have shown a monotonically declining distribution of channel open times (OT) with mean time constant of 0.2 ms. The distribution of channel closed times (CT) could be fitted by double-exponential curve with time constants 1 and 12 ms. We established that 5-HT acts on Ca(2+) channel activity indirectly via cytoplasm. 5-HT prolonged the OT (up to 0.3 ms) and shortened the CT proportionally for both constants to 0.4 and 6 ms correspondingly. A conclusion is made that enhancement of Ca(2+) macro-current by 5-HT is determined by kinetic changes, increase of the number of active channels, and increase of the probability of OT. At the same time the transmitter did not affect the unitary channel conductance.     

Control of Volume and Turgor in Stomatal Guard Cells

Water loss from plants is determined by the aperture of stomatal pores in the leaf epidermis, set by the level of vacuolar accumulation of potassium salt, and hence volume and turgor, of a pair of guard cells. Regulation of ion fluxes across the tonoplast, the key to regulation of stomatal aperture, can only be studied by tracer flux measurements. There are two transport systems in the tonoplast. The first is a Ca²⁺-activated channel, inhibited by phenylarsine oxide (PAO), responsible for the release of vacuolar K⁺(Rb⁺) in response to the “drought” hormone, abscisic acid (ABA). This channel is sensitive to pressure, down-regulated at low turgor and up-regulated at high turgor, providing a system for turgor regulation. ABA induces a transient stimulation of vacuolar ion efflux, during which the flux tracks the ion content (volume, turgor), suggesting ABA reduces the set-point of a control system. The second system, which is PAO-insensitive, is responsible for an ion flux from vacuole to cytoplasm associated with inward water flow following a hypo-osmotic transfer. It is suggested that this involves an aquaporin as sensor, and perhaps also as responder; deformation of the aquaporin may render it ion-permeable, or, alternatively, the deformed aquaporin may signal to an associated ion channel, activating it. Treatment with inhibitors of aquaporins, HgCl₂ or silver sulfadiazine, produces a large transient increase in ion release from the vacuole, also PAO-insensitive. It is suggested that this involves the same aquaporin, either rendered directly ion-permeable, or signalling to activate an associated ion channel.     

A structural overview of the ion channels of the TRPM family

The TRPM (transient receptor potential melastatin) family belongs to the superfamily of TRP cation channels. The TRPM subfamily is composed of eight members that are involved in diverse biological functions such as temperature sensing, inflammation, insulin secretion, and redox sensing. Since the first cloning of TRPM1 in 1998, tremendous progress has been made uncovering the function, structure, and pharmacology of this family. Complete structures of TRPM2, TRPM4, and TRPM8, as well as a partial structure of TRPM7, have been determined by cryo-EM, providing insights into their channel assembly, ion permeation, gating mechanisms, and structural pharmacology. Here we summarize the current knowledge about channel structure, emphasizing general features and principles of the structure of TRPM channels discovered since 2017. We also discuss some of the key unresolved issues in the field, including the molecular mechanisms underlying voltage and temperature dependence, as well as the functions of the TRPM channels’ C-terminal domains.     

Location of a constriction in the lumen of a transmembrane pore by targeted covalent attachment of polymer molecules

Few methods exist for obtaining the internal dimensions of transmembrane pores for which 3-D structures are lacking or for showing that structures determined by crystallography reflect the internal dimensions of pores in lipid bilayers. Several approaches, involving polymer penetration and transport, have revealed limiting diameters for various pores. But, in general, these approaches do not indicate the locations of constrictions in the channel lumen. Here, we combine cysteine mutagenesis and chemical modification with sulfhydryl-reactive polymers to locate the constriction in the lumen of the staphylococcal alpha-hemolysin pore, a model protein of known structure. The rates of reaction of each of four polymeric reagents (MePEG-OPSS) of different masses towards individual single cysteine mutants, comprising a set with cysteines distributed over the length of the lumen of the pore, were determined by macroscopic current recording. The rates for the three larger polymers (1.8, 2.5, and 5.0 kD) were normalized with respect to the rates of reaction with a 1.0-kD polymer for each of the seven positions in the lumen. The rate of reaction of the 5.0-kD polymer dropped dramatically at the centrally located Cys-111 residue and positions distal to Cys-111, whether the reagent was applied from the trans or the cis side of the bilayer. This semi-quantitative analysis sufficed to demonstrate that a constriction is located at the midpoint of the pore lumen, as predicted by the crystal structure, and although the constriction allows a 2.5-kD polymer to pass, transport of a 5.0-kD molecule is greatly restricted. In addition, PEG chains gave greater reductions in pore conductance when covalently attached to the narrower regions of the lumen, permitting further definition of the interior of the pore. The procedures described here should be applicable to other pores and to related structures such as the vestibules of ion channels.     

Proton pathways and H+/Cl- stoichiometry in bacterial chloride transporters

H+/Cl- antiport behavior has recently been observed in bacterial chloride channel homologs and eukaryotic CLC-family proteins. The detailed molecular-level mechanism driving the stoichiometric exchange is unknown. In the bacterial structure, experiments and modeling studies have identified two acidic residues, E148 and E203, as key sites along the proton pathway. The E148 residue is a major component of the fast gate, and it occupies a site crucial for both H+ and Cl- transport. E203 is located on the intracellular side of the protein; it is vital for H+, but not Cl-, transport. This suggests two independent ion transit pathways for H+ and Cl- on the intracellular side of the transporter. Previously, we utilized a new pore-searching algorithm, TransPath, to predict Cl- and H+ ion pathways in the bacterial ClC channel homolog, focusing on proton access from the extracellular solution. Here we employ the TransPath method and molecular dynamics simulations to explore H+ pathways linking E148 and E203 in the presence of Cl- ions located at the experimentally observed binding sites in the pore. A conclusion is that Cl- ions are required at both the intracellular (S(int)) and central (S(cen)) binding sites in order to create an electrostatically favorable H+ pathway linking E148 and E203; this electrostatic coupling is likely related to the observed 1H+/2Cl- stoichiometry of the antiporter. In addition, we suggest that a tyrosine residue side chain (Y445), located near the Cl- ion binding site at S(cen), is involved in proton transport between E148 and E203.     

Exploring blocker binding to a homology model of the open hERG K+ channel using docking and molecular dynamics methods

Binding of blockers to the human voltage-gated hERG potassium channel is studied using a combination of homology modelling, automated docking calculations and molecular dynamics simulations, where binding affinities are evaluated using the linear interaction energy method. A homology model was constructed based on the available crystal structure of the bacterial KvAP channel and the affinities of a series of sertindole analogues predicted using this model. The calculations reproduce the relative binding affinities of these compounds very well and indicate that both polar interactions near the intracellular opening of the selectivity filter as well as hydrophobic complementarity in the region around F656 are important for blocker binding. These results are consistent with recent alanine scanning mutation experiments on the blocking of the hERG channel by other compounds.     

Brownian dynamics simulations of the recognition of the scorpion toxin P05 with the small-conductance calcium-activated potassium channels

The recognition of the scorpion toxin P05 and the small-conductance, calcium-activated potassium (SK) channels, rsk1, rsk2, and rsk3, has been studied by means of the Brownian dynamics (BD) method. All of the 25 available structures of P05 in the RCSB Protein Data Bank determined by NMR were considered during the simulation, which indicated that the conformation of P05 affects both the recognition and binding between the two proteins significantly. Comparing the top four high-frequency structures of P05 binding to the SK channels, we found that the rsk2 channel, with high frequencies and lowest electrostatic interaction energies (E (int)(ES)), is the most favorable for P05 binding, while rsk3 is intermediate, and rsk1 is the least favorable. Among the 25 structures of P05, the 13th structure docks into the binding site of the rsk2 channel with the highest probability and most favorable electrostatic interactions. From the P05-rsk2 channel binding model, we identified the residues critical for the recognition of these two proteins through triplet contact analyses. P05 locates around the extracellular mouth of the SK channels and contacts the SK channels using its alpha-helix rather than beta-sheets. The critical triplet contacts for recognition between P05 and the rsk2 channel are Arg6 (P05)-Asp364 (SK), Arg7 (P05)-Asn368 (SK), and Arg13 (P05)-Asp341 (SK). The structure of the P05-rsk2 complex with the most favorable electrostatic interaction energy was further refined by molecular mechanics, showing that six hydrogen bonding interactions exist between P05 and the rsk2 channel: one hydrogen bond is formed between Arg6 (P05) and Asp364(D) (rsk2); Arg7 (P05) forms three hydrogen bonds with Asp341(B) (rsk2)) and Asp364(C) (rsk2); two hydrogen bonds are formed by Arg13 (P05) with Asp341(A) (rsk2) and Asp364(B) (rsk2). The simulation results are in good agreement with the previous molecular biological experiments and can explain the binding phenomena between P05 and SK channels at the level of molecular structure. The consistency between the results of the BD simulations and the experimental data indicated that our 3D model of the P05-rsk2 channel complex is reasonable and can be employed in further biological studies, such as rational design of the novel therapeutic agents blocking the small-conductance, calcium-activated and apamin-sensitive potassium channels, and for mutagenesis studies in both toxins and SK channels. In particular, both the BD simulations and the molecular mechanics refinements indicate that residue Asp364 of the rsk2 channel is critical for its recognition and binding functionality towards P05. This phenomenon has not been appreciated in the previous mutagenesis experiments, indicating that this might be a new clue for further functional study of SK channels.     

Bacterial ion channels and their eukaryotic homologues.

Due to the relative ease of obtaining their crystal structures, bacterial ion channels provide a unique opportunity to analyse structure and function of their eukaryotic homologues. This review describes prokaryotic channels whose structures have been determined. These channels are KcsA, a bacterial homologue of eukaryotic potassium channels, MscL, a bacterial mechanosensitive ion channel and ClC0, a prokaryotic homologue of the eukaryotic ClC family of anion-selective channels. General features of their structure and function are described with a special emphasis on the advantages that these channels offer for understanding the properties of their eukaryotic homologues. We present amino-acid sequences of eukaryotic proteins related in their primary sequences to bacterial mechanosensitive channels. The usefulness of bacterial mechanosensitive channels for the studies on general principles of mechanosensation is discussed.