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.

Ion channels: a first view of K+ channels in atomic glory

Crystal structures have been solved for the transbilayer pore domain of a bacterial K+ channel and the tetramerisation domain of voltage-gated K+ channel. These provide our first real structural insights into possible mechanisms of ion selectivity and permeation for K+ channels.

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

The ion channel of the nicotinic acetylcholine receptor (nAChR) is believed to be lined by transmembrane M2 helices. A "4-8-12" sequence motif, comprising serine (S) or threonine (T) residues at positions 4, 8 and 12 of M2, is conserved between different members, anion and cation selective, of the nAChR superfamily. Parallel bundles of 4-8-12 motif-containing helices are considered as simplified models of ion channels. The relationship between S and T sidechain conformations and channel-ion interactions is explored via evaluation of interaction energies of K+ and of Cl- ions with channel models. Energy calculations are used to determine optimal chi 2 (C alpha-C beta-O gamma-H gamma) values in the presence of K+ or Cl- ions. 4-8-12 motif-containing bundles may form favourable interactions with either cations or anions, dependent upon the chi 2 values adopted. Parallel-helix and tilted-helix bundles are considered, as are heteromeric models designed to mimic the Torpedo nAChR. The main conclusion of the study is that conformational flexibility at chi 2 enables both S and T residues to form favourable interactions with anions or cations. Consequently, there is apparently no difference between S and T residues in their interactions with permeant ions, which suggests that the presence of T vs. S residues within the 4-8-12 motif is not a major mechanism whereby anion/cation selectivity may be generated. The implications of these studies with respect to more elaborate models of nAChR and related receptors are considered.

Channels formed by the transmembrane helix of phospholamban: a simulation study

Phospholamban is a small membrane protein which can form cation selective ion channels in lipid bilayers. Each subunit contains a single, largely hydrophobic transmembrane helix. The helices are thought to assemble as a pentameric and approximately parallel bundle surrounding a central pore. A model of this assembly (PDB code IPSL) has been used as the starting point for molecular dynamics (MD) simulations of a system consisting of the pentameric helix bundle, plus 217 water molecules located within and at either mouth of the pore. Interhelix distance restraints were employed to maintain the integrity of the helix bundle during a 500 ps MD simulation. Water molecules within the pore exhibited reduced diffusional and rotational mobility. Interactions between the alpha-helix dipoles and the water dipoles, the latter aligned anti-parallel to the former, contribute to the stability of the system. Analysis of the potential energy of interaction of a K+ ion as it was moved through the pore suggested that unfavourable interactions of the cation with the aligned helix dipoles at the N-terminal mouth were overcome by favourable ion-water interactions. Comparable analysis for a Cl ion revealed that the ion-(pore + water) interactions were unfavourable along the whole of the pore, increasingly so from the N- to the C-terminal mouth. Overall, the interaction energy profiles were consistent with a pore selective for cations over anions. Pore radius profiles were used to predict a channel conductance of 50 to 70 ps in 0.2 M KCl, which compares well with an experimental value of 100 ps.

Ion channel stability and hydrogen bonding. Molecular modelling of channels formed by synthetic alamethicin analogues

Several analogues of the channel-forming peptaibol alamethicin have been demonstrated to exhibit faster switching between channel substates than does unmodified alamethicin. Molecular modelling studies are used to explore the possible molecular basis of these differences. Models of channels formed by alamethicin analogues were generated by restrained molecular dynamics in vacuo and refined by short molecular dynamics simulations with water molecules within and at either mouth of the channel. A decrease in backbone solvation was found to correlate with a decrease in open channel stability between alamethicin and an analogue in which all alpha-amino-isobutyric acid residues of alamethicin were replaced by leucine. A decrease in the extent of hydrogen-bonding at residue 7 correlates with lower open channel stabilities of analogues in which the glutamine at position 7 was replaced by smaller polar sidechains. These two observations indicate the importance of alamethicin/water H-bonds in stabilizing the open channel.

Single-channel data and missed events: analysis of a two-state Markov model

