Single channel kinetics of a glutamate receptor

The glutamate receptor-channel of locust muscle membrane was studied using the patch-clamp technique. Muscles were pretreated with concanavalin A to block receptor-channel desensitization, thus facilitating analysis of receptor-channel gating kinetics. Single channel kinetics were analyzed to aid in identification of the molecular basis of channel gating. Channel dwell-time distributions and dwell-time autocorrelation functions were calculated from single channel data recorded in the presence of 10(-4) M glutamate. Analysis of the dwell time distributions in terms of mixtures of exponential functions revealed there to be at least three open states of the receptor-channel and at least four closed states. Autocorrelation function analysis showed there to be at least three pathways linking the open states with the closed. This results in a minimal scheme for gating of the glutamate receptor-channel, which is suggestive of allosteric models of receptor-channel gating.

Glutamate receptor channel kinetics: the effect of glutamate concentration

Single channel recordings from the locust muscle D-glutamate receptor channel were obtained using glutamate concentrations ranging from 10(-6) to 10(-2) M. Channel kinetics were analyzed to aid in the development of a model for the gating mechanism. Analysis of channel dwell time histograms demonstrated that the channel possessed multiple open and closed states at concentrations of glutamate between 10(-5) and 10(-2) M. Correlations between successive dwell times showed that the gating mechanism was nonlinear (i.e., branched or cyclic) over the same glutamate concentration range. The glutamate concentration dependence of the channel open probability, and of the event frequency, was used to explore two possible allosteric gating mechanisms in more detail.

Simulation approaches to ion channel structure-function relationships

1. Introduction 475

1.1 Ion channels 475

1.1.1 Gramicidin 476

1.1.2 Helix bundle channels 477

1.1.3 K channels 480

1.1.4 Porins 483

1.1.5 Nicotinic acetylcholine receptor 483

1.1.6 Physiological properties 483

1.2 Simulations 484

1.2.1 Atomistic versus mean-field simulations 484

2. Atomistic simulations 485

2.1 Modelling of ion-interaction parameters 485

2.1.1 Interatomic distances and the problem of ionic radii 486

2.1.2 Solvation energy 487

2.1.3 Hydration shells and coordination numbers 489

2.1.4 Parameters in common use and transferability 491

2.1.5 Summary 491

2.2 Water in pores versus bulk 491

2.2.1 Simple pore models 494

2.2.2 gA 495

2.2.3 Alm 496

2.2.4 LS36 (and LS24) 496

2.2.5 Nicotinic receptor M2δ5 497

2.2.6 Influenza A M2 497

2.2.7 K channels 497

2.2.8 nAChR 498

2.2.9 Porins 498

2.2.10 Relevance 499

2.2.11 Problems with simulations 501

2.3 Dynamics of ions in pores 503

2.3.1 Simple pore models 503

2.3.2 Helix bundles 504

2.3.3 gA and KcsA 505

2.4 Energetics of permeation and ion selectivity 509

2.4.1 Potential and free energy profiles 509

2.4.2 gA 510

2.4.3 α-Helix bundles 511

2.4.4 KcsA 512

2.4.5 Ion selectivity 514

2.4.6 Problems of estimating energetic profiles 515

2.5 Conformational changes 516

2.5.1 gA 516

2.5.2 Alm and LS3 516

2.5.3 KcsA 517

2.6 Protonation states 523

3. Coarse-grained simulations 524

3.1 Introduction 524

3.1.1 Predicting conductance magnitudes 525

3.2 Electro-diffusion: the Nernst–Planck approach 526

3.2.1 Calculating the potential profile from Poisson and PB theory 528

3.2.2 Calculating the potential profile from BD simulations 530

3.2.3 Combining Nernst–Planck and Poisson: PNP 530

3.3 Beyond PNP 532

3.4 BD simulations 532

3.4.1 Basic theory in ion channels 532

3.4.2 Incorporating the environment 533

3.5 Applications 535

3.5.1 Model systems 535

3.5.1.1 Solving the Poisson and PB equation for channel-like geometries 535

3.5.1.2 Comparing PB, PNP and BD 536

3.5.2 Applications to known structures 537

3.5.2.1 gA 537

3.5.2.2 Porin 539

3.5.2.3 LS3 540

3.5.2.4 Alm 542

3.5.2.5 nAChR 542

3.5.2.6 KcsA 543

3.6 pKa calculations 543

3.7 Selectivity 544

3.7.1 Anion/cation selectivity 545

3.7.2 Monovalent/divalent ion selectivity 545

4. Problems 546

4.1 Atomistic simulations 546

4.1.1 Problems 546

4.1.2 Parameters 548

4.2 BD 549

4.3 Mean-field simulations 549

5. Conclusions 550

5.1 Progress 550

5.2 The future 550

6. Acknowledgements 551

7. References 551

Ion channels are proteins that form ‘holes’ in membranes through which selected ions move passively down their electrochemical gradients. The ions move quickly, at (nearly) diffusion limited rates (ca. 107 ions s−1 per channel). Ion channels are central to many properties of cell membranes. Traditionally they have been the concern of neuroscientists, as they control the electrical properties of the membranes of excitable cells (neurones, muscle; Hille, 1992). However, it is evident that ion channels are present in many types of cell, not all of which are electrically excitable, from diverse organisms, including plants, bacteria and viruses (where they are involved in functions such as cell homeostasis) in addition to animals. Thus ion channels are of general cell biological importance. They are also of biomedical interest, as several dizeases (‘channelopathies’) have been described which are caused by changes in properties of a specific ion channel (Ashcroft, 2000). Moreover, passive diffusion channels for substances other than ions are common (porins, aquaporins), as are active membrane transport processes coupled to ion gradients or ATP hydrolysis. An understanding of ion channels may also provide a gateway to understanding these processes.

