Molecular dynamics studies of ion permeation in VDAC.
Ontology highlight
ABSTRACT:
SUBMITTER: Rui H
PROVIDER: S-EPMC3030152 | biostudies-literature | 2011 Feb
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The voltage-dependent anion channel (VDAC) forms the major pore in the outer mitochondrial membrane. Its high conducting open state features a moderate anion selectivity. There is some evidence indicating that the electrophysiological properties of VDAC vary with the salt concentration. Using a theoretical approach the molecular basis for this concentration dependence was investigated. Molecular dynamics simulations and continuum electrostatic calculations performed on the mouse VDAC1 isoform clearly demonstrate that the distribution of fixed charges in the channel creates an electric field, which determines the anion preference of VDAC at low salt concentration. Increasing the salt concentration in the bulk results in a higher concentration of ions in the VDAC wide pore. This event induces a large electrostatic screening of the charged residues promoting a less anion selective channel. Residues that are responsible for the electrostatic pattern of the channel were identified using the molecular dynamics trajectories. Some of these residues are found to be conserved suggesting that ion permeation between different VDAC species occurs through a common mechanism. This inference is buttressed by electrophysiological experiments performed on bean VDAC32 protein akin to mouse VDAC.
| S-EPMC3229507 | biostudies-literature
Project description:In the exchange of metabolites and ions between the mitochondrion and the cytosol, the voltage-dependent anion channel (VDAC) is a key element, as it forms the major transport pathway for these compounds through the mitochondrial outer membrane. Numerous experimental studies have promoted the idea that VDAC acts as a regulator of essential mitochondrial functions. In this study, using a combination of molecular dynamics simulations, free-energy calculations, and electrophysiological measurements, we investigated the transport of ions through VDAC, with a focus on phosphate ions and metabolites. We showed that selectivity of VDAC towards small anions including monovalent phosphates arises from short-lived interactions with positively charged residues scattered throughout the pore. In dramatic contrast, permeation of divalent phosphate ions and phosphate metabolites (AMP and ATP) involves binding sites along a specific translocation pathway. This permeation mechanism offers an explanation for the decrease in VDAC conductance measured in the presence of ATP or AMP at physiological salt concentration. The binding sites occur at similar locations for the divalent phosphate ions, AMP and ATP, and contain identical basic residues. ATP features a marked affinity for a central region of the pore lined by two lysines and one arginine of the N-terminal helix. This cluster of residues together with a few other basic amino acids forms a "charged brush" which facilitates the passage of the anionic metabolites through the pore. All of this reveals that VDAC controls the transport of the inorganic phosphates and phosphate metabolites studied here through two different mechanisms.
| S-EPMC4393092 | biostudies-literature
Project description:The recent determination of cryo-EM structures of voltage-gated sodium (Nav) channels has revealed many details of these proteins. However, knowledge of ionic permeation through the Nav pore remains limited. In this work, we performed atomistic molecular dynamics (MD) simulations to study the structural features of various neuronal Nav channels based on homology modeling of the cryo-EM structure of the human Nav1.4 channel and, in addition, on the recently resolved configuration for Nav1.2. In particular, single Na+ permeation events during standard MD runs suggest that the ion resides in the inner part of the Nav selectivity filter (SF). On-the-fly free energy parametrization (OTFP) temperature-accelerated molecular dynamics (TAMD) was also used to calculate two-dimensional free energy surfaces (FESs) related to single/double Na+ translocation through the SF of the homology-based Nav1.2 model and the cryo-EM Nav1.2 structure, with different realizations of the DEKA filter domain. These additional simulations revealed distinct mechanisms for single and double Na+ permeation through the wild-type SF, which has a charged lysine in the DEKA ring. Moreover, the configurations of the ions in the SF corresponding to the metastable states of the FESs are specific for each SF motif. Overall, the description of these mechanisms gives us new insights into ion conduction in human Nav cryo-EM-based and cryo-EM configurations that could advance understanding of these systems and how they differ from potassium and bacterial Nav channels.
| S-EPMC10210251 | biostudies-literature
Project description: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.
| S-EPMC3030170 | biostudies-literature
Project description:Bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC) is activated to cation permeation upon lowering the solution pH. Its function can be modulated by anesthetic halothane. In the present work, we integrate molecular dynamics (MD) and Brownian dynamics (BD) simulations to elucidate the ion conduction, charge selectivity, and halothane modulation mechanisms in GLIC, based on recently resolved X-ray crystal structures of the open-channel GLIC. MD calculations of the potential of mean force (PMF) for a Na(+) revealed two energy barriers in the extracellular domain (R109 and K38) and at the hydrophobic gate of transmembrane domain (I233), respectively. An energy well for Na(+) was near the intracellular entrance: the depth of this energy well was modulated strongly by the protonation state of E222. The energy barrier for Cl(-) was found to be 3-4 times higher than that for Na(+). Ion permeation characteristics were determined through BD simulations using a hybrid MD/continuum electrostatics approach to evaluate the energy profiles governing the ion movement. The resultant channel conductance and a near-zero permeability ratio (P(Cl)/P(Na)) were comparable to experimental data. On the basis of these calculations, we suggest that a ring of five E222 residues may act as an electrostatic gate. In addition, the hydrophobic gate region may play a role in charge selectivity due to a higher dehydration energy barrier for Cl(-) ions. The effect of halothane on the Na(+) PMF was also evaluated. Halothane was found to perturb salt bridges in GLIC that may be crucial for channel gating and open-channel stability, but had no significant impact on the single ion PMF profiles.
| S-EPMC3071019 | biostudies-literature
Project description: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.
| S-EPMC1300661 | biostudies-other
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.
| S-EPMC1304654 | biostudies-literature
Project description:OmpF, a multiionic porin from Escherichia coli, is a useful protypical model system for addressing general questions about electrostatic interactions in the confinement of an aqueous molecular pore. Here, favorable anion locations in the OmpF pore were mapped by anomalous X-ray scattering of Br(–) ions from four different crystal structures and compared with Mg(2+) sites and Rb(+) sites from a previous anomalous diffraction study to provide a complete picture of cation and anion transfer paths along the OmpF channel. By comparing structures with various crystallization conditions, we find that anions bind in discrete clusters along the entire length of the OmpF pore, whereas cations find conserved binding sites with the extracellular, surface-exposed loops. Results from molecular dynamics simulations are consistent with the experimental data and help highlight the critical residues that preferentially contact either cations or anions during permeation. Analysis of these results provides new insights into the molecular mechanisms that determine ion selectivity in OmpF porin.
| S-EPMC3868647 | biostudies-literature
Project description:A signal property of connexin channels is the ability to mediate selective diffusive movement of molecules through plasma membrane(s), but the energetics and determinants of molecular movement through these channels have yet to be understood. Different connexin channels have distinct molecular selectivities that cannot be explained simply on the basis of size or charge of the permeants. To gain insight into the forces and interactions that underlie selective molecular permeation, we investigated the energetics of two uncharged derivatized sugars, one permeable and one impermeable, through a validated connexin26 (Cx26) channel structural model, using molecular dynamics and associated analytic tools. The system is a Cx26 channel equilibrated in explicit membrane/solvent, shown by Brownian dynamics to reproduce key conductance characteristics of the native channel. The results are consistent with the known difference in permeability to each molecule. The energetic barriers extend through most of the pore length, rather than being highly localized as in ion-specific channels. There is little evidence for binding within the pore. Force decomposition reveals how, for each tested molecule, interactions with water and the Cx26 protein vary over the length of the pore and reveals a significant contribution from hydrogen bonding and interaction with K(+). The flexibility of the pore width varies along its length, and the tested molecules have differential effects on pore width as they pass through. Potential sites of interaction within the pore are defined for each molecule. The results suggest that for the tested molecules, differences in hydrogen bonding and entropic factors arising from permeant flexibility substantially contribute to the energetics of permeation. This work highlights factors involved in selective molecular permeation that differ from those that define selectivity among atomic ions.
| S-EPMC4744205 | biostudies-literature
Project description:Fluoride ion channels of the Fluc family combat toxicity arising from accumulation of environmental F-. Although crystal structures are known, the densely packed pore region has precluded delineation of the ion pathway. Here we chart out the Fluc pore and characterize its chemical requirements for transport. A ladder of H-bond donating residues creates a 'polar track' demarking the ion-conduction pathway. Surprisingly, while track polarity is well conserved, polarity is nonetheless functionally dispensable at several positions. A threonine at one end of the pore engages in vital interactions through its ?-branched methyl group. Two critical central phenylalanines that directly coordinate F- through a quadrupolar-ion interaction cannot be functionally substituted by aromatic, non-polar, or polar sidechains. The only functional replacement is methionine, which coordinates F- through its partially positive ?-methylene in mimicry of phenylalanine's quadrupolar interaction. These results demonstrate the unusual chemical requirements for selectively transporting the strongly H-bonding F- anion.
| S-EPMC5636608 | biostudies-literature