A comparison between two prokaryotic potassium channels (KirBac1.1 and KcsA) in a molecular dynamics (MD) simulation study

The two potassium ion channels KirBac1.1 and KcsA are compared in a Molecular Dynamics (MD) simulation study. The location and motion of the potassium ions observed in the simulations are compared to those in the X-ray structures and previous simulations. In our simulations several of the crystallography resolved ion sites in KirBac1.1 are occupied by ions. In addition to this, two in KirBac1.1 unresolved sites where occupied by ions at sites that are in close correspondence to sites found in KcsA. There is every reason to believe that the conserved alignment of the selectivity filter in the potassium ion channel family corresponds to a very similar mechanism for ion transport across the filter. The gate residues, Phe146 in KirBac1.1 and Ala111 in KcsA acted in the simulations as effective barriers which never were passed by ions nor water molecules.     

Efficacy of external tetraethylammonium block of the KcsA potassium channel: Molecular and Brownian dynamics studies

Blockade of the KcsA potassium channel by externally applied tetraethylammonium is investigated using molecular dynamics calculations and Brownian dynamics simulations. In KcsA, the aromatic rings of four tyrosine residues located just external to the selectivity filter create an attractive energy well or a binding cage for a tetraethylammonium molecule. We first investigate the effects of re-orienting the four tyrosine residues such that the centers of the aromatic rings face the tetraethylammonium molecule directly. Then, we systematically move the residues inward in both orientations so that the radius of the binding cage formed by them becomes smaller. For each configuration, we construct a one-dimensional free energy profile by bringing in a tetraethylammonium molecule from the external reservoir toward the selectivity filter. The free energy profile is then converted to a one-dimensional potential energy profile, taking the available space between the tyrosine residues and the tetraethylammonium molecule into account. Incorporating this potential energy profile into the Brownian dynamics algorithm, we determine the conductance properties of the channel under various conditions, construct the current-tetraethylammonium-concentration curve and compare it with the experimentally determined inhibitory constant k_i for externally applied tetraethylammonium. We show that the experimentally determined binding affinity for externally applied tetraethylammonium can be replicated when each of the four tyrosine residues is moved inward by about 0.7 Å, irrespective of orientation of their aromatic rings.     

Electrostatic state of the cytoplasmic domain influences inactivation at the selectivity filter of the KcsA potassium channel

KcsA is a proton-activated K⁺ channel that is regulated at two gates: an activation gate located in the inner entrance of the pore and an inactivation gate at the selectivity filter. Previously, we revealed that the cytoplasmic domain (CPD) of KcsA senses proton and that electrostatic changes of the CPD influences the opening and closing of the activation gate. However, our previous studies did not reveal the effect of CPD on the inactivation gate because we used a non-inactivating mutant (E71A). In the present study, we used mutants that did not harbor the E71A mutation, and showed that the electrostatic state of the CPD influences the inactivation gate. Three novel CPD mutants were generated in which some negatively charged amino acids were replaced with neutral amino acids. These CPD mutants conducted K⁺, but showed various inactivation properties. Mutants carrying the D149N mutation showed high open probability and slow inactivation, whereas those without the D149N mutation showed low open probability and fast inactivation, similar to wild-type KcsA. In addition, mutants with D149N showed poor K⁺ selectivity, and permitted Na⁺ to flow. These results indicated that electrostatic changes in the CPD by D149N mutation triggered the loss of fast inactivation and changes in the conformation of selectivity filter. Additionally, the loss of fast inactivation induced by D149N was reversed by R153A mutation, suggesting that not only the electrostatic state of D149, but also that of R153 affects inactivation.     

Functional evidence for a supramolecular structure for the Streptomyces lividans potassium channel KcsA

Here we present functional evidence for involvement of poly-(R)-3-hydroxybutyrate (PHB) and inorganic polyphosphate (polyP) in ion conduction and selection at the intracellular side of the Streptomyces lividans potassium channel, KcsA. At < or = 25 degrees C, KcsA forms channels in planar bilayers that display signal characteristics of PHB/polyP channels at the intracellular side; i.e., a preference for divalent Mg(2+) cations at pH 7.2, and a preference for monovalent K+ cations at pH 6.8. Between 25 and 26 degrees C, KcsA undergoes a transition to a new conformation in which the channel exhibits high selectivity for K+, regardless of solution pH. This suggests that basic residues of the C-terminal polypeptides have moved closer to the polyP end unit, reducing its negative charge. The data support a supramolecular structure for KcsA in which influx of ions is prevented by the selectivity pore, whereas efflux of K+ is governed by a conductive core of PHB/polyP in partnership with the C-terminal polypeptide strands.     

Modelling the effect of a GHz electric field on the dynamics of K+ ions in KcsA potassium channel

The dynamics of potassium ions in a KcsA channel, located within a stochastically fluctuating medium, is modelled via the application of the molecular dynamics simulation method. We investigate the effect of presence and absence of an applied electric field, either constant or periodic, on the dynamics of the channel. It is found that the ions undergo a hopping motion when the channel is exposed to a constant electric field of strength 0.03 V/nm. Furthermore, an alternating electric field in the GHz range, normally present in the daily environment and encountered by most biological systems, is applied to the channel, showing that in this frequency range, the rigidity of the atomic bonds of the filter is increased, which in turn disturbs the ionic passage rate through the filter. Consequently, in this frequency range, the application of electric fields may affect the function of such channels.     

Computer simulation of the KvAP voltage-gated potassium channel: steered molecular dynamics of the voltage sensor

The recent crystal structures of the voltage-gated potassium channel KvAP and its isolated voltage-sensing 'paddle' (composed of segments S1-S4) challenge existing models of voltage gating and raise a number of questions about the structure of the physiologically relevant state. We investigate a possible gating mechanism based on the crystal structures in a 10 ns steered molecular dynamics simulation of KvAP in a membrane-mimetic octane layer. The structure of the full KvAP protein has been modified by restraining the S2-S4 domain to the conformation of the isolated high-resolution paddle structure. After an initial relaxation, the paddle tips are pulled through the membrane from the intracellular to the extracellular side, corresponding to a putative change from closed to open. We describe the effect of this large-scale motion on the central pore domain, which remains largely unchanged, on the protein hydrogen-bonding network and on solvent. We analyze the motion of the S3b-S4 portion of the protein and propose a possible coupling mechanism between the paddle motion and the opening of the channel. Interactions between the arginine residues in S4, solvent and chloride ions are likely to play a role in the gating charge.     

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

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

Construction of a cyclic nucleotide-gated KcsA K+ channel

The ability of an ion channel to open in response to a defined stimulus is central to its function. In ligand-gated channels, pore opening is conferred through transduction of a conformational change in a gating domain to the helices of the pore. Here, we present the construction of a designed cyclic nucleotide-gated (CNG) channel, named KcsA-CNG, by addition of a prokaryotic cyclic nucleotide-binding domain to a KcsA-derived K+ channel. This channel is functional in lipid bilayers at physiological pH and has the combined properties of both of its parent-derived components. It conducts K+ and is blocked by the K+ channel inhibitors Na+ and agitoxin-2. Channel open times are increased by about two orders of magnitude compared to wild-type KcsA. The average number of open channels increases by approximately 50% upon addition of cAMP. Although the absolute open probabilities are somewhat variable from one channel to the next, the property of cyclic nucleotide sensitivity is very reproducible. An apparent Kd value of approximately 90 nM was estimated. The successful construction of a cyclic nucleotide-gated KcsA K+ channel suggests that it should be possible to produce channels that will respond to novel ligands.     

Determinants of frequency-dependent regulation of Kv1.2-containing potassium channels

Voltage-gated potassium channels are important regulators of electrical excitation in many tissues, with Kv1.2 standing out as an essential contributor in the CNS. Genetic deletion of Kv1.2 invariably leads to early lethality in mice. In humans, mutations affecting Kv1.2 function are linked to epileptic encephalopathy and movement disorders. We have demonstrated that Kv1.2 is subject to a unique regulatory mechanism in which repetitive stimulation leads to dramatic potentiation of current. In this study, we explore the properties and molecular determinants of this use-dependent potentiation/activation. First, we examine how alterations in duty cycle (depolarization and repolarization/recovery times) affect the onset and extent of use-dependent activation. Also, we use trains of repetitive depolarizations to test the effects of a variety of Thr252 (S2-S3 linker) mutations on use-dependent activation. Substitutions of Thr with some sterically similar amino acids (Ser, Val, and Met, but not Cys) retain use-dependent activation, while bulky or charged amino acid substitutions eliminate use-dependence. Introduction of Thr at the equivalent position in other Kv1 channels (1.1, 1.3, 1.4), was not sufficient to transfer the phenotype. We hypothesize that use-dependent activation of Kv1.2 channels is mediated by an extrinsic regulator that binds preferentially to the channel closed state, with Thr252 being necessary but not sufficient for this interaction to alter channel function. These findings extend the conclusions of our recent demonstration of use-dependent activation of Kv1.2-containing channels in hippocampal neurons, by adding new details about the molecular mechanism underlying this effect.     

Quantum mechanical calculations of charge effects on gating the KcsA channel

A series of ab initio (density functional) calculations were carried out on side chains of a set of amino acids, plus water, from the (intracellular) gating region of the KcsA K⁺ channel. Their atomic coordinates, except hydrogen, are known from X-ray structures [D.A. Doyle, J.M. Cabral, R.A. Pfuetzner, A. Kuo, J.M. Gulbis, S.L. Cohen, B.T. Chait, R. MacKinnon, The structure of the potassium channel: molecular basis of K⁺ conduction and selectivity, Science 280 (1998) 69-77; R. MacKinnon, S.L. Cohen, A. Kuo, A. Lee, B.T. Chait, Structural conservation in prokaryotic and eukaryotic potassium channels, Science 280 (1998) 106-109; Y. Jiang, A. Lee, J. Chen, M. Cadene, B.T. Chait, R. MacKinnon, The open pore conformation of potassium channels. Nature 417 (2001) 523-526], as are the coordinates of some water oxygen atoms. The 1k4c structure is used for the starting coordinates. Quantum mechanical optimization, in spite of the starting configuration, places the atoms in positions much closer to the 1j95, more tightly closed, configuration. This state shows four water molecules forming a “basket” under the Q119 side chains, blocking the channel. When a hydrated K⁺ approaches this “basket”, the optimized system shows a strong set of hydrogen bonds with the K⁺ at defined positions, preventing further approach of the K⁺ to the basket. This optimized structure with hydrated K⁺ added shows an ice-like 12 molecule nanocrystal of water. If the water molecules exchange, unless they do it as a group, the channel will remain blocked. The “basket” itself appears to be very stable, although it is possible that the K⁺ with its hydrating water molecules may be more mobile, capable of withdrawing from the gate. It is also not surprising that water essentially freezes, or forms a kind of glue, in a nanometer space; this agrees with experimental results on a rather different, but similarly sized (nm dimensions) system [K.B. Jinesh, J.W.M. Frenken, Capillary condensation in atomic scale friction: how water acts like a glue, Phys. Rev. Lett. 96 (2006) 166103/1-4]. It also agrees qualitatively with simulations on channels [A. Anishkin, S. Sukharev, Water dynamics and dewetting transitions in the small mechanosensitive channel MscS, Biophys. J. 86 (2004) 2883-2895; O. Beckstein, M.S.P. Sansom, Liquid-vapor oscillations of water in hydrophobic nanopores, Proc. Natl Acad. Sci. U. S. A. 100 (2003) 7063-7068] and on featureless channel-like systems [J. Lu, M.E. Green, Simulation of water in a pore with charges: application to a gating mechanism for ion channels, Prog. Colloid Polym. Sci. 103 (1997) 121-129], in that it forms a boundary on water that is not obvious from the liquid state. The idea that a structure is stable, even if individual molecules exchange, is well known, for example from the hydration shell of ions. We show that when charges are added in the form of protons to the domains (one proton per domain), the optimized structure is open. No stable water hydrogen bonds hold it together; an opening of 11.0 Å appears, measured diagonally between non-neighboring domains as glutamine 119 carbonyl O-O distance. This is comparable to the opening in the MthK potassium channel structure that is generally agreed to be open. The appearance of the opening is in rather good agreement with that found by Perozo and coworkers. In contrast, in the uncharged structure this diagonal distance is 6.5 Å, and the water “basket” constricts the uncharged opening still further, with the ice-like structure that couples the K⁺ ion to the gating region freezing the entrance to the channel. Comparison with our earlier model for voltage gated channels suggests that a similar mechanism may apply in those channels.     

Effect of external electric fields on the potential energy profile of K+ ions in selective filter of the KcsA potassium channel