Amantadine blocks channel activity of the transmembrane segment of the NB protein from influenza B

NB is short auxiliary protein with ca. 100 amino acids, encoded in the viral genome of influenza B. It is believed to be similar to M2 from influenza A and Vpu from HIV-1 in that it demonstrates ion channel activity. Channels formed by the protein can be blocked by amantadine. We have synthesized the putative transmembrane segment of NB (IRG S20 IIITICVSL I30 VILIVFGCI A40 KIFI (NB, Lee)). Reconstituted in a lipid bilayer, the peptide shows channel activity. The addition of amantadine leads to dose-dependent loss of channel activity. Channel blocking is reversible. Channel behaviour of the peptide in the presence of amantadine is in accordance with findings for the intact channel. Thus, the synthetic transmembrane peptide captures the ion channel activity of the intact NB protein.

Proline-induced hinges in transmembrane helices: possible roles in ion channel gating

A number of ion channels contain transmembrane (TM) alpha-helices that contain proline-induced molecular hinges. These TM helices include the channel-forming peptide alamethicin (Alm), the S6 helix from voltage-gated potassium (Kv) channels, and the D5 helix from voltage-gated chloride (CLC) channels. For both Alm and KvS6, experimental data implicate hinge-bending motions of the helix in an aspect of channel gating. We have compared the hinge-bending motions of these TM helices in bilayer-like environments by multi-nanosecond MD simulations in an attempt to describe motions of these helices that may underlie possible modes of channel gating. Alm is an alpha-helical channel-forming peptide, which contains a central kink associated with a Gly-x-x-Pro motif in its sequence. Simulations of Alm in a TM orientation for 10 ns in an octane slab indicate that the Gly-x-x-Pro motif acts as a molecular hinge. The S6 helix from Shaker Kv channels contains a Pro-Val-Pro motif. Modeling studies and recent experimental data suggest that the KvS6 helix may be kinked in the vicinity of this motif. Simulations (10 ns) of an isolated KvS6 helix in an octane slab and in a POPC bilayer reveal hinge-bending motions. A pattern-matching approach was used to search for possible hinge-bending motifs in the TM helices of other ion channel proteins. This uncovered a conserved Gly-x-Pro motif in TM helix D5 of CLC channels. MD simulations of a model of hCLC1-D5 spanning an octane slab suggest that this channel also contains a TM helix that undergoes hinge-bending motion. In conclusion, our simulations suggest a model in which hinge-bending motions of TM helices may play a functional role in the gating mechanisms of several different families of ion channels.

Electrostatics studies and molecular dynamics simulations of a homology model of the Shaker K+ channel pore

A homology model of the pore domain of the Shaker K+ channel has been constructed using a bacterial K+ channel, KcsA, as a template structure. The model is in agreement with mutagenesis and sequence variability data. A number of structural features are conserved between the two channels, including a ring of tryptophan sidechains on the outer surface of the pore domain at the extracellular end of the helix bundle, and rings of acidic sidechains close to the extracellular mouth of the channel. One of these rings, that formed by four Asp447 sidechains at the mouth of the Shaker pore, is shown by pK(A) calculations to be incompletely ionized at neutral pH. The potential energy profile for a K+ ion moved along the central axis of the Shaker pore domain model selectivity filter reveals a shallow well, the depth of which is modulated by the ionization state of the Asp447 ring. This is more consistent with the high cation flux exhibited by the channel in its conductance value of 19 pS.

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 structure of the HIV-1 Vpu ion channel: modelling and simulation studies

Vpu is an 81 amino acid auxiliary protein in HIV-1 which exhibits channel activity. We used two homo-pentameric bundles with the helical transmembrane segments derived from FTIR spectroscopy in combination with a global molecular dynamics search protocol: (i) tryptophans (W) pointing into the pore, and (ii) W facing the lipids. Two equivalent bundles have been generated using a simulated annealing via a restrained molecular dynamics simulations (SA/MD) protocol. A fifth model was generated via SA/MD with all serines facing the pore. The latter model adopts a very stable structure during the 2 ns of simulation. The stability of the models with W facing the pore depends on the starting structure. A possible gating mechanism is outlined.

Amino acid distributions in integral membrane protein structures

Advances in structure determination of membrane proteins enable analysis of the propensities of amino acids in extramembrane versus transmembrane locations to be performed on the basis of structure rather than of sequence and predicted topology. Using 29 available structures of integral membrane proteins with resolutions better than 4 A the distributions of amino acids in the transmembrane domains were calculated. The results were compared to analysis based on just the sequences of the same transmembrane alpha-helices and significant differences were found. The distribution of residues between transmembrane alpha-helices and beta-strands was also compared. Large hydrophobic (Phe, Leu, Ile, Val) residues showed a clear preference for the protein surfaces facing the lipids for beta-barrels, but in alpha-helical proteins no such preference was seen, with these residues equally distributed between the interior and the surface of the protein. A notable exception to this was alanine, which showed a slight preference for the interior of alpha-helical membrane proteins. Aromatic residues were found to follow saddle-like distributions preferring to be located in the lipid/water interfaces. The resultant 'aromatic belts' were spaced more closely for beta-barrel than for alpha-helical membrane proteins. Charged residues could be shown to generally avoid surfaces facing the bilayer although they were found to occur frequently in the transmembrane region of beta-barrels. Indeed detailed comparison between alpha-helical and beta-barrel proteins showed many qualitative differences in residue distributions. This suggests that there may be subtle differences in the factors stabilising beta-barrels in bacterial outer membranes and alpha-helix bundles in all other membranes.

