Side-chain ionization states in a potassium channel

The catalytic effect of H2O on the hydrolysis of CO32−in hydrated clusters and its implication in the humidity driven CO2air capture

Incomplete hydration shells facilitate the hydrolysis of CO₃²⁻to be not only thermodynamically favorable but also kinetically favorable.

Unraveling the mechanism of selective ion transport in hydrophobic subnanometer channels

Recently reported synthetic organic nanopore (SONP) can mimic a key feature of natural ion channels, i.e., selective ion transport. However, the physical mechanism underlying the K(+)/Na(+) selectivity for the SONPs is dramatically different from that of natural ion channels. To achieve a better understanding of the selective ion transport in hydrophobic subnanometer channels in general and SONPs in particular, we perform a series of ab initio molecular dynamics simulations to investigate the diffusivity of aqua Na(+) and K(+) ions in two prototype hydrophobic nanochannels: (i) an SONP with radius of 3.2 Å, and (ii) single-walled carbon nanotubes (CNTs) with radii of 3-5 Å (these radii are comparable to those of the biological potassium K(+) channels). We find that the hydration shell of aqua Na(+) ion is smaller than that of aqua K(+) ion but notably more structured and less yielding. The aqua ions do not lower the diffusivity of water molecules in CNTs, but in SONP the diffusivity of aqua ions (Na(+) in particular) is strongly suppressed due to the rugged inner surface. Moreover, the aqua Na(+) ion requires higher formation energy than aqua K(+) ion in the hydrophobic nanochannels. As such, we find that the ion (K(+) vs. Na(+)) selectivity of the (8, 8) CNT is ∼20× higher than that of SONP. Hence, the (8, 8) CNT is likely the most efficient artificial K(+) channel due in part to its special interior environment in which Na(+) can be fully solvated, whereas K(+) cannot. This work provides deeper insights into the physical chemistry behind selective ion transport in nanochannels.

Molecular dynamics simulations of membrane proteins

Membrane protein structures are underrepresented in the Protein Data Bank (PDB) due to difficulties associated with expression and crystallization. As such, it is one area where computational studies, particularly Molecular Dynamics (MD) simulations, can provide useful additional information. Recently, there has been substantial progress in the simulation of lipid bilayers and membrane proteins embedded within them. Initial efforts at simulating membrane proteins embedded within a lipid bilayer were relatively slow and interactive processes, but recent advances now mean that the setup and running of membrane protein simulations is somewhat more straightforward, though not without its problems. In this chapter, we outline practical methods for setting up and running MD simulations of a membrane protein embedded within a lipid bilayer and discuss methodologies that are likely to contribute future improvements.

Energetics of Multi-Ion Conduction Pathways in Potassium Ion Channels

Potassium ion channels form pores in cell membranes, allowing potassium ions through while preventing the passage of sodium ions. Despite numerous high-resolution structures, it is not yet possible to relate their structure to their single molecule function other than at a qualitative level. Over the past decade, there has been a concerted effort using molecular dynamics to capture the thermodynamics and kinetics of conduction by calculating potentials of mean force (PMF). These can be used, in conjunction with the electro-diffusion theory, to predict the conductance of a specific ion channel. Here, we calculate seven independent PMFs, thereby studying the differences between two potassium ion channels, the effect of the CHARMM CMAP forcefield correction, and the sensitivity and reproducibility of the method. Thermodynamically stable ion-water configurations of the selectivity filter can be identified from all the free energy landscapes, but the heights of the kinetic barriers for potassium ions to move through the selectivity filter are, in nearly all cases, too high to predict conductances in line with experiment. This implies it is not currently feasible to predict the conductance of potassium ion channels, but other simpler channels may be more tractable.

Conductance properties of the inwardly rectifying channel, Kir3.2: molecular and Brownian dynamics study

Using the recently unveiled crystal structure, and molecular and Brownian dynamics simulations, we elucidate several conductance properties of the inwardly rectifying potassium channel, Kir3.2, which is implicated in cardiac and neurological disorders. We show that the pore is closed by a hydrophobic gating mechanism similar to that observed in Kv1.2. Once open, potassium ions move into, but not out of, the cell. The asymmetrical current-voltage relationship arises from the lack of negatively charged residues at the narrow intracellular mouth of the channel. When four phenylalanine residues guarding the intracellular gate are mutated to glutamate residues, the channel no longer shows inward rectification. Inward rectification is restored in the mutant Kir3.2 when it becomes blocked by intracellular Mg(2+). Tertiapin, a polypeptide toxin isolated from the honey bee, is known to block several subtypes of the inwardly rectifying channels with differing affinities. We identify critical residues in the toxin and Kir3.2 for the formation of the stable complex. A lysine residue of tertiapin protrudes into the selectivity filter of Kir3.2, while two other basic residues of the toxin form hydrogen bonds with acidic residues located just outside the channel entrance. The depth of the potential of mean force encountered by tertiapin is -16.1kT, thus indicating that the channel will be half-blocked by 0.4μM of the toxin.

Protonation state of E71 in KcsA and its role for channel collapse and inactivation