In this study, the potential energy profile of potassium ions in the selective filter part of a KcsA channel was investigated via the application of the molecular simulation method. For this purpose, using the molecular dynamics simulation, the effect of an applied electric field, either constant or oscillating, was studied on the dynamics of K ions in the filter. It was found that when the channel is exposed to a constant electric field of strength 0.03 V/nm, the ions experience a hopping motion. Furthermore, it was shown that the application of oscillating electric fields of 1 and 2.5 GHz, can increase the rigidity of the filter atomic bonds. By computing the potential energy of K ion in the filter, it was shown that the depth of the potential wells, corresponding to the filter sites, increased when an alternative field was applied. Therefore, exposing the channel to the GHz oscillating electric field could disturb the passing rate of ions through the filter, which in turn may affect the operation of these kinds of channels.     

Uncoupling charge movement from channel opening in voltage-gated potassium channels by ruthenium complexes

The Kv2.1 channel generates a delayed-rectifier current in neurons and is responsible for modulation of neuronal spike frequency and membrane repolarization in pancreatic β-cells and cardiomyocytes. As with other tetrameric voltage-activated K(+)-channels, it has been proposed that each of the four Kv2.1 voltage-sensing domains activates independently upon depolarization, leading to a final concerted transition that causes channel opening. The mechanism by which voltage-sensor activation is coupled to the gating of the pore is still not understood. Here we show that the carbon-monoxide releasing molecule 2 (CORM-2) is an allosteric inhibitor of the Kv2.1 channel and that its inhibitory properties derive from the CORM-2 ability to largely reduce the voltage dependence of the opening transition, uncoupling voltage-sensor activation from the concerted opening transition. We additionally demonstrate that CORM-2 modulates Shaker K(+)-channels in a similar manner. Our data suggest that the mechanism of inhibition by CORM-2 may be common to voltage-activated channels and that this compound should be a useful tool for understanding the mechanisms of electromechanical coupling.     

Investigating the size and dynamics of voltage-gated sodium channel fenestrations: A molecular dynamics study

Eukaryotic voltage-gated sodium channels (VGSCs) are essential for the initiation and propagation of action potentials in electrically excitable cells, and are important pharmaceutical targets for the treatment of neurological disorders such as epilepsy, cardiac arrhythmias, and chronic pain. Evidence suggests that small, hydrophobic, VGSC-blocking drugs can gain access to binding residues within the central cavity of these channels by passing through lateral, lipid-filled “fenestrations” which run between the exterior of the protein and its central pore. Here, we use molecular dynamics simulations to investigate how the size and shape of fenestrations change over time in several bacterial VGSC models and a homology model of Nav1.4. We show that over the course of the simulations, the size of the fenestrations is primarily influenced by rapid protein motions, such as amino acid side-chain rotation, and highlight that differences between fenestration bottleneck-contributing residues are the primary cause of variations in fenestration size between the 6 bacterial models. In the eukaryotic channel model, 2 fenestrations are wide, but 2 are narrow due to differences in the amino acid sequence in the 4 domains. Lipid molecules are found to influence the size of the fenestrations by protruding acyl chains into the fenestrations and displacing amino acid side-chains. Together, the results suggest that fenestrations provide viable pathways for small, flexible, hydrophobic drugs.     

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.     

Conducting-state properties of the KcsA potassium channel from molecular and Brownian dynamics simulations.

The mechanisms underlying transport of ions across the potassium channel are examined using electrostatic calculations and three-dimensional Brownian dynamics simulations. We first build open-state configurations of the channel with molecular dynamics simulations, by pulling the transmembrane helices outward until the channel attains the desired interior radius. To gain insights into ion permeation, we construct potential energy profiles experienced by an ion traversing the channel in the presence of other resident ions. These profiles reveal that in the absence of an applied field the channel accommodates three potassium ions in a stable equilibrium, two in the selectivity filter and one in the central cavity. In the presence of a driving potential, this three-ion state becomes unstable, and ion permeation across the channel is observed. These qualitative explanations are confirmed by the results of three-dimensional Brownian dynamics simulations. We find that the channel conducts when the ionizable residues near the extracellular entrance are fully charged and those near the intracellular side are partially charged. The conductance increases steeply as the radius of the intracellular mouth of the channel is increased from 2 A to 5 A. Our simulation results reproduce several experimental observations, including the current-voltage curves, conductance-concentration relationships, and outward rectification of currents.     

A new gating site in human aquaporin-4: Insights from molecular dynamics simulations

