Two rings of negative charges in the cytosolic vestibule of type-1 ryanodine receptor modulate ion fluxes

The tetrameric ryanodine receptor calcium release channels (RyRs) are cation-selective channels that have pore architecture similar to that of K+ channels. We recently identified, in close proximity to the selectivity filter motif GGGIG, a conserved lumenal DE motif that has a critical role in RyR ion permeation and selectivity. Here, we substituted three aspartate residues (D4938, D4945, D4953) with asparagine and four glutamate residues (E4942, E4948, E4952, E4955) with glutamine hypothesized to line the cytosolic vestibule of the skeletal muscle RyR (RyR1). Mutant single channel properties were determined using the planar lipid bilayer method. Two mutants (D4938N, D4945N) showed a reduced K+ ion conductance, with D4938N also exhibiting a reduced selectivity for Ca2+ compared to K+. The cytosolic location of D4938 and D4945 was confirmed using the polycation neomycin. Both D4938N and D4945N exhibited an attenuated block by neomycin to a greater extent from the cytosolic than lumenal side. By comparison, charge neutralization of lumenal loop residues (D4899Q, E4900N) eliminated the block from the lumenal but not the cytosolic side. The results suggest that, in addition to negatively charged residues on the lumenal side, rings of four negative charges formed by D4938 and D4945 in the cytosolic vestibule determine RyR ion fluxes.     

Permeation of large tetra-alkylammonium cations through mutant and wild-type voltage-gated sodium channels as revealed by relief of block at high voltage

Many large organic cations are potent blockers of K(+) channels and other cation-selective channels belonging to the P-region superfamily. However, the mechanism by which large hydrophobic cations enter and exit the narrow pores of these proteins is obscure. Previous work has shown that a conserved Lys residue in the DEKA locus of voltage-gated Na(+) channels is an important determinant of Na(+)/K(+) discrimination, exclusion of Ca(2+), and molecular sieving of organic cations. In this study, we sought to determine whether the Lys(III) residue of the DEKA locus interacts with internal tetra-alkylammonium cations (TAA(+)) that block Na(+) channels in a voltage-dependent fashion. We investigated block by a series of TAA(+) cations of the wild-type rat muscle Na(+) channel (DEKA) and two different mutants of the DEKA locus, DEAA and DERA, using whole-cell recording. TEA(+) and larger TAA(+) cations block both wild-type and DEAA channels. However, DEAA exhibits dramatic relief of block by large TAA(+) cations as revealed by a positive inflection in the macroscopic I-V curve at voltages greater than +140 mV. Paradoxically, relief of block at high positive voltage is observed for large (e.g., tetrapentylammonium) but not small (e.g., TEA(+)) symmetrical TAA(+) cations. The DEKA wild-type channel and the DERA mutant exhibit a similar relief-of-block phenomenon superimposed on background current rectification. The results indicate: (a) hydrophobic TAA(+) cations with a molecular diameter as large as 15 A can permeate Na(+) channels from inside to outside when driven by high positive voltage, and (b) the Lys(III) residue of the DEKA locus is an important determinant of inward rectification and internal block in Na(+) channels. From these observations, we suggest that hydrophobic interfaces between subunits, pseudosubunits, or packed helices of P-region channel proteins may function in facilitating blocker access to the pore, and may thus play an important role in the blocking and permeation behavior of large TAA(+) cations and potentially other kinds of local anesthetic molecules.     

Pore Dynamics and Conductance of RyR1 Transmembrane Domain

Ryanodine receptors (RyR) are calcium release channels, playing a major role in the regulation of muscular contraction. Mutations in skeletal muscle RyR (RyR1) are associated with congenital diseases such as malignant hyperthermia and central core disease (CCD). The absence of high-resolution structures of RyR1 has limited our understanding of channel function and disease mechanisms at the molecular level. Previously, we have reported a hypothetical structure of the RyR1 pore-forming region, obtained by homology modeling and supported by mutational scans, electrophysiological measurements, and cryo-electron microscopy. Here, we utilize the expanded model encompassing six transmembrane helices to calculate the RyR1 pore region conductance, to analyze its structural stability, and to hypothesize the mechanism of the Ile4897 CCD-associated mutation. The calculated conductance of the wild-type RyR1 suggests that the proposed pore structure can sustain ion currents measured in single-channel experiments. We observe a stable pore structure on timescales of 0.2 μs, with multiple cations occupying the selectivity filter and cytosolic vestibule, but not the inner chamber. We further suggest that stability of the selectivity filter critically depends on the interactions between the I4897 residue and several hydrophobic residues of the neighboring subunit. Loss of these interactions in the case of polar substitution I4897T results in destabilization of the selectivity filter, a possible cause of the CCD-specific reduced Ca²⁺ conductance.     

State-dependent block of BK channels by synthesized shaker ball peptides

Crystal structures of potassium channels have strongly corroborated an earlier hypothetical picture based on functional studies, in which the channel gate was located on the cytoplasmic side of the pore. However, accessibility studies on several types of ligand-sensitive K(+) channels have suggested that their activation gates may be located near or within the selectivity filter instead. It remains to be determined to what extent the physical location of the gate is conserved across the large K(+) channel family. Direct evidence about the location of the gate in large conductance calcium-activated K(+) (BK) channels, which are gated by both voltage and ligand (calcium), has been scarce. Our earlier kinetic measurements of the block of BK channels by internal quaternary ammonium ions have raised the possibility that they may lack a cytoplasmic gate. We show in this study that a synthesized Shaker ball peptide (ShBP) homologue acts as a state-dependent blocker for BK channels when applied internally, suggesting a widening at the intracellular end of the channel pore upon gating. This is consistent with a gating-related conformational change at the cytoplasmic end of the pore-lining helices, as suggested by previous functional and structural studies on other K(+) channels. Furthermore, our results from two BK channel mutations demonstrate that similar types of interactions between ball peptides and channels are shared by BK and other K(+) channel types.     

Effect of sequence hydrophobicity and bilayer width upon the minimum length required for the formation of transmembrane helices in membranes

The minimum hydrophobic length necessary to form a transmembrane (TM) helix in membranes was investigated using model membrane-inserted hydrophobic helices. The fluorescence of a Trp at the center of the sequence and its sensitivity to quenching were used to ascertain helix position within the membrane. Peptides with hydrophobic cores composed of poly(Leu) were compared to sequences containing a poly 1:1 Leu:Ala core (which have a hydrophobicity typical of natural TM helices). Studies varying bilayer width revealed that the poly(Leu) core peptides predominately formed a TM state when the bilayer width exceeded hydrophobic sequence length by (i.e. when negative mismatch was) up to ∼ 11-12 Å (e.g. the case of a 11-12 residue hydrophobic sequence in bilayers with a biologically relevant width, i.e. dioleoylphosphatidylcholine (DOPC) bilayers), while poly(LeuAla) core peptides formed predominantly TM state with negative mismatch of up to 9 Å (a 13 residue hydrophobic sequence in DOPC bilayers). This indicates that minimum length necessary to form a predominating amount of a TM state (minimum TM length) is only modestly hydrophobicity-dependent for the sequences studied here, and a formula that defines the minimum TM length as a function of hydrophobicity for moderately-to-highly hydrophobic sequences was derived. The minimum length able to form a stable TM helix for alternating LeuAla sequences, and that for sequences with a Leu block followed by an Ala block, was similar, suggesting that a hydrophobicity gradient along the sequence may not be an important factor in TM stability. TM stability was also similar for sequences flanked by different charged ionizable residues (Lys, His, Asp). However, ionizable flanking residues destabilized the TM configuration much more when charged than when uncharged. The ability of short hydrophobic sequences to form TM helices in membranes in the presence of substantial negative mismatch implies that lipid bilayers have a considerable ability to adjust to negative mismatch, and that short TM helices may be more common than generally believed. Factors that modulate the ability of bilayers to adjust to mismatch may strongly affect the configuration of short hydrophobic helices.     

Pore Dynamics and Conductance of RyR1 Transmembrane Domain

Ryanodine receptors (RyR) are calcium release channels, playing a major role in the regulation of muscular contraction. Mutations in skeletal muscle RyR (RyR1) are associated with congenital diseases such as malignant hyperthermia and central core disease (CCD). The absence of high-resolution structures of RyR1 has limited our understanding of channel function and disease mechanisms at the molecular level. Previously, we have reported a hypothetical structure of the RyR1 pore-forming region, obtained by homology modeling and supported by mutational scans, electrophysiological measurements, and cryo-electron microscopy. Here, we utilize the expanded model encompassing six transmembrane helices to calculate the RyR1 pore region conductance, to analyze its structural stability, and to hypothesize the mechanism of the Ile4897 CCD-associated mutation. The calculated conductance of the wild-type RyR1 suggests that the proposed pore structure can sustain ion currents measured in single-channel experiments. We observe a stable pore structure on timescales of 0.2 μs, with multiple cations occupying the selectivity filter and cytosolic vestibule, but not the inner chamber. We further suggest that stability of the selectivity filter critically depends on the interactions between the I4897 residue and several hydrophobic residues of the neighboring subunit. Loss of these interactions in the case of polar substitution I4897T results in destabilization of the selectivity filter, a possible cause of the CCD-specific reduced Ca²⁺ conductance.     

Open-state conformation of the KcsA K+ channel: Monte Carlo normal mode following simulations

Potassium channels fluctuate between closed and open states. The detailed mechanism of the conformational changes opening the intracellular pore in the K+ channel from Streptomyces lividans (KcsA) is unknown. Applying Monte Carlo normal mode following, we find that gating involves rotation and unwinding of the TM2 bundle, lateral movement of the TM2 helices away from the channel axis, and disappearance of the TM2 bundle. The open-state conformation of KcsA exhibits a very wide inner vestibule, with a radius approximately 5-7 A and inner helices bent at the A98-G99 hinge. Computed conformational changes demonstrate that spin labeling and X-ray experiments illuminate different stages in gating: transition begins with clockwise rotation of the TM2 helices ending at a final state with the TM2 bend hinged near residues A98-G99. The concordance between the computational and experimental results provides atomic-level insights into the structural rearrangements of the channel's inner pore.     

A prokaryotic glutamate receptor: homology modelling and molecular dynamics simulations of GluR0

GluR0 is a prokaryotic homologue of mammalian glutamate receptors that forms glutamate-activated, potassium-selective ion channels. The topology of its transmembrane (TM) domain is similar to that of simple potassium channels such as KcsA. Two plausible alignments of the sequence of the TM domain of GluR0 with KcsA are possible, differing in the region of the P helix. We have constructed homology models based on both alignments and evaluated them using 6 ns duration molecular dynamics simulations in a membrane-mimetic environment. One model, in which an insertion in GluR0 relative to KcsA is located in the loop between the M1 and P helices, is preferred on the basis of lower structural drift and maintenance of the P helix conformation during simulation. This model also exhibits inter-subunit salt bridges that help to stabilise the TM domain tetramer. During the simulation, concerted K(+) ion-water movement along the selectivity filter is observed, as is the case in simulations of KcsA. K(+) ion exit from the central cavity is associated with opening of the hydrophobic gate formed by the C-termini of the M2 helices. In the intact receptor the opening of this gate will be controlled by interactions with the extramembranous ligand-binding domains.     

Pharmacology and Surface Electrostatics of the K Channel Outer Pore Vestibule

In spite of a generally well-conserved outer vestibule and pore structure, there is considerable diversity in the pharmacology of K channels. We have investigated the role of specific outer vestibule charged residues in the pharmacology of K channels using tetraethylammonium (TEA) and a trivalent TEA analog, gallamine. Similar to Shaker K channels, gallamine block of Kv3.1 channels was more sensitive to solution ionic strength than was TEA block, a result consistent with a contribution from an electrostatic potential near the blocking site. In contrast, TEA block of another type of K channel (Kv2.1) was insensitive to solution ionic strength and these channels were resistant to block by gallamine. Neutralizing either of two lysine residues in the outer vestibule of these Kv2.1 channels conferred ionic strength sensitivity to TEA block. Kv2.1 channels with both lysines neutralized were sensitive to block by gallamine, and the ionic strength dependence of this block was greater than that for TEA. These results demonstrate that Kv3.1 (like Shaker) channels contain negatively charged residues in the outer vestibule of the pore that influence quaternary ammonium pharmacology. The presence of specific lysine residues in wild-type Kv2.1 channels produces an outer vestibule with little or no net charge, with important consequences for quaternary ammonium block. Neutralizing these key lysines results in a negatively charged vestibule with pharmacological properties approaching those of other types of K channels.     

Hydrophobic organization of α-helix membrane bundle in bacteriorhodopsin

The hydrophobic organization of the intramembraneα-helical bundle in bacteriorhodopsin (BRh) was assessed based on a new approach to characterization of spatial hydrophobic properties of transmembrane (TM)α-helical peptides. The method employs two independent techniques: Monte Carlo simulations of nonpolar solvent around TM peptides and analysis of molecular hydrophobicity potential on their surfaces. The results obtained by the two methods agree with each other and permit precise hydrophobicity mapping of TM peptides. Superimposition of such data on the experimentally derived spatial model of the membrane moiety together with 2D maps of hydrophobic hydrophilic contacts provide considerable insight into the hydrophobic organization of BRh. The helix bundle is stabilized to a large extent by hydrophobic interactions between helices—neighbors in the sequence of BRh, by long-range interactions in helix pairs C-E, C-F, and C-G, and by nonpolar contracts between retinal and helices C, D, E, F. Unlike globular proteins, no polar contacts between residues distantly separated in the sequence of BRh were found in the bundle. One of the most striking results of this study is the finding that the hydrophobic organization of BRh is significantly different from those in bacterial photoreaction centers. Thus, TMα-helices in BRh expose their most nonpolar sides to the bilayer as well as to the neighboring helices and to the interior of the bundle. Some of them contact lipids with their relatively hydrophilic surfaces. No correlation was found between disposition of the most hydrophobic and the most variable sides of the TM helices.     

KcsA crystal structure as framework for a molecular model of the Na(+) channel pore

The crystal structure of the pore-forming part of the KcsA bacterial K(+)-selective channel suggests a possible motif for related voltage-gated channels. We examined the hypothesis that the spacial orientation of the KcsA M1 and M2 alpha-helices also predicts the backbone location of S5 and S6 helices of the voltage-gated Na(+) channel. That channel's P region structure is expected to be different because selectivity is determined by side-chain interactions rather than by main-chain carbonyls, and its outer vestibule accommodates relatively large toxin molecules, tetrodotoxin (TTX) and saxitoxin (STX), which interact with selectivity ring residues. The Na(+) channel P loop was well-modeled by the alpha-helix-turn-beta-strand motif, which preserves the relationships for toxin interaction with the Na(+) channel found experimentally. This outer vestibule was docked into the extracellular part of the inverted teepee structure formed by the S5 and S6 helices that were spacially located by coordinates of the KcsA M1 and M2 helix main chains [Doyle et al. (1998) Science 280, 69-74], but populated with side chains of the respective S5 and S6 structures. van der Waals contacts were optimized with minimal adjustment of the S5, S6, and P loop structures, forming a densely packed pore structure. Nonregular external S5-P and P-S6 segments were not modeled here, except the P-S6 segment of domain II. The resulting selectivity region structure is consistent with Na(+) channel permeation properties, offering suggestions for the molecular processes involved in selectivity. The ability to construct a Na(+) channel pore model consistent with most of the available biophysical and mutational information suggests that the KcsA structural framework may be conserved in voltage-gated channels.     

Hydrophobic interface between two regulators of K+ conductance domains critical for calcium-dependent activation of large conductance Ca2+-activated K+ channels

It has been suggested that the large conductance Ca(2)+-activated K(+) channel contains one or more domains known as regulators of K(+) conductance (RCK) in its cytosolic C terminus. Here, we show that the second RCK domain (RCK2) is functionally important and that it forms a heterodimer with RCK1 via a hydrophobic interface. Mutant channels lacking RCK2 are nonfunctional despite their tetramerization and surface expression. The hydrophobic residues that are expected to form an interface between RCK1 and RCK2, based on the crystal structure of the bacterial MthK channel, are well conserved, and the interactions of these residues were confirmed by mutant cycle analysis. The hydrophobic interaction appears to be critical for the Ca(2+)-dependent gating of the large conductance Ca(2+)-activated K(+) channel.     

Functional Characterization of the Cardiac Ryanodine Receptor Pore-Forming Region

