Dataset Information

Permeation and block of the Kv1.2 channel examined using brownian and molecular dynamics.

ABSTRACT: Using both Brownian and molecular dynamics, we replicate many of the salient features of Kv1.2, including the current-voltage-concentration profiles and the binding affinity and binding mechanisms of charybdotoxin, a scorpion venom. We also elucidate how structural differences in the inner vestibule can give rise to significant differences in its permeation characteristics. Current-voltage-concentration profiles are constructed using Brownian dynamics simulations, based on the crystal structure 2A79. The results are compatible with experimental data, showing similar conductance, rectification, and saturation with current. Unlike KcsA, for example, the inner pore of Kv1.2 is mainly hydrophobic and neutral, and to explore the consequences of this, we investigate the effect of mutating neutral proline residues at the mouth of the inner vestibule to charged aspartate residues. We find an increased conductance, less inward rectification, and quicker saturation of the current-voltage profile. Our simulations use modifications to our Brownian dynamics program that extend the range of channels that can be usefully modeled. Using molecular dynamics, we investigate the binding of the charybdotoxin scorpion venom to the outer vestibule of the channel. A potential of mean force is derived using umbrella sampling, giving a dissociation constant within a factor of ?2 to experimentally derived constants. The residues involved in the toxin binding are in agreement with experimental mutagenesis studies. We thus show that the experimental observations on the voltage-gated channel, including the toxin-channel interaction, can reliably be replicated by using the two widely used computational tools.

SUBMITTER: Gordon D

PROVIDER: S-EPMC3297784 | biostudies-literature | 2011 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Permeation and block of the Kv1.2 channel examined using brownian and molecular dynamics.

Gordon Dan D Chung Shin-Ho SH

Biophysical journal 20111201 11

Using both Brownian and molecular dynamics, we replicate many of the salient features of Kv1.2, including the current-voltage-concentration profiles and the binding affinity and binding mechanisms of charybdotoxin, a scorpion venom. We also elucidate how structural differences in the inner vestibule can give rise to significant differences in its permeation characteristics. Current-voltage-concentration profiles are constructed using Brownian dynamics simulations, based on the crystal structure ...[more]

PMID: 22261055

Similar Datasets

Project description:Identification of the molecular interaction governing ion conduction through biological pores is one of the most important goals of modern electrophysiology. Grand canonical Monte Carlo Brownian dynamics (GCMC/BD) and three-dimensional Poisson-Nernst-Plank (3d-PNP) electrodiffusion algorithms offer powerful and general approaches to study of ion permeation through wide molecular pores. A detailed analysis of ion flows through the staphylococcal alpha-hemolysin channel based on series of simulations at different concentrations and transmembrane potentials is presented. The position-dependent diffusion coefficient is approximated on the basis of a hydrodynamic model. The channel conductance calculated by GCMC/BD is approximately 10% higher than (electrophysiologically measured) experimental values, whereas results from 3d-PNP are always 30-50% larger. Both methods are able to capture all important electrostatic interactions in equilibrium conditions. The asymmetric conductance upon the polarity of the transmembrane potential observed experimentally is reproduced by GCMC/BD and 3d-PNP. The separation of geometrical and energetic influence of the channel on ion conduction reveals that such asymmetries arise from the permanent charge distribution inside the pore. The major determinant of the asymmetry is unbalanced charge in the triad of polar residues D127, D128, and K131. The GCMC/BD or 3d-PNP calculations reproduce also experimental reversal potentials and permeability rations in asymmetric ionic solutions. The weak anionic selectivity of the channel results from the presence of the salt bridge between E111 and K147 in the constriction zone. The calculations also reproduce the experimentally derived dependence of the reversible potential to the direction of the salt gradient. The origin of such effect arises from the asymmetrical distribution of energetic barriers along the channel axis, which modulates the preferential ion passage in different directions.

Project description:Structure and dynamics of voltage-gated ion channels, in particular the motion of the S4 helix, is a highly interesting and hotly debated topic in current membrane protein research. It has critical implications for insertion and stabilization of membrane proteins as well as for finding how transitions occur in membrane proteins-not to mention numerous applications in drug design. Here, we present a full 1 micros atomic-detail molecular dynamics simulation of an integral Kv1.2 ion channel, comprising 120,000 atoms. By applying 0.052 V/nm of hyperpolarization, we observe structural rearrangements, including up to 120 degrees rotation of the S4 segment, changes in hydrogen-bonding patterns, but only low amounts of translation. A smaller rotation ( approximately 35 degrees ) of the extracellular end of all S4 segments is present also in a reference 0.5 micros simulation without applied field, which indicates that the crystal structure might be slightly different from the natural state of the voltage sensor. The conformation change upon hyperpolarization is closely coupled to an increase in 3(10) helix contents in S4, starting from the intracellular side. This could support a model for transition from the crystal structure where the hyperpolarization destabilizes S4-lipid hydrogen bonds, which leads to the helix rotating to keep the arginine side chains away from the hydrophobic phase, and the driving force for final relaxation by downward translation is partly entropic, which would explain the slow process. The coordinates of the transmembrane part of the simulated channel actually stay closer to the recently determined higher-resolution Kv1.2 chimera channel than the starting structure for the entire second half of the simulation (0.5-1 micros). Together with lipids binding in matching positions and significant thinning of the membrane also observed in experiments, this provides additional support for the predictive power of microsecond-scale membrane protein simulations.

Dataset Information

Permeation and block of the Kv1.2 channel examined using brownian and molecular dynamics.

Publications

Permeation and block of the Kv1.2 channel examined using brownian and molecular dynamics.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets