Electrostatics of the intracellular vestibule of K+ channels

Molecular Dynamics of Ion Conduction through the Selectivity Filter of the NaVAb Sodium Channel

The determination of the atomic structures of voltage-gated bacterial sodium channels using X-ray crystallography has provided a first view of this family of membrane proteins. Molecular dynamics simulations offer one approach to clarify the underlying mechanism of permeation and selectivity in these channels. However, it appears that the intracellular gate of the pore domain is either closed or only open partially in the available X-ray structures. The lack of structure with a fully open intracellular gate poses a special challenge to computational studies aimed at simulating ion conduction. To circumvent this problem, we simulated a model of the Na_VAb channel in which the transmembrane S5 and S6 helices of the pore domain have been truncated to provide direct open access to the intracellular entryway to the pore. Molecular dynamics simulations were carried out over a range of membrane potential and ion concentration of sodium and potassium. The simulations show that the Na_VAb selectivity filter is essentially a cationic pore supporting the conduction of ions at a rate comparable to aqueous diffusion with no significant selectivity for sodium. Conductance and selectivity vary as a function of ion concentration for both cations. Permeation occurs primarily via a knock-on mechanism for both sodium and potassium, although the ion ordering in single file along the pore is not strictly maintained. The character of the outward current appears quite different from the inward current, with a buildup on ions in the selectivity filter prior to escape toward the extracellular side, indicating the presence of a rectification effect that is overcome by nonphysiological applied voltages.

Mechanisms of Drug Binding to Voltage-Gated Sodium Channels

Voltage-gated sodium (Na⁺) channels are expressed in virtually all electrically excitable tissues and are essential for muscle contraction and the conduction of impulses within the peripheral and central nervous systems. Genetic disorders that disrupt the function of these channels produce an array of Na⁺ channelopathies resulting in neuronal impairment, chronic pain, neuromuscular pathologies, and cardiac arrhythmias. Because of their importance to the conduction of electrical signals, Na⁺ channels are the target of a wide variety of local anesthetic, antiarrhythmic, anticonvulsant, and antidepressant drugs. The voltage-gated family of Na⁺ channels is composed of α-subunits that encode for the voltage sensor domains and the Na⁺-selective permeation pore. In vivo, Na⁺ channel α-subunits are associated with one or more accessory β-subunits (β₁-β₄) that regulate gating properties, trafficking, and cell-surface expression of the channels. The permeation pore of Na⁺ channels is divided in two parts: the outer mouth of the pore is the site of the ion selectivity filter, while the inner cytoplasmic pore serves as the channel activation gate. The cytoplasmic lining of the permeation pore is formed by the S6 segments that include highly conserved aromatic amino acids important for drug binding. These residues are believed to undergo voltage-dependent conformational changes that alter drug binding as the channels cycle through the closed, open, and inactivated states. The purpose of this chapter is to broadly review the mechanisms of Na⁺ channel gating and the models used to describe drug binding and Na⁺ channel inhibition.

Energetics of Ion Permeation in an Open-Activated TRPV1 Channel

Ion channels enable diffusion of ions down physiological electrochemical gradients. Modulation of ion permeation is crucial for the physiological functioning of cells, and misregulation of ion channels is linked to a myriad of channelopathies. The ion permeation mechanism in the transient receptor potential (TRP) ion channel family is currently not understood at an atomistic level. In this work, we employed a simulation strategy for ion permeation (molecular-dynamics simulations with bias-exchange metadynamics) to study and compare monovalent (Na(+), K(+)) ion permeation in the open-activated TRP vanniloid-1 (TRPV1) ion channel. Using ∼3.6 μs of simulation trajectories, we obtained atomistic evidence for the nonselective nature of TRPV1. Our analysis shows that solvated monovalent ions permeate through the selectivity filter with comparable energetic barriers via a two-site mechanism. Finally, we confirmed that an intracellular binding site is located between the intracellular gate residues I679 and E684.

A Continuum Poisson-Boltzmann Model for Membrane Channel Proteins

Membrane proteins constitute a large portion of the human proteome and perform a variety of important functions as membrane receptors, transport proteins, enzymes, signaling proteins, and more. Computational studies of membrane proteins are usually much more complicated than those of globular proteins. Here, we propose a new continuum model for Poisson-Boltzmann calculations of membrane channel proteins. Major improvements over the existing continuum slab model are as follows: (1) The location and thickness of the slab model are fine-tuned based on explicit-solvent MD simulations. (2) The highly different accessibilities in the membrane and water regions are addressed with a two-step, two-probe grid-labeling procedure. (3) The water pores/channels are automatically identified. The new continuum membrane model is optimized (by adjusting the membrane probe, as well as the slab thickness and center) to best reproduce the distributions of buried water molecules in the membrane region as sampled in explicit water simulations. Our optimization also shows that the widely adopted water probe of 1.4 Å for globular proteins is a very reasonable default value for membrane protein simulations. It gives the best compromise in reproducing the explicit water distributions in membrane channel proteins, at least in the water accessible pore/channel regions. Finally, we validate the new membrane model by carrying out binding affinity calculations for a potassium channel, and we observe good agreement with the experimental results.

Pore size matters for potassium channel conductance

Ion channels are membrane proteins that mediate efficient ion transport across the hydrophobic core of cell membranes, an unlikely process in their absence. K⁺ channels discriminate K⁺ over cations with similar radii with extraordinary selectivity and display a wide diversity of ion transport rates, covering differences of two orders of magnitude in unitary conductance. The pore domains of large- and small-conductance K⁺ channels share a general architectural design comprising a conserved narrow selectivity filter, which forms intimate interactions with permeant ions, flanked by two wider vestibules toward the internal and external openings. In large-conductance K⁺ channels, the inner vestibule is wide, whereas in small-conductance channels it is narrow. Here we raise the idea that the physical dimensions of the hydrophobic internal vestibule limit ion transport in K⁺ channels, accounting for their diversity in unitary conductance.

Pore dimensions and the role of occupancy in unitary conductance of Shaker K channels

The resistance of the inner vestibule limits Shaker’s conductance.

K channels mediate the selective passage of K⁺ across the plasma membrane by means of intimate interactions with ions at the pore selectivity filter located near the external face. Despite high conservation of the selectivity filter, the K⁺ transport properties of different K channels vary widely, with the unitary conductance spanning a range of over two orders of magnitude. Mutation of Pro475, a residue located at the cytoplasmic entrance of the pore of the small-intermediate conductance K channel Shaker (Pro475Asp (P475D) or Pro475Gln (P475Q)), increases Shaker’s reported ∼20-pS conductance by approximately six- and approximately threefold, respectively, without any detectable effect on its selectivity. These findings suggest that the structural determinants underlying the diversity of K channel conductance are distinct from the selectivity filter, making P475D and P475Q excellent probes to identify key determinants of the K channel unitary conductance. By measuring diffusion-limited unitary outward currents after unilateral addition of 2 M sucrose to the internal solution to increase its viscosity, we estimated a pore internal radius of capture of ∼0.82 Å for all three Shaker variants (wild type, P475D, and P475Q). This estimate is consistent with the internal entrance of the Kv1.2/2.1 structure if the effective radius of hydrated K⁺ is set to ∼4 Å. Unilateral exposure to sucrose allowed us to estimate the internal and external access resistances together with that of the inner pore. We determined that Shaker resistance resides mainly in the inner cavity, whereas only ∼8% resides in the selectivity filter. To reduce the inner resistance, we introduced additional aspartate residues into the internal vestibule to favor ion occupancy. No aspartate addition raised the maximum unitary conductance, measured at saturating [K⁺], beyond that of P475D, suggesting an ∼200-pS conductance ceiling for Shaker. This value is approximately one third of the maximum conductance of the large conductance K (BK) channel (the K channel of highest conductance), reducing the energy gap between their K⁺ transport rates to ∼1 kT. Thus, although Shaker’s pore sustains ion translocation as the BK channel’s does, higher energetic costs of ion stabilization or higher friction with the ion’s rigid hydration cage in its narrower aqueous cavity may entail higher resistance.

Voltage-Gated Sodium Channels: Mechanistic Insights From Atomistic Molecular Dynamics Simulations

The permeation of ions and other molecules across biological membranes is an inherent requirement of all cellular organisms. Ion channels, in particular, are responsible for the conduction of charged species, hence modulating the propagation of electrical signals. Despite the universal physiological implications of this property, the molecular functioning of ion channels remains ambiguous. The combination of atomistic structural data with computational methodologies, such as molecular dynamics (MD) simulations, is now considered routine to investigate structure-function relationships in biological systems. A fuller understanding of conduction, selectivity, and gating, therefore, is steadily emerging due to the applicability of these techniques to ion channels. However, because their structure is known at atomic resolution, studies have consistently been biased toward K(+) channels, thus the molecular determinants of ionic selectivity, activation, and drug blockage in Na(+) channels are often overlooked. The recent increase of available crystallographic data has eminently encouraged the investigation of voltage-gated sodium (NaV) channels via computational methods. Here, we present an overview of simulation studies that have contributed to our understanding of key principles that underlie ionic conduction and selectivity in Na(+) channels, in comparison to the K(+) channel analogs.

Mechanism for attenuated outward conductance induced by mutations in the cytoplasmic pore of Kir2.1 channels

Outward currents through Kir2.1 channels regulate the electrical properties of excitable cells. These currents are subject to voltage-dependent attenuation by the binding of polyamines to high- and low-affinity sites, which leads to inward rectification, thereby controlling cell excitability. To examine the effects of positive charges at the low-affinity site in the cytoplasmic pore on inward rectification, we studied a mutant Kir channel (E224K/H226E) and measured single-channel currents and streaming potentials (V_stream), the latter provide the ratio of water to ions queued in a single-file permeation process in the selectivity filter. The water-ion coupling ratio was near one at a high K⁺ concentration ([K⁺]) for the wild-type channel and increased substantially as [K⁺] decreased. On the other hand, fewer ions occupied the selectivity filter in the mutant at all [K⁺]. A model for the Kir channel involving a K⁺ binding site in the wide pore was introduced. Model analyses revealed that the rate constants associated with the binding and release to and from the wide-pore K⁺ binding site was modified in the mutant. These effects lead to the reduced contribution of a conventional two-ion permeation mode to total conductance, especially at positive potentials, thereby inward rectification.

Extension of the free energy workflow FEW towards implicit solvent/implicit membrane MM-PBSA calculations

BACKGROUND

The number of high-resolution structures of pharmacologically relevant membrane proteins has been strongly increasing. This makes computing relative affinities of chemically similar compounds binding to a membrane protein possible in order to guide decision making in drug design. However, the preparation step of such calculations is time-consuming and complex.

METHODS

We extended the free energy workflow tool FEW, available in AMBER, towards facilitating the setup of molecular dynamics simulations with explicit membrane, and the setup and execution of effective binding energy calculations according to a 1-trajectory implicit solvent/implicit membrane MM-PBSA approach for multiple ligands binding to the same membrane protein.

RESULTS

We validated the implemented protocol initially on two model systems, a sodium ion in the presence of an implicit membrane slab and a proton traversing the M2 proton-channel of the influenza A virus. For the latter, we found a good agreement for several important events along the proton pathway with those obtained in a recent computational study. Finally, we performed a case study on effective binding energy calculations for a set of inhibitors binding to the M2 proton-channel.

CONCLUSIONS

From the case study, we estimate a considerable speed up in the setup and analysis times for implicit solvent/implicit membrane MM-PBSA calculations by the extended version of FEW compared to a manual preparation.

GENERAL SIGNIFICANCE

Together with the overall runtime and the analysis results, this suggests that such type of calculations can be valuable in later stages of drug design projects on membrane proteins. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

The voltage-dependent gate in MthK potassium channels is located at the selectivity filter

Understanding how ion channels open and close their pores is crucial for understanding their physiological roles. We used intracellular quaternary ammonium blockers to locate the voltage-dependent gate in MthK potassium channels from Methanobacterium thermoautotrophicum with electrophysiology and X-ray crystallography. Blockers bind in an aqueous cavity between two putative gates, an intracellular gate and the selectivity filter. Thus, these blockers directly probe gate location: an intracellular gate will prevent binding when closed, whereas a selectivity filter gate will always allow binding. A kinetic analysis of tetrabutylammonium block of single MthK channels combined with X-ray crystallographic analysis of the pore with tetrabutylantimony unequivocally determined that the voltage-dependent gate, like the C-type inactivation gate in eukaryotic channels, is located at the selectivity filter. State-dependent binding kinetics suggests that MthK inactivation leads to conformational changes within the cavity and intracellular pore entrance.

Crystal structure of a voltage-gated K+ channel pore module in a closed state in lipid membranes

Voltage-gated K(+) channels underlie the electrical excitability of cells. Each subunit of the functional tetramer consists of the tandem fusion of two modules, an N-terminal voltage-sensor and a C-terminal pore. To investigate how sensor coupling to the pore generates voltage-dependent channel opening, we solved the crystal structure and characterized the function of a voltage-gated K(+) channel pore in a lipid membrane. The structure of a functional channel in a membrane environment at 3.1 Å resolution establishes an unprecedented connection between channel structure and function. The structure is unique in delineating an ion-occupied ready to conduct selectivity filter, a confined aqueous cavity, and a closed activation gate, embodying a dynamic entity trapped in an unstable closed state.

Voltage-gated sodium channels: biophysics, pharmacology, and related channelopathies

Voltage-gated sodium channels (VGSC) are multi-molecular protein complexes expressed in both excitable and non-excitable cells. They are primarily formed by a pore-forming multi-spanning integral membrane glycoprotein (α-subunit) that can be associated with one or more regulatory β-subunits. The latter are single-span integral membrane proteins that modulate the sodium current (I_Na) and can also function as cell adhesion molecules. In vitro some of the cell-adhesive functions of the β-subunits may play important physiological roles independently of the α-subunits. Other endogenous regulatory proteins named “channel partners” or “channel interacting proteins” (ChiPs) like caveolin-3 and calmodulin/calmodulin kinase II (CaMKII) can also interact and modulate the expression and/or function of VGSC. In addition to their physiological roles in cell excitability and cell adhesion, VGSC are the site of action of toxins (like tetrodotoxin and saxitoxin), and pharmacologic agents (like antiarrhythmic drugs, local anesthetics, antiepileptic drugs, and newly developed analgesics). Mutations in genes that encode α- and/or β-subunits as well as the ChiPs can affect the structure and biophysical properties of VGSC, leading to the development of diseases termed sodium “channelopathies”.  This review will outline the structure, function, and biophysical properties of VGSC as well as their pharmacology and associated channelopathies and highlight some of the recent advances in this field.

Permeation properties of the hair cell mechanotransducer channel provide insight into its molecular structure

Mechanoelectric transducer (MET) channels, located near stereocilia tips, are opened by deflecting the hair bundle of sensory hair cells. Defects in this process result in deafness. Despite this critical function, the molecular identity of MET channels remains a mystery. Inherent channel properties, particularly those associated with permeation, provide the backbone for the molecular identification of ion channels. Here, a novel channel rectification mechanism is identified, resulting in a reduced pore size at positive potentials. The apparent difference in pore dimensions results from Ca(2+) binding within the pore, occluding permeation. Driving force for permeation at hyperpolarized potentials is increased because Ca(2+) can more easily be removed from binding within the pore due to the presence of an electronegative external vestibule that dehydrates and concentrates permeating ions. Alterations in Ca(2+) binding may underlie tonotopic and Ca(2+)-dependent variations in channel conductance. This Ca(2+)-dependent rectification provides targets for identifying the molecular components of the MET channel.

Concentration dependent ion selectivity in VDAC: a molecular dynamics simulation study

The voltage-dependent anion channel (VDAC) forms the major pore in the outer mitochondrial membrane. Its high conducting open state features a moderate anion selectivity. There is some evidence indicating that the electrophysiological properties of VDAC vary with the salt concentration. Using a theoretical approach the molecular basis for this concentration dependence was investigated. Molecular dynamics simulations and continuum electrostatic calculations performed on the mouse VDAC1 isoform clearly demonstrate that the distribution of fixed charges in the channel creates an electric field, which determines the anion preference of VDAC at low salt concentration. Increasing the salt concentration in the bulk results in a higher concentration of ions in the VDAC wide pore. This event induces a large electrostatic screening of the charged residues promoting a less anion selective channel. Residues that are responsible for the electrostatic pattern of the channel were identified using the molecular dynamics trajectories. Some of these residues are found to be conserved suggesting that ion permeation between different VDAC species occurs through a common mechanism. This inference is buttressed by electrophysiological experiments performed on bean VDAC32 protein akin to mouse VDAC.

Constant electric field simulations of the membrane potential illustrated with simple systems

Advances in modern computational methods and technology make it possible to carry out extensive molecular dynamics simulations of complex membrane proteins based on detailed atomic models. The ultimate goal of such detailed simulations is to produce trajectories in which the behavior of the system is as realistic as possible. A critical aspect that requires consideration in the case of biological membrane systems is the existence of a net electric potential difference across the membrane. For meaningful computations, it is important to have well validated methodologies for incorporating the latter in molecular dynamics simulations. A widely used treatment of the membrane potential in molecular dynamics consists of applying an external uniform electric field E perpendicular to the membrane. The field acts on all charged particles throughout the simulated system, and the resulting applied membrane potential V is equal to the applied electric field times the length of the periodic cell in the direction perpendicular to the membrane. A series of test simulations based on simple membrane-slab models are carried out to clarify the consequences of the applied field. These illustrative tests demonstrate that the constant-field method is a simple and valid approach for accounting for the membrane potential in molecular dynamics studies of biomolecular systems. This article is part of a Special Issue entitled: Membrane protein structure and function.

The crystal structure of a voltage-gated sodium channel

Voltage-gated sodium (Na(V)) channels initiate electrical signalling in excitable cells and are the molecular targets for drugs and disease mutations, but the structural basis for their voltage-dependent activation, ion selectivity and drug block is unknown. Here we report the crystal structure of a voltage-gated Na(+) channel from Arcobacter butzleri (NavAb) captured in a closed-pore conformation with four activated voltage sensors at 2.7 Å resolution. The arginine gating charges make multiple hydrophilic interactions within the voltage sensor, including unanticipated hydrogen bonds to the protein backbone. Comparisons to previous open-pore potassium channel structures indicate that the voltage-sensor domains and the S4-S5 linkers dilate the central pore by pivoting together around a hinge at the base of the pore module. The NavAb selectivity filter is short, ∼4.6 Å wide, and water filled, with four acidic side chains surrounding the narrowest part of the ion conduction pathway. This unique structure presents a high-field-strength anionic coordination site, which confers Na(+) selectivity through partial dehydration via direct interaction with glutamate side chains. Fenestrations in the sides of the pore module are unexpectedly penetrated by fatty acyl chains that extend into the central cavity, and these portals are large enough for the entry of small, hydrophobic pore-blocking drugs. This structure provides the template for understanding electrical signalling in excitable cells and the actions of drugs used for pain, epilepsy and cardiac arrhythmia at the atomic level.

The contribution of individual subunits to the coupling of the voltage sensor to pore opening in Shaker K channels: effect of ILT mutations in heterotetramers

Voltage-gated ion channels couple conformational change(s) of the voltage-sensing domain to those of the opening of an intracellular gate to allow ionic conduction. Much larger positive potentials are required to couple these conformational changes to the opening of the gate of Shaker K(+) channels with the concurrent mutations V369I, I372L, and S376T (ILT) at the N-terminal end of the S4 segment. We used cut-open oocyte voltage clamp to study the biophysical and thermodynamical properties of heterotetrameric concatemerized channels with different stoichiometries of ILT mutations. The voltage-sensing domains of ILT mutant channels require smaller depolarization to activate but their intracellular gate does not immediately follow the movement of the voltage-sensing domain, requiring larger depolarization to open. Our results demonstrate that each subunit contributes equally to the rightward shift of the conductance-voltage relationship and that a single ILT-containing subunit is sufficient to induce a large enthalpic and entropic barrier, limiting opening of the intracellular gate.

Voltage profile along the permeation pathway of an open channel

For ion channels, the transmembrane potential plays a critical role by acting as a driving force for permeant ions. At the microscopic level, the transmembrane potential is thought to decay nonlinearly across the ion permeation pathway because of the irregular three-dimensional shape of the channel's pore. By taking advantage of the current structural and functional understanding of cyclic nucleotide-gated channels, in this study we experimentally explore the transmembrane potential's distribution across the open pore. As a readout for the voltage drop, we engineered cysteine residues along the selectivity filter and scanned the sensitivity of their modification rates by Ag(+) to the transmembrane potential. The experimental data, which indicate that the majority of the electric field drops across the selectivity filter, are in good agreement with continuum electrostatic calculations using a homology model of an open CNG channel. By focusing the transmembrane potential across the selectivity filter, the electromotive driving force is coupled with the movement of permeant ions in the filter, maximizing the efficiency of this process.

APBSmem: a graphical interface for electrostatic calculations at the membrane

Electrostatic forces are one of the primary determinants of molecular interactions. They help guide the folding of proteins, increase the binding of one protein to another and facilitate protein-DNA and protein-ligand binding. A popular method for computing the electrostatic properties of biological systems is to numerically solve the Poisson-Boltzmann (PB) equation, and there are several easy-to-use software packages available that solve the PB equation for soluble proteins. Here we present a freely available program, called APBSmem, for carrying out these calculations in the presence of a membrane. The Adaptive Poisson-Boltzmann Solver (APBS) is used as a back-end for solving the PB equation, and a Java-based graphical user interface (GUI) coordinates a set of routines that introduce the influence of the membrane, determine its placement relative to the protein, and set the membrane potential. The software Jmol is embedded in the GUI to visualize the protein inserted in the membrane before the calculation and the electrostatic potential after completing the computation. We expect that the ease with which the GUI allows one to carry out these calculations will make this software a useful resource for experimenters and computational researchers alike. Three examples of membrane protein electrostatic calculations are carried out to illustrate how to use APBSmem and to highlight the different quantities of interest that can be calculated.

An iris-like mechanism of pore dilation in the CorA magnesium transport system

Magnesium translocation across cell membranes is essential for numerous physiological processes. Three recently reported crystal structures of the CorA magnesium transport system revealed a surprising architecture, with a bundle of giant alpha-helices forming a 60-A-long pore that extends beyond the membrane before widening into a funnel-shaped cytosolic domain. The presence of divalent cations in putative intracellular regulation sites suggests that these structures correspond to the closed conformation of CorA. To examine the nature of the conduction pathway, we performed 110-ns molecular-dynamics simulations of two of these structures in a lipid bilayer with and without regulatory ions. The results show that a 15-A-long hydrophobic constriction straddling the membrane-cytosol interface constitutes a steric bottleneck whose location coincides with an electrostatic barrier opposing cation translocation. In one of the simulations, structural relaxation after the removal of regulatory ions led to concerted changes in the tilt of the pore helices, resulting in iris-like dilation and spontaneous hydration of the hydrophobic neck. This simple and robust mechanism is consistent with the regulation of pore opening by intracellular magnesium concentration, and explains the unusual architecture of CorA.

Voltage-dependent gating in a "voltage sensor-less" ion channel

The voltage sensitivity of voltage-gated cation channels is primarily attributed to conformational changes of a four transmembrane segment voltage-sensing domain, conserved across many levels of biological complexity. We have identified a remarkable point mutation that confers significant voltage dependence to Kir6.2, a ligand-gated channel that lacks any canonical voltage-sensing domain. Similar to voltage-dependent Kv channels, the Kir6.2[L157E] mutant exhibits time-dependent activation upon membrane depolarization, resulting in an outwardly rectifying current-voltage relationship. This voltage dependence is convergent with the intrinsic ligand-dependent gating mechanisms of Kir6.2, since increasing the membrane PIP2 content saturates Po and eliminates voltage dependence, whereas voltage activation is more dramatic when channel Po is reduced by application of ATP or poly-lysine. These experiments thus demonstrate an inherent voltage dependence of gating in a "ligand-gated" K+ channel, and thereby provide a new view of voltage-dependent gating mechanisms in ion channels. Most interestingly, the voltage- and ligand-dependent gating of Kir6.2[L157E] is highly sensitive to intracellular [K+], indicating an interaction between ion permeation and gating. While these two key features of channel function are classically dealt with separately, the results provide a framework for understanding their interaction, which is likely to be a general, if latent, feature of the superfamily of cation channels.

PIP(2)-binding site in Kir channels: definition by multiscale biomolecular simulations

Phosphatidylinositol bisphosphate (PIP(2)) is an activator of mammalian inwardly rectifying potassium (Kir) channels. Multiscale simulations, via a sequential combination of coarse-grained and atomistic molecular dynamics, enabled exploration of the interactions of PIP(2) molecules within the inner leaflet of a lipid bilayer membrane with possible binding sites on Kir channels. Three Kir channel structures were investigated: X-ray structures of KirBac1.1 and of a Kir3.1-KirBac1.3 chimera and a homology model of Kir6.2. Coarse-grained simulations of the Kir channels in PIP(2)-containing lipid bilayers identified the PIP(2)-binding site on each channel. These models of the PIP(2)-channel complexes were refined by conversion to an atomistic representation followed by molecular dynamics simulation in a lipid bilayer. All three channels were revealed to contain a conserved binding site at the N-terminal end of the slide (M0) helix, at the interface between adjacent subunits of the channel. This binding site agrees with mutagenesis data and is in the proximity of the site occupied by a detergent molecule in the Kir chimera channel crystal. Polar contacts in the coarse-grained simulations corresponded to long-lived electrostatic and H-bonding interactions between the channel and PIP(2) in the atomistic simulations, enabling identification of key side chains.

Structure and electrostatic property of cytoplasmic domain of ZntB transporter

ZntB is the distant homolog of CorA Mg(2+) transporter within the metal ion transporter superfamily. It was early reported that the ZntB from Salmonella typhimurium facilitated efflux of Zn(2+) and Cd(2+), but not Mg(2+). Here, we report the 1.90 A crystal structure of the intracellular domain of ZntB from Vibrio parahemolyticus. The domain forms a funnel-shaped homopentamer that is similar to the full-length CorA from Thermatoga maritima, but differs from two previously reported dimeric structures of truncated CorA intracellular domains. However, no Zn(2+) or Cd(2+) binding sites were identified in the high-resolution structure. Instead, 25 well-defined Cl(-) ions were observed and some of these binding sites are highly conserved within the ZntB family. Continuum electrostatics calculations suggest that the central pore of the funnel is highly attractive for cations, especially divalents. The presence of the bound Cl(-) ions increases the stability of cations along the pore suggesting they could be important in enhancing cation transport.

Atypical mechanism of conduction in potassium channels

Potassium channels can conduct passively K+ ions with rates of up to approximately 10(8) ions per second at physiological conditions, and they are selective to these species by a factor of 10(4) over Na+ ions. Ion conduction has been proposed to involve transitions between 2 main states, with 2 or 3 K+ ions occupying the selectivity filter separated by an intervening water molecule. The largest free energy barrier of such a process was reported to be of the order of 2-3 kcal mol(-1). Here, we present an alternative mechanism for conduction of K+ in potassium channels where site vacancies are involved, and we propose that coexistence of several ion permeation mechanisms is energetically possible. Conduction can be described as a more anarchic phenomenon than previously characterized by the concerted translocations of K+-water-K+.

A structural model for K2P potassium channels based on 23 pairs of interacting sites and continuum electrostatics

K(2P)Ø, the two-pore domain potassium background channel that determines cardiac rhythm in Drosophila melanogaster, and its homologues that establish excitable membrane activity in mammals are of unknown structure. K(2P) subunits have two pore domains flanked by transmembrane (TM) spans: TM1-P1-TM2-TM3-P2-TM4. To establish spatial relationships in K(2P)Ø, we identified pairs of sites that display electrostatic compensation. Channels silenced by the addition of a charge in pore loop 1 (P1) or P2 were restored to function by countercharges at specific second sites. A three-dimensional homology model was determined using the crystal structure of K(V)1.2, effects of K(2P)Ø mutations to establish alignment, and compensatory charge-charge pairs. The model was refined and validated by continuum electrostatic free energy calculations and covalent linkage of introduced cysteines. K(2P) channels use two subunits arranged so that the P1 and P2 loops contribute to one pore, identical P loops face each other diagonally across the pore, and the channel complex has bilateral symmetry with a fourfold symmetric selectivity filter.

Long-pore electrostatics in inward-rectifier potassium channels