Ryanodine receptors are homotetrameric intracellular calcium release channels. The efficiency of these channels is underpinned by exceptional rates of cation translocation through the open channel and this is achieved at the expense of the high degree of selectivity characteristic of many other types of channel. Crystallization of prokaryotic potassium channels has provided insights into the structures and mechanisms responsible for ion selection and movement in these channels, however no equivalent structural detail is currently available for ryanodine receptors. Nevertheless both molecular modeling and cryo-electron microscopy have identified the probable pore-forming region (PFR) of the ryanodine receptor (RyR) and suggest that this region contains structural elements equivalent to those of the PFRs of potassium-selective channels. The aim of the current study was to establish if the isolated putative cardiac RyR (RyR2) PFR could form a functional ion channel. We have expressed and purified the RyR2 PFR and shown that function is retained following reconstitution into planar phospholipid bilayers. Our data provide the first direct experimental evidence to support the proposal that the conduction pathway of RyR2 is formed by structural elements equivalent to those of the potassium channel PFR.     

Structural understanding of the transmembrane domains of inositol triphosphate receptors and ryanodine receptors towards calcium channeling.

Inositol 1,4,5-triphosphate receptors (Insp(3)Rs) and ryanodine receptors (ryRs) act as cationic channels transporting calcium ions from the endoplasmic reticulum to cytosol by forming tetramers and are proteins localized to the endoplasmic reticulum (ER). Despite the absence of classical calcium-binding motifs, calcium channeling occurs at the transmembrane domain. We have investigated putative calcium binding motifs in these sequences. Prediction methods indicate the presence of six transmembrane helices in the C-terminal domain, one of the three domains conserved between Insp(3)R and ryR receptors. The recently identified crystal structure of the K(+) channel, which also forms tetramers, revealed that two transmembrane helices, an additional pore helix and a selectivity filter are responsible for selective K(+) ion channeling. The last three TM helices of Insp(3)R and ryR are particularly well conserved and we found analogous pore helix and selectivity filter motif in these sequences. We obtained a three-dimensional structural model for the transmembrane tetramer by extrapolating the distant structural similarity to the K(+) channels.     

Ionic blockade of the rat connexin40 gap junction channel by large tetraalkylammonium ions.

The rat connexin40 gap junction channel is permeable to monovalent cations including tetramethylammonium and tetraethylammonium ions. Larger tetraalkyammonium (TAA(+)) ions beginning with tetrabutylammonium (TBA(+)) reduced KCl junctional currents disproportionately. Ionic blockade by tetrapentylammonium (TPeA(+)) and tetrahexylammonium (THxA(+)) ions were concentration- and voltage-dependent and occurred only when TAA(+) ions were on the same side as net K(+) efflux across the junction, indicative of block of the ionic permeation pathway. The voltage-dependent dissociation constants (K(m)(V(j))) were lower for THxA(+) than TPeA(+), consistent with steric effects within the pore. The K(m)-V(j) relationships for TPeA(+) and THxA(+) were fit with different reaction rate models for a symmetrical (homotypic) connexin gap junction channel and were described by either a one- or two-site model that assumed each ion traversed the entire V(j) field. Bilateral addition of TPeA(+) ions confirmed a common site of interaction within the pore that possessed identical K(m)(V(j)) values for cis-trans concentrations of TPeA(+) ions as indicated by the modeled I-V relations and rapid channel block that precluded unitary current measurements. The TAA(+) block of K(+) currents and bilateral TPeA(+) interactions did not alter V(j)-gating of Cx40 gap junctions. N-octyl-tributylammonium and -triethylammonium also blocked rCx40 channels with higher affinity and faster kinetics than TBA(+) or TPeA(+), indicative of a hydrophobic site within the pore near the site of block.     

Molecular docking of the scorpion toxin Tc1 to the structural model of the voltage-gated potassium channel Kv1.1 from human Homo sapiens

In this study, structural model of the pore loop region of the voltage-gated potassium channel Kv1.1 from human Homo sapiens was constructed based on the crystallographic structure of KcsA by structural homology. The pore loop region of Kv1.1 exhibits similar folds as that of KcsA. The structural feature of the selectivity filter of Kv1.1 is nearly identical to that of KcsA, whereas most of the structural variations occur in the turret as well as in the inner and outer helices. Molecular docking experiments of the scorpion toxin Tc1 from Tityus cambridgei to the outer vestibule of KcsA as well as Kv1.1 were subsequently performed with various initial Tc1 orientations. Tc1 was found to form the most stable complexes with these two K+ channels when the side chain of Lys14 occupies the pore of the selectivity filter through electrostatic interaction. Tc1 binds preferentially towards Kv1.1 than KcsA due to stronger hydrophobic and electrostatic interactions formed between the toxin and the selectivity filter and outer vestibule of Kv1.1. Furthermore, surface complementarity of the outer vestibules of the channels to the Tc1 spatial conformations also plays an important role in stabilizing both the Tc1/KcsA and Tc1/Kv1.1 complexes.     

