Dataset Information

Molecular Dynamics Simulations of KirBac1.1 Mutants Reveal Global Gating Changes of Kir Channels.

ABSTRACT: Prokaryotic inwardly rectifying (KirBac) potassium channels are homologous to mammalian Kir channels. Their activity is controlled by dynamical conformational changes that regulate ion flow through a central pore. Understanding the dynamical rearrangements of Kir channels during gating requires high-resolution structure information from channels crystallized in different conformations and insight into the transition steps, which are difficult to access experimentally. In this study, we use MD simulations on wild type KirBac1.1 and an activatory mutant to investigate activation gating of KirBac channels. Full atomistic MD simulations revealed that introducing glutamate in position 143 causes significant widening at the helix bundle crossing gate, enabling water flux into the cavity. Further, global rearrangements including a twisting motion as well as local rearrangements at the subunit interface in the cytoplasmic domain were observed. These structural rearrangements are similar to recently reported KirBac3.1 crystal structures in closed and open conformation, suggesting that our simulations capture major conformational changes during KirBac1.1 opening. In addition, an important role of protein-lipid interactions during gating was observed. Slide-helix and C-linker interactions with lipids were strengthened during activation gating.

SUBMITTER: Linder T

PROVIDER: S-EPMC4415035 | biostudies-literature | 2015 Apr

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Molecular Dynamics Simulations of KirBac1.1 Mutants Reveal Global Gating Changes of Kir Channels.

Linder Tobias T Wang Shizhen S Zangerl-Plessl Eva-Maria EM Nichols Colin G CG Stary-Weinzinger Anna A

Journal of chemical information and modeling 20150403 4

Prokaryotic inwardly rectifying (KirBac) potassium channels are homologous to mammalian Kir channels. Their activity is controlled by dynamical conformational changes that regulate ion flow through a central pore. Understanding the dynamical rearrangements of Kir channels during gating requires high-resolution structure information from channels crystallized in different conformations and insight into the transition steps, which are difficult to access experimentally. In this study, we use MD si ...[more]

PMID: 25794351

Similar Datasets

Project description:Inward rectifier (Kir) potassium channels are characterized by two transmembrane helices per subunit, plus an intracellular C-terminal domain that controls channel gating in response to changes in concentration of various ligands. Based on the crystal structure of the tetrameric C-terminal domain of Kir3.1, it is possible to build a homology model of the ATP-binding C-terminal domain of Kir6.2. Molecular dynamics simulations have been used to probe the dynamics of Kir C-terminal domains and to explore the relationship between their dynamics and possible mechanisms of channel gating. Multiple simulations, each of 10 ns duration, have been performed for Kir3.1 (crystal structure) and Kir6.2 (homology model), in both their monomeric and tetrameric forms. The Kir6.2 simulations were performed with and without bound ATP. The results of the simulations reveal comparable conformational stability for the crystal structure and the homology model. There is some decrease in conformational flexibility when comparing the monomers with the tetramers, corresponding mainly to the subunit interfaces in the tetramer. The beta-phosphate of ATP interacts with the side chain of K185 in the Kir6.2 model and simulations. The flexibility of the Kir6.2 tetramer is not changed greatly by the presence of bound ATP, other than in two loop regions. Principal components analysis of the simulated dynamics suggests loss of symmetry in both the Kir3.1 and Kir6.2 tetramers, consistent with "dimer-of-dimers" motion of subunits in C-terminal domains of the corresponding Kir channels. This is suggestive of a gating model in which a transition between exact tetrameric symmetry and dimer-of-dimers symmetry is associated with a change in transmembrane helix packing coupled to gating of the channel. Dimer-of-dimers motion of the C-terminal domain tetramer is also supported by coarse-grained (anisotropic network model) calculations. It is of interest that loss of exact rotational symmetry has also been suggested to play a role in gating in the bacterial Kir homolog, KirBac1.1, and in the nicotinic acetylcholine receptor channel.

Project description:The recently determined crystal structure of a chimeric Kv1.2-Kv2.1 Kv channel at 2.4 A resolution motivated this molecular-dynamics simulation study of the chimeric channel and its mutants embedded in a DPPC membrane. For the channel protein, we used two types of C-terminus: E+ and Eo. E+ contains, and Eo lacks, the EGEE residue quartet located distal to the S6 helix. For both E+ and Eo, the following trend was observed: When S4 helices were restrained at the same position as in the x-ray structure (S4high), the S6 gate remained open for 12 ns. The results were similar when the S4 helices were pulled downward 7 A (S4low). However, S4middle (or S4low) facilitated the S6 gate-narrowing for the following mutated channels (shown in order of increasing effect): 1), E395W; 2), E395W-F401A-F402A; and 3), E395W-F401A-F402A-V478W. The amino acid numbering system is that used for the Shaker channel. Even though all four subunits were set at S4low, S6 gate-narrowing was often brought about by movements of only two opposing S6 helices toward the central axis of the pore, resulting in a twofold symmetry-like structure. A free-energy profile analysis over the ion conduction pathway shows that the two opposing S6 helices whose peptide backbones are approximately 10.4 A distant from each other lead to an energetic barrier of approximately 25 kJ/mol. S6 movement was coupled with translocation of the S4-S5 linker toward the central axis of the same subunit, and the coupling was mediated by salt bridges formed between the inner (intracellular side) end of S4 and that of S6. Simulations in which S4 of only one subunit was pulled down to S4low showed that a weak intersubunit coordination is present for S5 movement, whereas the coupling between the S4-S5 linker and S6 is largely an intrasubunit one. In general, whereas subunit-based behavior appears to be dominant and to permit heteromeric conformations of the pore domain, direct intersubunit coupling of S5 or S6 is weak. Therefore, the "concerted transition" of the pore domain that has been predicted based on electrophysiological analyses is likely to be mediated mainly by the dual effects of S4 and the S4-S5 linker; these segments of one subunit can interact with both S5 of the same subunit and that of the adjacent subunit.

Dataset Information

Molecular Dynamics Simulations of KirBac1.1 Mutants Reveal Global Gating Changes of Kir Channels.

Publications

Molecular Dynamics Simulations of KirBac1.1 Mutants Reveal Global Gating Changes of Kir Channels.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure