Dataset Information

Brownian dynamics simulations of ion transport through the VDAC.

ABSTRACT: It is important to gain a physical understanding of ion transport through the voltage-dependent anion channel (VDAC) because this channel provides primary permeation pathways for metabolites and electrolytes between the cytosol and mitochondria. We performed grand canonical Monte Carlo/Brownian dynamics (GCMC/BD) simulations to explore the ion transport properties of human VDAC isoform 1 (hVDAC1; PDB:2K4T) embedded in an implicit membrane. When the MD-derived, space-dependent diffusion constant was used in the GCMC/BD simulations, the current-voltage characteristics and ion number profiles inside the pore showed excellent agreement with those calculated from all-atom molecular-dynamics (MD) simulations, thereby validating the GCMC/BD approach. Of the 20 NMR models of hVDAC1 currently available, the third one (NMR03) best reproduces both experimental single-channel conductance and ion selectivity (i.e., the reversal potential). In addition, detailed analyses of the ion trajectories, one-dimensional multi-ion potential of mean force, and protein charge distribution reveal that electrostatic interactions play an important role in the channel structure and ion transport relationship. Finally, the GCMC/BD simulations of various mutants based on NMR03 show good agreement with experimental ion selectivity. The difference in ion selectivity between the wild-type and the mutants is the result of altered potential of mean force profiles that are dominated by the electrostatic interactions.

SUBMITTER: Lee KI

PROVIDER: S-EPMC3030170 | biostudies-literature | 2011 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Brownian dynamics simulations of ion transport through the VDAC.

Lee Kyu Il KI Rui Huan H Pastor Richard W RW Im Wonpil W

Biophysical journal 20110201 3

It is important to gain a physical understanding of ion transport through the voltage-dependent anion channel (VDAC) because this channel provides primary permeation pathways for metabolites and electrolytes between the cytosol and mitochondria. We performed grand canonical Monte Carlo/Brownian dynamics (GCMC/BD) simulations to explore the ion transport properties of human VDAC isoform 1 (hVDAC1; PDB:2K4T) embedded in an implicit membrane. When the MD-derived, space-dependent diffusion constant ...[more]

PMID: 21281575

Similar Datasets

Project description:Voltage-dependent anion channel (VDAC), the major channel of the mitochondrial outer membrane, serves as a principal pathway for ATP, ADP, and other respiratory substrates across this membrane. Using umbrella-sampling simulations, we established the thermodynamic and kinetic components governing ATP transport across the VDAC1 channel. We found that there are several low-affinity binding sites for ATP along the translocation pathway and that the main barrier for ATP transport is located around the center of the channel and is formed predominantly by residues in the N-terminus. The binding affinity of ATP to an open channel was found to be in the millimolar to micromolar range. However, we show that this weak binding increases the ATP translocation probability by about 10-fold compared with the VDAC pore in which attractive interactions were artificially removed. Recently, it was found that free dimeric tubulin induces a highly efficient, reversible blockage of VDAC reconstituted into planar lipid membranes. It was proposed that by blocking VDAC permeability for ATP/ADP and other mitochondrial respiratory substrates tubulin controls mitochondrial respiration. Using the Rosetta protein-protein docking algorithm, we established a tentative structure of the VDAC-tubulin complex. An extensive set of equilibrium and nonequilibrium (under applied electric field) molecular dynamics (MD) simulations was used to establish the conductance of the open and blocked channel. It was found that the presence of the unstructured C-terminal tail of tubulin in the VDAC pore decreases its conductance by more than 40% and switches its selectivity from anionic to cationic. The subsequent 1D potential of mean force (PMF) computations for the VDAC-tubulin complex show that the state renders ATP transport virtually impossible. A number of residues pivotal for tubulin binding to the channel were identified that help to clarify the molecular details of VDAC-tubulin interaction and to provide new insight into the mechanism of the control of mitochondria respiration by VDAC.

