Molecular dynamics of the KcsA K(+) channel in a bilayer membrane

Probing the energy landscape of activation gating of the bacterial potassium channel KcsA

The bacterial potassium channel KcsA, which has been crystallized in several conformations, offers an ideal model to investigate activation gating of ion channels. In this study, essential dynamics simulations are applied to obtain insights into the transition pathways and the energy profile of KcsA pore gating. In agreement with previous hypotheses, our simulations reveal a two phasic activation gating process. In the first phase, local structural rearrangements in TM2 are observed leading to an intermediate channel conformation, followed by large structural rearrangements leading to full opening of KcsA. Conformational changes of a highly conserved phenylalanine, F114, at the bundle crossing region are crucial for the transition from a closed to an intermediate state. 3.9 µs umbrella sampling calculations reveal that there are two well-defined energy barriers dividing closed, intermediate, and open channel states. In agreement with mutational studies, the closed state was found to be energetically more favorable compared to the open state. Further, the simulations provide new insights into the dynamical coupling effects of F103 between the activation gate and the selectivity filter. Investigations on individual subunits support cooperativity of subunits during activation gating.

Physics of ion channels

We review the basic physics involved in transport of ions across membrane channels in cells. Electrochemical forces that control the diffusion of ions are discussed both from microscopic and macroscopic perspectives. A case is made for use of Brownian dynamics as the minimal phenomenological model that provides a bridge between experiments and more fundamental theoretical approaches. Application of Brownian and molecular dynamics methods to channels with known molecular structures is discussed.

The open gate structure of the membrane-embedded KcsA potassium channel viewed from the cytoplasmic side

Crystallographic studies of channel proteins have provided insight into the molecular mechanisms of ion channels, even though these structures are obtained in the absence of the membrane and some structural portions have remained unsolved. Here we report the gating structure of the membrane-embedded KcsA potassium channel using atomic force microscopy (AFM). The activation gate of the KcsA channel is located on the intracellular side, and the cytoplasmic domain was truncated to clear the view of this location. Once opened, the individual subunits in the tetramer were resolved with the pore open at the center. Furthermore, AFM was able to capture the previously unsolved bulge helix at the entrance. A molecular dynamics simulation revealed that the bulge helices fluctuated dramatically at the open entryway. This dynamic behavior was observed as vigorous open-channel noise in the single-channel current recordings. The role of the bulge helices in the open gate structure is discussed.

High-risk long QT syndrome mutations in the Kv7.1 (KCNQ1) pore disrupt the molecular basis for rapid K(+) permeation

Type 1 long QT syndrome (LQT1) is caused by loss-of-function mutations in the KCNQ1 gene, which encodes the K(+) channel (Kv7.1) that underlies the slowly activating delayed rectifier K(+) current in the heart. Intragenic risk stratification suggests LQT1 mutations that disrupt conserved amino acid residues in the pore are an independent risk factor for LQT1-related cardiac events. The purpose of this study is to determine possible molecular mechanisms that underlie the loss of function for these high-risk mutations. Extensive genotype-phenotype analyses of LQT1 patients showed that T322M-, T322A-, or G325R-Kv7.1 confers a high risk for LQT1-related cardiac events. Heterologous expression of these mutations with KCNE1 revealed they generated nonfunctional channels and caused dominant negative suppression of WT-Kv7.1 current. Molecular dynamics simulations of analogous mutations in KcsA (T85M-, T85A-, and G88R-KcsA) demonstrated that they disrupted the symmetrical distribution of the carbonyl oxygen atoms in the selectivity filter, which upset the balance between the strong attractive and K(+)-K(+) repulsive forces required for rapid K(+) permeation. We conclude high-risk LQT1 mutations in the pore likely disrupt the architectural and physical properties of the K(+) channel selectivity filter.

Mechanism of Cd2+ coordination during slow inactivation in potassium channels

In K+ channels, rearrangements of the pore outer vestibule have been associated with C-type inactivation gating. Paradoxically, the crystal structure of Open/C-type inactivated KcsA suggests these movements to be modest in magnitude. In this study, we show that under physiological conditions, the KcsA outer vestibule undergoes relatively large dynamic rearrangements upon inactivation. External Cd2+ enhances the rate of C-type inactivation in an cysteine mutant (Y82C) via metal-bridge formation. This effect is not present in a non-inactivating mutant (E71A/Y82C). Tandem dimer and tandem tetramer constructs of equivalent cysteine mutants in KcsA and Shaker K+ channels demonstrate that these Cd2+ metal bridges are formed only between adjacent subunits. This is well supported by molecular dynamics simulations. Based on the crystal structure of Cd2+ -bound Y82C-KcsA in the closed state, together with electron paramagnetic resonance distance measurements in the KcsA outer vestibule, we suggest that subunits must dynamically come in close proximity as the channels undergo inactivation.

Thermodynamic coupling between activation and inactivation gating in potassium channels revealed by free energy molecular dynamics simulations

