The roles of serine and threonine sidechains in ion channels: a modelling study

Interaction of insecticides with heteromeric GABA-gated chloride channels from zebrafish Danio rerio (Hamilton)

The ionotropic GABA_A receptor (GABA_AR) is the main fast inhibitory post-synaptic receptor and is also an important insecticidal target. Effect of insecticides on fish has attracted intensive attention. However, no systematic study on heteromeric zebrafish GABA_AR expressed in oocytes has been reported to date. In this study, the α1 subunit, the β2S subunit lacking 47 amino acid residues compared with the β2L subunit, and the γ2 subunit having five transmembrane domains were isolated from zebrafish Danio rerio. The responses of homomeric and heteromeric (α1, β2S and γ2) channels to agonists and GABA_AR-targeted compounds were detected with two-electrode voltage clamp. Dose-dependent responses were observed in heteromeric α1β2S, β2Sγ2, and α1β2Sγ2 GABAR channels with EC₅₀ values at 21.75, 6291, and 33.69 μM for GABA-induced current and 3.28, 155.5, and 3.79 mM for β-alanine-induced current, respectively. However, no response was induced by benzamidine in all GABAR channels. Abamectin, dieldrin, fluralaner and fipronil could strongly inhibited GABA-induced inward current ≥50% at 10^-6 M, while α-endosulfan, flufiprole and ethiprole only inhibited GABA-induced current <50%. This study has clarified the interaction of insecticides with the heteromeric GABA_AR channel, which could help us further explore the potential function and toxicological importance of GABA_ARs from D. rerio.

An amino-acid substitution in the influenza-B NB protein affects ion-channel gating

The effects of site-directed mutations in NB, a protein encoded by the influenza B virus that has been shown to form cation-selective ion channels at pH 6.0, were studied on ion channel characteristics in artificial lipid bilayers. It was thought that the residues in the hydrophobic region of NB we selected for mutation might be involved in the transport of cations across the channel and that changes in these residues might affect channel properties such as gating and ion-selectivity. Serine residues at positions 20 and 28, threonine at position 24 and cysteine at position 26 were replaced by alanine. We found that the mutation S20A gave channels that did not gate and that remained open most of the time. Proton permeability of NB channels, as detected by fluorescence quenching, was also altered by the mutation S20A: channels were no longer proton-permeable. The other mutations, S28A, T24A and C26A, did not have any detectable effect on the activity or proton permeability of channels formed by NB. The results indicate that serine 20 may have an important role in normal function of NB channels.

Simulations of ion permeation through a potassium channel: molecular dynamics of KcsA in a phospholipid bilayer

Potassium channels enable K(+) ions to move passively across biological membranes. Multiple nanosecond-duration molecular dynamics simulations (total simulation time 5 ns) of a bacterial potassium channel (KcsA) embedded in a phospholipid bilayer reveal motions of ions, water, and protein. Comparison of simulations with and without K(+) ions indicate that the absence of ions destabilizes the structure of the selectivity filter. Within the selectivity filter, K(+) ions interact with the backbone (carbonyl) oxygens, and with the side-chain oxygen of T75. Concerted single-file motions of water molecules and K(+) ions within the selectivity filter of the channel occur on a 100-ps time scale. In a simulation with three K(+) ions (initially two in the filter and one in the cavity), the ion within the central cavity leaves the channel via its intracellular mouth after approximately 900 ps; within the cavity this ion interacts with the Ogamma atoms of two T107 side chains, revealing a favorable site within the otherwise hydrophobically lined cavity. Exit of this ion from the channel is enabled by a transient increase in the diameter of the intracellular mouth. Such "breathing" motions may form the molecular basis of channel gating.

Permeability of wild-type and mutant cystic fibrosis transmembrane conductance regulator chloride channels to polyatomic anions

Permeability of the cystic fibrosis transmembrane conductance regulator (CFTR) chloride channel to polyatomic anions of known dimensions was studied in stably transfected Chinese hamster ovary cells by using the patch clamp technique. Biionic reversal potentials measured with external polyatomic anions gave the permeability ratio (P/P) sequence NO > Cl > HCO > formate > acetate. The same selectivity sequence but somewhat higher permeability ratios were obtained when anions were tested from the cytoplasmic side. Pyruvate, propanoate, methane sulfonate, ethane sulfonate, and gluconate were not measurably permeant (P/P < 0.06) from either side of the membrane. The relationship between permeability ratios from the outside and ionic diameters suggests a minimum functional pore diameter of approximately 5.3 A. Permeability ratios also followed a lyotropic sequence, suggesting that permeability is dependent on ionic hydration energies. Site-directed mutagenesis of two adjacent threonines in TM6 to smaller, less polar alanines led to a significant (24%) increase in single channel conductance and elevated permeability to several large anions, suggesting that these residues do not strongly bind permeating anions, but may contribute to the narrowest part of the pore.

Structural features of isolated M2 helices of nicotinic receptors. Simulated annealing via molecular dynamics studies

The nicotinic acetylcholine receptor is an integral membrane protein and a ligand-gated cation channel. It has stoichiometry alpha 2 beta gamma delta, the subunits arranged symmetrically around an approximate five-fold axis. Five M2 helices, one from each subunit, form a parallel helix bundle surrounding a central pore. Simulated annealing via restrained molecular dynamics (SA/MD) has been employed to generate ensembles of isolated M2 transmembrane helices. Four ensembles of two different M2 helix sequences, M2 delta and M2 gamma, have been generated by SA/MD. The ensembles differed in their treatment of electrostatic interactions. Analysis of the simulated structures showed that intra-helical H-bonds were more strongly conserved in the C-terminal (and more hydrophobic) segment of M2 helices. Conformations of polar sidechains have been analyzed, placing particular emphasis on EK (and QK) pairs at the N-termini of M2 delta (and M2 gamma) helices. Conformations of EK sidechain pairs were obtained for the high resolution structures in the protein database in order to guide our analysis of simulated structures. Serine and threonine sidechain conformations in the M2 models also have been determined. Implications of studies of isolated M2 helices for models of the intact pore region of the nicotinic receptor are discussed.

Kinked structures of isolated nicotinic receptor M2 helices: a molecular dynamics study

The pore-lining M2 helix of the nicotinic acetylcholine receptor exhibits a pronounced kink when the corresponding ion channel is in a closed conformation [N. Unwin (1993) Journal of Molecular Biology, Vol. 229, pp. 1101-1124]. We have performed molecular dynamics simulations of isolated 22-residue M2 helices in order to identify a possible molecular origin of this kink. In order to sample a wide range of conformational space, a simulated annealing protocol was used to generate five initial M2 helix structures, each of which was subsequently used as the basis of 300 ps MD simulations. Two helix sequences (M2 alpha and M2 delta) were studied in this manner, resulting in a total of ten 300 ps trajectories. Kinked helices present in the trajectories were identified and energy minimized to yield a total of five different stable kinked structures. For comparison, a similar molecular dynamics simulation of a Leu23 helix yielded no stable kinked structures. In four of the five kinked helices, the kink was stabilized by H bonds between the helix backbone and polar side-chain atoms. Comparison with data from the literature on site-directed mutagenesis of M2 residues suggests that such polar side-chain to main-chain H bonds may also contribute to kinking of M2 helices in the intact channel protein.

Parallel helix bundles and ion channels: molecular modeling via simulated annealing and restrained molecular dynamics

A parallel bundle of transmembrane (TM) alpha-helices surrounding a central pore is present in several classes of ion channel, including the nicotinic acetylcholine receptor (nAChR). We have modeled bundles of hydrophobic and of amphipathic helices using simulated annealing via restrained molecular dynamics. Bundles of Ala20 helices, with N = 4, 5, or 6 helices/bundle were generated. For all three N values the helices formed left-handed coiled coils, with pitches ranging from 160 A (N = 4) to 240 A (N = 6). Pore radius profiles revealed constrictions at residues 3, 6, 10, 13, and 17. A left-handed coiled coil and a similar pattern of pore constrictions were observed for N = 5 bundles of Leu20. In contrast, N = 5 bundles of Ile20 formed right-handed coiled coils, reflecting loosened packing of helices containing beta-branched side chains. Bundles formed by each of two classes of amphipathic helices were examined: (a) M2a, M2b, and M2c derived from sequences of M2 helices of nAChR; and (b) (LSSLLSL)3, a synthetic channel-forming peptide. Both classes of amphipathic helix formed left-handed coiled coils. For (LSSLLSL)3 the pitch of the coil increased as N increased from 4 to 6. The M2c N = 5 helix bundle is discussed in the context of possible models of the pore domain of nAChR.

