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.

Potassium channel structure: domain by domain.

Since the determination of the structure of a bacterial potassium channel, the ion channel community has managed to gain momentum in the quest for a complete picture. The information is coming at a steady flow, on a domain by domain basis. Recent discoveries are starting to reveal clues to the complex manner in which potassium channels show enormous diversity of function and also to their methods of regulation. Currently, the structures of four domains are known, with the most recent addition being the Kvbeta structure. As efforts continue in the study of the transmembrane domains, especially the voltage-sensing apparatus, there has been a new realization with respect to the identification and role of the cytoplasmic domains in protein-protein interactions in particular. An additional discovery, considerably aided by recent genomic analysis, is that potassium channels comprising subunits with two pore regions and four transmembrane helices combined in a dimeric fashion are abundant and are probable targets for local anesthetics.

Simulation studies of the interaction of antimicrobial peptides and lipid bilayers.

Experimental studies of a number of antimicrobial peptides are sufficiently detailed to allow computer simulations to make a significant contribution to understanding their mechanisms of action at an atomic level. In this review we focus on simulation studies of alamethicin, melittin, dermaseptin and related antimicrobial, membrane-active peptides. All of these peptides form amphipathic alpha-helices. Simulations allow us to explore the interactions of such peptides with lipid bilayers, and to understand the effects of such interactions on the conformational dynamics of the peptides. Mean field methods employ an empirical energy function, such as a simple hydrophobicity potential, to provide an approximation to the membrane. Mean field approaches allow us to predict the optimal orientation of a peptide helix relative to a bilayer. Molecular dynamics simulations that include an atomistic model of the bilayer and surrounding solvent provide a more detailed insight into peptide-bilayer interactions. In the case of alamethicin, all-atom simulations have allowed us to explore several steps along the route from binding to the membrane surface to formation of transbilayer ion channels. For those antimicrobial peptides such as dermaseptin which prefer to remain at the surface of a bilayer, molecular dynamics simulations allow us to explore the favourable interactions between the peptide helix sidechains and the phospholipid headgroups.

Alamethicin channels - modelling via restrained molecular dynamics simulations.

Alamethicin channels have been modelled as approximately parallel bundles of transbilayer helices containing between N = 4 and 8 helices per bundle. Initial models were generated by in vacuo restrained molecular dynamics (MD) simulations, and were refined by 60 ps MD simulations with water molecules present within and at the mouths of the central pore. The helix bundles were stabilized by networks of H-bonds between intra-pore water molecules and Gln-7 side-chains. Channel conductances were predicted on the basis of pore radius profiles, and suggested that the N = 4 bundle formed an occluded pore, whereas pores with N > or = 5 helices per bundle were open. Continuum electrostatics calculations suggested that the N = 6 pore is cation-selective, whereas pores with N > or = 7 helices per bundle were predicted to be somewhat less ion-selective.

Potassium and sodium ions in a potassium channel studied by molecular dynamics simulations.

We have performed simulations of both a single potassium ion and a single sodium ion within the pore of the bacterial potassium channel KcsA. For both ions there is a dehydration energy barrier at the cytoplasmic mouth suggesting that the crystal structure is a closed conformation of the channel. There is a potential energy barrier for a sodium ion in the selectivity filter that is not seen for potassium. Radial distribution functions for both ions with the carbonyl oxygens of the selectivity filter indicate that sodium may interact more tightly with the filter than does potassium. This suggests that the key to the ion selectivity of KcsA is the greater dehydration energy of Na(+) ions, and helps to explain the block of KcsA by internal Na(+) ions.

Structure and orientation of the mammalian antibacterial peptide cecropin P1 within phospholipid membranes.

Cecropins are positively charged antibacterial peptides that act by permeating the membrane of susceptible bacteria. To gain insight into the mechanism of membrane permeation, the secondary structure and the orientation within phospholipid membranes of the mammalian cecropin P1 (CecP) was studied using attenuated total reflectance Fourier-transform infrared (ATR-FTIR) spectroscopy and molecular dynamics simulations. The shape and frequency of the amide I and II absorption peaks of CecP within acidic PE/PG multibilayers (phosphatidylethanolamine/phosphatidylglycerol) in a 7:3 (w/w) ratio (a phospholipid composition similar to that of many bacterial membranes), indicated that the peptide is predominantly alpha-helical. Polarized ATR-FTIR spectroscopy was used to determine the orientation of the peptide relative to the bilayer normal of phospholipid multibilayers. The ATR dichroic ratio of the amide I band of CecP peptide reconstituted into oriented PE/PG phospholipid membranes indicated that the peptide is preferentially oriented nearly parallel to the surface of the lipid membranes. A similar secondary structure and orientation were found when zwitterionic phosphatidylcholine phospholipids were used. The incorporation of CecP did not significantly change the order parameters of the acyl chains of the multibilayer, further suggesting that CecP does not penetrate the hydrocarbon core of the membranes. Molecular dynamics simulations were used to gain insight into possible effects of transmembrane potential on the orientation of CecP relative to the membrane. The simulations appear to confirm that CecP adopts an orientation parallel to the membrane surface and does not insert into the bilayer in response to a cis positive transmembrane voltage difference. Taken together, the results further support a "carpet-like" mechanism, rather than the formation of transmembrane pores, as the mode of action of CecP. According to this model, formation of a layer of peptide monomers on the membrane surface destablizes the phospholipid packing of the membrane leading to its eventual disintegration.

Interactions of alpha-helices with lipid bilayers: a review of simulation studies.

Membrane proteins, of which the majority seem to contain one or more alpha-helix, constitute approx. 30% of most genomes. A complete understanding of the nature of helix/bilayer interactions is necessary for an understanding of the structural principles underlying membrane proteins. This review describes computer simulation studies of helix/bilayer interactions. Key experimental studies of the interactions of alpha-helices and lipid bilayers are briefly reviewed. Surface associated helices are found in some membrane-bound enzymes (e.g. prostaglandin synthase), and as stages in the mechanisms of antimicrobial peptides and of pore-forming bacterial toxins. Transmembrane alpha-helices are found in most integral membrane proteins, and also in channels formed by amphipathic peptides or by bacterial toxins. Mean field simulations, in which the lipid bilayer is approximated as a hydrophobic continuum, have been used in studies of membrane-active peptides (e.g. alamethicin, melittin, magainin and dermaseptin) and of simple membrane proteins (e.g. phage Pf1 coat protein). All atom molecular dynamics simulations of fully solvated bilayers with transmembrane helices have been applied to: the constituent helices of bacteriorhodopsin; peptide-16 (a simple model TM helix); and a number of pore-lining helices from ion channels. Surface associated helices (e.g. melittin and dermaseptin) have been simulated, as have alpha-helical bundles such as bacteriorhodopsin and alamethicin. From comparison of the results from the two classes of simulation, it emerges that a major theoretical challenge is to exploit the results of all atom simulations in order to improve the mean field approach.

Simulation of voltage-dependent interactions of alpha-helical peptides with lipid bilayers.

Pore formation in lipid bilayers by channel-forming peptides and toxins is thought to follow voltage-dependent insertion of amphipathic alpha-helices into lipid bilayers. We have developed an approximate potential for use within the CHARMm molecular mechanics program which enables one to simulate voltage-dependent interaction of such helices with a lipid bilayer. Two classes of helical peptides which interact with lipid bilayers have been studied: (a) delta-toxin, a 26 residue channel-forming peptide from Staphylococcus aureus; and (b) synthetic peptides corresponding to the alpha 5 and alpha 7 helices of the pore-forming domain of Bacillus thuringiensis CryIIIA delta-endotoxin. Analysis of delta-toxin molecular dynamics (MD) simulations suggested that the presence of a transbilayer voltage stabilized the inserted location of delta-toxin helices, but did not cause insertion per se. A series of simulations for the alpha 5 and alpha 7 peptides revealed dynamic switching of the alpha 5 helix between a membrane-associated and a membrane-inserted state in response to a transbilayer voltage. In contrast the alpha 7 helix did not exhibit such switching but instead retained a membrane associated state. These results are in agreement with recent experimental studies of the interactions of synthetic alpha 5 and alpha 7 peptides with lipid bilayers.

Edit, cut and paste in the nicotinic acetylcholine receptor gene family of Drosophila melanogaster.

Nicotinic acetylcholine receptors (nAChRs) are important for fast synaptic cholinergic transmission. They are targets of drugs/chemicals for human and animal health as well as for pest control. With the advent of genome sequencing, entire nAChR gene families have now been described for vertebrates and invertebrates. Mostly, these are extensive with a large number of distinct subunits, making possible many nAChR subtypes differing in transmitter affinity, channel conductance, ion selectivity, desensitization, modulation and pharmacology. The smallest nAChR gene family to date is that of the fruit fly, Drosophila melanogaster, with only 10 members. This apparently compact family belies its true diversity as 4 of the 10 subunits show alternative splicing. Also, using Drosophila, A-to-I pre-mRNA editing has been demonstrated for the first time in nAChRs. Such is the extent of this variation, that one subunit alone (Dalpha6) can potentially generate far more isoforms than seen in entire gene families from other species. We present here three-dimensional models constructed for insect nAChRs, which show that many variations introduced by alternative splicing and RNA editing may influence receptor function.