Patch-clamp recording permits investigation of the gating kinetics of single ion channels. Careful statistical analysis of kinetic data can yield clues as to the molecular events underlying channel gating. However, it is important that such analysis should take full account of the limitations that arise from the finite time resolution of patch-clamp recording techniques. Single-ion-channel data are generally interpreted in terms of Markov process models of channel gating mechanisms. Experimental channel records suffer from time interval omission, i.e. failure to detect brief channel openings and closings. This leads to an identifiability problem when analysing single-channel data, i.e. different gating mechanisms provide equally convincing descriptions of the same experimental data. We consider a two-state Markov model of receptor-channel gating in which the channel opening rate is proportional to the agonist concentration, C in equilibrium with OA. By using computer-simulated data, the approximate likelihood of the data is maximized to yield parameter estimates for the model. At a single agonist concentration there is an identifiability problem in that two pairs of parameter estimates are obtained. The 'true' parameter estimates cannot be distinguished from the 'false' ones. By considering data corresponding to a range of agonist concentrations one may identify the 'true' parameter estimates as those that do not change as the agonist concentration is increased. Alternatively, one may identify the 'true' parameter estimates directly by maximizing a global likelihood, the latter being obtained by simultaneous consideration of data obtained at several different agonist concentrations.(ABSTRACT TRUNCATED AT 250 WORDS)

Channel gating: twist to open

Recent studies of the bacterial mechanosensitive channel MscL have combined a number of different approaches to come up with a model for the channel gating mechanism.

The pore-lining region of shaker voltage-gated potassium channels: comparison of beta-barrel and alpha-helix bundle models

Although there is a large body of site-directed mutagenesis data that identify the pore-lining sequence of the voltage-gated potassium channel, the structure of this region remains unknown. We have interpreted the available biochemical data as a set of topological and orientational restraints and employed these restraints to produce molecular models of the potassium channel pore region, H5. The H5 sequence has been modeled either as a tetramer of membrane-spanning beta-hairpins, thus producing an eight-stranded beta-barrel, or as a tetramer of incompletely membrane-spanning alpha-helical hairpins, thus producing an eight-staved alpha-helix bundle. In total, restraints-directed modeling has produced 40 different configurations of the beta-barrel model, each configuration comprising an ensemble of 20 structures, and 24 different configurations of the alpha-helix bundle model, each comprising an ensemble of 24 structures. Thus, over 1300 model structures for H5 have been generated. Configurations have been ranked on the basis of their predicted pore properties and on the extent of their agreement with the biochemical data. This ranking is employed to identify particular configurations of H5 that may be explored further as models of the pore-lining region of the voltage-gated potassium channel pore.

Modelling packing interactions in parallel helix bundles: pentameric bundles of nicotinic receptor M2 helices

The transbilayer pore of the nicotinic acetylcholine receptor (nAChR) is formed by a pentameric bundle of M2 helices. Models of pentameric bundles of M2 helices have been generated using simulated annealing via restrained molecular dynamics. The influence of: (a) the initial C alpha template; and (b) screening of sidechain electrostatic interactions on the geometry of the resultant M2 helix bundles is explored. Parallel M2 helices, in the absence of sidechain electrostatic interactions, pack in accordance with simple ridges-in-grooves considerations. This results in a helix crossing angle of ca. +12 degrees, corresponding to a left-handed coiled coil structure for the bundle as a whole. Tilting of M2 helices away from the central pore axis at their C-termini and/or inclusion of sidechain electrostatic interactions may perturb such ridges-in-grooves packing. In the most extreme cases right-handed coiled coils are formed. An interplay between inter-helix H-bonding and helix bundle geometry is revealed. The effects of changes in electrostatic screening on the dimensions of the pore mouth are described and the significance of these changes in the context of models for the nAChR pore domain is discussed.

Protein-water-ion interactions in a model of the pore domain of a potassium channel: a simulation study

A model of the selectivity filter of a voltage-gated K+ (Kv) channel formed by an eight-stranded beta-barrel is compared with physiological properties of the channel. Continuum electrostatic calculations suggest that only two of the eight Asp sidechains at the extracellular mouth of the pore will ionise. A ring of four Tyr sidechains forms the narrowest region of the pore. Molecular dynamic simulations of the potential energy of a K+ ion as translated along the model pore indicate that the two ionised Asp sidechains and the hydroxyl groups of the Tyr sidechains stabilise the partially desolvated ion as it passes through the narrowest region.

Peptide-surfactant interactions: A combined spectroscopic and molecular dynamics simulation approach

In the present contribution, we report a combined spectroscopic and computational approach aiming to unravel at atomic resolution the effect of the anionic SDS detergent on the structure of two model peptides, the α-helix TrpCage and the β-stranded TrpZip. A detailed characterization of the specific amino acids involved is performed. Monomeric (single molecules) and micellar SDS species differently interact with the α-helix and β-stranded peptides, emphasizing the different mechanisms occurring below and above the critical aggregation concentration (CAC). Below the CAC, the α-helix peptide is fully unfolded, losing its hydrophobic core and its Asp-Arg salt bridge, while the β-stranded peptide keeps its native structure with its four Trp well oriented. Above the CAC, the SDS micelles have the same effect on both peptides, that is, destabilizing the tertiary structure while keeping their secondary structure. Our studies will be helpful to deepen our understanding of the action of the denaturant SDS on peptides and proteins.