Inward-rectifier potassium (Kir) channels differ from the canonical K⁺ channel structure in that they possess a long extended pore (∼85 Å) for ion conduction that reaches deeply into the cytoplasm. This unique structural feature is presumably involved in regulating functional properties specific to Kir channels, such as conductance, rectification block, and ligand-dependent gating. To elucidate the underpinnings of these functional roles, we examine the electrostatics of an ion along this extended pore. Homology models are constructed based on the open-state model of KirBac1.1 for four mammalian Kir channels: Kir1.1/ROMK, Kir2.1/IRK, Kir3.1/GIRK, and Kir6.2/KATP. By solving the Poisson-Boltzmann equation, the electrostatic free energy of a K⁺ ion is determined along each pore, revealing that mammalian Kir channels provide a favorable environment for cations and suggesting the existence of high-density regions in the cytoplasmic domain and cavity. The contribution from the reaction field (the self-energy arising from the dielectric polarization induced by the ion's charge in the complex geometry of the pore) is unfavorable inside the long pore. However, this is well compensated by the electrostatic interaction with the static field arising from the protein charges and shielded by the dielectric surrounding. Decomposition of the static field provides a list of residues that display remarkable correspondence with existing mutagenesis data identifying amino acids that affect conduction and rectification. Many of these residues demonstrate interactions with the ion over long distances, up to 40 Å, suggesting that mutations potentially affect ion or blocker energetics over the entire pore. These results provide a foundation for understanding ion interactions in Kir channels and extend to the study of ion permeation, block, and gating in long, cation-specific pores.

Trans-channel interactions in batrachotoxin-modified skeletal muscle sodium channels: voltage-dependent block by cytoplasmic amines, and the influence of mu-conotoxin GIIIA derivatives and permeant ions

External mu-conotoxins and internal amine blockers inhibit each other's block of voltage-gated sodium channels. We explore the basis of this interaction by measuring the shifts in voltage-dependence of channel inhibition by internal amines induced by two mu-conotoxin derivatives with different charge distributions and net charges. Charge changes on the toxin were made at residue 13, which is thought to penetrate most deeply into the channel, making it likely to have the strongest individual interaction with an internal charged ligand. When an R13Q or R13E molecule was bound to the channel, the voltage dependence of diethylammonium (DEA)-block shifted toward more depolarized potentials (23 mV for R13Q, and 16 mV for R13E). An electrostatic model of the repulsion between DEA and the toxin simulated these data, with a distance between residue 13 of the mu-conotoxin and the DEA-binding site of approximately 15 A. Surprisingly, for tetrapropylammonium, the shifts were only 9 mV for R13Q, and 7 mV for R13E. The smaller shifts associated with R13E, the toxin with a smaller net charge, are generally consistent with an electrostatic interaction. However, the smaller shifts observed for tetrapropylammonium than for DEA suggest that other factors must be involved. Two observations indicate that the coupling of permeant ion occupancy of the channel to blocker binding may contribute to the overall amine-toxin interaction: 1), R13Q binding decreases the apparent affinity of sodium for the conducting pore by approximately 4-fold; and 2), increasing external [Na(+)] decreases block by DEA at constant voltage. Thus, even though a number of studies suggest that sodium channels are occupied by no more than one ion most of the time, measurable coupling occurs between permeant ions and toxin or amine blockers. Such interactions likely determine, in part, the strength of trans-channel, amine-conotoxin interactions.

Ion channel gates: comparative analysis of energy barriers

The energetic profile of an ion translated along the axis of an ion channel should reveal whether the structure corresponds to a functionally open or closed state of the channel. In this study, we explore the combined use of Poisson-Boltzmann electrostatic calculations and evaluation of van der Waals interactions between ion and pore to provide an initial appraisal of the gating state of a channel. This approach is exemplified by its application to the bacterial inward rectifier potassium channel KirBac3.1, where it reveals the closed gate to be formed by a ring of leucine (L124) side chains. We have extended this analysis to a comparative survey of gating profiles, including model hydrophobic nanopores, the nicotinic acetylcholine receptor, and a number of potassium channel structures and models. This enables us to identify three gating regimes, and to show the limitation of this computationally inexpensive method. For a (closed) gate radius of 0.4 nm < R < 0.8 nm, a hydrophobic gate may be present. For a gate radius of 0.2 nm < R < 0.4 nm, both electrostatic and van der Waals interactions will contribute to the barrier height. Below R = 0.2 nm, repulsive van der Waals interactions are likely to dominate, resulting in a sterically occluded gate. In general, the method is more useful when the channel is wider; for narrower channels, the flexibility of the protein may allow otherwise-unsurmountable energetic barriers to be overcome.

Generation, comparison, and merging of pathways between protein conformations: gating in K-channels

We present a general framework for the generation, alignment, comparison, and hybridization of motion pathways between two known protein conformations. The framework, which is rooted in probabilistic motion-planning techniques in robotics, allows for the efficient generation of collision-free motion pathways, while considering a wide range of degrees of freedom involved in the motion. Within the framework, we provide the means to hybridize pathways, thus producing, the motion pathway of the lowest energy barrier out of the many pathways proposed by our algorithm. This method for comparing and hybridizing pathways is modular, and may be used within the context of molecular dynamics and Monte Carlo simulations. The framework was implemented within the Rosetta software suite, where the protein is represented in atomic detail. The K-channels switch between open and closed conformations, and we used the overall framework to investigate this transition. Our analysis suggests that channel-opening may follow a three-phase pathway. First, the channel unlocks itself from the closed state; second, it opens; and third, it locks itself in the open conformation. A movie that depicts the proposed pathway is available in the Supplementary Material (Movie S1) and at http://www.cs.tau.ac.il/~angela/SuppKcsA.html.

Double bilayers and transmembrane gradients: a molecular dynamics study of a highly charged peptide

The position and extent of movement of a charged peptide within a membrane bilayer provides much controversy. In our study, we have examined the nature of the highly charged helix-turn-helix motif (S3b and S4) to address how a highly charged peptide is stabilized within a bilayer in the presence of various transmembrane electrical potentials. Our double-bilayer simulation results show how the variation of the salt concentrations between the inner and outer bath establishes a transmembrane potential. Our results also show that important features of the peptide affected by changes in electrical potential are the center of mass depth, the swivel/kink degrees of conformation, and the hydrogen-bonding patterns. As the voltage gradient across the bilayer increased, the center of mass of the peptide shifted in a direction toward the outer bath. The peptide also has a higher percent helical content and the swivel/kink conformation is more rigid for nonpolarized systems where no voltage drop occurred between salt baths. Our results also provide some suggestions for how this domain may be affected by environmental changes as part of the voltage sensor in a K-channel.

Molecular mechanism of pH sensing in KcsA potassium channels

The bacterial potassium channel KcsA is gated by high concentrations of intracellular protons, allowing the channel to open at pH < 5.5. Despite prior attempts to determine the mechanism responsible for pH gating, the proton sensor has remained elusive. We have constructed a KcsA channel mutant that remains open up to pH 9.0 by replacing key ionizable residues from the N and C termini of KcsA with residues mimicking their protonated counterparts with respect to charge. A series of individual and combined mutations were investigated by using single-channel recordings in lipid bilayers. We propose that these residues are the proton-binding sites and at neutral pH they form a complex network of inter- and intrasubunit salt bridges and hydrogen bonds near the bundle crossing that greatly stabilize the closed state. In our model, these residues change their ionization state at acidic pH, thereby disrupting this network, modifying the electrostatic landscape near the channel gate, and favoring channel opening.

Intrinsic electrostatic potential in the BK channel pore: role in determining single channel conductance and block

The internal vestibule of large-conductance Ca²⁺ voltage-activated K⁺ (BK) channels contains a ring of eight negative charges not present in K⁺ channels of lower conductance (Glu386 and Glu389 in hSlo) that modulates channel conductance through an electrostatic mechanism (Brelidze, T.I., X. Niu, and K.L. Magleby. 2003. Proc. Natl. Acad. Sci. USA. 100:9017–9022). In BK channels there are also two acidic amino acid residues in an extracellular loop (Asp326 and Glu329 in hSlo). To determine the electrostatic influence of these charges on channel conductance, we expressed wild-type BK channels and mutants E386N/E389N, D326N, E329Q, and D326N/E329Q channels on Xenopus laevis oocytes, and measured the expressed currents under patch clamp. Contribution of E329 to the conductance is negligible and single channel conductance of D326N/E329Q channels measured at 0 mV in symmetrical 110 mM K⁺ was 18% lower than the control. Current–voltage curves displayed weak outward rectification for D326N and the double mutant. The conductance differences between the mutants and wild-type BK were caused by an electrostatic effect since they were enhanced at low K⁺ (30 mM) and vanished at high K⁺ (1 M K⁺). We determine the electrostatic potential change, Δφ, caused by the charge neutralization using TEA⁺ block for the extracellular charges and Ba²⁺ for intracellular charges. We measured 13 ± 2 mV for Δφ at the TEA⁺ site when turning off the extracellular charges, and 17 ± 2 mV for the Δφ at the Ba²⁺ site when the intracellular charges were turned off. To understand the electrostatic effect of charge neutralizations, we determined Δφ using a BK channel molecular model embedded in a lipid bilayer and solving the Poisson-Boltzmann equation. The model explains the experimental results adequately and, in particular, gives an economical explanation to the differential effect on the conductance of the neutralization of charges D326 and E329.

A tyrosine substitution in the cavity wall of a k channel induces an inverted inactivation

Ion permeation and gating kinetics of voltage-gated K channels critically depend on the amino-acid composition of the cavity wall. Residue 470 in the Shaker K channel is an isoleucine, making the cavity volume in a closed channel insufficiently large for a hydrated K(+) ion. In the cardiac human ether-a-go-go-related gene channel, which exhibits slow activation and fast inactivation, the corresponding residue is tyrosine. To explore the role of a tyrosine at this position in the Shaker channel, we studied I470Y. The activation became slower, and the inactivation faster and more complex. At +60 mV the channel inactivated with two distinct rates (tau(1) = 20 ms, tau(2) = 400 ms). Experiments with tetraethylammonium and high K(+) concentrations suggest that the slower component was of the P/C-type. In addition, an inactivation component with inverted voltage dependence was introduced. A step to -40 mV inactivates the channel with a time constant of 500 ms. Negative voltage steps do not cause the channel to recover from this inactivated state (tau >> 10 min), whereas positive voltage steps quickly do (tau = 2 ms at +60 mV). The experimental findings can be explained by a simple branched kinetic model with two inactivation pathways from the open state.

Evaluation of ion binding to DNA duplexes using a size-modified Poisson-Boltzmann theory

Poisson-Boltzmann (PB) theory is among the most widely applied electrostatic theories in biological and chemical science. Despite its reasonable success in explaining a wide variety of phenomena, it fails to incorporate two basic physical effects, ion size and ion-ion correlations, into its theoretical treatment. Recent experimental work has shown significant deviations from PB theory in competitive monovalent and divalent ion binding to a DNA duplex. The experimental data for monovalent binding are consistent with a hypothesis that attributes these deviations to counterion size. To model the observed differences, we have generalized an existing size-modified Poisson-Boltzmann (SMPB) theory and developed a new numerical implementation that solves the generalized theory around complex, atomistic representations of biological molecules. The results of our analysis show that good agreement to data at monovalent ion concentrations up to approximately 150 mM can be attained by adjusting the ion-size parameters in the new size-modified theory. SMPB calculations employing calibrated ion-size parameters predict experimental observations for other nucleic acid structures and salt conditions, demonstrating that the theory is predictive. We are, however, unable to model the observed deviations in the divalent competition data with a theory that only accounts for size but neglects ion-ion correlations, highlighting the need for theoretical descriptions that further incorporate ion-ion correlations. The accompanying numerical solver has been released publicly, providing the general scientific community the ability to compute SMPB solutions around a variety of different biological structures with only modest computational resources.

Accessibility of four arginine residues on the S4 segment of the Bacillus halodurans sodium channel

The voltage-gated Na(+) channel of Bacillus halodurans (NaChBac) is composed of six transmembrane segments (S1-S6), with a pore-forming region composed of segments S5 and S6 and a voltage-sensing domain composed of segments S1-S4. The S4 segment forms the core of the voltage sensor. We explored the accessibility of four arginine residues on the S4 segment of NaChBac, which are positioned at every third position from each other. These arginine residues on the S4 segment were replaced with cysteines using site-directed mutagenesis. Na(+) currents were recorded using the whole-cell configuration of the patch-clamp technique. We tested the effect of the sulfhydryl reagents applied from inside and outside the cellular space in the open and closed conformations. Structural models of the voltage sensor of NaChBac were constructed based on the recently crystallized KvAP and Kv1.2 K(+) channels to visualize arginine residue accessibility. Our results suggest that arginine accessibility did not change significantly between the open and closed conformations, supporting the idea of a small movement of the S4 segment during gating. Molecular modeling of the closed conformation also supported a small movement of S4, which is mainly characterized by a rotation and a tilt along the periphery of the pore. Interestingly, the second arginine residue of the S4 segment (R114) was accessible to sulfhydryl reagents from both sides of the membrane in the closed conformation and, based on our model, seemed to be at the junction of the intracellular and extracellular water crevices.

Ion conductance vs. pore gating and selectivity in KcsA channel: Modeling achievements and perspectives

KcsA potassium channel belongs to a wide family of allosteric proteins that switch between closed and open states conformations in response to a stimulus, and act as a regulator of cation activity in living cells. The gating mechanism and cation selectivity of such channels have been extensively studied in the literature, with a revival emphasis these latter years, due to the publication of the crystallized structure of KcsA. Despite the increasing number of research and review papers on these topics, quantitative interpretation of these processes at the atomic scale is far from achieved. On the basis of available experimental and theoretical data, and by including our recent results, we review the progresses in this field of activity and discuss the weaknesses that should be corrected. In this spirit, we partition the channel into the filter, cavity, extra and intracellular media, in order to analyze separately the specificity of each region. Special emphasis is brought to the study of an open state for the channel and to the different properties generated by the opening. The influence of water as a structural and dynamical component of the channel properties in closed and open states, as well as in the sequential motions of the cations, is analyzed using molecular dynamics simulations and ab initio calculations. The polarization and charge transfer effects on the ions' dynamics and kinetics are discussed in terms of partial charge models.

Molecular Dynamics Simulation of the Cytosolic Mouth in Kcv-Type Potassium Channels

The functional effect of mutations near the intracellular mouth of the short viral Kcv potassium channel was studied by molecular dynamics simulations. As a model system we used the analogously mutated and truncated KirBac1.1, a channel with known crystal structure that shares genuine local sequence motifs with Kcv. By a novel simulated annealing methodology for structural averaging, information about the structure and dynamics of the intracellular mouth was extracted and complemented by Poisson-Boltzmann and 3D-RISM (reference interaction site model) integral equation theory for the determination of the K+ free energy surface. Besides the wild-type analogue of Kcv with its experimental reference activity (truncated KirBac1.1), two variants were studied: a deletion mutant where the N-terminus is further truncated by eight amino acids, showing inactivity in the Kcv reference system, and a point mutant where the kink-forming proline at position 13 is substituted by alanine, resulting in hyperactivity. The computations reveal that the change of activity is closely related to a hydrophilic intracellular constriction formed by the C-terminal residues of the monomers. Hyperactivity of the point mutant is correlated with both sterical and electrostatic factors, while inactivity of the deletion mutant is related to a loss of specific salt bridge patterns between the C- and N-terminus at the constriction and to the consequences for ion passage barriers, as revealed by integral equation theory. The cytosolic gate, however, is probably formed by the N-terminal segment up to the proline kink and not by the constriction. The results are compared with design principles found for other channels.

Charge at the lidocaine binding site residue Phe-1759 affects permeation in human cardiac voltage-gated sodium channels

Our homology molecular model of the open/inactivated state of the Na(+) channel pore predicts, based on extensive mutagenesis data, that the local anaesthetic lidocaine docks eccentrically below the selectivity filter, such that physical occlusion is incomplete. Electrostatic field calculations suggest that the drug's positively charged amine produces an electrostatic barrier to permeation. To test the effect of charge at this pore level on permeation in hNa(V)1.5 we replaced Phe-1759 of domain IVS6, the putative binding site for lidocaine's alkylamino end, with positively and negatively charged residues as well as the neutral cysteine and alanine. These mutations eliminated use-dependent lidocaine block with no effect on tonic/rested state block. Mutant whole cell currents were kinetically similar to wild type (WT). Single channel conductance (gamma) was reduced from WT in both F1759K (by 38%) and F1759R (by 18%). The negatively charged mutant F1759E increased gamma by 14%, as expected if the charge effect were electrostatic, although F1759D was like WT. None of the charged mutations affected Na(+)/K(+) selectivity. Calculation of difference electrostatic fields in the pore model predicted that lidocaine produced the largest positive electrostatic barrier, followed by lysine and arginine, respectively. Negatively charged glutamate and aspartate both lowered the barrier, with glutamate being more effective. Experimental data were in rank order agreement with the predicted changes in the energy profile. These results demonstrate that permeation rate is sensitive to the inner pore electrostatic field, and they are consistent with creation of an electrostatic barrier to ion permeation by lidocaine's charge.

Neutralization of a single arginine residue gates open a two-pore domain, alkali-activated K+ channel

Potassium channels share a common selectivity filter that determines the conduction characteristics of the pore. Diversity in K+ channels is given by how they are gated open. TASK-2, TALK-1, and TALK-2 are two-pore region (2P) KCNK K+ channels gated open by extracellular alkalinization. We have explored the mechanism for this alkalinization-dependent gating using molecular simulation and site-directed mutagenesis followed by functional assay. We show that the side chain of a single arginine residue (R224) near the pore senses pH in TASK-2 with an unusual pKa of 8.0, a shift likely due to its hydrophobic environment. R224 would block the channel through an electrostatic effect on the pore, a situation relieved by its deprotonation by alkalinization. A lysine residue in TALK-2 fulfills the same role but with a largely unchanged pKa, which correlates with an environment that stabilizes its positive charge. In addition to suggesting unified alkaline pH-gating mechanisms within the TALK subfamily of channels, our results illustrate in a physiological context the principle that hydrophobic environment can drastically modulate the pKa of charged amino acids within a protein.

Mechanism of intracellular block of the KcsA K+ channel by tetrabutylammonium: insights from X-ray crystallography, electrophysiology and replica-exchange molecular dynamics simulations

The mechanism of intracellular blockade of the KcsA potassium channel by tetrabutylammonium (TBA) is investigated through functional, structural and computational studies. Using planar-membrane electrophysiological recordings, we characterize the binding kinetics as well as the dependence on the transmembrane voltage and the concentration of the blocker. It is found that the apparent affinity of the complex is significantly greater than that of any of the eukaryotic K(+) channels studied previously, and that the off-rate increases with the applied transmembrane voltage. In addition, we report a crystal structure of the KcsA-TBA complex at 2.9 A resolution, with TBA bound inside the large hydrophobic cavity located at the center of the channel, consistent with the results of previous functional and structural studies. Of particular interest is the observation that the presence of TBA has a negligible effect on the channel structure and on the position of the potassium ions occupying the selectivity filter. Inspection of the electron density corresponding to TBA suggests that the ligand may adopt more than one conformation in the complex, though the moderate resolution of the data precludes a definitive interpretation on the basis of the crystallographic refinement methods alone. To provide a rationale for these observations, we carry out an extensive conformational sampling of an atomic model of TBA bound in the central cavity of KcsA, using the Hamiltonian replica-exchange molecular dynamics simulation method. Comparison of the simulated and experimental density maps indicates that the latter does reflect at least two distinct binding orientations of TBA. The simulations show also that the relative population of these binding modes is dependent on the ion configuration occupying the selectivity filter, thus providing a clue to the nature of the voltage-dependence of the binding kinetics.

Application of the Poisson-Nernst-Planck theory with space-dependent diffusion coefficients to KcsA

The Poisson-Nernst-Planck electrodiffusion theory serves to compute charge fluxes and is here applied to the ion current through a protein channel. KcsA was selected as an example because of the abundance of experimental and theoretical data. The potassium channels MthK and KvAP were used as templates to define two open channel models for KcsA. Channel boundary surfaces and protein charge distributions were defined according to atomic radii and partial atomic charges. To establish the sensitivity of the results to these parameters, two different sets were used. Assigning the potassium diffusion coefficients equal to the value for free-diffusion in water (1.96 x 10(-9) m(2)/s), the computed currents overestimated the experimental data. Ion distributions inside the channel suggest that the overestimate is not due to an excess of charge shielding. A good agreement with the experimental data was achieved by reducing the potassium diffusion coefficient inside the channel to 1.96 x 10(-10) m(2)/s, a value of substantial motility but nonetheless in accord with the intuitive notion that the channel has a high affinity for the ions and therefore slows them down. These results are independent of the open channel model and the parameterization adopted for atomic radii and partial atomic charges. The method offers a reliable estimate of the channel current with low computational effort.