Structure-function relationships in helix-bundle channels probed via total chemical synthesis of alamethicin dimers: effects of a Gln7 to Asn7 mutation.

Alamethicin channels are prototypical helix bundles that may serve as tractable models for more complex protein ion channels. Solid-phase peptide synthesis of alamethicin analogues using FMOC-amino acid fluorides followed by chemical dimerization of these peptides facilitates structure-function studies of particular channel states in bilayer membranes. State 3 in particular, tentatively assigned to a hexameric helix bundle, is sufficiently long-lived that current-voltage measurements can be made during the lifetime of an individual channel opening. Molecular models of hexameric helix bundles, generated using restrained molecular dynamics with simulated annealing, indicate that a Gln7-->Asn7 (Q7-->N7) mutation will increase channel diameter locally. Experimentally, the conductance of state 3 of the N7-alm channel is found to be larger than that of the Q7-alm channel when ion flow is in the usual direction (cations entering the C-terminal end of the channel). When ion flow is in the opposite direction, no difference in the conductances of state 3 of Q7 and state 3 of N7 channels is observed. These results indicate that the effect of a change in pore diameter at position 7 is dependent on the magnitude of other barriers to permeation and that these barriers are voltage-dependent.

The dielectric properties of water within model transbilayer pores.

Ion channels contain extended columns of water molecules within their transbilayer pores. The dynamic properties of such intrapore water have been shown to differ from those of water in its bulk state. In previous molecular dynamics simulations of two classes of model pore (parallel bundles of Ala20 alpha-helices and antiparallel barrels of Ala10 beta-strands), a substantially reduced translational and rotational mobility of waters was observed within the pore relative to bulk water. Molecular dynamics simulations in the presence of a transpore electrostatic field (i.e., a voltage drop along the pore axis) have been used to estimate the resultant polarization (due to reorientation) of the intrapore water, and hence to determine the local dielectric behavior within the pore. It is shown that the local dielectric constant of water within a pore is reduced for models formed by parallel alpha-helix bundles, but not by those formed by beta-barrels. This result is discussed in the context of electrostatics calculations of ion permeation through channels, and the effect of the local dielectric of water within a helix bundle pore is illustrated with a simple Poisson-Boltzmann calculation.

Simulation studies of alamethicin-bilayer interactions.

Alamethicin is an alpha-helical peptide that forms voltage-activated ion channels. Experimental data suggest that channel formation occurs via voltage-dependent insertion of alamethicin helices into lipid bilayers, followed by self-assembly of inserted helices to form a parallel helix bundle. Changes in the kink angle of the alamethicin helix about its central proline residue have also been suggested to play a role in channel gating. Alamethicin helices generated by simulated annealing and restrained molecular dynamics adopt a kink angle similar to that in the x-ray crystal structure, even if such simulations start with an idealized unkinked helix. This suggests that the kinked helix represents a stable conformation of the molecule. Molecular dynamics simulations in the presence of a simple bilayer model and a transbilayer voltage difference are used to explore possible mechanisms of helix insertion. The bilayer is represented by a hydrophobicity potential. An alamethicin helix inserts spontaneously in the absence of a transbilayer voltage. Application of a cis positive voltage decreases the time to insertion. The helix kink angle fluctuates during the simulations. Insertion of the helix is associated with a decrease in the mean kink angle, thus helping the alamethicin molecule to span the bilayer. The simulation results are discussed in terms of models of alamethicin channel gating.

Intrinsic rectification of ion flux in alamethicin channels: studies with an alamethicin dimer.

Covalent dimers of alamethicin form conducting structures with gating properties that permit measurement of current-voltage (I-V) relationships during the lifetime of a single channel. These I-V curves demonstrate that the alamethicin channel is a rectifier that passes current preferentially, with voltages of the same sign as that of the voltage that induced opening of the channel. The degree of rectification depends on the salt concentration; single-channel I-V relationships become almost linear in 3 M potassium chloride. These properties may be qualitatively understood by using Poisson-Nernst-Planck theory and a modeled structure of the alamethicin pore.