The prototypical prokaryotic potassium channel KcsA alters its pore depending on the ambient potassium; at high potassium, it exists in a conductive form, and at low potassium, it collapses into a nonconductive structure with reduced ion occupancy. We present solid-state NMR studies of KcsA in which we test the hypothesis that an important channel-inactivation process, known as C-type inactivation, proceeds via a state similar to this collapsed state. We test this using an inactivation-resistant mutant E71A, and show that E71A is unable to collapse its pore at both low potassium and low pH, suggesting that the collapsed state is structurally similar to the inactivated state. We also show that E71A has a disordered selectivity filter. Using site-specific K(+) titrations, we detect a local change at E71 that is coupled to channel collapse at low K(+). To gain more insight into this change, we site specifically measure the chemical shift tensors of the side-chain carboxyls of E71 and its hydrogen bond partner D80, and use the tensors to assign protonation states to E71 and D80 at high K(+) and neutral pH. Our measurements show that E71 is protonated at pH 7.5 and must have an unusually perturbed pK(a) (> 7.5) suggesting that the change at E71 is a structural rearrangement rather than a protonation event. The results offer new mechanistic insights into why the widely used mutant KcsA-E71A does not inactivate and establish the ambient K(+) level as a means to populate the inactivated state of KcsA in a controlled way.

Synthetic cation-selective nanotube: permeant cations chaperoned by anions

The ability to design ion-selective, synthetic nanotubes which mimic biological ion channels may have significant implications for the future treatment of bacteria, diseases, and as ultrasensitive biosensors. We present the design of a synthetic nanotube made from carbon atoms that selectively allows monovalent cations to move across and rejects all anions. The cation-selective nanotube mimics some of the salient properties of biological ion channels. Before practical nanodevices are successfully fabricated it is vital that proof-of-concept computational studies are performed. With this in mind we use molecular and stochastic dynamics simulations to characterize the dynamics of ion permeation across a single-walled (10, 10), 36 Å long, carbon nanotube terminated with carboxylic acid with an effective radius of 5.08 Å. Although cations encounter a high energy barrier of 7 kT, its height is drastically reduced by a chloride ion in the nanotube. The presence of a chloride ion near the pore entrance thus enables a cation to enter the pore and, once in the pore, it is chaperoned by the resident counterion across the narrow pore. The moment the chaperoned cation transits the pore, the counterion moves back to the entrance to ferry another ion. The synthetic nanotube has a high sodium conductance of 124 pS and shows linear current-voltage and current-concentration profiles. The cation-anion selectivity ratio ranges from 8 to 25, depending on the ionic concentrations in the reservoirs.

Symmetry-restrained molecular dynamics simulations improve homology models of potassium channels

Most crystallized homo-oligomeric ion channels are highly symmetric, which dramatically decreases conformational space and facilitates building homology models (HMs). However, in molecular dynamics (MD) simulations channels deviate from ideal symmetry and accumulate thermal defects, which complicate the refinement of HMs using MD. In this work we evaluate the ability of symmetry constrained MD simulations to improve HMs accuracy, using an approach conceptually similar to Critical Assessment of techniques for protein Structure Prediction (CASP) competition: build HMs of channels with known structure and evaluate the efficiency of proposed methods in improving HMs accuracy (measured as deviation from experimental structure). Results indicate that unrestrained MD does not improve the accuracy of HMs, instantaneous symmetrization improves accuracy but not stability of HMs during subsequent unrestrained MD, while gradually imposing symmetry constraints improves both accuracy (by 5-50%) and stability of HMs. Moreover, accuracy and stability are strongly correlated, making stability a reliable criterion in predicting the accuracy of new HMs. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

Proton inhibition of unitary currents of vanilloid receptors

Protons, which are released during inflammation and injury, regulate many receptors and ion channels involved in pain transduction, including capsaicin channels (transient receptor potential vanilloid receptors 1). Whereas extracellular acidification both sensitizes and directly activates the channel, it also causes concomitant reduction of the unitary current amplitudes. Here, we investigate the mechanisms and molecular basis of this inhibitory effect of protons on channel conductance. Single-channel recordings showed that the unitary current amplitudes decreased with extracellular pH in a dose-dependent manner, consistent with a model in which protons bind to a site within the channel with an apparent pKa of ∼6. The inhibition was voltage dependent, ∼65% at −60 mV and 37% at +60 mV when pH was reduced from 7.4 to 5.5. The unitary current amplitudes reached saturation at [K⁺] ≥ 1 M, and notably the maximum amplitudes did not converge with different pHs, inconsistent with a blockade model based on surface charge screening or competitive inhibition of permeating ions. Mutagenesis experiments uncovered two acidic residues critical for proton inhibition, one located at the pore entrance and the other on the pore helix. Based on homology to the KcsA structure, the two acidic residues, along with another basic residue also on the pore helix, could form a triad interacting with each other through extensive hydrogen bonds and electrostatic contacts, suggesting that protons may mediate the interactions between the selectivity filter and pore helix, thereby altering the local structure in the filter region and consequently the conductance of the channel.

Molecular modeling and simulation studies of ion channel structures, dynamics and mechanisms

Ion channels are integral membrane proteins that enable selected ions to flow passively across membranes. Channel proteins have been the focus of computational approaches to relate their three-dimensional (3D) structure to their physiological function. We describe a number of computational tools to model ion channels. Homology modeling may be used to construct structural models of channels based on available X-ray structures. Electrostatics calculations enable an approximate evaluation of the energy profile of an ion passing through a channel. Molecular dynamics simulations and free-energy calculations provide information on the thermodynamics and kinetics of channel function.

Monte Carlo simulation for statistical mechanics model of ion-channel cooperativity in cell membranes

Voltage-gated ion channels are key molecules for the generation and propagation of electrical signals in excitable cell membranes. The voltage-dependent switching of these channels between conducting and nonconducting states is a major factor in controlling the transmembrane voltage. In this study, a statistical mechanics model of these molecules has been discussed on the basis of a two-dimensional spin model. A new Hamiltonian and a new Monte Carlo simulation algorithm are introduced to simulate such a model. It was shown that the results well match the experimental data obtained from batrachotoxin-modified sodium channels in the squid giant axon using the cut-open axon technique.

Intrinsic aqueduct orifices facilitate K+ channel gating

The ion-conducting pore of potassium channels, which can open and close to regulate ion passage, was at long thought to be a one-dimensional pore structure with a water-filled central cavity. Here, we find four orifices in the KcsA potassium channel, which are perpendicular to the pore and stretch out from the cavity. Equilibrium molecular dynamics simulations show that water molecules can flow between the cavity and orifices. Targeted molecular dynamics simulations show that during the opening process, water molecules can move into the cavity through the orifices to facilitate channel gating, whereas blocking the aqueduct orifices makes the channel difficult to open.

TransPath: a computational method for locating ion transit pathways through membrane proteins

The finely tuned structures of membrane channel proteins allow selective passage of ions through the available aqueous pores. To understand channel function, it is crucial to locate the pores and study their physical and chemical properties. Here, we propose a new pore-searching algorithm (TransPath), which uses the Configurational Bias Monte Carlo (CBMC) method to generate transmembrane trajectories driven by both geometric and electrostatic features. The trajectories are binned into groups determined by a vector distance criterion. From each group, a representative trajectory is selected based on the Rosenbluth weight, and the geometrically optimal path is obtained by simulated annealing. Candidate ion pathways can then be determined by analysis of the radius and potential profiles. The proposed method and its implementation are illustrated using the bacterial KcsA potassium channel as an example. The procedure is then applied to the more complex structures of the bacterial E. coli chloride channel homolog and a homology model of the ClC-0 channel.

Conformational changes in the selectivity filter of the open-state KcsA channel: an energy minimization study

Potassium channels switch between closed and open conformations and selectively conduct K(+) ions. There are at least two gates. The TM2 bundle at the intracellular site is the primary gate of KcsA, and rearrangements at the selectivity filter (SF) act as the second gate. The SF blocks ion flow via an inactivation process similar to C-type inactivation of voltage-gated K(+) channels. We recently generated the open-state conformation of the KcsA channel. We found no major, possibly inactivating, structural changes in the SF associated with this massive inner-pore rearrangement, which suggests that the gates might act independently. Here we energy-minimize the open state of wild-type and mutant KcsA, validating in silico structures of energy-minimized SFs by comparison with crystallographic structures, and use these data to gain insight into how mutation, ion depletion, and K(+) to Na(+) substitution influence SF conformation. Both E71 or D80 protonations/mutations and the presence/absence of protein-buried water molecule(s) modify the H-bonding network stabilizing the P-loops, spawning numerous SF conformations. We find that the inactivated state corresponds to conformations with a partially unoccupied or an entirely empty SF. These structures, involving modifications in all four P-loops, are stabilized by H-bonds between amide H and carbonyl O atoms from adjacent P-loops, which block ion passage. The inner portions of the P-loops are more rigid than the outer parts. Changes are localized to the outer binding sites, with innermost site S4 persisting in the inactivated state. Strong binding by Na(+) locally contracts the SF around Na(+), releasing ligands that do not participate in Na(+) coordination, and occluding the permeation pathway. K(+) selectivity primarily appears to arise from the inability of the SF to completely dehydrate Na(+) ions due to basic structural differences between liquid water and the "quasi-liquid" SF matrix.

Conformational changes and gating at the selectivity filter of potassium channels

The translocation of ions and water across cell membranes is a prerequisite for many of life's processes. K(+) channels are a diverse family of integral membrane proteins through which K(+) can pass selectively. There is an ongoing debate about the nature of conformational changes associated with the opening and closing and conductive and nonconductive states of potassium (K(+)) channels. These changes depend on the membrane potential, the K(+) concentration gradient, and large scale motions of transmembrane helices and associated residues. Experiments also suggest that local structural changes in the selectivity filter may act as the dominant gate referred to as C-type inactivation. Herein we present an extensive computational study on KirBac, which supports the existence of a physical gate or constriction in the selectivity filter (SF) of K(+) channels. Our computations identify a new selectivity filter structure, which is likely associated with C-type inactivation. Specifically, the four peptide chains that comprise the filter adopt an unusual structure in which their dihedrals alternate between left- and right-handed Ramachandran angles, which also justifies the need for conservation of glycine in the K(+) selectivity filter, since it is the only residue able to play this bifunctional role.

Molecular dynamics simulations of membrane proteins

Membrane protein structures are underrepresented in the Protein Data Bank (PDB) because of difficulties associated with expression and crystallization. As such, it is one area in which computational studies, particularly molecular dynamics (MD), can provide useful additional information. Recently, there has been substantial progress in the simulation of lipid bilayers and membrane proteins embedded within them. Initial efforts at simulating membrane proteins embedded within a lipid bilayer were relatively slow and interactive processes, but recent advances now mean that the setup and running of membrane protein simulations is somewhat more straightforward, although not without its problems. In this chapter, we outline practical methods for setting up and running MD simulations of a membrane protein embedded within a lipid bilayer and discuss methodologies that are likely to contribute future improvements.

Cooperative transport in a potassium ion channel

Our current understanding of ion permeation through the selectivity filter of the KcsA potassium channel is based on the concept of a multi-ion transport mechanism. The details of this concerted movement, however, are not well understood. In the present paper we report on molecular dynamics simulations which provides new insights. It is shown that ion translocation is based on the collective hopping of ions and water molecules which is mediated by the flexible charged carbonyl groups lining the backbone of the pore. In particular, there is strong evidence for pairwise translocations where one ion and one water molecule form a bound state. We suggest a physical explanation of the observed phenomena employing a simple lattice model. It is argued that the water molecules can act as rectifiers during the hopping of ion-water pairs.