Structure and function of channel-forming peptaibols

Transport of ions through channels is fundamental to a number of physiological processes, especially the electrical properties of excitable cells (Hille, 1992). To understand this process at a molecular level requires atomic resolution structures of channel proteins.

Hydrophilic surface maps of channel-forming peptides: analysis of amphipathic helices

Ion channels may be formed by bundles of amphipathic alpha-helices aligned parallel to one another and spanning a lipid bilayer membrane, with the hydrophilic faces of the helices lining a central pore. In order to provide insight into the packing of such helices in bundles, a method has been developed to evaluate hydrophilic surface maps of amphipathic alpha-helices and to display these surfaces in a readily interpretable form. The procedure is based upon empirical energy calculations of interactions of a water molecule with an amphipathic alpha-helix. The method has been applied to three channel-forming peptides: Staphylococcal delta-toxin; alamethicin; and a synthetic leucine- and serine-containing peptide. Particular emphasis is placed upon the effects of sidechain conformational flexibility on hydrophilic surface maps. A family of models of the delta-toxin helix is generated by a simulated annealing procedure. The results of hydrophilic surface map analyses provide more exact definition of the centre of the hydrophilic face of amphipathic helices, and of the variation of the position of the centre in response to changes in sidechain conformation. This information is used to define families of preliminary models for a given ion channel, as is illustrated for delta-toxin.

Alamethicin and related peptaibols--model ion channels

Peptaibols are considered as models of those ion channels which consist of a bundle of transbilayer helices surrounding a central pore. X-Ray diffraction and NMR studies have yielded high resolution structures for several peptaibols. In conjunction with other spectroscopic investigations and molecular dynamics simulations, these studies suggest that peptaibols form proline-kinked alpha-helices, and that there may be "hinge-bending" movement of the helix in the region of the central proline residue. The amphipathicity of peptaibol helices is analyzed in relation to their channel-forming properties. Studies of the interactions of peptaibols with lipid bilayers suggest that they are helical when in a membrane-like environment, and that the helix orientation relative to the bilayer is sensitive to the peptaibol:lipid ratio, and to the degree of hydration of the bilayer. Electrical studies reveal that many peptaibols form multiple-conductance level channels in a voltage-dependent fashion. Analysis of conductance levels provides support for the "barrel stave" model of channel formation, whereby different conductance levels correspond to different numbers of monomers in a helix bundle. Alternative models for voltage-activation are discussed, and the roles of molecular dipoles and of hinge-bending in this process are considered. Two molecular models for an N = 6 bundle of alamethicin helices are presented and their electrostatic properties analyzed. The relevance of studies of peptaibols to channel and transport proteins in general is considered.

Ion channel formation by zervamicin-IIB. A molecular modelling study

Zervamicin-IIB (Zrv-IIB) is a 16 residue peptaibol which forms voltage-activated, multiple conductance level channels in planar lipid bilayers. A molecular model of Zrv-IIB channels is presented. The structure of monomeric Zrv-IIB is based upon the crystal structure of Zervamicin-Leu. The helical backbone is kinked by a hydroxyproline residue at position 10. Zrv-IIB channels are modelled as helix bundles of from 4 to 8 parallel helices surrounding a central pore. The monomers are packed with their C-terminal helical segments in close contact, and the bundles are stabilized by hydrogen bonds between glutamine 11 and hydroxyproline 10 of adjacent helices. Interaction energy profiles for movement of three different probes species (K+, Cl- and water) through the central pore are analyzed. The conformations of: (a) the sidechain of glutamine 3; (b) the hydroxyl group of hydroxyproline 10; and (c) the C-terminal hydroxyl group are "optimized" in order to maximize favourable interactions between the channel and the probes, resulting in favourable interaction energy profiles for all three. This suggests that conformational flexibility of polar sidechains enables the channel lining to mimic an aqueous environment.