Project description:The association of the scorpion toxin Lq2 and a potassium ion (K(+)) channel has been studied using the Brownian dynamics (BD) simulation method. All of the 22 available structures of Lq2 in the Brookhaven Protein Data Bank (PDB) determined by NMR were considered during the simulation, which indicated that the conformation of Lq2 affects the binding between the two proteins significantly. Among the 22 structures of Lq2, only 4 structures dock in the binding site of the K(+) channel with a high probability and favorable electrostatic interactions. From the 4 candidates of the Lq2-K(+) channel binding models, we identified a good three-dimensional model of Lq2-K(+) channel complex through triplet contact analysis, electrostatic interaction energy estimation by BD simulation and structural refinement by molecular mechanics. Lq2 locates around the extracellular mouth of the K(+) channel and contacts the K(+) channel using its beta-sheet rather than its alpha-helix. Lys27, a conserved amino acid in the scorpion toxins, plugs the pore of the K(+) channel and forms three hydrogen bonds with the conserved residues Tyr78(A-C) and two hydrophobic contacts with Gly79 of the K(+) channel. In addition, eight hydrogen-bonds are formed between residues Arg25, Cys28, Lys31, Arg34 and Tyr36 of Lq2 and residues Pro55, Tyr78, Gly79, Asp80, and Tyr82 of K(+) channel. Many of them are formed by side chains of residues of Lq2 and backbone atoms of the K(+) channel. Thirteen hydrophobic contacts exist between residues Met29, Asn30, Lys31 and Tyr36 of Lq2 and residues Pro55, Ala58, Gly79, Asp80 and Tyr82 of the K(+) channel. These favorable interactions stabilize the association between the two proteins. These observations are in good agreement with the experimental results and can explain the binding phenomena between scorpion toxins and K(+) channels at the level of molecular structure. The consistency between the BD simulation and the experimental data indicates that our three-dimensional model of Lq2-K(+) channel complex is reasonable and can be used in further biological studies such as rational design of blocking agents of K(+) channels and mutagenesis in both toxins and K(+) channels.

Project description:The recognition of the scorpion toxin maurotoxin (MTX) by the voltage-gated potassium (Kv1) channels, Kv1.1, Kv1.2, and Kv1.3, has been studied by means of Brownian dynamics (BD) simulations. All of the 35 available structures of MTX in the Protein Data Bank (http://www.rcsb.org/pdb) determined by nuclear magnetic resonance were considered during the simulations, which indicated that the conformation of MTX significantly affected both the recognition and the binding between MTX and the Kv1 channels. Comparing the top five highest-frequency structures of MTX binding to the Kv1 channels, we found that the Kv1.2 channel, with the highest docking frequencies and the lowest electrostatic interaction energies, was the most favorable for MTX binding, whereas Kv1.1 was intermediate, and Kv1.3 was the least favorable one. Among the 35 structures of MTX, the 10th structure docked into the binding site of the Kv1.2 channel with the highest probability and the most favorable electrostatic interactions. From the MTX-Kv1.2 binding model, we identified the critical residues for the recognition of these two proteins through triplet contact analyses. MTX locates around the extracellular mouth of the Kv1 channels, making contacts with its beta-sheets. Lys23, a conserved amino acid in the scorpion toxins, protrudes into the pore of the Kv1.2 channel and forms two hydrogen bonds with the conserved residues Gly401(D) and Tyr400(C) and one hydrophobic contact with Gly401(C) of the Kv1.2 channel. The critical triplet contacts for recognition between MTX and the Kv1.2 channel are Lys23(MTX)-Asp402(C)(Kv1), Lys27(MTX)-Asp378(D)(Kv1), and Lys30(MTX)-Asp402(A)(Kv1). In addition, six hydrogen-bonding interactions are formed between residues Lys23, Lys27, Lys30, and Tyr32 of MTX and residues Gly401, Tyr400, Asp402, Asp378, and Thr406 of Kv1.2. Many of them are formed by side chains of residues of MTX and backbone atoms of the Kv1.2 channel. Five hydrophobic contacts exist between residues Pro20, Lys23, Lys30 and Tyr32 of MTX and residues Asp402, Val404, Gly401, and Arg377 of the Kv1.2 channel. The simulation results are in agreement with the previous molecular biology experiments and explain the binding phenomena between MTX and Kv1 channels at the molecular level. The consistency between the results of the BD simulations and the experimental data indicated that our three-dimensional model of the MTX-Kv1.2 channel complex is reasonable and can be used in additional biological studies, such as rational design of novel therapeutic agents blocking the voltage-gated channels and in mutagenesis studies in both the toxins and the Kv1 channels. In particular, both the BD simulations and the molecular mechanics refinements indicate that residue Asp378 of the Kv1.2 channel is critical for its recognition and binding functionality toward MTX. This phenomenon has not been appreciated in the previous mutagenesis experiments, indicating this might be a new clue for additional functional study of Kv1 channels.

Dataset Information

Brownian dynamics simulations of ion transport through the VDAC.

Publications

Brownian dynamics simulations of ion transport through the VDAC.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets