The Voltage-Dependent Deactivation of the KvAP Channel Involves the Breakage of Its S4 Helix

Voltage-gated potassium channels (Kv) allow ion permeation upon changes of the membrane electrostatic potential (Vm). Each subunit of these tetrameric channels is composed of six transmembrane helices, of which the anti-parallel helix bundle S1-S4 constitutes the voltage-sensor domain (VSD) and S5-S6 forms the pore domain. Here, using 82 molecular dynamics (MD) simulations involving 266 replicated VSDs, we report novel responses of the archaebacterial potassium channel KvAP to membrane polarization. We show that the S4 α-helix, which is straight in the experimental crystal structure solved under depolarized conditions (Vm ∼ 0), breaks into two segments when the cell membrane is hyperpolarized (Vm < 0), and reversibly forms a single straight helix following depolarization (Vm = 0). The outermost segment of S4 translates along the normal to the membrane, bringing new perspective to previously paradoxical accessibility experiments that were initially thought to imply the displacement of the whole VSD across the membrane. The novel model is applied through steered and unbiased MD simulations to the recently solved whole structure of KvAP. The simulations show that the resting state involves a re-orientation of the S5 α-helix by ∼ 5-6 degrees in respect to the normal of the bilayer, which could result in the constriction and closure of the selectivity filter. Our findings support the idea that the breakage of S4 under (hyper)polarization is a general feature of Kv channels with a non-swapped topology

The Voltage-Dependent Deactivation of the KvAP Channel Involves the Breakage of Its S4 Helix

We show that the response of the KvAP voltage-sensor domain to membrane polarization consists in two co-occurring conformational changes: the rupture of the Asp62-Arg133 salt-bridge and the breakage of the S4 helix. This finding is consistent with the hypothesis that the breakage of S4 is a general feature of non-swapped Kv channels. Upon polarization, the N-ter segment of the S4 helix translates toward the intracellular side of the membrane, which allows to solve previously paradoxical avidin accessibility measurements. Simulations involving the recently solved whole structure of KvAP show how these conformational changes induce a reorientation of the S5 and pore helices, putatively leading to the constriction and closure of the selectivity filter. This study also links mechanistic insights in prokaryotic potassium channel membrane potential sensing with the perspective of finding new antibiotic targets

Unidirectional Transport Mechanism in an ATP Dependent Exporter

ATP-binding cassette (ABC) transporters use the energy of ATP binding and hydrolysis to move a large variety of compounds across biological membranes. P-glycoprotein, involved in multidrug resistance, is the most investigated eukaryotic family member. Although a large number of biochemical and structural approaches have provided important information, the conformational dynamics underlying the coupling between ATP binding/hydrolysis and allocrite transport remains elusive. To tackle this issue, we performed molecular dynamic simulations for different nucleotide occupancy states of Sav1866, a prokaryotic P-glycoprotein homologue. The simulations reveal an outward-closed conformation of the transmembrane domain that is stabilized by the binding of two ATP molecules. The hydrolysis of a single ATP leads the X-loop, a key motif of the ATP binding cassette, to interfere with the transmembrane domain and favor its outward-open conformation. Our findings provide a structural basis for the unidirectionality of transport in ABC exporters and suggest a ratio of one ATP hydrolyzed per transport cycle

Generalized solvent boundary potential for computer simulations

Copyright 2001 American Institute of Physics. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics. The following article appeared in The Journal of Chemical Physics and may be found at http://dx.doi.org/10.1063/1.1336570.A general approach has been developed to allow accurate simulations of a small region part of a large macromolecular system while incorporating the influence of the remaining distant atoms with an effective boundary potential. The method is called the Generalized Solvent Boundary Potential (GSBP). By representing the surrounding solvent as a continuum dielectric, both the solvent-shielded static field from the distant atoms of the macromolecule and the reaction field from the dielectricsolvent acting on the atoms in the region of interest are included. The static field is calculated once, using the finite-difference Poisson–Boltzmann (PB) equation, and the result is stored on a discrete grid for efficient simulations. The solventreaction field is developed using a basis-set expansion whose coefficients correspond to generalized electrostatic multipoles. A matrix representing the reaction field Green’s function between those generalized multipoles is calculated only once using the PB equation and stored for efficient simulations. In the present work, the formalism is applied to both spherical and orthorhombic simulation regions for which orthonormal basis-sets exist based on spherical harmonics or cartesian Legendre polynomials. The GSBP method is also tested and illustrated with simple model systems and two detailed atomic systems: the active site region of aspartyl-tRNA synthetase (spherical region) and the interior of the KcsA potassium channel (orthorhombic region). Comparison with numerical finite-difference PB calculations shows that GSBP can accurately describe all long-range electrostatic interactions and remain computationally inexpensive

Calcium regulates acid-sensing ion channel 3 activation by competing with protons in the channel pore and at an allosteric binding site

The extracellular Ca2+ concentration changes locally under certain physiological and pathological conditions. Such variations affect the function of ion channels of the nervous system, and consequently also neuronal signalling. We investigated here the mechanisms by which Ca2+ controls the activity of acid-sensing ion channel (ASIC) 3. ASICs are neuronal, H+-gated Na+ channels involved in several physiological and pathological processes, including the expression of fear, learning, pain sensation and neurodegeneration after ischemic stroke. It was previously shown that Ca2+ negatively modulates the ASIC pH dependence. While protons are default activators of ASIC3, this channel can also be activated at pH7.4 by removal of the extracellular Ca2+. Two previous studies concluded that low pH opens ASIC3 by displacing Ca2+ ions that block the channel pore at physiological pH. We show here that an acidic residue, distant from the pore, controls, together with pore residues, the modulation of ASIC3 by Ca2+. Our study identifies a new regulatory site in ASIC3 and demonstrates that ASIC3 activation involves an allosteric mechanism together with Ca2+ unbinding from the channel pore. We provide a molecular analysis of a regulatory mechanism found in many ion channels

Extracellular Blockade of K+ Channels by Tea: Results from Molecular Dynamics Simulations of the Kcsa Channel

TEA is a classical blocker of K+ channels. From mutagenesis studies, it has been shown that external blockade by TEA is strongly dependent upon the presence of aromatic residue at Shaker position 449 which is located near the extracellular entrance to the pore (Heginbotham, L., and R. MacKinnon. 1992. Neuron. 8:483–491). The data suggest that TEA interacts simultaneously with the aromatic residues of the four monomers. The determination of the 3-D structure of the KcsA channel using X-ray crystallography (Doyle, D.A., J.M. Cabral, R.A. Pfuetzner, A. Kuo, J.M. Gulbis, S.L. Cohen, B.T. Chait, and R. MacKinnon. 1998. Science. 280:69–77) has raised some issues that remain currently unresolved concerning the interpretation of these observations. In particular, the center of the Tyr82 side chains in KcsA (corresponding to position 449 in Shaker) forms a square of 11.8-Å side, a distance which is too large to allow simultaneous interactions of a TEA molecule with the four aromatic side chains. In this paper, the external blockade by TEA is explored by molecular dynamics simulations of an atomic model of KcsA in an explicit phospholipid bilayer with aqueous salt solution. It is observed, in qualitative accord with the experimental results, that TEA is stable when bound to the external side of the wild-type KcsA channel (with Tyr82), but is unstable when bound to a mutant channel in which the tyrosine residue has been substituted by a threonine. The free energy profile of TEA relative to the pore is calculated using umbrella sampling simulations to characterize quantitatively the extracellular blockade. It is found, in remarkable agreement with the experiment, that the TEA is more stably bound by 2.3 kcal/mol to the channel with four tyrosine residues. In the case of the wild-type KcsA channel, TEA (which has the shape of a flattened oblate spheroid) acts as an ideal plug blocking the pore. In contrast, it is considerably more off-centered and tilted in the case of the mutant channel. The enhanced stability conferred by the tyrosine residues does not arise from Π–cation interactions, but appears to be due to differences in the hydration structure of the TEA. Finally, it is shown that the experimentally observed voltage dependence of TEA block, which is traditionally interpreted in terms of the physical position of the TEA along the axis of the pore, must arise indirectly via coupling with the ions in the pore

A microscopic view of ion conduction through the K+ channel

Recent results from x-ray crystallography and molecular dynamics free-energy simulations have revealed the existence of a number of specific cation-binding sites disposed along the narrow pore of the K(+) channel from Streptomyces lividans (KcsA), suggesting that K(+) ions might literally “hop” in single file from one binding site to the next as permeation proceeds. In support of this view, it was found that the ion configurations correspond to energy wells of similar depth and that ion translocation is opposed only by small energy barriers. Although such features of the multiion potential energy surface are certainly essential for achieving a high throughput rate, diffusional and dissipative dynamical factors must also be taken into consideration to understand how rapid conduction of K(+) is possible. To elucidate the mechanism of ion conduction, we established a framework theory enabling the direct simulation of nonequilibrium fluxes by extending the results of molecular dynamics over macroscopically long times. In good accord with experimental measurements, the simulated maximum conductance of the channel at saturating concentration is on the order of 550 and 360 pS for outward and inward ions flux, respectively, with a unidirectional flux-ratio exponent of 3. Analysis of the ion-conduction process reveals a lack of equivalence between the cation-binding sites in the selectivity filter