Side-chain ionization states in a potassium channel

KcsA is a bacterial K+ channel that is gated by pH. Continuum dielectric calculations on the crystal structure of the channel protein embedded in a low dielectric slab suggest that side chains E71 and D80 of each subunit, which lie adjacent to the selectivity filter region of the channel, form a proton-sharing pair in which E71 is neutral (protonated) and D80 is negatively charged at pH 7. When K+ ions are introduced into the system at their crystallographic positions the pattern of proton sharing is altered. The largest perturbation is for a K+ ion at site S3, i.e., interacting with the carbonyls of T75 and V76. The presence of multiple K+ ions in the filter increases the probability of E71 being ionized and of D80 remaining neutral (i.e., protonated). The ionization states of the protein side chains influence the potential energy profile experienced by a K+ ion as it is translated along the pore axis. In particular, the ionization state of the E71-D80 proton-sharing pair modulates the shape of the potential profile in the vicinity of the selectivity filter. Such reciprocal effects of ion occupancy on side-chain ionization states, and of side-chain ionization states on ion potential energy profiles will complicate molecular dynamics simulations and related studies designed to calculate ion permeation energetics.

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.

Membrane proteins: Aquaporins--channels without ions

Recently determined structures have shed new light on the way that aquaporins act as passive, but selective, pores for the transport of small molecules--such as water or glycerol--across membranes.

Voltage-dependent insertion of alamethicin at phospholipid/water and octane/water interfaces

Understanding the binding and insertion of peptides in lipid bilayers is a prerequisite for understanding phenomena such as antimicrobial activity and membrane-protein folding. We describe molecular dynamics simulations of the antimicrobial peptide alamethicin in lipid/water and octane/water environments, taking into account an external electric field to mimic the membrane potential. At cis-positive potentials, alamethicin does not insert into a phospholipid bilayer in 10 ns of simulation, due to the slow dynamics of the peptide and lipids. However, in octane N-terminal insertion occurs at field strengths from 0.33 V/nm and higher, in simulations of up to 100 ns duration. Insertion of alamethicin occurs in two steps, corresponding to desolvation of the Gln7 side chain, and the backbone of Aib10 and Gly11. The proline induced helix kink angle does not change significantly during insertion. Polyalanine and alamethicin form stable helices both when inserted in octane and at the water/octane interface, where they partition in the same location. In water, both polyalanine and alamethicin partially unfold in multiple simulations. We present a detailed analysis of the insertion of alamethicin into the octane slab and the influence of the external field on the peptide structure. Our findings give new insight into the mechanism of channel formation by alamethicin and the structure and dynamics of membrane-associated helices.

Hinges, swivels and switches: the role of prolines in signalling via transmembrane alpha-helices

Extracellular signals are transduced across membranes via conformational changes in the transmembrane domains (TMs) of ion channels and G-protein-coupled receptors (GPCRs). Experimental and simulation studies indicate that such conformational switches in transmembrane (alpha-helices can be generated by proline-containing motifs that form molecular hinges. Computational approaches tested on model channel-forming peptides (e.g. alamethicin) reveal functional mechanisms in gap-junction proteins (such as connexin) and voltage-gated K+ channels. Similarly, functionally important roles for proline-based switches in TM6 and TM7 were identified in GPCRs. However, hinges in transmembrane helices are not confined to proline-containing sequence motifs, as evidenced by a non-proline hinge in the M2 helix of the nicotinic acetylcholine receptor. This helix lines the pore and plays a key role in the gating of this channel.

Transmembrane peptide NB of influenza B: a simulation, structure, and conductance study

The putative transmembrane segment of the ion channel forming peptide NB from influenza B was synthesized by standard solid-phase peptide synthesis. Insertion into the planar lipid bilayer revealed ion channel activity with conductance levels of 20, 61, 107, and 142 pS in a 0.5 M KCl buffer solution. In addition, levels at -100 mV show conductances of 251 and 413 pS. A linear current-voltage relation reveals a voltage-independent channel formation. In methanol and in vesicles the peptide appears to adopt an alpha-helical-like structure. Computational models of alpha-helix bundles using N = 4, 5, and 6 NB peptides per bundle revealed water-filled pores after 1 ns of MD simulation in a solvated lipid bilayer. Calculated conductance values [using HOLE (Smart et al. (1997) Biophys. J. 72, 1109-1126)] of ca. 20, 60, and 90 pS, respectively, suggested that the multiple conductance levels seen experimentally must correspond to different degrees of oligomerization of the peptide to form channels.

Simulations of ion channels--watching ions and water move

Ion channels mediate electrical excitability in neurons and muscle. Three-dimensional structures for model peptide channels and for a potassium (K+) channel have been combined with computer simulations to permit rigorous exploration of structure-function relations of channels. Water molecules and ions within transbilayer pores tend to diffuse more slowly than in bulk solutions. In the narrow selectivity filter of the bacterial K+ channel (i.e. the region of the channel that discriminates between different species of ions) a column of water molecules and K+ ions moves in a concerted fashion. By combining atomistic simulations (in which all atoms of the channel molecule, water and ions are treated explicitly) with continuum methods (in which the description of the channel system is considerably simplified) it is possible to simulate some of the physiological properties of channels.

The nicotinic acetylcholine receptor: from molecular model to single-channel conductance

The nicotinic acetylcholine receptor (nAChR) is the archetypal ligand-gated ion channel. A model of the alpha7 homopentameric nAChR is described in which the pore-lining M2 helix bundle is treated atomistically and the remainder of the molecule is treated as a "low resolution" cylinder. The surface charge on the cylinder is derived from the distribution of charged amino acids in the amino acid sequence (excluding the M2 segments). This model is explored in terms of its predicted single-channel properties. Based on electrostatic potential profiles derived from the model, the one-dimensional Poisson-Nernst-Planck equation is used to calculate single-channel current/voltage curves. The predicted single-channel conductance is three times higher (ca. 150 pS) than that measured experimentally, and the predicted ion selectivity agrees with the observed cation selectivity of nAChR. Molecular dynamics (MD) simulations are used to estimate the self-diffusion coefficients (D) of water molecules within the channel. In the narrowest region of the pore, D is reduced ca. threefold relative to that of bulk water. Assuming that the diffusion of ions scales with that of water, this yields a revised prediction of the single-channel conductance (ca. 50 pS) in good agreement with the experimental value. We conclude that combining atomistic (MD) and continuum electrostatics calculations is a promising approach to bridging the gap between structure and physiology of ion channels.

Molecular contacts in the transmembrane c-subunit oligomer of F-ATPases identified by tryptophan substitution mutagenesis

When isolated in its monomeric form, subunit c of the proton transporting ATP synthase of Escherichia coli was shown to fold in a hairpin-like structure consisting of two hydrophobic membrane spanning helices and a short connecting hydrophilic loop. In the plasma membrane of Escherichia coli, however, about 9-12 c-subunit monomers form an oligomeric complex that functions in transmembrane proton conduction and in energy transduction to the catalytic F1 domain. The arrangement of the monomers and the molecular architecture of the complex were studied by tryptophan scanning mutagenesis and restrained MD simulations. Residues 12-24 of the N-terminal transmembrane segment of subunit c were individually substituted by the large and moderately hydrophobic tryptophan side chain. Effects on the activity of the mutant proteins were studied in selective growth experiments and various ATP synthase specific activity assays. The results identify potential intersubunit contacts and structurally non-distorted, accessible residues in the c-oligomer and add constraints to the arrangement of monomers in the oligomeric complex. Results from our mutagenesis experiments were interpreted in structural models of the c-oligomer that have been obtained by restrained MD simulations. Different stoichiometries and monomer orientations were applied in these calculations. A cylindrical complex consisting of 10 monomers that are arranged in two concentric rings with the N-terminal helices of the monomers located at the periphery shows the best match with the experimental data.

Structure-based prediction of the conductance properties of ion channels

The HOLE procedure allows the prediction of the absolute conductance of an ion channel model from its structure. The original prediction method uses an empirically corrected Ohmic method. It is most successful, with predictions being reliable to within a factor of two. A new modification of the procedure is presented in which the self-diffusion coefficients of water molecules from molecular dynamics simulation are used to replace the empirical correction factor. A "prediction" of the conductance for the porin OmpF by the new method is made and shown to be very close to the experimental value. HOLE also allows the prediction of the effect that the addition of non-electrolyte polymers will have on channel conductance. The method has great potential to yield structural information from data provided by single channel recordings but needs further validation by making measurements on channels of known structure. Preliminary results are given of single channel records establishing the effects of non-electrolytes on the conductance of gramicidin D channels. As an example of the potential uses of the procedure application is made to examine the oligomerization of alpha-toxin (alpha-hemolysin) channels. A model for the alpha-toxin hexamer, based on the crystal structure for the heptamer, is generated using molecular mechanics methods. The compatibility of the structures with single channel conductance data is assessed using HOLE.

Alamethicin channels in a membrane: molecular dynamics simulations

Alamethicin (Alm) is a 20 residue peptide which forms a kinked alpha-helix in membrane and membrane-mimetic environments. Ion channels formed by intramembraneous aggregates of Alm are thought to be formed by bundles of approximately parallel Alm helices surrounding a central bilayer pore. Different channel conductance levels correspond to different numbers of helices per bundle, ranging from N = 5 to N > 8. Calculation of the predicted pKA values of the ring of Glu18 sidechains at the C-terminal mouth of the pore suggests that at neutral pH most or all of these sidechains will remain protonated. Nanosecond molecular dynamics (MD) simulations of N = 5, 6, 7 and 8 bundles of Alm helices in a POPC bilayer have been run, corresponding to a total simulation time of 4 ns. These simulations explore the stability and conformational dynamics of these helix bundle channels when embedded in a full phospholipid bilayer in an aqueous environment. The structural and dynamic properties of water in these model channels are examined. As in earlier in vacuo simulations (J. Breed, R. Sankararamakrishnan, I. D. Kerr and M. S. P. Sansom, Biophys. J., 1996, 70, 1643) the dipole moments of water molecules within the pores are aligned antiparallel to the helix dipoles. This helps to contribute to the stability of the helix bundles.

Homology modeling and molecular dynamics simulation studies of an inward rectifier potassium channel

A homology model has been generated for the pore-forming domain of Kir6.2, a component of an ATP-sensitive K channel, based on the x-ray structure of the bacterial channel KcsA. Analysis of the lipid-exposed and pore-lining surfaces of the model reveals them to be compatible with the known features of membrane proteins and Kir channels, respectively. The Kir6.2 homology model was used as the starting point for nanosecond-duration molecular dynamics simulations in a solvated phospholipid bilayer. The overall drift from the model structure was comparable to that seen for KcsA in previous similar simulations. Preliminary analysis of the interactions of the Kir6.2 channel model with K(+) ions and water molecules during these simulations suggests that concerted single-file motion of K(+) ions and water through the selectivity filter occurs. This is similar to such motion observed in simulations of KcsA. This suggests that a single-filing mechanism is conserved between different K channel structures and may be robust to changes in simulation details. Comparison of Kir6.2 and KcsA suggests some degree of flexibility in the filter, thus complicating models of ion selectivity based upon a rigid filter.

Transmembrane domains of viral ion channel proteins: a molecular dynamics simulation study

Nanosecond molecular dynamics simulations in a fully solvated phospholipid bilayer have been performed on single transmembrane alpha-helices from three putative ion channel proteins encoded by viruses: NB (from influenza B), CM2 (from influenza C), and Vpu (from HIV-1). alpha-Helix stability is maintained within a core region of ca. 28 residues for each protein. Helix perturbations are due either to unfavorable interactions of hydrophobic residues with the lipid headgroups or to the need of the termini of short helices to extend into the surrounding interfacial environment in order to form H-bonds. The requirement of both ends of a helix to form favorable interactions with lipid headgroups and/or water may also lead to tilting and/or kinking of a transmembrane alpha-helix. Residues that are generally viewed as poor helix formers in aqueous solution (e.g., Gly, Ile, Val) do not destabilize helices, if located within a helix that spans a lipid bilayer. However, helix/bilayer mismatch such that a helix ends abruptly within the bilayer core destabilizes the end of the helix, especially in the presence of Gly and Ala residues. Hydrogen bonding of polar side-chains with the peptide backbone and with one another occurs when such residues are present within the bilayer core, thus minimizing the energetic cost of burying such side-chains.

Simulation studies on bacteriorhodopsin bundle of transmembrane alpha segments

Bacteriorhodopsin (BR) is a membrane protein which pumps protons through the plasma membrane. Seven transmembrane BR helical segments are subjected to simulation studies in order to investigate the packing process of transmembrane helices. A Monte Carlo simulated annealing protocol is employed to optimize the helix bundle system. Helix packing is optimized according to a semi-empirical potential mainly composed of six components: a bilayer potential, a crossing angle potential, a helix dipole potential, a helix-helix distance potential, a helix orientation potential and a helix-helix distance restraint potential (a loop potential). Necessary parameters are derived from theoretical studies and statistical analysis of experimentally determined protein structures. The structures from the simulations are compared with the experimentally determined structures in terms of geometry. The structures generated show similar shapes to the experimentally suggested structure even without the helix-helix distance restraint potential. However, the relative locations of individual helices were reproduced only when the helix-helix distance restraint potential was used with restraint conditions. Our results suggest that transmembrane helix bundles resembling those observed experimentally 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.

Structure and dynamics of the pore-lining helix of the nicotinic receptor: MD simulations in water, lipid bilayers, and transbilayer bundles

Multiple nanosecond duration molecular dynamics simulations on the pore-lining M2 helix of the nicotinic acetylcholine receptor reveal how its structure and dynamics change as a function of environment. In water, the M2 helix partially unfolds to form a molecular hinge in the vicinity of a central Leu residue that has been implicated in the mechanism of ion channel gating. In a phospholipid bilayer, either as a single transmembrane helix, or as part of a pentameric helix bundle, the M2 helix shows less flexibility, but still exhibits a kink in the vicinity of the central Leu. The single M2 helix tilts relative to the bilayer normal by 12 degrees, in agreement with recent solid state NMR data (Opella et al., Nat Struct Biol 6:374-379, 1999). The pentameric helix bundle, a model for the pore domain of the nicotinic receptor and for channels formed by M2 peptides in a bilayer, is remarkably stable over a 2-ns MD simulation in a bilayer, provided one adjusts the pK(A)s of ionizable residues to their calculated values (when taking their environment into account) before starting the simulation. The resultant transbilayer pore shows fluctuations at either mouth which transiently close the channel. Proteins 2000;39:47-55.

Membrane simulations: bigger and better?

Molecular dynamics simulations of biological membranes have come of age. Simulations of pure lipid bilayers are extending our understanding of both optimal simulation procedures and the detailed structural dynamics of lipids in these systems. Simulation methods established using simple bilayer-embedded peptides are being extended to a wide range of membrane proteins and membrane protein models, and are beginning to reveal some of the complexities of membrane protein structural dynamics and their relationship to biological function.

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.

Protonation of lysine residues inverts cation/anion selectivity in a model channel

A dimeric alamethicin analog with lysine at position 18 in the sequence (alm-K18) was previously shown to form stable anion-selective channels in membranes at pH 7.0 [Starostin, A. V., R. Butan, V. Borisenko, D. A. James, H. Wenschuh, M. S. Sansom, and G. A. Woolley. 1999. Biochemistry. 38:6144-6150]. To probe the charge state of the conducting channel and how this might influence cation versus anion selectivity, we performed a series of single-channel selectivity measurements at different pH values. At pH 7.0 and below, only anion-selective channels were found with P(K(+))/P(Cl(-)) = 0. 25. From pH 8-10, a mixture of anion-selective, non-selective, and cation-selective channels was found. At pH > 11 only cation-selective channels were found with P(K(+))/P(Cl(-)) = 4. In contrast, native alamethicin-Q18 channels (with Gln in place of Lys at position 18) were cation-selective (P(K(+))/P(Cl(-)) = 4) at all pH values. Continuum electrostatics calculations were then carried out using an octameric model of the alm-K18 channel embedded in a low dielectric slab to simulate a membrane. Although the calculations can account for the apparent pK(a) of the channel, they fail to correctly predict the degree of selectivity. Although a switch from cation- to anion-selectivity as the channel becomes protonated is indicated, the degree of anion-selectivity is severely overestimated, suggesting that the continuum approach does not adequately represent some aspect of the electrostatics of permeation in these channels. Side-chain conformational changes upon protonation, conformational changes, and deprotonation caused by permeating cations and counterion binding by lysine residues upon protonation are considered as possible sources of the overestimation.

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.

Exploring models of the influenza A M2 channel: MD simulations in a phospholipid bilayer

The M2 protein of influenza A virus forms homotetrameric helix bundles, which function as proton-selective channels. The native form of the protein is 97 residues long, although peptides representing the transmembrane section display ion channel activity, which (like the native channel) is blocked by the antiviral drug amantadine. As a small ion channel, M2 may provide useful insights into more complex channel systems. Models of tetrameric bundles of helices containing either 18 or 22 residues have been simulated while embedded in a fully hydrated 1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphatidylcholine bilayer. Several different starting models have been used. These suggest that the simulation results, at least on a nanosecond time scale, are sensitive to the exact starting structure. Electrostatics calculations carried out on a ring of four ionizable aspartate residues at the N-terminal mouth of the channel suggest that at any one time, only one will be in a charged state. Helix bundle models were mostly stable over the duration of the simulation, and their helices remained tilted relative to the bilayer normal. The M2 helix bundles form closed channels that undergo breathing motions, alternating between a tetramer and a dimer-of-dimers structure. Under these conditions either the channel forms a pocket of trapped waters or it contains a column of waters broken predominantly at the C-terminal mouth of the pore. These waters exhibit restricted motion in the pore and are effectively "frozen" in a way similar to those seen in previous simulations of a proton channel formed by a four-helix bundle of a synthetic leucine-serine peptide (, Biophys. J. 77:2400-2410).

Structure and dynamics of K channel pore-lining helices: a comparative simulation study

Isolated pore-lining helices derived from three types of K-channel have been analyzed in terms of their structural and dynamic features in nanosecond molecular dynamics (MD) simulations while spanning a lipid bilayer. The helices were 1) M1 and M2 from the bacterial channel KcsA (Streptomyces lividans), 2) S5 and S6 from the voltage-gated (Kv) channel Shaker (Drosophila melanogaster), and 3) M1 and M2 from the inward rectifier channel Kir6.2 (human). In the case of the Kv and Kir channels, for which x-ray structures are not known, both short and long models of each helix were considered. Each helix was incorporated into a lipid bilayer containing 127 palmitoyloleoylphosphatidylcholine molecules, which was solvated with approximately 4000 water molecules, yielding approximately 20, 000 atoms in each system. Nanosecond MD simulations were used to aid the definition of optimal lengths for the helix models from Kv and Kir. Thus the study corresponds to a total simulation time of 10 ns. The inner pore-lining helices (M2 in KcsA and Kir, S6 in Shaker) appear to be slightly more flexible than the outer pore-lining helices. In particular, the Pro-Val-Pro motif of S6 results in flexibility about a molecular hinge, as was suggested by previous in vacuo simulations (, Biopolymers. 39:503-515). Such flexibility may be related to gating in the corresponding intact channel protein molecules. Analysis of H-bonds revealed interactions with both water and lipid molecules in the water/bilayer interfacial region. Such H-bonding interactions may lock the helices in place in the bilayer during the folding of the channel protein (as is implicit in the two-stage model of membrane protein folding). Aromatic residues at the extremities of the helices underwent complex motions on both short (<10 ps) and long (>100 ps) time scales.

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.

Molecular dynamics of synthetic leucine-serine ion channels in a phospholipid membrane

Molecular dynamics calculations were carried out on models of two synthetic leucine-serine ion channels: a tetrameric bundle with sequence (LSLLLSL)(3)NH(2) and a hexameric bundle with sequence (LSSLLSL)(3)NH(2). Each protein bundle is inserted in a palmitoyloleoylphosphatidylcholine bilayer membrane and solvated by simple point charge water molecules inside the pore and at both mouths. Both systems appear to be stable in the absence of an electric field during the 4 ns of molecular dynamics simulation. The water motion in the narrow pore of the four-helix bundle is highly restricted and may provide suitable conditions for proton transfer via a water wire mechanism. In the wider hexameric pore, the water diffuses much more slowly than in bulk but is still mobile. This, along with the dimensions of the pore, supports the observation that this peptide is selective for monovalent cations. Reasonable agreement of predicted conductances with experimentally determined values lends support to the validity of the simulations.

Putting the parts together

The structures of several different domains and subunits of potassium channels have recently been solved. Reassembling these fragments into a working model of an intact voltage-gated channel will be a major challenge.

The mechanism of channel formation by alamethicin as viewed by molecular dynamics simulations

Alamethicin is a 20-residue channel-forming peptide that forms a stable amphipathic alpha-helix in membrane and membrane-mimetic environments. This helix contains a kink induced by a central Gly-X-X-Pro sequence motif. Alamethicin channels are activated by a cis positive transbilayer voltage. Channel activation is suggested to correspond to voltage-induced insertion of alamethicin helices in the bilayer. Alamethicin forms multi-conductance channels in lipid bilayers. These channels are formed by parallel bundles of transmembrane helices surrounding a central pore. A change in the number of helices per bundle switches the single channel conductance level. Molecular dynamics simulations of alamethicin in a number of different environments have been used to explore its channel-forming properties. These simulations include: (i) alamethicin in solution in water and in methanol; (ii) a single alamethicin helix at the surface of a phosphatidylcholine bilayer; (iii) single alamethicin helices spanning a phosphatidylcholine bilayer; and (iv) channels formed by bundles of 5, 6, 7 or 8 alamethicin helices spanning a phosphatidylcholine bilayer. The total simulation time is c. 30 ns. Thus, these simulations provide a set of dynamic snapshots of a possible mechanism of channel formation by this peptide.

Engineering charge selectivity in alamethicin channels

The peptide alamethicin provides a system for engineering ion channel charge selectivity. To define alamethicin charge selectivity experimentally, we measured single-channel current-voltage relationships in KCl gradients using covalently linked peptide dimers. Two factors were found to contribute to the charge selectivity of these channels: (i) the ionic strength of the surrounding solutions; and (ii) the distribution of fixed charge on the peptide. Native alamethicin channels exhibited either cation selectivity or anion selectivity depending on which end of the channel was at the low salt side of the membrane. When the glutamine residue at position 18 in the sequence was replaced with a lysine residue, an anion-selective channel was obtained regardless of which end of the channel was at the low salt side of the membrane.

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.

Effective diffusion coefficients of K+ and Cl- ions in ion channel models

We have used molecular dynamics simulations, corresponding to a total simulation time of 11 ns, to investigate the effective short-time local diffusion coefficient of potassium and chloride ions in a series of model ion channels. These models, which include channels formed by the fungal peptide alamethicin, by a synthetic leucine-serine peptide, and by the pore-lining M2 helix bundle of the nicotinic acetylcholine receptor, have a range of different secondary structures, diameters and hydrophobicities. We find that the diffusion coefficients of both ions are appreciably reduced in the narrower channels, the extent of the reduction being similar for both the anionic and cationic species. This suggests that a difference in mobility cannot be the source of the ion selectivity exhibited by some of the channels (for example, the leucine-serine peptide). We find no evidence for a reduction in mobility of either ion in the nAChR model. These results are broadly in line with a previous similar study of Na+ ions, and may be useful in Poisson-Nernst-Planck, Eyring rate theory or Brownian dynamics calculations of channel conductance.