Aquaporin-4 (AQP4) is the predominant water channel in different organs and tissues. An alteration of its physiological functioning is responsible for several disorders of water regulation and, thus, is considered an attractive target with a promising therapeutic and diagnostic potential. Molecular dynamics (MD) simulations performed on the AQP4 tetramer embedded in a bilayer of lipid molecules allowed us to analyze the role of spontaneous fluctuations occurring inside the pore. Following the approach by Hashido et al. [Hashido M, Kidera A, Ikeguchi M (2007) Biophys J 93: 373–385], our analysis on 200 ns trajectory discloses three domains inside the pore as key elements for water permeation. Herein, we describe the gating mechanism associated with the well-known selectivity filter on the extracellular side of the pore and the crucial regulation ensured by the NPA motifs (asparagine, proline, alanine). Notably, on the cytoplasmic side, we find a putative gate formed by two residues, namely, a cysteine belonging to the loop D (C178) and a histidine from loop B (H95). We observed that the spontaneous reorientation of the imidazole ring of H95 acts as a molecular switch enabling H-bond interaction with C178. The occurrence of such local interaction seems to be responsible for the narrowing of the pore and thus of a remarkable decrease in water flux rate. Our results are in agreement with recent experimental observations and may represent a promising starting point to pave the way for the discovery of chemical modulators of AQP4 water permeability.     

Dissipative particle dynamics simulation of a gold nanoparticle system

Dissipative particle dynamics (DPD) was carried out to study systems containing gold atoms, organic ether (oligohydroquinonyl ether terminated with a thiol group) and organic solvents. The components in the simulated system are very different in size and chemical nature. Our simulation showed that the reproduction of the macroscopic experimental phase separation, properly dividing the polymeric molecule into beads, selecting the size of gold bead, and choosing the appropriate interaction parameters between beads are crucial. In addition, the solvent effect was the dominant factor for the formation of spherical aggregates of Au atoms and organic ether molecules. We report the interaction strengths between the solvent and gold clusters. Our work has demonstrated that DPD methods can be applied to the study of complex meso-scale systems.     

Local and global structure of the monomeric subunit of the potassium channel KcsA probed by NMR

KcsA is a homotetrameric 68-kDa membrane-associated potassium channel which selectively gates the flux of potassium ions across the membrane. The channel is known to undergo a pH-dependent open-to-closed transition. Here we describe an NMR study of the monomeric subunit of the channel (KcsA^M), solubilized in SDS micelles. Chemical shift, solvent exchange, backbone ¹⁵N relaxation and residual dipolar coupling (RDC) data show the TM1 helix to remain intact, but the TM2 helix contains a distinct kink, which is subject to concentration-independent but pH-dependent conformational exchange on a microsecond time scale. The kink region, centered at G99, was previously implicated in the gating of the tetrameric KcsA channel. An RDC-based model of KcsA^M at acidic pH orients TM1 and the two helical segments of the kinked TM2 in a configuration reminiscent of the open conformation of the channel. Thus, the transition between states appears to be an inherent capability of the monomer, with the tetrameric assembly exerting a modulatory effect upon the transition which gives the channel its physiological gating profile.     

Simulations of ion permeation through a potassium channel: molecular dynamics of KcsA in a phospholipid bilayer

Potassium channels enable K(+) ions to move passively across biological membranes. Multiple nanosecond-duration molecular dynamics simulations (total simulation time 5 ns) of a bacterial potassium channel (KcsA) embedded in a phospholipid bilayer reveal motions of ions, water, and protein. Comparison of simulations with and without K(+) ions indicate that the absence of ions destabilizes the structure of the selectivity filter. Within the selectivity filter, K(+) ions interact with the backbone (carbonyl) oxygens, and with the side-chain oxygen of T75. Concerted single-file motions of water molecules and K(+) ions within the selectivity filter of the channel occur on a 100-ps time scale. In a simulation with three K(+) ions (initially two in the filter and one in the cavity), the ion within the central cavity leaves the channel via its intracellular mouth after approximately 900 ps; within the cavity this ion interacts with the Ogamma atoms of two T107 side chains, revealing a favorable site within the otherwise hydrophobically lined cavity. Exit of this ion from the channel is enabled by a transient increase in the diameter of the intracellular mouth. Such "breathing" motions may form the molecular basis of channel gating.