The amount of ionic current flowing through K⁺ channels is determined by the interplay between two separate time-dependent processes: activation and inactivation gating. Activation is concerned with the stimulus-dependent opening of the main intracellular gate, whereas inactivation is a spontaneous conformational transition of the selectivity filter toward a nonconductive state occurring on a variety of timescales. A recent analysis of multiple x-ray structures of open and partially open KcsA channels revealed the mechanism by which movements of the inner activation gate, formed by the inner helices from the four subunits of the pore domain, bias the conformational changes at the selectivity filter toward a nonconductive inactivated state. This analysis highlighted the important role of Phe103, a residue located along the inner helix, near the hinge position associated with the opening of the intracellular gate. In the present study, we use free energy perturbation molecular dynamics simulations (FEP/MD) to quantitatively elucidate the thermodynamic basis for the coupling between the intracellular gate and the selectivity filter. The results of the FEP/MD calculations are in good agreement with experiments, and further analysis of the repulsive, van der Waals dispersive, and electrostatic free energy contributions reveals that the energetic basis underlying the absence of inactivation in the F103A mutation in KcsA is the absence of the unfavorable steric interaction occurring with the large Ile100 side chain in a neighboring subunit when the intracellular gate is open and the selectivity filter is in a conductive conformation. Macroscopic current analysis shows that the I100A mutant indeed relieves inactivation in KcsA, but to a lesser extent than the F103A mutant.

Ion solvation and structural stability in a sodium channel investigated by molecular dynamics calculations

The stability and ion binding properties of the homo-tetrameric pore domain of a prokaryotic, voltage-gated sodium channel are studied by extensive all-atom molecular dynamics simulations, with the channel protein being embedded in a fully hydrated lipid bilayer. It is found that Na(+) ion presents in a mostly hydrated state inside the wide pore of the selectivity filter of the sodium channel, in sharp contrast to the nearly fully dehydrated state for K(+) ions in potassium channels. Our results also indicate that Na(+) ions make contact with only one or two out of the four polypeptide chains forming the selectivity filter, and surprisingly, the selectivity filter exhibits robust stability for various initial ion configurations even in the absence of ions. These findings are quite different from those in potassium channels. Furthermore, an electric field above 0.5V/nm is suggested to be able to induce Na(+) permeation through the selectivity filter.

Hinge-bending motions in the pore domain of a bacterial voltage-gated sodium channel

Computational methods and experimental data are used to provide structural models for NaChBac, the homo-tetrameric voltage-gated sodium channel from the bacterium Bacillus halodurans, with a closed and partially open pore domain. Molecular dynamic (MD) simulations on membrane-bound homo-tetrameric NaChBac structures, each comprising six helical transmembrane segments (labeled S1 through S6), reveal that the shape of the lumen, which is defined by the bundle of four alpha-helical S6 segments, is modulated by hinge bending motions around the S6 glycine residues. Mutation of these glycine residues into proline and alanine affects, respectively, the structure and conformational flexibility of the S6 bundle. In the closed channel conformation, a cluster of stacked phenylalanine residues from the four S6 helices hinders diffusion of water molecules and Na(+) ions. Activation of the voltage sensor domains causes destabilization of the aforementioned cluster of phenylalanines, leading to a more open structure. The conformational change involving the phenylalanine cluster promotes a kink in S6, suggesting that channel gating likely results from the combined action of hinge-bending motions of the S6 bundle and concerted reorientation of the aromatic phenylalanine side-chains.

Charybdotoxin and margatoxin acting on the human voltage-gated potassium channel hKv1.3 and its H399N mutant: an experimental and computational comparison

The effect of the pore-blocking peptides charybdotoxin and margatoxin, both scorpion toxins, on currents through human voltage-gated hK(v)1.3 wild-type and hK(v)1.3_H399N mutant potassium channels was characterized by the whole-cell patch clamp technique. In the mutant channels, both toxins hardly blocked current through the channels, although they did prevent C-type inactivation by slowing down the current decay during depolarization. Molecular dynamics simulations suggested that the fast current decay in the mutant channel was a consequence of amino acid reorientations behind the selectivity filter and indicated that the rigidity-flexibility in that region played a key role in its interactions with scorpion toxins. A channel with a slightly more flexible selectivity filter region exhibits distinct interactions with scorpion toxins. Our studies suggest that the toxin-channel interactions might partially restore rigidity in the selectivity filter and thereby prevent the structural rearrangements associated with C-type inactivation.

Computer Simulations of Voltage-Gated Cation Channels

The relentless growth in computational power has seen increasing applications of molecular dynamics (MD) simulation to the study of membrane proteins in realistic membrane environments, which include explicit membrane lipids, water and ions. The concomitant increasing availability of membrane protein structures for ion channels, and transporters -- to name just two examples -- has stimulated many of these MD studies. In the case of voltage-gated cation channels (VGCCs) recent computational works have focused on ion-conduction and gating mechanisms, along with their regulation by agonist/antagonist ligands. The information garnered from these computational studies is largely inaccessible to experiment and is crucial for understanding the interplay between the structure and function as well as providing new directions for experiments. This article highlights recent advances in probing the structure and function of potassium channels and offers a perspective on the challenges likely to arise in making analogous progress in characterizing sodium channels.

Kinetics of the terminal electron transfer step in cytochrome c oxidase

Cytochrome c oxidase (cco) catalyzes the oxygen reduction reaction in most aerobically respiring organisms. Decades of research have uncovered many aspects relating to structure and function of this enzyme. However, the origin of the unusually fast terminal electron transfer step from heme a to heme a(3) in cco has been the subject of intense discussions over recent years. Yet, no satisfactory consensus has been achieved. Carrying out large-scale molecular dynamics simulation of the protein embedded in a solvated membrane, we obtain a reorganization free energy λ = 0.57 eV. Evaluation of the quantized single-mode rate equation using the experimental rate and the computed reorganization free energy gives a value of 1.5 meV for the average electronic coupling (H(ab)) between heme a and heme a(3). Thus, according to our calculations, the nanosecond electron transfer (ET) is due to a small but significant activation barrier (ΔG(‡) = 0.12 eV) in combination with effective electronic coupling between the two cofactors. The activation free energy is caused predominantly by collective reorganization of protein residues. We show that our results are consistent with the weak temperature dependence observed in experiment if one allows for very minor variations in the donor-acceptor distance as the temperature changes.

Oligomerization at the membrane: potassium channel structure and function

Cell membranes present a naturally impervious barrier to aqueous solutes, such that the physiochemical environment on either side of the lipid bilayer can substantially differ. Integral membrane proteins are embedded in this heterogeneous lipid environment, wherein the juxtaposition of apolar and polar molecular surfaces defines factors such as transverse orientation, the surface area available for oligomerisation and the symmetry of resultant assemblies. This chapter focuses on potassium channels -representative molecular pores that play a critical role in electrical signalling by enabling selective transport of K(+) ions across cell membranes. Oligomerization is central to K(+) channel action; individual subunits are nonfunctional and conduction, selectivity and gating involve manipulation of the common subunit interface of the tetramer. Regulation of channel activity can be viewed from the perspective that the pore of K(+) channels has coopted other proteins, utilizing a process of hetero-oligomerisation to absorb new functions that both enable the pore to respond to extrinsic signals and provide an electrical signature.

Gating at the selectivity filter of ion channels that conduct Na+ and K+ ions

The NaK channel is a cation selective channel with similar permeability for K(+) and Na(+). The available crystallographic structure of wild-type (WT) NaK is usually associated with a conductive state of the channel. Here, potential of mean force for complete conduction events of Na(+) and K(+) ions through NaK show that: i), large energy barriers prevent the passage of ions through the WT NaK structure, ii), the barriers are correlated to the presence of a hydrogen bond between Asp-66 and Asn-68, and iii), the structure of NaK mutated to mimic cyclic nucleotide-gated channels conducts Na(+) and K(+). These results support the hypothesis that the filter of cation selective channels can adopt at least two different structures: a conductive one, represented by the x-ray structures of the NaK-CNG chimeras, and a closed one, represented by the x-ray structures of the WT NaK.

Structural correlates of selectivity and inactivation in potassium channels

Potassium channels are involved in a tremendously diverse range of physiological applications requiring distinctly different functional properties. Not surprisingly, the amino acid sequences for these proteins are diverse as well, except for the region that has been ordained the "selectivity filter". The goal of this review is to examine our current understanding of the role of the selectivity filter and regions adjacent to it in specifying selectivity as well as its role in gating/inactivation and possible mechanisms by which these processes are coupled. Our working hypothesis is that an amino acid network behind the filter modulates selectivity in channels with the same signature sequence while at the same time affecting channel inactivation properties. This article is part of a Special Issue entitled: Membrane protein structure and function.

Charybdotoxin unbinding from the mKv1.3 potassium channel: a combined computational and experimental study

Charybdotoxin, belonging to the group of so-called scorpion toxins, is a short peptide able to block many voltage-gated potassium channels, such as mKv1.3, with high affinity. We use a reliable homology model based on the high-resolution crystal structure of the 94% sequence identical homologue Kv1.2 for charybdotoxin docking followed by molecular dynamics simulations to investigate the mechanism and energetics of unbinding, tracing the behavior of the channel protein and charybdotoxin during umbrella-sampling simulations as charybdotoxin is moved away from the binding site. The potential of mean force is constructed from the umbrella sampling simulations and combined with K(d) and free energy values gained experimentally using the patch-clamp technique to study the free energy of binding at different ion concentrations and the mechanism of the charybdotoxin-mKv1.3 binding process. A possible charybdotoxin binding mechanism is deduced that includes an initial hydrophobic contact followed by stepwise electrostatic interactions and finally optimization of hydrogen bonds and salt bridges.

A multipoint hydrogen-bond network underlying KcsA C-type inactivation

In the prokaryotic potassium channel KcsA activation gating at the inner bundle gate is followed by C-type inactivation at the selectivity filter. Entry into the C-type inactivated state has been directly linked to the strength of the H-bond interaction between residues Glu-71 and Asp-80 behind the filter, and is allosterically triggered by the rearrangement of the inner bundle gate. Here, we show that H-bond pairing between residues Trp-67 and Asp-80, conserved in most K⁺ channels, constitutes another critical interaction that determines the rate and extent of KcsA C-type inactivation. Disruption of the equivalent interaction in Shaker (Trp-434-Asp-447) and Kv1.2 (Trp-366-Asp-379) leads also to modulation of the inactivation process, suggesting that these residues also play an analogous role in the inactivation gating of Kv channels. The present results show that in KcsA C-type inactivation gating is governed by a multipoint hydrogen-bond network formed by the triad Trp-67-Glu71-Asp-80. This triad exerts a critical role in the dynamics and conformational stability of the selectivity filter and might serve as a general modulator of selectivity filter gating in other members of the K⁺ channel family.

A permeation theory for single-file ion channels: one- and two-step models

How many steps are required to model permeation through ion channels? This question is investigated by comparing one- and two-step models of permeation with experiment and MD simulation for the first time. In recent MD simulations, the observed permeation mechanism was identified as resembling a Hodgkin and Keynes knock-on mechanism with one voltage-dependent rate-determining step [Jensen et al., PNAS 107, 5833 (2010)]. These previously published simulation data are fitted to a one-step knock-on model that successfully explains the highly non-Ohmic current-voltage curve observed in the simulation. However, these predictions (and the simulations upon which they are based) are not representative of real channel behavior, which is typically Ohmic at low voltages. A two-step association/dissociation (A/D) model is then compared with experiment for the first time. This two-parameter model is shown to be remarkably consistent with previously published permeation experiments through the MaxiK potassium channel over a wide range of concentrations and positive voltages. The A/D model also provides a first-order explanation of permeation through the Shaker potassium channel, but it does not explain the asymmetry observed experimentally. To address this, a new asymmetric variant of the A/D model is developed using the present theoretical framework. It includes a third parameter that represents the value of the "permeation coordinate" (fractional electric potential energy) corresponding to the triply occupied state n of the channel. This asymmetric A/D model is fitted to published permeation data through the Shaker potassium channel at physiological concentrations, and it successfully predicts qualitative changes in the negative current-voltage data (including a transition to super-Ohmic behavior) based solely on a fit to positive-voltage data (that appear linear). The A/D model appears to be qualitatively consistent with a large group of published MD simulations, but no quantitative comparison has yet been made. The A/D model makes a network of predictions for how the elementary steps and the channel occupancy vary with both concentration and voltage. In addition, the proposed theoretical framework suggests a new way of plotting the energetics of the simulated system using a one-dimensional permeation coordinate that uses electric potential energy as a metric for the net fractional progress through the permeation mechanism. This approach has the potential to provide a quantitative connection between atomistic simulations and permeation experiments for the first time.

Coarse grained model for exploring voltage dependent ion channels

The relationship between the membrane voltage and the gating of voltage activated ion channels and other systems have been a problem of great current interest. Unfortunately, reliable molecular simulations of external voltage effects present a major challenge, since meaningful converging microscopic simulations are not yet available and macroscopic treatments involve major uncertainties in terms of the dielectric used and other key features. This work extends our coarse grained (CG) model to simulations of membrane/protein systems under external potential. Special attention is devoted to a consistent modeling of the effect of external potential due to the electrodes, emphasizing semimacroscopic description of the electrolytes in the solution regions between the membranes and the electrodes, as well as the coupling between the combined potential from the electrodes plus the electrolytes and the protein ionized groups. We also provide a clear connection to microscopic treatment of the electrolytes and thus can explore possible conceptual problems that are hard to resolve by other current approaches. For example, we obtain a clear description of the charge distribution in the entire electrolyte system, including near the electrodes in membrane/electrodes systems (where continuum models do not seem to provide the relevant results). Furthermore, the present treatment provides an insight on the distribution of the electrolyte charges before and after equilibration across the membrane, and thus on the nature of the gating charge. The different aspects of the model have been carefully validated by considering problems ranging for the simple Debye-Huckel, and the Gouy-Chapman models to the evaluation of the electrolyte distribution between two electrodes, as well as the effect of extending the simulation system by periodic replicas. Overall the clear connection to microscopic descriptions combined with the power of the CG modeling seems to offer a powerful tool for exploring the balance between the protein conformational energy and the interaction with the external potential in voltage activated channels. To illustrate these features we present a preliminary study of the gating charge in the voltage activated Kv1.2 channel, using the actual change in the electrolyte charge distribution rather than the conventional macroscopic estimate. We also discuss other special features of the model, which include the ability to capture the effect of changes in the protonation states of the protein residues during the close to open voltage induced transition. This article is part of a Special Issue entitled: Membrane protein structure and function.

K2P potassium channels, mysterious and paradoxically exciting

New evidence reveals that the common electrolyte disorder hypokalemia can induce K2P1 channels that are normally selective for K+ to break the rules and conduct Na+. This defiant behavior leads to paradoxical depolarization of many cells in the heart, increasing the risk for lethal arrhythmia. The new research resolves a mystery uncovered 50 years ago and bestows an array of new riddles. Here, I discuss how K2P1 might achieve this alchemy--through stable residence of the K+ selectivity filter in a Na+-conductive state between its open and C-inactive configurations--and predict that other K+ channels and environmental stimuli will be discovered to produce the same excitatory misconduct.

Probing the thermodynamics of competitive ion binding using minimum energy structures

Ion binding is known to affect the properties of biomolecules and is directly involved in many biochemical pathways. Because of the highly polar environments where ions are found, a quantum-mechanical treatment is preferable for understanding the energetics of competitive ion binding. Due to computational cost, a quantum mechanical treatment may involve several approximations, however, whose validity can be difficult to determine. Using thermodynamic cycles, we show how intuitive models for complicated ion binding reactions can be built up from simplified, isolated ion-ligand binding site geometries suitable for quantum mechanical treatment. First, the ion binding free energies of individual, minimum energy structures determine their intrinsic ion selectivities. Next, the relative propensity for each minimum energy structure is determined locally from the balance of ion-ligand and ligand-ligand interaction energies. Finally, the environment external to the binding site exerts its influence both through long-ranged dispersive and electrostatic interactions with the binding site as well as indirectly through shifting the binding site compositional and structural preferences. The resulting picture unifies field-strength, topological control, and phase activation viewpoints into a single theory that explicitly indicates the important role of solute coordination state on overall reaction energetics. As an example, we show that the Na(+) → K(+) selectivities can be recovered by correctly considering the conformational contribution to the selectivity. This can be done even when constraining configuration space to the neighborhood around a single, arbitrarily chosen, minimum energy structure. Structural regions around minima for K(+)- and Na(+)-water clusters are exhibited that display both rigid/mechanical and disordered/entropic selectivity mechanisms for both Na(+) and K(+). Thermodynamic consequences of the theory are discussed with an emphasis on the role of coordination structure in determining experimental properties of ions in complex biological environments.

The role of bulk protein in local models of ion-binding to proteins: comparative study of KcsA, its semisynthetic analog with a locked-in binding site, and valinomycin

In studying ion-selectivity in biomaterials, it is common to study ion-protein interactions within a local neighborhood around the ion. This local system analysis for the S(2) site of KcsA, its semisynthetic analog, and valinomycin yields the free energy change in exchanging K(+) with Na(+) in quantitative agreement with the value obtained by considering ion-interactions with the entire system. But the energetics of ion binding in the local system and in the entire system differ significantly and lead to different conclusions regarding the physical basis of ion selectivity. For configurations sampled from an all-atom simulation, we show that the selectivity free energy can be decomposed into a contribution arising from interactions of the ion with its local neighborhood, ΔW(local), and a term arising from the field imposed on the ion and the binding site by the rest of the medium, ΔW(ϕ). The local contribution ΔW(local) is numerically close to the actual free energy difference because the field contribution is small. The field contribution is small because of cancellation of inversely related ion-medium and site-medium interactions. Our analysis presents a rigorous foundation for the numerical success of the local system analysis and shows that its implications do not always hold for the entire protein.

Dynamics may significantly influence the estimation of interatomic distances in biomolecular X-ray structures

Atomic positions obtained by X-ray crystallography are time and space averages over many molecules in the crystal. Importantly, interatomic distances, calculated between such average positions and frequently used in structural and mechanistic analyses, can be substantially different from the more appropriate time-average and ensemble-average interatomic distances. Using crystallographic B-factors, one can deduce corrections, which have so far been applied exclusively to small molecules, to obtain correct average distances as a function of the type of atomic motion. Here, using 4774 high-quality protein X-ray structures, we study the significance of such corrections for different types of atomic motion. Importantly, we show that for distances shorter than 5 Å, corrections greater than 0.5 Å may apply, especially for noncorrelated or anticorrelated motion. For example, 14% of the studied structures have at least one pair of atoms with a correction of ≥ 0.5 Å in the case of noncorrelated motion. Using molecular dynamics simulations of villin headpiece, ubiquitin, and SH3 domain unit cells, we demonstrate that the majority of average interatomic distances in these proteins agree with noncorrelated corrections, suggesting that such deviations may be truly relevant. Importantly, we demonstrate that the corrections do not significantly affect stereochemistry and the overall quality of final refined X-ray structures, but can provide marked improvements in starting unrefined models obtained from low-resolution X-ray data. Finally, we illustrate the potential mechanistic and biological significance of the calculated corrections for KcsA ion channel and show that they provide indirect evidence that motions in its selectivity filter are highly correlated.

Energetics of double-ion occupancy in the gramicidin A channel