The protonation state of the Glu-71/Asp-80 residues in the KcsA potassium channel: a first-principles QM/MM molecular dynamics study

Although a few x-ray structures of the KcsA K(+) channel have been crystallized several issues concerning the mechanisms of the ionic permeation and the protonation state of the selectivity filter ionizable side chains are still open. Using a first-principles quantum mechanical/molecular mechanical simulation approach, we have investigated the protonation state of Glu-71 and Asp-80, two important residues located in the vicinity of the selectivity filter. Results from the dynamics show that a proton is shared between the two residues, with a slight preference for Glu-71. The proton is found to exchange on the picosecond timescale, an interesting phenomenon that cannot be observed in classical molecular dynamics. Simulations of different ionic loading states of the filter show that the probability for the proton transfer is correlated with the filter occupancy. In addition, the Glu-71/Asp-80 pair is able to modulate the potential energy profile experienced by a K(+) ion as it translates along the pore axis. These theoretical predictions, along with recent experimental results, suggest that changes of the filter structure could be associated with a shift in the Glu-Asp protonation state, which in turn would influence the ion translocation.

Channel-mediated potassium uptake in Helicobacter pylori is essential for gastric colonization

To date, the biological role of prokaryotic K(+) channels remains unknown. Helicobacter pylori contains a gene encoding a putative K(+) channel (HpKchA) of the two-transmembrane RCK (regulation of K(+) conductance) domain family, but lacks known bacterial K(+) uptake systems. A H. pylori DeltahpKchA mutant presented a strong growth defect at low K(+) concentration, which was compensated by KCl addition. The role of the separate RCK domain was investigated in H. pylori by mutagenesis of its internal start codon, which led to a K(+)-dependent intermediate growth phenotype, consistent with RCK activating channel function. Tagging HpKchA C-terminally, we detected a 1:1 stoichiometry of the full-length HpKchA and the separate RCK domain. We constructed single amino-acid exchanges within the unusual selectivity filter of HpKchA (ATGFGA) in H. pylori and observed complete loss (G74A), a slight defect (G76A or F75G) or wild-type (A77D) channel function. HpKchA was essential for colonization of the murine stomach. These data show, for the first time, a biological function for a prokaryotic K(+) channel, as a K(+) uptake system, essential for the persistence of H. pylori in the gastric environment.

Ion hydration in nanopores and the molecular basis of selectivity

Using a simple model, it is shown that the cost of constraining a hydrated potassium ion inside a narrow pore is smaller than the cost of constraining hydrated sodium or lithium ions in pores of radius around 1.5 A. The opposite is true for pores of radius around 2.5 A. The reason for the selectivity in the first region is that the potassium ion allows for a greater distortion of its hydration shell and can therefore maintain a better coordination, and the reason for the reverse selectivity in the second region is that the smaller ions retain their hydration shells in these pores. This is relevant to the molecular basis of ion selective channels, and since this mechanism does not depend on the molecular details of the pore, it could also operate in all sorts of nanotubes.

Voltage-dependent gating at the KcsA selectivity filter

The prokaryotic K(+) channel KcsA, although lacking a 'standard' voltage-sensing domain, shows voltage-dependent gating that leads to an increase in steady-state open probability of almost two orders of magnitude between +150 and -150 mV. Here we show that voltage-dependent gating in KcsA is associated with the movement of approximately 0.7 equivalent electronic charges. This charge movement produces an increase in the rate of entry into a long-lived inactivated state and seems to be independent of the proton-activation mechanism. Charge neutralization at position 71 renders the channel essentially voltage-independent by preventing entry into the inactivated state. A mechanism for voltage-dependent gating at the selectivity filter is proposed that is based on the reorientation of the carboxylic moiety of Glu71 and its influence in the conformational dynamics of the selectivity filter.

Conformational dynamics of M2 helices in KirBac channels: helix flexibility in relation to gating via molecular dynamics simulations

KirBac1.1 and 3.1 are bacterial homologues of mammalian inward rectifier K channels. We have performed extended molecular dynamics simulations (five simulations, each of >20 ns duration) of the transmembrane domain of KirBac in two membrane environments, a palmitoyl oleoyl phosphatidylcholine bilayer and an octane slab. Analysis of these simulations has focused on the conformational dynamics of the pore-lining M2 helices, which form the cytoplasmic hydrophobic gate of the channel. Principal components analysis reveals bending of M2, with a molecular hinge at the conserved glycine (Gly134 in KirBac1.1, Gly120 in KirBac3.1). More detailed analysis reveals a dimer-of-dimers type motion. The first two eigenvectors describing the motions of M2 correspond to helix kink and swivel motions. The conformational flexibility of M2 seen in these simulations correlates with differences in M2 conformation between that seen in the X-ray structures of closed channels (KcsA and KirBac) in which the helix is undistorted, and in open channels (e.g. MthK) in which the M2 helix is kinked. Thus, the simulations, albeit on a time scale substantially shorter than that required for channel gating, suggest a gating model in which the intrinsic flexibility of M2 about a molecular hinge is coupled to conformational transitions of an intracellular 'gatekeeper' domain, the latter changing conformation in response to ligand binding.

Anionic phospholipid interactions with the potassium channel KcsA: simulation studies