Surface binding of alamethicin stabilizes its helical structure: molecular dynamics simulations

Alamethicin is an amphipathic alpha-helical peptide that forms ion channels. An early event in channel formation is believed to be the binding of alamethicin to the surface of a lipid bilayer. Molecular dynamics simulations are used to compare the structural and dynamic properties of alamethicin in water and alamethicin bound to the surface of a phosphatidylcholine bilayer. The bilayer surface simulation corresponded to a loosely bound alamethicin molecule that interacted with lipid headgroups but did not penetrate the hydrophobic core of the bilayer. Both simulations started with the peptide molecule in an alpha-helical conformation and lasted 2 ns. In water, the helix started to unfold after approximately 300 ps and by the end of the simulation only the N-terminal region of the peptide remained alpha-helical and the molecule had collapsed into a more compact form. At the surface of the bilayer, loss of helicity was restricted to the C-terminal third of the molecule and the rod-shaped structure of the peptide was retained. In the surface simulation about 10% of the peptide/water H-bonds were replaced by peptide/lipid H-bonds. These simulations suggest that some degree of stabilization of an amphipathic alpha-helix occurs at a bilayer surface even without interactions between hydrophobic side chains and the acyl chain core of the bilayer.

An anion-selective analogue of the channel-forming peptide alamethicin

The peptide alamethicin self-assembles to form helix bundle ion channels in membranes. Previous macroscopic measurements have shown that these channels are mildly cation-selective. Models indicate that a source of cation selectivity is a zone of partial negative charge toward the C-terminal end of the peptide. We synthesized an alamethicin derivative with a lysine in this zone (replacing the glutamine at position 18 in the sequence). Microscopic (single-channel) measurements demonstrate that dimeric alamethicin-lysine18 (alm-K18) forms mildly anion-selective channels under conditions where channels formed by the parent peptide are cation-selective. Long-range electrostatic interactions can explain the inversion of ion selectivity and the conductance properties of alamethicin channels.

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.

An extensive and diverse gene family of nicotinic acetylcholine receptor alpha subunits in Caenorhabditis elegans

Using reverse transcription-polymerase chain reactions the transcription of eight novel candidate nicotinic acetylcholine receptor (nAChR) alpha subunit genes has been demonstrated in the nematode Caenorhabditis elegans. Together with five other alpha subunit genes described elsewhere by ourselves (unc-38) and other workers (deg-3, acr-4, Ce21 and acr-6), this is now the largest known family of nAChR alpha subunit genes in a single species. By homology we have identified four groups of alpha subunits: DEG-3-like; ACR-16[Ce21]-like; UNC-38-like and ACR-8-like. Five C. elegans nAChR alpha subunits contain a modification in loop C of the ACh binding site in which the normally conserved Tyr-x-Cys-Cys, is replaced by a distinct motif (Tyr-x-x-Cys-Cys). Variation is also found in the channel lining M2 regions, including the replacement in four subunits of the highly conserved leucine at the 9' position by valine and most notably, the replacement in all ACR-8-like subunits of the highly conserved glutamic acid at the -1' position by histidine. Restrained molecular dynamics simulations have been used to generate homo-pentameric M2 helix bundle models for alpha subunits and possible functional implications examined. The calculated electrostatic potential energy profile for the M2 region of ACR-8 differs strikingly from that of ACR-16[Ce21] largely due to the presence of histidine at the -1' position, suggesting a possible perturbation of nAChR channel action permeability in the presence of this subunit type.

Defining the transmembrane helix of M2 protein from influenza A by molecular dynamics simulations in a lipid bilayer

Integral membrane proteins containing at least one transmembrane (TM) alpha-helix are believed to account for between 20% and 30% of most genomes. There are several algorithms that accurately predict the number and position of TM helices within a membrane protein sequence. However, these methods tend to disagree over the beginning and end residues of TM helices, posing problems for subsequent modeling and simulation studies. Molecular dynamics (MD) simulations in an explicit lipid and water environment are used to help define the TM helix of the M2 protein from influenza A virus. Based on a comparison of the results of five different secondary structure prediction algorithms, three different helix lengths (an 18mer, a 26mer, and a 34mer) were simulated. Each simulation system contained 127 POPC molecules plus approximately 3500-4700 waters, giving a total of approximately 18,000-21,000 atoms. Two simulations, each of 2 ns duration, were run for the 18mer and 26mer, and five separate simulations were run for the 34mer, using different starting models generated by restrained in vacuo MD simulations. The total simulation time amounted to 11 ns. Analysis of the time-dependent secondary structure of the TM segments was used to define the regions that adopted a stable alpha-helical conformation throughout the simulation. This analysis indicates a core TM region of approximately 20 residues (from residue 22 to residue 43) that remained in an alpha-helical conformation. Analysis of atomic density profiles suggested that the 18mer helix revealed a local perturbation of the lipid bilayer. Polar side chains on either side of this region form relatively long-lived H-bonds to lipid headgroups and water molecules.

An alamethicin channel in a lipid bilayer: molecular dynamics simulations

We present the results of 2-ns molecular dynamics (MD) simulations of a hexameric bundle of Alm helices in a 1-palmitoyl-2-oleoylphosphatidylcholine bilayer. These simulations explore the dynamic properties of a model of a helix bundle channel in a complete phospholipid bilayer in an aqueous environment. We explore the stability and conformational dynamics of the bundle in a phospholipid bilayer. We also investigate the effect on bundle stability of the ionization state of the ring of Glu18 side chains. If all of the Glu18 side chains are ionised, the bundle is unstable; if none of the Glu18 side chains are ionized, the bundle is stable. pKA calculations suggest that either zero or one ionized Glu18 is present at neutral pH, correlating with the stable form of the helix bundle. The structural and dynamic properties of water in this model channel were examined. As in earlier in vacuo simulations (Breed et al., 1996 .Biophys. J. 70:1643-1661), the dipole moments of water molecules within the pore were aligned antiparallel to the helix dipoles. This contributes to the stability of the helix bundle.

Ion channels: structure of a molecular brake

A combination of crystallographic and mutagenesis studies on the HERG K+ channel, a key determinant of cardiac excitability, has suggested how the protein's extramembraneous amino-terminal domain might act as a 'molecular brake' that slows down channel deactivation.

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.

Peptide-bilayer interactions: simulations of dermaseptin B, an antimicrobial peptide

Dermaseptins, a family of antimicrobial peptides, are believed to act by forming amphipathic alpha-helices which associate with the cell membrane, leading to its permeabilisation and disruption. A simple mean field method is described for simulation of the interactions of peptides with lipid bilayers which includes an approximate representation of the electrostatic effects of the head-group region of the bilayer. Starting from an atomistic model of a PC phospholipid bilayer we calculate an average electrostatic potential along the bilayer normal. By combining the interaction of the peptide with this electrostatic potential and with the hydrophobic core of the membrane we arrive at a more complete description of peptide-bilayer energetics than would be obtained using sidechain hydrophobicities alone. Using this interaction potential in MD simulations of the frog skin peptide dermaseptin B reveals that the lipid bilayer stabilises the alpha-helical conformation of the peptide. This is in agreement with FTIR data. A surface associated orientation thus appears to be the most stable arrangement of the peptide, at least at zero ionic strength and without taking account of possible peptide-peptide interactions.