To understand the energetics of double-ion occupancy in gramicidin A (gA) channels, the 2D potential of mean force (PMF) is calculated for two ions at different positions along the channel axis. The cross sections of this 2D PMF are compared with available one-ion PMFs to highlight the effect of one ion on the permeation dynamics of the other. It is found that, if the first ion stays on one side in the channel, the second ion has to pass over an additional barrier to move into the outer binding site. At the same time, both outer and inner binding sites for the second ion become shallower than those in the one-ion PMF. The calculated ion-ion repulsion for a doubly occupied channel is about 2 kcal/mol, in good agreement with previous experimental estimates. The number of water molecules inside the channel and their dipole moment are calculated to interpret the energetics of double-ion occupancy. As the first ion moves into the outer binding site and then further into the channel, the oxygen atoms of the single-file water column in the channel are oriented to point toward the ion. The observed dipole moment distribution of a singly occupied channel has only one sharp peak, and the water alignment is essentially perfect once the ion is in the inner binding site. For this reason, there is an energy penalty to accommodate a second ion at the opposite end of the channel.

On the structural basis of modal gating behavior in K(+) channels

Modal-gating shifts represent an effective regulatory mechanism by which ion channels control the extent and time course of ionic fluxes. Under steady-state conditions, the K(+) channel KcsA shows three distinct gating modes, high-P(o), low-P(o) and a high-frequency flicker mode, each with about an order of magnitude difference in their mean open times. Here we show that in the absence of C-type inactivation, mutations at the pore-helix position Glu71 unmask a series of kinetically distinct modes of gating in a side chain-specific way. These gating modes mirror those seen in wild-type channels and suggest that specific interactions in the side chain network surrounding the selectivity filter, in concert with ion occupancy, alter the relative stability of pre-existing conformational states of the pore. The present results highlight the key role of the selectivity filter in regulating modal gating behavior in K(+) channels.

Insights on the permeability of wide protein channels: measurement and interpretation of ion selectivity

Ion channels are hollow proteins that have evolved to exhibit discrimination between charged solutes. This property, known as ion selectivity is critical for several biological functions. By using the bacterial porin OmpF as a model system of wide protein channels, we demonstrate that significant insights can be gained when selectivity measurements are combined with electrodiffusion continuum models and simulations based on the atomic structure. A correct interpretation of the mechanisms ruling the many sources of channel discrimination is a first, indispensable step for the understanding of the controlled movement of ions into or out of cells characteristic of many physiological processes. We conclude that the scattered information gathered from several independent approaches should be appropriately merged to provide a unified and coherent picture of the channel selectivity.

Approximation of super-ions for single-file diffusion of multiple ions through narrow pores

The general theory of the single-file multiparticle diffusion in the narrow pores could be greatly simplified in the case of inverted bell-like shape of the single-particle energy profile, which is often observed in biological ion channels. There is a narrow and deep groove in the energy landscape of multiple interacting ions in such profiles, which corresponds to the pre-defined optimal conduction pathway in the configurational space. If such groove exists, the motion of multiple ions can be reduced to the motion of single quasiparticle, called the superion, which moves in one-dimensional effective potential. The concept of the superions dramatically reduces the computational complexity of the problem and provides very clear physical interpretation of conduction phenomena in the narrow pores.

Exploring the ion selectivity properties of a large number of simplified binding site models

The ability to discriminate between different cations efficiently is essential for the proper physiological functioning of many membrane transport proteins. One obvious mechanism of ion selectivity is when a binding site is structurally constrained by the protein architecture and its geometry is precisely adapted to fit an ion of a given size. This mechanism is not effective in the case of flexible protein binding sites that are able to deform structurally or to adapt to a bound ion. In this study, the concept of nontrivial ion selectivity arising in a highly flexible protein binding site conceptually represented as a microdroplet of ligands confined to a small volume is explored. The environment imposed by the spatial confinement is a critical feature of the reduced models. A large number of reduced binding site models (1077) comprising typical ion-coordinating ligands (carbonyl, hydroxyl, carboxylate, water) are constructed and characterized for Na(+)/K(+) and Ca(2+)/Ba(2+) size selectivity using free energy perturbation molecular dynamics simulations. Free energies are highly correlated with the sum of ion-ligand and ligand-ligand mean interactions, but the relative balance of those two contributions is different for K(+)-selective and Na(+)-selective binding sites. The analysis indicates that both the number and the type of ligands are important factors in ion selectivity.

Structural basis for the coupling between activation and inactivation gates in K(+) channels