Molecular dynamics (MD) simulations have been used to unmask details of specific interactions of anionic phospholipids with intersubunit binding sites on the surface of the bacterial potassium channel KcsA. Crystallographic data on a diacyl glycerol fragment at this site were used to model phosphatidylethanolamine (PE), or phosphatidylglycerol (PG), or phosphatidic acid (PA) at the intersubunit binding sites. Each of these models of a KcsA-lipid complex was embedded in phosphatidyl choline bilayer and explored in a 20 ns MD simulation. H-bond analysis revealed that in terms of lipid-protein interactions PA > PG >> PE and revealed how anionic lipids (PG and PA) bind to a site provided by two key arginine residues (R(64) and R(89)) at the interface between adjacent subunits. A 27 ns simulation was performed in which KcsA (without any lipids initially modeled at the R(64)/R(89) sites) was embedded in a PE/PG bilayer. There was a progressive specific increase over the course of the simulation in the number of H-bonds of PG with KcsA. Furthermore, two specific PG binding events at R(64)/R(89) sites were observed. The phosphate oxygen atoms of bound PG formed H-bonds to the guanidinium group of R(89), whereas the terminal glycerol H-bonded to R(64). Overall, this study suggests that simulations can help identify and characterize sites for specific lipid interactions on a membrane protein surface.

Through the channel and around the channel: Validating and comparing microscopic approaches for the evaluation of free energy profiles for ion penetration through ion channels

Microscopic calculations of free energy profiles for ion transport through biological ion channels present a very serious challenge to modern simulation approaches. The main problem is due to the major convergence problems associated with the heterogeneous landscape of the electrostatic environment in ion channels and with the need to evaluate the profile associated with the transfer of the ion from bulk water to the channel environment. This problem is compounded by the lack of reliable and relevant benchmarks that can discriminate between alternative approaches. The present study is aimed at reducing the above problems by defining benchmarks that are directly relevant to ion channels and can also give converging results. This is done by constructing a series of models of a truncated gramicidin channel with different numbers of water molecules and by comparing the profiles for going around the channel and through the channel. These discriminating models are then used to validate and compare the adiabatic charging free energy perturbation (FEP) approach combined with an umbrella sampling approach (Warshel, A. J. Phys. Chem. 1982, 86, 2218) and the potential of mean force (PMF) approach used frequently in studies of ion channels. It is found that both approaches work quite well until one moves to the case of the fully solvated channel. In this limit, the PMF approach may give different results for the overall work of going through the channel and around the channel, while the FEP approach gives physically consistent results. The present benchmark also indicates that the weighted histogram analysis method (WHAM) approach does not offer a significant advantage over earlier approaches at least as much as studies of ion channels are concerned. Finally, it is concluded that the FEP approach may be more useful in evaluating the overall barrier for moving ions from water to ion channels and that in some cases it might be beneficial to use the FEP approach for selective points along the channel and then to connect these points by PMF calculations.

Ion conduction and selectivity in K(+) channels

Potassium (K(+)) channels are tetrameric membrane-spanning proteins that provide a selective pore for the conductance of K(+) across the cell membranes. These channels are most remarkable in their ability to discriminate K(+) from Na(+) by more than a thousandfold and conduct at a throughput rate near diffusion limit. The recent progress in the structural characterization of K(+) channel provides us with a unique opportunity to understand their function at the atomic level. With their ability to go beyond static structures, molecular dynamics simulations based on atomic models can play an important role in shaping our view of how ion channels carry out their function. The purpose of this review is to summarize the most important findings from experiments and computations and to highlight a number of fundamental mechanistic questions about ion conduction and selectivity that will require further work.

Ion transit pathways and gating in ClC chloride channels

ClC chloride channels possess a homodimeric structure in which each monomer contains an independent chloride ion pathway. ClC channel gating is regulated by chloride ion concentration, pH and voltage. Based on structural and physiological evidence, it has been proposed that a glutamate residue on the extracellular end of the selectivity filter acts as a fast gate. We utilized a new search algorithm that incorporates electrostatic information to explore the ion transit pathways through wild-type and mutant bacterial ClC channels. Examination of the chloride ion permeation pathways supports the importance of the glutamate residue in gating. An external chloride binding site previously postulated in physiological experiments is located near a conserved basic residue adjacent to the gate. In addition, access pathways are found for proton migration to the gate, enabling pH control at hyperpolarized membrane potentials. A chloride ion in the selectivity filter is required for the pH-dependent gating mechanism.

In silico activation of KcsA K+ channel by lateral forces applied to the C-termini of inner helices

Crystallographic studies of K(+) channels in the closed (KcsA) and open (MthK) states suggest that Gly(99) (KcsA numbering) in the inner helices serves as a gating hinge during channel activation. However, some P-loop channels have larger residues in the corresponding position. The comparison of x-ray structures of KcsA and MthK shows that channel activation alters backbone torsions and helical H-bonds in residues 95-105. Importantly, the changes in Gly(99) are not the largest ones. This raises questions about the mechanism of conformational changes upon channel gating. In this work, we have built a model of the open KcsA using MthK as a template and simulated opening and closing of KcsA by constraining C-ends of the inner helices at a gradually changing distance from the pore axis without restraining mobility of the helices along the axis. At each imposed distance, the energy was Monte Carlo-minimized. The channel-opening and channel-closing trajectories arrived to the structures in which the backbone geometry was close to that seen in MthK and KcsA, respectively. In the channel-opening trajectory, the constraints-induced lateral forces caused kinks at midpoints of the inner helices between Val(97) and Gly(104) but did not destroy interdomain contacts, the pore helices, and the selectivity filter. The simulated activation of the Gly(99)Ala mutant yielded essentially similar results. Analysis of interresidue energies shows that the N-terminal parts of the inner helices form strong attractive contacts with the pore helices and the outer helices. The lateral forces induce kinks at the position where the helix-breaking torque is maximal and the intersegment contacts vanish. This mechanism may be conserved in different P-loop channels.

Ions and blockers in potassium channels: insights from free energy simulations

Potassium ion channels enable efficient and selective permeation of K+ ions across nonpolar biological membranes. Here we review the results of recent free energy calculations related to the permeation of monovalent cations through K+ channels and to the channel inhibition by blocker compounds. In particular, the progress in computational studies of the bacterial KcsA channel is discussed.

Filter flexibility and distortion in a bacterial inward rectifier K+ channel: simulation studies of KirBac1.1

The bacterial channel KirBac1.1 provides a structural homolog of mammalian inward rectifier potassium (Kir) channels. The conformational dynamics of the selectivity filter of Kir channels are of some interest in the context of possible permeation and gating mechanisms for this channel. Molecular dynamics simulations of KirBac have been performed on a 10-ns timescale, i.e., comparable to that of ion permeation. The results of five simulations (total simulation time 50 ns) based on three different initial ion configurations and two different model membranes are reported. These simulation data provide evidence for limited (<0.1 nm) filter flexibility during the concerted motion of ions and water molecules within the filter, such local changes in conformation occurring on an approximately 1-ns timescale. In the absence of K(+) ions, the KirBac selectivity filter undergoes more substantial distortions. These resemble those seen in comparable simulations of other channels (e.g., KcsA and KcsA-based homology models) and are likely to lead to functional closure of the channel. This suggests filter distortions may provide a mechanism of K-channel gating in addition to changes in the hydrophobic gate formed at the intracellular crossing point of the M2 helices. The simulation data also provide evidence for interactions of the "slide" (pre-M1) helix of KirBac with phospholipid headgroups.

Functional influence of the pore helix glutamate in the KcsA K+ channel

The E71 residue is buried near the selectivity filter in the KcsA K+ channel and forms a carboxyl-carboxylate bridge with D80. We have investigated the importance of E71 by examining neutralization mutants at this position using biochemical and electrophysiological methods. E71 mutations differentially destabilize the detergent-solubilized tetramer; among them, the E71V neutralization mutant has a relatively subtle effect. The E71V channel displays electrical activity when reconstituted into planar lipid bilayers. In single channel recordings, the mutant retains K+/Na+ selectivity, and its conductance in the outward direction is unaltered. Some conduction properties are changed: inward conductance is increased. Our results show that that the E71 side chain is not a primary determinant of ion selectivity or conduction in the wild-type channel, either directly or through the E71:D80 carboxyl-carboxylate bridge.

Minimalist molecular model for nanopore selectivity

Using a simple model it is shown that the cost of constraining a hydrated potassium ion inside a narrow nanopore is smaller than the cost of constraining the smaller hydrated sodium ion. The former allows for a greater distortion of its hydration shell and can therefore maintain a better coordination. We propose that in this way the larger ion can go through narrow pores more easily. This is relevant to the molecular basis of ion selective nanopores and since this mechanism does not depend on the molecular details of the pore, it could also operate in all sorts of nanotubes, from biological to synthetic.

Lipid-protein interactions of integral membrane proteins: a comparative simulation study

The interactions between membrane proteins and their lipid bilayer environment play important roles in the stability and function of such proteins. Extended (15-20 ns) molecular dynamics simulations have been used to explore the interactions of two membrane proteins with phosphatidylcholine bilayers. One protein (KcsA) is an alpha-helix bundle and embedded in a palmitoyl oleoyl phosphatidylcholine bilayer; the other (OmpA) is a beta-barrel outer-membrane protein and is in a dimyristoyl phosphatidylcholine bilayer. The simulations enable analysis in detail of a number of aspects of lipid-protein interactions. In particular, the interactions of aromatic amphipathic side chains (i.e., Trp, Tyr) with lipid headgroups, and "snorkeling" interactions of basic side chains (i.e., Lys, Arg) with phosphate groups are explored. Analysis of the number of contacts and of H-bonds reveal fluctuations on an approximately 1- to 5-ns timescale. There are two clear bands of interacting residues on the surface of KcsA, whereas there are three such bands on OmpA. A large number of Arg-phosphate interactions are seen for KcsA; for OmpA, the number of basic-phosphate interactions is smaller and shows more marked fluctuations with respect to time. Both classes of interaction occur in clearly defined interfacial regions of width approximately 1 nm. Analysis of lateral diffusion of lipid molecules reveals that "boundary" lipid molecules diffuse at about half the rate of bulk lipid. Overall, these simulations present a dynamic picture of lipid-protein interactions: there are a number of more specific interactions but even these fluctuate on an approximately 1- to 5-ns timescale.

Electrostatics of Ion Stabilization in a ClC Chloride Channel Homologue from Escherichia coli

The structural determinants of electrostatics of ion stabilization within EcClC, a ClC-type chloride channel homologue from Escherichia coli, are studied using a continuum dielectric approximation. Specifically, the ion occupancy is investigated in the wild-type protein and a mutant thereof, and the contribution to the electrostatic binding free energy of local and non-local interactions is characterized at the single-residue level. This analysis shows that, in spite of the desolvation cost and the strong ion-ion repulsion, all previously reported binding sites can be occupied simultaneously. The stabilizing effect of the protein arises from hydrogen bonding as well as from longer-range favorable interactions, such as with the strictly conserved Lys131 side-chain. The latter is involved in the stabilization of the conserved GSGIP motif that delimits two of the binding sites. Interestingly, an additional low-affinity binding site, mediated by a structurally analogous motif including the side-chain of Arg340, can be identified on the extracellular side of the permeation pathway. Finally, it is shown that, in contrast to K-channels, and in analogy to the SBP/PBP sulfate/phosphate-binding proteins, the contribution of helix macrodipoles to chloride binding in EcClC is only marginal.