Light at the end of the Ca(2+)-release channel tunnel: structures and mechanisms involved in ion translocation in ryanodine receptor channels

RyR and InsP3R are Ca(2+)-release channels. When induced to open by the appropriate stimulus, these channels allow Ca2+ to leave intracellular storage organelles at an astonishing rate. Investigations of the ion-handling properties of isolated RyR channels have demonstrated that, at least in comparison to voltage-gated channels of surface membranes, these channels display limited powers of discrimination between physiologically relevant cations and this relative lack of selectivity is likely to contribute to the ability of Ca(2+)-release channels to maintain high rates of cation translocation without compromising function. A range of ion-handling properties in RyR are consistent with the proposal that this channel functions as a single-ion channel and theoretical considerations indicate that the high rates of ion translocation monitored for RyR would require the pore of such a structure to be short and possess a large capture radius. Measurements of the dimensions of regions of RyR involved in ion conduction and discrimination indicate that this is likely to be the case. In each monomer of RyR/InsP3R, residues making up the last two trans-membrane spanning domains and a luminal loop linking these two helices contribute to the formation of the channel pore. The luminal loops of both RyR and InsP3R contain amino acid sequences similar to those known to form the selectivity filter of K+ channels. In addition the luminal loops of both Ca(2+)-release channels contain sequences that are likely to form helices that may be analogous to the pore helix visualised in KcsA. The correlation in structural elements of the luminal loops of RyR/InsP3R and KcsA has prompted us to speculate on the tertiary arrangement for this region of the Ca(2+)-release channels using the established structure of KcsA as a framework.     

Polar/Ionizable Residues in Transmembrane Segments: Effects on Helix-Helix Packing

The vast majority of membrane proteins are anchored to biological membranes through hydrophobic α-helices. Sequence analysis of high-resolution membrane protein structures show that ionizable amino acid residues are present in transmembrane (TM) helices, often with a functional and/or structural role. Here, using as scaffold the hydrophobic TM domain of the model membrane protein glycophorin A (GpA), we address the consequences of replacing specific residues by ionizable amino acids on TM helix insertion and packing, both in detergent micelles and in biological membranes. Our findings demonstrate that ionizable residues are stably inserted in hydrophobic environments, and tolerated in the dimerization process when oriented toward the lipid face, emphasizing the complexity of protein-lipid interactions in biological membranes.     

Local hydrophobicity stabilizes secondary structures in proteins

The probability of occurrence of helix and β-sheet residues in 47 globular proteins was determined as a function of local hydrophobicity, which was defined by the sum of the Nozaki-Tanford transfer free energies at two nearest-neighbors on both sides of the amino acid sequence. In general, hydrophilic amino acids favor neither helix nor β-sheet formations when neighbor residues are also hydrophilic but favor helix formation at higher local hydrophobicity. On the other hand, some hydrophobic amino acids such as Met, Leu, and Ile favor helix formation when neighbor residues are hydrophilic. None of the hydrophobic amino acids favor β-sheet formation with hydrophilic neighbors, but most of them strongly favor β-sheet formation at high local hydrophobicity. When the average of 20 amino acids is taken, both helix and β-sheet residue probabilities are higher at higher local hydrophobicity, although the increase is steeper for β-sheets. Therefore, β-sheet formation is more influenced by local hydrophobicity than helix formation. Generally, helices are nearer the surface and tend to have hydrophilic and hydrophobic faces at opposite sides. The tendency of alternating regions of hydrophilic and hydrophobic residues in a helical sequence was revealed by calculating the correlation of the Nozaki-Tanford values. Such amphipathic helices may be important in protein–protein and protein–lipid interactions and in forming hydrophilic channels in the membrane. The choice of 30 nonhomologous proteins as the data set did not alter the above results.