The coupled interplay between activation and inactivation gating is a functional hallmark of K(+) channels. This coupling has been experimentally demonstrated through ion interaction effects and cysteine accessibility, and is associated with a well defined boundary of energetically coupled residues. The structure of the K(+) channel KcsA in its fully open conformation, in addition to four other partial channel openings, richly illustrates the structural basis of activation-inactivation gating. Here, we identify the mechanistic principles by which movements on the inner bundle gate trigger conformational changes at the selectivity filter, leading to the non-conductive C-type inactivated state. Analysis of a series of KcsA open structures suggests that, as a consequence of the hinge-bending and rotation of the TM2 helix, the aromatic ring of Phe 103 tilts towards residues Thr 74 and Thr 75 in the pore-helix and towards Ile 100 in the neighbouring subunit. This allows the network of hydrogen bonds among residues Trp 67, Glu 71 and Asp 80 to destabilize the selectivity filter, allowing entry to its non-conductive conformation. Mutations at position 103 have a size-dependent effect on gating kinetics: small side-chain substitutions F103A and F103C severely impair inactivation kinetics, whereas larger side chains such as F103W have more subtle effects. This suggests that the allosteric coupling between the inner helical bundle and the selectivity filter might rely on straightforward mechanical deformation propagated through a network of steric contacts. Average interactions calculated from molecular dynamics simulations show favourable open-state interaction-energies between Phe 103 and the surrounding residues. We probed similar interactions in the Shaker K(+) channel where inactivation was impaired in the mutant I470A. We propose that side-chain rearrangements at position 103 mechanically couple activation and inactivation in KcsA and a variety of other K(+) channels.

Ion selectivity of the KcsA channel: a perspective from multi-ion free energy landscapes

Potassium (K(+)) channels are specialized membrane proteins that are able to facilitate and regulate the conduction of K(+) through cell membranes. Comprising five specific cation binding sites (S(0)-S(4)) formed by the backbone carbonyl groups of conserved residues common to all K(+) channels, the narrow selectivity filter allows fast conduction of K(+) while being highly selective for K(+) over Na(+). To extend our knowledge of the microscopic mechanism underlying selectivity in K(+) channels, we characterize the free energy landscapes governing the entry and translocation of a Na(+) or a K(+) from the extracellular side into the selectivity filter of KcsA. The entry process of an extracellular ion is examined in the presence of two additional K(+) in the pore, and the three-ion potential of mean force is computed using extensive all-atom umbrella sampling molecular dynamics simulations. A comparison of the potentials of mean force yields a number of important results. First, the free energy minima corresponding to configurations with extracellular K(+) or Na(+) in binding site S(0) or S(1) are similar in depth, suggesting that the thermodynamic selectivity governed by the free energy minima for those two binding sites is insignificant. Second, the free energy barriers between stable multi-ion configurations are generally higher for Na(+) than for K(+), implying that the kinetics of ion conduction is slower when a Na(+) enters the pore. Third, the region corresponding to binding site S(2) near the center of the narrow pore emerges as the most selective for K(+) over Na(+). In particular, while there is a stable minimum for K(+) in site S(2), Na(+) faces a steep free energy increase with no local free energy well in this region. Lastly, analysis shows that selectivity is not correlated with the overall coordination number of the ion entering the pore, but is predominantly affected by changes in the type of coordinating ligands (carbonyls versus water molecules). These results further highlight the importance of the central region near binding site S(2) in the selectivity filter of K(+) channels.

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.

Theory of single-file multiparticle diffusion in narrow pores

Single-file diffusion of multiple strongly interacting particles in a one-dimensional pore is described within a general analytical framework. The theory accounts for nonequilibrium conditions, explicit particle-particle interactions, external potential acting on the particles and the fluctuations of the number of particles due to their exchange with external equilibrium reservoirs. It is shown that the problem can be reduced to a closed hierarchical set of partial differential equations of increasing dimensionality, which can be solved numerically. Our framework allows computing any macroscopic characteristic of multiparticle diffusion in the pore. It is shown that the pore occupancy probabilities and the current are rational functions of external concentrations in the steady state. The theory is tested on a simplified model of the narrow rigid pore inspired by the selectivity filter of biological ion channel. Perspectives and limitations of the theory are discussed.

Conformational dynamics of the inner pore helix of voltage-gated potassium channels

Voltage-gated potassium (Kv) channels control the electrical excitability of neurons and muscles. Despite this key role, how these channels open and close or gate is not fully understood. Gating is usually attributed to the bending and straightening of pore-lining helices at glycine and proline residues. In this work we focused on the role of proline in the Pro-Val-Pro (PVP) motif of the inner S6 helix in the Kv1.2 channel. We started by developing a simple hinged-rod model to fully explore the configurational space of bent helices and we related these configurations to the degree of pore opening. We then carried out fully atomistic simulations of the S6 helices and compared these simulations to the hinged-rod model. Both methods suggest that Kv1 channels are not tightly closed when the inner helices are straight, unlike what is seen in the non-PVP containing channels KcsA and KirBac. These results invite the possibility that the S6 helices may be kinked when Kv1 channels are closed. Our simulations indicate that the wild-type helix adopts multiple spatially distinct configurations, which is consistent with its role in adopting a closed state and an open state. The two most dominant configurational basins correspond to a 6 A movement of the helix tail accompanied by the PVP region undergoing a local alpha-helix to 3(10)-helix transition. We explored how single point mutations affect the propensity of the S6 helix to adopt particular configurations. Interestingly, mutating the first proline, P405 (P473 in Shaker), to alanine completely removed the bistable nature of the S6 helix possibly explaining why this mutation compromises the channel. Next, we considered four other mutations in the area known to affect channel gating and we saw similarly dramatic changes to the helix's dynamics and range of motion. Our results suggest a possible mechanism of helix pore closure and they suggest differences in the closed state of glycine-only channels, like KcsA, and PVP containing channels.

Determination of the charge profile in the KcsA selectivity filter using ab initio calculations and molecular dynamics simulations

The charge profile of K+ and Na+ ions moving in a single file across the filter of the KcsA channel is determined using both molecular dynamics simulations and ab initio calculations. We show a strong correlation between the charge variation and the ion location resulting in a saw-tooth profile, which provides additional information on the influence of charge transfer on the permeation and selectivity of the channel.

Protein dynamics detected in a membrane-embedded potassium channel using two-dimensional solid-state NMR spectroscopy

We report longitudinal (15)N relaxation rates derived from two-dimensional ((15)N, (13)C) chemical shift correlation experiments obtained under magic angle spinning for the potassium channel KcsA-Kv1.3 reconstituted in multilamellar vesicles. Thus, we demonstrate that solid-state NMR can be used to probe residue-specific backbone dynamics in a membrane-embedded protein. Enhanced backbone mobility was detected for two glycine residues within the selectivity filter that are highly conserved in potassium channels and that are of core relevance to the filter structure and ion selectivity.

Ion selectivity in the KcsA potassium channel from the perspective of the ion binding site

To understand the thermodynamic exclusion of Na(+) relative to K(+) from the S(2) site of the selectivity filter, the distribution P(X)(epsilon) (X = K(+) or Na(+)) of the binding energy (epsilon) of the ion with the channel is analyzed using the potential distribution theorem. By expressing the excess chemical potential of the ion as a sum of mean-field epsilon and fluctuation mu(flux,X)(ex) contributions, we find that selectivity arises from a higher value of mu(flux,Na(+))(ex) relative to mu(flux,K(+))(ex). To understand the role of site-site interactions on mu(ex)(flux,X), we decompose P(X)(epsilon) into n-dependent distributions, where n is the number of ion-coordinating ligands within a distance lambda from the ion. For lambda comparable to typical ion-oxygen bond distances, investigations building on this multistate model reveal an inverse correlation between favorable ion-site and site-site interactions: the ion-coordination states that most influence the thermodynamics of the ion are also those for which the binding site is energetically less strained and vice versa. This correlation motivates understanding entropic effects in ion binding to the site and leads to the finding that mu(flux,X)(ex) is directly proportional to the average site-site interaction energy, a quantity that is sensitive to the chemical type of the ligand coordinating the ion. Increasing the coordination number around Na(+) only partially accounts for the observed magnitude of selectivity; acknowledging the chemical type of the ion-coordinating ligand is essential.

Nucleotide recognition by CopA, a Cu+-transporting P-type ATPase

Heavy metal pumps constitute a large subgroup in P-type ion-transporting ATPases. One of the outstanding features is that the nucleotide binding N-domain lacks residues critical for ATP binding in other well-studied P-type ATPases. Instead, they possess an HP-motif and a Gly-rich sequence in the N-domain, and their mutations impair ATP binding. Here, we describe 1.85 A resolution crystal structures of the P- and N-domains of CopA, an archaeal Cu(+)-transporting ATPase, with bound nucleotides. These crystal structures show that CopA recognises the adenine ring completely differently from other P-type ATPases. The crystal structure of the His462Gln mutant, in the HP-motif, a disease-causing mutation in human Cu(+)-ATPases, shows that the Gln side chain mimics the imidazole ring, but only partially, explaining the reduction in ATPase activity. These crystal structures lead us to propose a role of the His and a mechanism for removing Mg(2+) from ATP before phosphoryl transfer.

Dynamics, energetics, and selectivity of the low-K+ KcsA channel structure

Potassium channels are a diverse family of integral membrane proteins through which K(+) can pass selectively. There is ongoing debate about the nature of conformational changes associated with the opening/closing and conductive/nonconductive states of potassium channels. The channels partly exert their function by varying their conductance through a mechanism known as C-type inactivation. Shortly after the activation of K(+) channels, their selectivity filter stops conducting ions at a rate that depends on various stimuli. The molecular mechanism of C-type inactivation has not been fully understood yet. However, the X-ray structure of the KcsA channel obtained in the presence of low K(+) concentration is thought to be representative of a K(+) channel in the C-type inactivated state. Here, extensive, fully atomistic molecular dynamics and free-energy simulations of the low-K(+) KcsA structure in an explicit lipid bilayer are performed to evaluate the stability of this structure and the selectivity of its binding sites. We find that the low-K(+) KcsA structure is stable on the timescale of the molecular dynamics simulations performed, and that ions preferably remain in S1 and S4. In the absence of ions, the selectivity filter evolves toward an asymmetric architecture, as already observed in other computations of the high-K(+) structure of KcsA and KirBac. The low-K(+) KcsA structure is not permeable by Na(+), K(+), or Rb(+), and the selectivity of its binding sites is different from that of the high-K(+) structure.