Ion channel gating: insights via molecular simulations

Ion channels are gated, i.e. they can switch conformation between a closed and an open state. Molecular dynamics simulations may be used to study the conformational dynamics of ion channels and of simple channel models. Simulations on model nanopores reveal that a narrow (<4 A) hydrophobic region can form a functionally closed gate in the channel and can be opened by either a small (approximately 1 A) increase in pore radius or an increase in polarity. Modelling and simulation studies confirm the importance of hydrophobic gating in K channels, and support a model in which hinge-bending of the pore-lining M2 (or S6 in Kv channels) helices underlies channel gating. Simulations of a simple outer membrane protein, OmpA, indicate that a gate may also be formed by interactions of charged side chains within a pore, as is also the case in ClC channels.

Potassium channel, ions, and water: simulation studies based on the high resolution X-ray structure of KcsA

Interactions of Na(+), K(+), Rb(+), and Cs(+) ions within the selectivity filter of a potassium channel have been investigated via multiple molecular dynamics simulations (total simulation time, 48 ns) based on the high resolution structure of KcsA, embedded in a phospholipid bilayer. As in simulations based on a lower resolution structure of KcsA, concerted motions of ions and water within the filter are seen. Despite the use of a higher resolution structure and the inclusion of four buried water molecules thought to stabilize the filter, this region exhibits a significant degree of flexibility. In particular, pronounced distortion of filter occurs if no ions are present within it. The two most readily permeant ions, K(+) and Rb(+), are similar in their interactions with the selectivity filter. In contrast, Na(+) ions tend to distort the filter by binding to a ring of four carbonyl oxygens. The larger Cs(+) ions result in a small degree of expansion of the filter relative to the x-ray structure. Cs(+) ions also appear to interact differently with the gate region of the channel, showing some tendency to bind within a predominantly hydrophobic pocket. The four water molecules buried between the back of the selectivity filter and the remainder of the protein show comparable mobility to the surrounding protein and do not exchange with water molecules within the filter or the central cavity. A preliminary comparison of the use of particle mesh Ewald versus cutoff protocols for the treatment of long-range electrostatics suggests some difference in the kinetics of ion translocation within the filter.

KcsA closed and open: modelling and simulation studies

Bacterial homologues of mammalian potassium channels provide structures of two states of a gated K channel. Thus, the crystal structure of KcsA represents a closed state whilst that of MthK represents an open state. Using homology modelling and molecular dynamics simulations we have built a model of the transmembrane domain of KcsA in an open state and have compared its conformational stability with that of the same domain of KcsA in a closed state. Approximate Born energy calculations of monovalent cations within the two KcsA channel states suggest that the intracellular hydrophobic gate in the closed state provides a barrier of height approximately 5 kT to ion permeation, whilst in the open state the barrier is absent. Simulations (10 ns duration) in an octane slab (a simple membrane mimetic) suggest that closed- and open-state models are of comparable conformational stability, both exhibiting conformational drifts of approximately 3.3 A Calpha RMSD relative to the respective starting models. Substantial conformational fluctuations are observed in the intracellular gate region during both simulations (closed state and open state). In the simulation of open-state KcsA, rapid ( < 5 ns) exit of all three K+ ions occurs through the intracellular mouth of the channel. Helix kink and swivel motion is observed at the molecular hinge formed by residue G99 of the M2 helix. This motion is more substantial for the open- than for the closed-state model of the channel.

Dielectric self-energy in Poisson-Boltzmann and Poisson-Nernst-Planck models of ion channels

We demonstrated previously that the two continuum theories widely used in modeling biological ion channels give unreliable results when the radius of the conduit is less than two Debye lengths. The reason for this failure is the neglect of surface charges on the protein wall induced by permeating ions. Here we attempt to improve the accuracy of the Poisson-Boltzmann and Poisson-Nernst-Planck theories, when applied to channel-like environments, by including a specific dielectric self-energy term to overcome spurious shielding effects inherent in these theories. By comparing results with Brownian dynamics simulations, we show that the inclusion of an additional term in the equations yields significant qualitative improvements. The modified theories perform well in very wide and very narrow channels, but are less successful at intermediate sizes. The situation is worse in multi-ion channels because of the inability of the continuum theories to handle the ion-to-ion interactions correctly. Thus, further work is required if these continuum theories are to be reliably salvaged for quantitative studies of biological ion channels in all situations.

Modeling permeation energetics in the KcsA potassium channel

The thermodynamics of cation permeation through the KcsA K(+) channel selectivity filter is studied from the perspective of a physically transparent semimicroscopic model using Monte Carlo free energy integration. The computational approach chosen permits dissection of the separate contributions to ionic stabilization arising from different parts of the channel (selectivity filter carbonyls, single-file water, cavity water, reaction field of bulk water, inner helices, ionizable residues). All features play important roles; their relative significance varies with the ion's position in the filter. The cavity appears to act as an electrostatic buffer, shielding filter ions from structural changes in the inner pore. The model exhibits K(+) vs. Na(+) selectivity, and roughly isoenergetic profiles for K(+) and Rb(+), and discriminates against Cs(+), all in agreement with experimental data. It also indicates that Ba(2+) and Na(+) compete effectively with permeant ions at a site near the boundary between the filter and the cavity, in the vicinity of the barium blocker site.

Recent advances in ion channel research

The field of ion channels has entered into a rapid phase of development in the last few years, partly due to the breakthroughs in determination of the crystal structures of membrane proteins and advances in computer simulations of biomolecules. These advances have finally enabled the long-dreamed goal of relating function of a channel to its underlying molecular structure. Here we present simplified accounts of the competing permeation theories and then discuss their application to the potassium, gramicidin A and calcium channels.

K(+) versus Na(+) ions in a K channel selectivity filter: a simulation study

Molecular dynamics simulations of a bacterial potassium channel (KcsA) embedded in a phospholipid bilayer reveal significant differences in interactions of the selectivity filter with K(+) compared with Na(+) ions. K(+) ions and water molecules within the filter undergo concerted single-file motion in which they translocate between adjacent sites within the filter on a nanosecond timescale. In contrast, Na(+) ions remain bound to sites within the filter and do not exhibit translocation on a nanosecond timescale. Furthermore, entry of a K(+) ion into the filter from the extracellular mouth is observed, whereas this does not occur for a Na(+) ion. Whereas K(+) ions prefer to sit within a cage of eight oxygen atoms of the filter, Na(+) ions prefer to interact with a ring of four oxygen atoms plus two water molecules. These differences in interactions in the selectivity filter may contribute to the selectivity of KcsA for K(+) ions (in addition to the differences in dehydration energy between K(+) and Na(+)) and the block of KcsA by internal Na(+) ions. In our simulations the selectivity filter exhibits significant flexibility in response to changes in ion/protein interactions, with a somewhat greater distortion induced by Na(+) than by K(+) ions.

Modeling diverse range of potassium channels with Brownian dynamics

Using the experimentally determined KcsA structure as a template, we propose a plausible explanation for the diversity of potassium channels seen in nature. A simplified model of KcsA is constructed from its atomic resolution structure by smoothing out the protein-water boundary and representing the atoms forming the channel protein as a homogeneous, low dielectric medium. The properties of the simplified and atomic-detail models, deduced from electrostatic calculations and Brownian dynamics simulations, are shown to be qualitatively similar. We then study how the current flowing across the simplified model channel changes as the shape of the intrapore region is modified. This is achieved by increasing the radius of the intracellular pore systematically from 1.5 to 5 A while leaving the dimensions of the selectivity filter and inner chamber unaltered. The strengths of the dipoles located near the entrances of the channel, the carbonyl groups lining the selectivity filter, and the helix macrodipoles are kept constant. The channel conductance increases steadily as the radius of the intracellular pore is increased. The rate-limiting step for both the outward and inward current is the time it takes for an ion to cross the residual energy barrier located in the intrapore region. The current-voltage relationship obtained with symmetrical solutions is linear when the applied potential is less than approximately 100 mV but deviates slightly from Ohm's law at higher applied potentials. The nonlinearity in the current-voltage curve becomes less pronounced as the radius of the intracellular pore is increased. When the strengths of the dipoles near the intracellular entrance are reduced, the channel shows a pronounced inward rectification. Finally, the conductance exhibits the saturation property observed experimentally. We discuss the implications of these findings on the transport of ions across the potassium channels and membrane channels in general.

Molecular dynamics simulations of biological reactions

This review considers the author perspective on the emergence of molecular dynamics (MD) simulations of biological processes. It starts with the 1976 simulation of the primary event in rhodopsin, moves to the earliest simulations of enzymatic reactions and electron transfer reactions and ends up with recent simulations of proton translocations and ion transport in proteins. The emphasis is placed on our progress in simulations of actual biological reactions and functional properties, rather than on studies of general properties such as structure and thermal motions. In most cases it has been possible to develop special strategies that capture the relevant dynamics of the given biological process. The predictive power of our early simulations of fast biological process (e.g. vision and photosynthesis) and the insight obtained from these studies is pointed out. Critical examinations of dynamical effects in different biological processes is reviewed. This includes the finding that dynamical effects are unlikely to contribute significantly to enzyme catalysis or to other processes with significant activation barriers. Even in the case of ion channels it is found that the most important effects are associated with energetics rather than dynamics. Nevertheless, MD simulations provide what is probably the most realistic description of the actual reactive events. The resulting insight is crucial in studies of fast photobiological reactions and instructive in cases of slower processes.

Simulations of ion current in realistic models of ion channels: the KcsA potassium channel

Realistic studies of ion current in biologic channels present a major challenge for computer simulation approaches. All-atom molecular dynamics simulations involve serious time limitations that prevent their use in direct evaluation of ion current in channels with significant barriers. The alternative use of Brownian dynamics (BD) simulations can provide the current for simplified macroscopic models. However, the time needed for accurate calculations of electrostatic energies can make BD simulations of ion current expensive. The present work develops an approach that overcomes some of the above challenges and allows one to simulate ion currents in models of biologic channels. Our method provides a fast and reliable estimate of the energetics of the system by combining semimacroscopic calculations of the self-energy of each ion and an implicit treatment of the interactions between the ions, as well as the interactions between the ions and the protein-ionizable groups. This treatment involves the use of the semimacroscopic version of the protein dipole Langevin dipole (PDLD/S) model in its linear response approximation (LRA) implementation, which reduces the uncertainties about the value of the protein "dielectric constant." The resulting free energy surface is used to generate the forces for on-the-fly BD simulations of the corresponding ion currents. Our model is examined in a preliminary simulation of the ion current in the KcsA potassium channel. The complete free energy profile for a single ion transport reflects reasonable energetics and captures the effect of the protein-ionized groups. This calculated profile indicates that we are dealing with the channel in its closed state. Reducing the barrier at the gate region allows us to simulate the ion current in a reasonable computational time. Several limiting cases are examined, including those that reproduce the observed current, and the nature of the productive trajectories is considered. The ability to simulate the current in realistic models of ion channels should provide a powerful tool for studies of the biologic function of such systems, including the analysis of the effect of mutations, pH, and electric potentials.