Potassium channels: watching a voltage-sensor tilt and twist

The fourth transmembrane helix (S4) is the primary voltage-sensor of voltage-gated ion channels. Recent studies have used fluorescence resonance energy transfer as a spectroscopic ruler to determine the nature and magnitude of the voltage-induced movement of S4 that leads to channel opening.

Analysis of post-perturbation gating kinetics of single ion channels

Analysis of mean dwell-times as a function of the number of channel openings elapsed since a stepwise perturbation in ion-channel kinetics is shown to provide information concerning the topology of the underlying gating mechanism. The difference between the post-perturbation mean dwell-time and the corresponding equilibrium mean is shown to decay as the sum of Ng-1 geometric terms in k, the number of openings since the perturbation, where Ng is the minimum number of gateway states in the channel gating mechanism. The method is illustrated by consideration of various simple gating schemes. A modification of the method accommodating the presence of channel inactivation or desensitization is described. Application of the method to a delayed-rectifier type K+ channel of NG108-15 cells reveals that Ng greater than or equal to 2, consistent with a branched gating mechanism.

Membrane proteins: A tale of barrels and corks

Crystal structures have been solved for two bacterial outer membrane proteins, FhuA and FepA, which mediate active transport of chelated iron. Analysis of ligand-induced changes in the structure of FhuA has provided our first structural insights into an active transport mechanism for a complex solute.

Simulation of the packing of idealized transmembrane alpha-helix bundles

The aim of this study is to investigate if the packing motifs of native transmembrane helices can be produced by simulations with simple potentials and to develop a method for the rapid generation of initial candidate models for integral membrane proteins composed of bundles of transmembrane helices. Constituent residues are mapped along the helix axis in order to maintain the amino acid sequence-dependent properties of the helix. Helix packing is optimized according to a semi-empirical potential mainly composed of four components: a bilayer potential, a crossing angle potential, a helix dipole potential and a helix-helix distance potential. A Monte Carlo simulated annealing protocol is employed to optimize the helix bundle system. Necessary parameters are derived from theoretical studies and statistical analysis of experimentally determined protein structures. Preliminary testing of the method has been conducted with idealized seven Ala(20) helix bundles. The structures generated show a high degree of compactness. It was observed that both bacteriorhodopsin-like and delta-endotoxin-like structures are generated in seven-helix bundle simulations, within which the composition varies dependent upon the cooling rate. The simulation method has also been employed to explore the packing of N = 4 and N = 12 transmembrane helix bundles. The results suggest that seven and 12 transmembrane helix bundles resembling those observed experimentally (e.g., bacteriorhodopsin, rhodopsin and cytochrome c oxidase subunit I) may be generated by simulations using simple potentials.

Simulation studies on bacteriorhodopsin alpha-helices

Bacteriorhodopsin (BR) is a membrane protein which pumps protons through the plasma membrane. Transmembrane BR helical segments are subjected to simulation studies in order to investigate the effect of bilayer environment in various simulation conditions. A bilayer potential is introduced to the system to mimic the lipid membrane. The structures from the simulations are compared with the experimentally determined structures in terms of geometrical properties. Electrostatic contribution to the helix packing is also investigated. The simulation results show that the packing geometry of the transmembrane helices is highly affected by the bilayer potential. The results obtained from the simulations may be used for further simulation studies and analysis in investigating transmembrane helix packing.

Molecular dynamics simulations of isolated transmembrane helices of potassium channels

In the middle of the S6 helix in voltage-gated potassium channels there is a highly conserved Pro-Val-Pro motif, while the equivalent M2 helix of inward rectifier potassium channels contains a conserved glycine residue in a comparable position. The structural implications of these conserved motifs are of interest given the evidence that S6 and M2 are components of the lining of their respective pores. Multiple sequence alignment and TM helix prediction methods were used to define consensus regions for S6 and M2. Ensembles of 50 structures for each helix were generated by simulated annealing and restrained molecular dynamics. Time-dependent fluctuations of S6 and M2 were investigated by long time scale molecular dynamics simulations on representative members of each ensemble carried out in vacuo in the presence and absence of a hydrophobic potential that mimics a lipid bilayer. The results are discussed in terms of the structural basis of the kink in S6 and M2 and of a putative functional role for flexible helices as "molecular swivels."