Brownian dynamics simulations of ion transport through the VDAC

Dynamics accelerate the kinetics of ion diffusion through channels: Continuous-time random walk models beyond the mean field approximation

Ion channels are proteins that play a significant role in physiological processes, including neuronal excitability and signal transduction. However, the precise mechanisms by which these proteins facilitate ion diffusion through cell membranes are not well understood. This is because experimental techniques to characterize ion channel activity operate on a time scale too large to understand the role of the various protein conformations on diffusion. Meanwhile, computational approaches operate on a time scale too short to rationalize the observed behavior at the microscopic scale. In this paper, we present a continuous-time random walk model that aims to bridge the scales between the atomistic models of ion channels and the experimental measurement of their conductance. We show how diffusion slows down in complex systems by using 3D lattices that map out the pore geometry of two channels: Nav1.7 and gramicidin. We also introduce spatial and dynamic site disorder to account for system heterogeneity beyond the mean field approximation. Computed diffusion coefficients show that an increase in spatial disorder slows down diffusion kinetics, while dynamic disorder has the opposite effect. Our results imply that microscopic or phenomenological models based on the potential of mean force data overlook the functional importance of protein dynamics on ion diffusion through channels.

Fast prediction of antibiotic permeability through membrane channels using Brownian dynamics

The efficient permeation across the Gram-negative bacterial membrane is an important step in the overall process of antibacterial action of a molecule and the one that has posed a significant hurdle on the way toward approved antibiotics. Predicting the permeability for a large library of molecules and assessing the effect of different molecular transformations on permeation rates of a given molecule is critical to the development of effective antibiotics. We present a computational approach for obtaining estimates of molecular permeability through a porin channel in a matter of hours using a Brownian dynamics approach. The fast sampling using a temperature acceleration scheme enables the approximate estimation of permeability using the inhomogeneous solubility diffusion model. Although the method is a significant approximation to similar all-atom approaches tested previously, we show that the present approach predicts permeabilities that correlate fairly well with the respective experimental permeation rates from liposome swelling experiments and accumulation rates from antibiotic accumulation assays, and is significantly, i.e., about 14 times, faster compared with a previously reported approach. The possible applications of the scheme in high-throughput screening for fast permeators are discussed.

The Single Residue K12 Governs the Exceptional Voltage Sensitivity of Mitochondrial Voltage-Dependent Anion Channel Gating

The voltage-dependent anion channel (VDAC) is a β-barrel channel of the mitochondrial outer membrane (MOM) that passively transports ions, metabolites, polypeptides, and single-stranded DNA. VDAC responds to a transmembrane potential by "gating," i.e. transitioning to one of a variety of low-conducting states of unknown structure. The gated state results in nearly complete suppression of multivalent mitochondrial metabolite (such as ATP and ADP) transport, while enhancing calcium transport. Voltage gating is a universal property of β-barrel channels, but VDAC gating is anomalously sensitive to transmembrane potential. Here, we show that a single residue in the pore interior, K12, is responsible for most of VDAC's voltage sensitivity. Using the analysis of over 40 μs of atomistic molecular dynamics (MD) simulations, we explore correlations between motions of charged residues inside the VDAC pore and geometric deformations of the β-barrel. Residue K12 is bistable; its motions between two widely separated positions along the pore axis enhance the fluctuations of the β-barrel and augment the likelihood of gating. Single channel electrophysiology of various K12 mutants reveals a dramatic reduction of the voltage-induced gating transitions. The crystal structure of the K12E mutant at a resolution of 2.6 Å indicates a similar architecture of the K12E mutant to the wild type; however, 60 μs of atomistic MD simulations using the K12E mutant show restricted motion of residue 12, due to enhanced connectivity with neighboring residues, and diminished amplitude of barrel motions. We conclude that β-barrel fluctuations, governed particularly by residue K12, drive VDAC gating transitions.

A Deep Dive into VDAC1 Conformational Diversity Using All-Atom Simulations Provides New Insights into the Structural Origin of the Closed States

The voltage-dependent anion channel 1 (VDAC1) is a crucial mitochondrial transporter that controls the flow of ions and respiratory metabolites entering or exiting mitochondria. As a voltage-gated channel, VDAC1 can switch between a high-conducting “open” state and a low-conducting “closed” state emerging at high transmembrane (TM) potentials. Although cell homeostasis depends on channel gating to regulate the transport of ions and metabolites, structural hallmarks characterizing the closed states remain unknown. Here, we performed microsecond accelerated molecular dynamics to highlight a vast region of VDAC1 conformational landscape accessible at typical voltages known to promote closure. Conformers exhibiting durable subconducting properties inherent to closed states were identified. In all cases, the low conductance was due to the particular positioning of an unfolded part of the N-terminus, which obstructed the channel pore. While the N-terminal tail was found to be sensitive to voltage orientation, our models suggest that stable low-conducting states of VDAC1 predominantly take place from disordered events and do not result from the displacement of a voltage sensor or a significant change in the pore. In addition, our results were consistent with conductance jumps observed experimentally and corroborated a recent study describing entropy as a key factor for VDAC gating.

How to Enter a Bacterium: Bacterial Porins and the Permeation of Antibiotics

Despite tremendous successes in the field of antibiotic discovery seen in the previous century, infectious diseases have remained a leading cause of death. More specifically, pathogenic Gram-negative bacteria have become a global threat due to their extraordinary ability to acquire resistance against any clinically available antibiotic, thus urging for the discovery of novel antibacterial agents. One major challenge is to design new antibiotics molecules able to rapidly penetrate Gram-negative bacteria in order to achieve a lethal intracellular drug accumulation. Protein channels in the outer membrane are known to form an entry route for many antibiotics into bacterial cells. Up until today, there has been a lack of simple experimental techniques to measure the antibiotic uptake and the local concentration in subcellular compartments. Hence, rules for translocation directly into the various Gram-negative bacteria via the outer membrane or via channels have remained elusive, hindering the design of new or the improvement of existing antibiotics. In this review, we will discuss the recent progress, both experimentally as well as computationally, in understanding the structure-function relationship of outer-membrane channels of Gram-negative pathogens, mainly focusing on the transport of antibiotics.

The Open State Selectivity of the Bean Seed VDAC Depends on Stigmasterol and Ion Concentration

The voltage-dependent anion channel (VDAC) is the major pathway for metabolites and ions transport through the mitochondrial outer membrane. It can regulate the flow of solutes by switching to a low conductance state correlated with a selectivity reversal, or by a selectivity inversion of its open state. The later one was observed in non-plant VDACs and is poorly characterized. We aim at investigating the selectivity inversion of the open state using plant VDAC purified from Phaseolus coccineus (PcVDAC) to evaluate its physiological role. Our main findings are: (1) The VDAC selectivity inversion of the open state occurs in PcVDAC, (2) Ion concentration and stigmasterol affect the occurrence of the open state selectivity inversion and stigmasterol appears to interact directly with PcVDAC. Interestingly, electrophysiological data concerning the selectivity inversion of the PcVDAC open state suggests that the phenomenon probably does not have a significant physiological effect in vivo.

Millisecond-Long Simulations of Antibiotics Transport through Outer Membrane Channels

To reach their target site inside Gram-negative bacteria, almost all antibiotics need to cross the outer membrane. Computational modeling of such processes can be numerically demanding due to the size of the systems and especially due to the timescales involved. Recently, a hybrid Brownian and molecular dynamics approach, i.e., Brownian dynamics including explicit atoms (BRODEA), has been developed and evaluated for studying the transport of monoatomic ions through membrane channels. Later on, this numerically efficient scheme has been applied to determine the free energy surfaces of the ciprofloxacin and enrofloxacin translocation through the porin OmpC using temperature-accelerated simulations. To improve the usability and accuracy of the approach, schemes to approximate the position-dependent diffusion constant of the molecule while traversing the pore had to be established. To this end, we have studied the translocation of the charged phosphonic acid antibiotic fosfomycin through the porin OmpF from Escherichia coli devising and benchmarking several diffusion models. To test the efficiency and sensitivity of these models, the effect of OmpF mutations on the permeation of fosfomycin was analyzed. Permeation events have been recorded over millisecond-long biased and unbiased simulations, from which thermodynamics and kinetics quantities of the translocation processes were determined. As a result, the use of the BRODEA approach, together with the appropriate diffusion model, was seen to accurately reproduce the findings observed in electrophysiology experiments and all-atom molecular dynamics simulations. These results suggest that the BRODEA approach can become a valuable tool for screening numerous compounds to evaluate their outer membrane permeability, a property important in the development of new antibiotics.

Free Energy Calculations for Protein-Ligand Binding Prediction

Computational prediction of protein-ligand binding involves initial determination of the binding mode and subsequent evaluation of the strength of the protein-ligand interactions, which directly correlates with ligand binding affinities. As a consequence of increasing computer power, rigorous approaches to calculate protein-ligand binding affinities, such as free energy perturbation (FEP) methods, are becoming an essential part of the toolbox of computer-aided drug design. In this chapter, we provide a general overview of these methods and introduce the QFEP modules, which are open-source API workflows based on our molecular dynamics (MD) package Q. The module QligFEP allows estimation of relative binding affinities along ligand series, while QresFEP is a module to estimate binding affinity shifts caused by single-point mutations of the protein. We herein provide guidelines for the use of each of these modules based on data extracted from ligand-design projects. While these modules are stand-alone, the combined use of the two workflows in a drug-design project yields complementary perspectives of the ligand binding problem, providing two sides of the same coin. The selected case studies illustrate how to use QFEP to approach the two key questions associated with ligand binding prediction: identifying the most favorable binding mode from different alternatives and establishing structure-affinity relationships that allow the rational optimization of hit compounds.

The inverse relation between mitochondrial transmembrane potential and proteins α-helix in neuronal-like cells under static magnetic field and the role of VDAC

In this study, a correlation between cell channel α-helices displacement and the mitochondrial transmembrane potential after exposure of 3, 7, 15 and 24 h of neuronal-like cells to a uniform magnetic field at the intensity of 2 mT was shown. Fourier Transform Infrared (FTIR) Spectroscopy and fluorescence techniques were used to analyze the secondary structure of protein content and mitochondrial transmembrane potential, respectively. The main result of this study was represented by a significant inverse relation between the mitochondrial transmembrane potential and the intensity of the Amide I band that can be associated with time exposure. Given that mitochondrial transmembrane potential should be related to the gating state of voltage-dependent anion channel (VDAC) in mitochondrial membrane, this result could have a relevant role in medicine. Indeed, VDAC's irregular behavior can be associated with several varieties of mitochondria-associated pathologies and various forms of cancer and neurodegeneration.

QligFEP: an automated workflow for small molecule free energy calculations in Q

The process of ligand binding to a biological target can be represented as the equilibrium between the relevant solvated and bound states of the ligand. This which is the basis of structure-based, rigorous methods such as the estimation of relative binding affinities by free energy perturbation (FEP). Despite the growing capacity of computing power and the development of more accurate force fields, a high throughput application of FEP is currently hampered due to the need, in the current schemes, of an expert user definition of the "alchemical" transformations between molecules in the series explored. Here, we present QligFEP, a solution to this problem using an automated workflow for FEP calculations based on a dual topology approach. In this scheme, the starting poses of each of the two ligands, for which the relative affinity is to be calculated, are explicitly present in the MD simulations associated with the (dual topology) FEP transformation, making the perturbation pathway between the two ligands univocal. We show that this generalized method can be applied to accurately estimate solvation free energies for amino acid sidechain mimics, as well as the binding affinity shifts due to the chemical changes typical of lead optimization processes. This is illustrated in a number of protein systems extracted from other FEP studies in the literature: inhibitors of CDK2 kinase and a series of A2A adenosine G protein-coupled receptor antagonists, where the results obtained with QligFEP are in excellent agreement with experimental data. In addition, our protocol allows for scaffold hopping perturbations to identify the binding affinities between different core scaffolds, which we illustrate with a series of Chk1 kinase inhibitors. QligFEP is implemented in the open-source MD package Q, and works with the most common family of force fields: OPLS, CHARMM and AMBER.

Ten Years of High Resolution Structural Research on the Voltage Dependent Anion Channel (VDAC)-Recent Developments and Future Directions

Mitochondria are evolutionarily related to Gram-negative bacteria and both comprise two membrane systems with strongly differing protein composition. The major protein in the outer membrane of mitochondria is the voltage-dependent anion channel (VDAC), which mediates signal transmission across the outer membrane but also the exchange of metabolites, most importantly ADP and ATP. More than 30 years after its discovery three identical high-resolution structures were determined in 2008. These structures show a 19-stranded anti-parallel beta-barrel with an N-terminal helix located inside. An odd number of beta-strands is also shared by Tom40, another member of the VDAC superfamily. This indicates that this superfamily is evolutionarily relatively young and that it has emerged in the context of mitochondrial evolution. New structural information obtained during the last decade on Tom40 can be used to cross-validate the structure of VDAC and vice versa. Interpretation of biochemical and biophysical studies on both protein channels now rests on a solid basis of structural data. Over the past 10 years, complementary structural and functional information on proteins of the VDAC superfamily has been collected from in-organello, in-vitro, and in silico studies. Most of these findings have confirmed the validity of the original structures. This short article briefly reviews the most important advances on the structure and function of VDAC superfamily members collected during the last decade and summarizes how they enhanced our understanding of the channel.

Phosphorylation of rat brain purified mitochondrial Voltage-Dependent Anion Channel by c-Jun N-terminal kinase-3 modifies open-channel noise

The drift kinetic energy of ionic flow through single ion channels cause vibrations of the pore walls which are observed as open-state current fluctuations (open-channel noise) during single-channel recordings. Vibration of the pore wall leads to transitions among different conformational sub-states of the channel protein in the open-state. Open-channel noise analysis can provide important information about the different conformational sub-state transitions and how biochemical modifications of ion channels would affect their transport properties. It has been shown that c-Jun N-terminal kinase-3 (JNK3) becomes activated by phosphorylation in various neurodegenerative diseases and phosphorylates outer mitochondrion associated proteins leading to neuronal apoptosis. In our earlier work, JNK3 has been reported to phosphorylate purified rat brain mitochondrial voltage-dependent anion channel (VDAC) in vitro and modify its conductance and opening probability. In this article we have compared the open-state noise profile of the native and the JNK3 phosphorylated VDAC using Power Spectral Density vs frequency plots. Power spectral density analysis of open-state noise indicated power law with average slope value α ≈1 for native VDAC at both positive and negative voltage whereas average α value < 0.5 for JNK3 phosphorylated VDAC at both positive and negative voltage. It is proposed that 1/f¹ power law in native VDAC open-state noise arises due to coupling of ionic transport and conformational sub-states transitions in open-state and this coupling is perturbed as a result of channel phosphorylation.

CHARMM-GUI 10 years for biomolecular modeling and simulation

CHARMM-GUI, http://www.charmm-gui.org, is a web-based graphical user interface that prepares complex biomolecular systems for molecular simulations. CHARMM-GUI creates input files for a number of programs including CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Since its original development in 2006, CHARMM-GUI has been widely adopted for various purposes and now contains a number of different modules designed to set up a broad range of simulations: (1) PDB Reader & Manipulator, Glycan Reader, and Ligand Reader & Modeler for reading and modifying molecules; (2) Quick MD Simulator, Membrane Builder, Nanodisc Builder, HMMM Builder, Monolayer Builder, Micelle Builder, and Hex Phase Builder for building all-atom simulation systems in various environments; (3) PACE CG Builder and Martini Maker for building coarse-grained simulation systems; (4) DEER Facilitator and MDFF/xMDFF Utilizer for experimentally guided simulations; (5) Implicit Solvent Modeler, PBEQ-Solver, and GCMC/BD Ion Simulator for implicit solvent related calculations; (6) Ligand Binder for ligand solvation and binding free energy simulations; and (7) Drude Prepper for preparation of simulations with the CHARMM Drude polarizable force field. Recently, new modules have been integrated into CHARMM-GUI, such as Glycolipid Modeler for generation of various glycolipid structures, and LPS Modeler for generation of lipopolysaccharide structures from various Gram-negative bacteria. These new features together with existing modules are expected to facilitate advanced molecular modeling and simulation thereby leading to an improved understanding of the structure and dynamics of complex biomolecular systems. Here, we briefly review these capabilities and discuss potential future directions in the CHARMM-GUI development project. © 2016 Wiley Periodicals, Inc.

Gating of Connexin Channels by transjunctional-voltage: Conformations and models of open and closed states

Voltage is an important physiologic regulator of channels formed by the connexin gene family. Connexins are unique among ion channels in that both plasma membrane inserted hemichannels (undocked hemichannels) and intercellular channels (aggregates of which form gap junctions) have important physiological roles. The hemichannel is the fundamental unit of gap junction voltage-gating. Each hemichannel displays two distinct voltage-gating mechanisms that are primarily sensitive to a voltage gradient formed along the length of the channel pore (the transjunctional voltage) rather than sensitivity to the absolute membrane potential (V_m or V_i-o). These transjunctional voltage dependent processes have been termed V_j- or fast-gating and loop- or slow-gating. Understanding the mechanism of voltage-gating, defined as the sequence of voltage-driven transitions that connect open and closed states, first and foremost requires atomic resolution models of the end states. Although ion channels formed by connexins were among the first to be characterized structurally by electron microscopy and x-ray diffraction in the early 1980's, subsequent progress has been slow. Much of the current understanding of the structure-function relations of connexin channels is based on two crystal structures of Cx26 gap junction channels. Refinement of crystal structure by all-atom molecular dynamics and incorporation of charge changing protein modifications has resulted in an atomic model of the open state that arguably corresponds to the physiologic open state. Obtaining validated atomic models of voltage-dependent closed states is more challenging, as there are currently no methods to solve protein structure while a stable voltage gradient is applied across the length of an oriented channel. It is widely believed that the best approach to solve the atomic structure of a voltage-gated closed ion channel is to apply different but complementary experimental and computational methods and to use the resulting information to derive a consensus atomic structure that is then subjected to rigorous validation. In this paper, we summarize our efforts to obtain and validate atomic models of the open and voltage-driven closed states of undocked connexin hemichannels. This article is part of a Special Issue entitled: Gap Junction Proteins edited by Jean Claude Herve.

Modulation of Asymmetric Flux in Heterotypic Gap Junctions by Pore Shape, Particle Size and Charge

Gap junction channels play a vital role in intercellular communication by connecting cytoplasm of adjoined cells through arrays of channel-pores formed at the common membrane junction. Their structure and properties vary depending on the connexin isoform(s) involved in forming the full gap junction channel. Lack of information on the molecular structure of gap junction channels has limited the development of computational tools for single channel studies. Currently, we rely on cumbersome experimental techniques that have limited capabilities. We have earlier reported a simplified Brownian dynamics gap junction pore model and demonstrated that variations in pore shape at the single channel level can explain some of the differences in permeability of heterotypic channels observed in in vitro experiments. Based on this computational model, we designed simulations to study the influence of pore shape, particle size and charge in homotypic and heterotypic pores. We simulated dye diffusion under whole cell voltage clamping. Our simulation studies with pore shape variations revealed a pore shape with maximal flux asymmetry in a heterotypic pore. We identified pore shape profiles that match the in silico flux asymmetry results to the in vitro results of homotypic and heterotypic gap junction formed out of Cx43 and Cx45. Our simulation results indicate that the channel's pore-shape established flux asymmetry and that flux asymmetry is primarily regulated by the sizes of the conical and/or cylindrical mouths at each end of the pore. Within the set range of particle size and charge, flux asymmetry was found to be independent of particle size and directly proportional to charge magnitude. While particle charge was vital to creating flux asymmetry, charge magnitude only scaled the observed flux asymmetry. Our studies identified the key factors that help predict asymmetry. Finally, we suggest the role of such flux asymmetry in creating concentration imbalances of messenger molecules in cardiomyocytes. We also assess the potency of fibroblasts in aggravating such imbalances through Cx43-Cx45 heterotypic channels in fibrotic heart tissue.

Voltage Dependence of Conformational Dynamics and Subconducting States of VDAC-1

The voltage-dependent anion channel 1 (VDAC-1) is an important protein of the outer mitochondrial membrane that transports energy metabolites and is involved in apoptosis. The available structures of VDAC proteins show a wide β-stranded barrel pore, with its N-terminal α-helix (N-α) bound to its interior. Electrophysiology experiments revealed that voltage, its polarity, and membrane composition modulate VDAC currents. Experiments with VDAC-1 mutants identified amino acids that regulate the gating process. However, the mechanisms for how these factors regulate VDAC-1, and which changes they trigger in the channel, are still unknown. In this study, molecular dynamics simulations and single-channel experiments of VDAC-1 show agreement for the current-voltage relationships of an "open" channel and they also show several subconducting transient states that are more cation selective in the simulations. We observed voltage-dependent asymmetric distortions of the VDAC-1 barrel and the displacement of particular charged amino acids. We constructed conformational models of the protein voltage response and the pore changes that consistently explain the protein conformations observed at opposite voltage polarities, either in phosphatidylethanolamine or phosphatidylcholine membranes. The submicrosecond VDAC-1 voltage response shows intrinsic structural changes that explain the role of key gating amino acids and support some of the current gating hypotheses. These voltage-dependent protein changes include asymmetric barrel distortion, its interaction with the membrane, and significant displacement of N-α amino acids.

Computational simulations of asymmetric fluxes of large molecules through gap junction channel pores

Gap junction channels are formed out of connexin isoforms, which enable molecule and ion selective diffusion amongst neighboring cells. HeLa cells expressing distinct connexins (Cx) allow the formation of heterotypic channels, where we observed a molecular charge-independent preferential flux of large fluorescent molecules in the Cx45 to Cx43 direction. We hypothesize that the pore's shape is a significant factor along-side charge and transjunctional voltages for this asymmetric flux. To test this hypothesis, we developed a 3D computational model simulating Brownian diffusion of large molecules in a gap junction channel pore. The basic pore contour was derived from x-ray crystallographic structures of Cx43 and Cx26 and approximated using basic geometric shapes. Lucifer yellow dye molecules and cesium counter-ions were modeled as spheres using their respective Stokes radii. Our simulation results from simple diffusion and constant concentration gradient experiments showed that only charged particles yield asymmetric fluxes in heterotypic pores. While increasing the inner mouth size resulted in a near-quadratic rise in flux, the rise was asymptotic for outer mouth radii increase. Probability maps and average force per particle per pore section explain the asymmetric flux with variation in pore shape. Furthermore, the simulation results are in agreement with our in vitro experimental results with HeLa cells in Cx43-Cx45 heterotypic configurations. The presence of asymmetric fluxes can help us to understand effects of the molecular structure of the pore and predict potential differences in vivo.

Current state of theoretical and experimental studies of the voltage-dependent anion channel (VDAC)

Voltage-dependent anion channel (VDAC), the major channel of the mitochondrial outer membrane provides a controlled pathway for respiratory metabolites in and out of the mitochondria. In spite of the wealth of experimental data from structural, biochemical, and biophysical investigations, the exact mechanisms governing selective ion and metabolite transport, especially the role of titratable charged residues and interactions with soluble cytosolic proteins, remain hotly debated in the field. The computational advances hold a promise to provide a much sought-after solution to many of the scientific disputes around solute and ion transport through VDAC and hence, across the mitochondrial outer membrane. In this review, we examine how Molecular Dynamics, Free Energy, and Brownian Dynamics simulations of the large β-barrel channel, VDAC, advanced our understanding. We will provide a short overview of non-conventional techniques and also discuss examples of how the modeling excursions into VDAC biophysics prospectively aid experimental efforts. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Simulations of outer membrane channels and their permeability

Channels in the outer membrane of Gram-negative bacteria provide essential pathways for the controlled and unidirectional transport of ions, nutrients and metabolites into the cell. At the same time the outer membrane serves as a physical barrier for the penetration of noxious substances such as antibiotics into the bacteria. Most antibiotics have to pass through these membrane channels to either reach cytoplasmic bound targets or to further cross the hydrophobic inner membrane. Considering the pharmaceutical significance of antibiotics, understanding the functional role and mechanism of these channels is of fundamental importance in developing strategies to design new drugs with enhanced permeation abilities. Due to the biological complexity of membrane channels and experimental limitations, computer simulations have proven to be a powerful tool to investigate the structure, dynamics and interactions of membrane channels. Considerable progress has been made in computer simulations of membrane channels during the last decade. The goal of this review is to provide an overview of the computational techniques and their roles in modeling the transport across outer membrane channels. A special emphasis is put on all-atom molecular dynamics simulations employed to better understand the transport of molecules. Moreover, recent molecular simulations of ion, substrate and antibiotics translocation through membrane pores are briefly summarized. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Dual mechanism of ion permeation through VDAC revealed with inorganic phosphate ions and phosphate metabolites

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.

Voltage regulation of connexin channel conductance

Voltage is an important parameter that regulates the conductance of both intercellular and plasma membrane channels (undocked hemichannels) formed by the 21 members of the mammalian connexin gene family. Connexin channels display two forms of voltage-dependence, rectification of ionic currents and voltage-dependent gating. Ionic rectification results either from asymmetries in the distribution of fixed charges due to heterotypic pairing of different hemichannels, or by channel block, arising from differences in the concentrations of divalent cations on opposite sides of the junctional plaque. This rectification likely underpins the electrical rectification observed in some electrical synapses. Both intercellular and undocked hemichannels also display two distinct forms of voltage-dependent gating, termed Vj (fast)-gating and loop (slow)-gating. This review summarizes our current understanding of the molecular determinants and mechanisms underlying these conformational changes derived from experimental, molecular-genetic, structural, and computational approaches.

Voltage-dependent anion channels: the wizard of the mitochondrial outer membrane

Voltage dependent anion channels (VDACs) are the most abundant proteins in the outer mitochondrial membrane. Although they are essential in metabolite exchange, cell defense and apoptosis, the molecular mechanism of these VDAC-mediated processes remains elusive. Here we review recent progress in terms of VDACs’ structure and regulation, with a special focus on the molecular aspects of gating and the interaction with effector proteins.

Computational investigation of cholesterol binding sites on mitochondrial VDAC

The mitochondrial voltage-dependent anion channel (VDAC) allows passage of ions and metabolites across the mitochondrial outer membrane. Cholesterol binds mammalian VDAC, and we investigated the effects of binding to human VDAC1 with atomistic molecular dynamics simulations that totaled 1.4 μs. We docked cholesterol to specific sites on VDAC that were previously identified with NMR, and we tested the reliability of multiple docking results in each site with simulations. The most favorable binding modes were used to build a VDAC model with cholesterol occupying five unique sites, and during multiple 100 ns simulations, cholesterol stably and reproducibly remained bound to the protein. For comparison, VDAC was simulated in systems with identical components but with cholesterol initially unbound. The dynamics of loops that connect adjacent β-strands were most affected by bound cholesterol, with the averaged root-mean-square fluctuation (RMSF) of multiple residues altered by 20–30%. Cholesterol binding also stabilized charged residues inside the channel and localized the surrounding electrostatic potentials. Despite this, ion diffusion through the channel was not significantly affected by bound cholesterol, as evidenced by multi-ion potential of mean force measurements. Although we observed modest effects of cholesterol on the open channel, our model will be particularly useful in experiments that investigate how cholesterol affects VDAC function under applied electrochemical forces and also how other ligands and proteins interact with the channel.

Charged residues distribution modulates selectivity of the open state of human isoforms of the voltage dependent anion-selective channel

Voltage Dependent Anion-selective Channels (VDACs) are pore-forming proteins located in the outer mitochondrial membrane. They are responsible for the access of ions and energetic metabolites into the inner membrane transport systems. Three VDAC isoforms exist in mammalian, but their specific role is unknown. In this work we have performed extensive (overall ∼5 µs) Molecular Dynamics (MD) simulations of the human VDAC isoforms to detect structural and conformational variations among them, possibly related to specific functional roles of these proteins. Secondary structure analysis of the N-terminal domain shows a high similarity among the three human isoforms of VDAC but with a different plasticity. In particular, the N-terminal domain of the hVDAC1 is characterized by a higher plasticity, with a ∼20% occurrence for the 'unstructured' conformation throughout the folded segment, while hVDAC2, containing a peculiar extension of 11 amino acids at the N-terminal end, presents an additional 310-helical folded portion comprising residues 10' to 3, adhering to the barrel wall. The N-terminal sequences of hVDAC isoforms are predicted to have a low flexibility, with possible consequences in the dynamics of the human VDACs. Clear differences were found between hVDAC1 and hVDAC3 against hVDAC2: a significantly modified dynamics with possible important consequence on the voltage-gating mechanism. Charge distribution inside and at the mouth of the pore is responsible for a different preferential localization of ions with opposite charge and provide a valuable rationale for hVDAC1 and hVDAC3 having a Cl-/K+ selectivity ratio of 1.8, whereas hVDAC2 of 1.4. Our conclusion is that hVDAC isoforms, despite sharing a similar scaffold, have modified working features and a biological work is now requested to give evidence to the described dissimilarities.

Conductance simulation of the purinergic P2X2, P2X4, and P2X7 ionic channels using a combined Brownian dynamics and molecular dynamics approach

This paper investigates the application of an original combined approach of molecular and Brownian dynamic methods with quantum chemistry calculations for modeling the process of conductance of ion channels using purinergic P2X family receptors P2X2, P2X4, and P2X7 as a case study. A simplified model of the ionic channel in the lipid bilayer has been developed. A high level of conductance (30 pS) of P2X2 ionic channel together with the key role of Asp349 in forming the selectivity filter of P2X2 has been shown by using this approach. Calculated P2X2 permeability to monovalent cations Li(+), Na(+), and K(+) conforms to the free diffusion coefficient of these ions, which shows the low selectivity of P2X2 ionic channel.

Molecular dynamics simulations of membrane proteins under asymmetric ionic concentrations

A computational method is developed to allow molecular dynamics simulations of biomembrane systems under realistic ionic gradients and asymmetric salt concentrations while maintaining the conventional periodic boundary conditions required to minimize finite-size effects in an all-atom explicit solvent representation. The method, which consists of introducing a nonperiodic energy step acting on the ionic species at the edge of the simulation cell, is first tested with illustrative applications to a simple membrane slab model and a phospholipid membrane bilayer. The nonperiodic energy-step method is then used to calculate the reversal potential of the bacterial porin OmpF, a large cation-specific β-barrel channel, by simulating the I-V curve under an asymmetric 10:1 KCl concentration gradient. The calculated reversal potential of 28.6 mV is found to be in excellent agreement with the values of 26–27 mV measured from lipid bilayer experiments, thereby demonstrating that the method allows realistic simulations of nonequilibrium membrane transport with quantitative accuracy. As a final example, the pore domain of Kv1.2, a highly selective voltage-activated K⁺ channel, is simulated in a lipid bilayer under conditions that recreate, for the first time, the physiological K⁺ and Na⁺ concentration gradients and the electrostatic potential difference of living cells.

Ion permeation in the NanC porin from Escherichia coli: free energy calculations along pathways identified by coarse-grain simulations

Using the X-ray structure of a recently discovered bacterial protein, the N-acetylneuraminic acid-inducible channel (NanC), we investigate computationally K(+) and Cl(-) ions' permeation. We identify ion permeation pathways that are likely to be populated using coarse-grain Monte Carlo simulations. Next, we use these pathways as reaction coordinates for umbrella sampling-based free energy simulations. We find distinct tubelike pathways connecting specific binding sites for K(+) and, more pronounced, for Cl(-) ions. Both ions permeate the porin preserving almost all of their first hydration shell. The calculated free energy barriers are G(#) ≈ 4 kJ/mol and G(#) ≈ 8 kJ/mol for Cl(-) and K(+), respectively. Within the approximations associated with these values, discussed in detail in this work, we suggest that the porin is slightly selective for Cl(-) versus K(+). Our suggestion is consistent with the experimentally observed weak Cl(-) over K(+) selectivity. A rationale for the latter is suggested by a comparison with previous calculations on strongly anion selective porins.

Identification of a cation transport pathway in Neisseria meningitidis PorB

Neisseria meningitidis is the main causative agent of bacterial meningitis. In its outer membrane, the trimeric Neisserial porin PorB is responsible for the diffusive transport of essential hydrophilic solutes across the bilayer. Previous molecular dynamics simulations based on the recent crystal structure of PorB have suggested the presence of distinct solute translocation pathways through this channel. Although PorB has been electrophysiologically characterized as anion-selective, cation translocation through nucleotide-bound PorB during pathogenesis is thought to be instrumental for host cell death. As a result, we were particularly interested in further characterizing cation transport through the pore. We combined a structural approach with additional computational analysis. Here, we present two crystal structures of PorB at 2.1 and 2.65 Å resolution. The new structures display additional electron densities around the protruding loop 3 (L3) inside the pore. We show that these electron densities can be identified as monovalent cations, in our case Cs(+), which are tightly bound to the inner channel. Molecular dynamics simulations reveal further ion interactions and the free energy landscape for ions inside PorB. Our results suggest that the crystallographically identified locations of Cs(+) form a cation transport pathway inside the pore. This finding suggests how positively charged ions are translocated through PorB when the channel is inserted into mitochondrial membranes during Neisserial infection, a process which is considered to dissipate the mitochondrial transmembrane potential gradient and thereby induce apoptosis.

Molecular origin of VDAC selectivity towards inorganic ions: a combined molecular and Brownian dynamics study

The voltage-dependent anion channel (VDAC) serves as the major pore for metabolites and electrolytes in the outer mitochondrial membrane. To refine our understanding of ion permeation through this channel we performed an extensive Brownian (BD) and molecular dynamics (MD) study on the mouse VDAC isoform 1 wild-type and mutants (K20E, D30K, K61E, E158K and K252E). The selectivity and the conductance of the wild-type and of the variant channels computed from the BD trajectories are in agreement with experimental data. The calculated selectivity is shown to be very sensitive to slight conformational changes which may have some bearing on the variability of the selectivity values measured on the VDAC open state. The MD and BD free energy profiles of the ion permeation suggest that the pore region comprising the N-terminal helix and the barrel band encircling it predominantly controls the ion transport across the channel. The overall 12μs BD and 0.9μs MD trajectories of the mouse VDAC isoform 1 wild-type and mutants feature no distinct pathways for ion diffusion and no long-lived ion-protein interactions. The dependence of ion distribution in the wild-type channel with the salt concentration can be explained by an ionic screening of the permanent charges of the protein arising from the pore. Altogether these results bolster the role of electrostatic features of the pore as the main determinant of VDAC selectivity towards inorganic anions.

Flexibility of the N-terminal mVDAC1 segment controls the channel's gating behavior

Since the solution of the molecular structures of members of the voltage dependent anion channels (VDACs), the N-terminal α-helix has been the main focus of attention, since its strategic location, in combination with its putative conformational flexibility, could define or control the channel’s gating characteristics. Through engineering of two double-cysteine mVDAC1 variants we achieved fixing of the N-terminal segment at the bottom and midpoint of the pore. Whilst cross-linking at the midpoint resulted in the channel remaining constitutively open, cross-linking at the base resulted in an “asymmetric” gating behavior, with closure only at one electric field´s orientation depending on the channel’s orientation in the lipid bilayer. Additionally, and while the native channel adopts several well-defined closed states (S1 and S2), the cross-linked variants showed upon closure a clear preference for the S2 state. With native-channel characteristics restored following reduction of the cysteines, it is evident that the conformational flexibility of the N-terminal segment plays indeed a major part in the control of the channel’s gating behavior.

β-Barrel mobility underlies closure of the voltage-dependent anion channel

The voltage-dependent anion channel (VDAC) is the major protein in the outer mitochondrial membrane, where it mediates transport of ATP and ADP. Changes in its permeability, induced by voltage or apoptosis-related proteins, have been implicated in apoptotic pathways. The three-dimensional structure of VDAC has recently been determined as a 19-stranded β-barrel with an in-lying N-terminal helix. However, its gating mechanism is still unclear. Using solid-state NMR spectroscopy, molecular dynamics simulations, and electrophysiology, we show that deletion of the rigid N-terminal helix sharply increases overall motion in VDAC's β-barrel, resulting in elliptic, semicollapsed barrel shapes. These states quantitatively reproduce conductance and selectivity of the closed VDAC conformation. Mutation of the N-terminal helix leads to a phenotype intermediate to the open and closed states. These data suggest that the N-terminal helix controls entry into elliptic β-barrel states which underlie VDAC closure. Our results also indicate that β-barrel channels are intrinsically flexible.

Voltage-dependent conformational changes in connexin channels

Channels formed by connexins display two distinct types of voltage-dependent gating, termed V(j)- or fast-gating and loop- or slow-gating. Recent studies, using metal bridge formation and chemical cross-linking have identified a region within the channel pore that contributes to the formation of the loop-gate permeability barrier. The conformational changes are remarkably large, reducing the channel pore diameter from 15 to 20Å to less than 4Å. Surprisingly, the largest conformational change occurs in the most stable region of the channel pore, the 3(10) or parahelix formed by amino acids in the 42-51 segment. The data provide a set of positional constraints that can be used to model the structure of the loop-gate closed state. Less is known about the conformation of the V(j)-gate closed state. There appear to be two different mechanisms; one in which conformational changes in channel structure are linked to a voltage sensor contained in the N-terminus of Cx26 and Cx32 and a second in which the C-terminus of Cx43 and Cx40 may act either as a gating particle to block the channel pore or alternatively to stabilize the closed state. The later mechanism utilizes the same domains as implicated in effecting pH gating of Cx43 channels. It is unclear if the two V(j)-gating mechanisms are related or if they represent different gating mechanisms that operate separately in different subsets of connexin channels. A model of the V(j)-closed state of Cx26 hemichannel that is based on the X-ray structure of Cx26 and electron crystallographic structures of a Cx26 mutation suggests that the permeability barrier for V(j)-gating is formed exclusively by the N-terminus, but recent information suggests that this conformation may not represent a voltage-closed state. Closed state models are considered from a thermodynamic perspective based on information from the 3.5Å Cx26 crystal structure and molecular dynamics (MD) simulations. The applications of computational and experimental methods to define the path of allosteric molecular transitions that link the open and closed states are discussed. This article is part of a Special Issue entitled: The Communicating junctions, composition, structure and characteristics.

Molecular dynamics simulations of the Cx26 hemichannel: insights into voltage-dependent loop-gating

Loop-gating is one of two voltage-dependent mechanisms that regulate the open probability of connexin channels. The loop-gate permeability barrier is formed by a segment of the first extracellular loop (E1) (the parahelix) and appears to be accompanied by straightening of the bend angle between E1 and the first transmembrane domain (TM1). Here, all-atom molecular dynamics simulations are used to identify and characterize interacting van der Waals and electrostatic networks that stabilize the parahelices and TM1/E1 bend angles of the open Cx26 hemichannel. Dynamic fluctuations in an electrostatic network in each subunit are directly linked to the stability of parahelix structure and TM1/E1 bend angle in adjacent subunits. The electrostatic network includes charged residues that are pore-lining and thus positioned to be voltage sensors. We propose that the transition to the closed state is initiated by voltage-driven disruption of the networks that stabilize the open-state parahelix configuration, allowing the parahelix to protrude into the channel pore to form the loop-gate barrier. Straightening of the TM1/E1 bend appears to be a consequence of the reorganization of the interacting networks that accompany the conformational change of the parahelix. The electrostatic network extends across subunit boundaries, suggesting a concerted gating mechanism.

Web interface for Brownian dynamics simulation of ion transport and its applications to beta-barrel pores

Brownian dynamics (BD) based on accurate potential of mean force is an efficient and accurate method for simulating ion transport through wide ion channels. Here, a web-based graphical user interface (GUI) is presented for carrying out grand canonical Monte Carlo (GCMC) BD simulations of channel proteins: http://www.charmm-gui.org/input/gcmcbd. The webserver is designed to help users avoid most of the technical difficulties and issues encountered in setting up and simulating complex pore systems. GCMC/BD simulation results for three proteins, the voltage dependent anion channel (VDAC), α-Hemolysin (α-HL), and the protective antigen pore of the anthrax toxin (PA), are presented to illustrate the system setup, input preparation, and typical output (conductance, ion density profile, ion selectivity, and ion asymmetry). Two models for the input diffusion constants for potassium and chloride ions in the pore are compared: scaling of the bulk diffusion constants by 0.5, as deduced from previous all-atom molecular dynamics simulations of VDAC, and a hydrodynamics based model (HD) of diffusion through a tube. The HD model yields excellent agreement with experimental conductances for VDAC and α-HL, while scaling bulk diffusion constants by 0.5 leads to underestimates of 10-20%. For PA, simulated ion conduction values overestimate experimental values by a factor of 1.5-7 (depending on His protonation state and the transmembrane potential), implying that the currently available computational model of this protein requires further structural refinement.

Plant VDAC: facts and speculations

The voltage-dependent anion-selective channel (VDAC) is the most abundant protein in the mitochondrial outer membrane and the major transport pathway for a large variety of compounds ranging from ions to large polymeric molecules such as DNA and tRNA. Plant VDACs feature a secondary structure content and electrophysiological properties akin to those of VDACs from other organisms. They however undergo a specific regulation. The general importance of VDAC in plant physiology has only recently emerged. Besides their role in metabolite transport, plant VDACs are also involved in the programmed cell death triggered in response to biotic and abiotic stresses. Moreover, their colocalization in non-mitochondrial membranes suggests a diversity of function. This review summarizes our current understanding of the structure and function of plant VDACs. This article is part of a Special Issue entitled: VDAC structure, function, and regulation of mitochondrial metabolism.

Concentration dependent ion selectivity in VDAC: a molecular dynamics simulation study

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.

Molecular dynamics simulations of the Cx26 hemichannel: evaluation of structural models with Brownian dynamics

The recently published crystal structure of the Cx26 gap junction channel provides a unique opportunity for elucidation of the structure of the conductive connexin pore and the molecular determinants of its ion permeation properties (conductance, current-voltage [I-V] relations, and charge selectivity). However, the crystal structure was incomplete, most notably lacking the coordinates of the N-terminal methionine residue, which resides within the pore, and also lacking two cytosolic domains. To allow computational studies for comparison with the known channel properties, we completed the structure. Grand canonical Monte Carlo Brownian dynamics (GCMC/BD) simulations of the completed and the published Cx26 hemichannel crystal structure indicate that the pore is too narrow to permit significant ion flux. The GCMC/BD simulations predict marked inward current rectification and almost perfect anion selectivity, both inconsistent with known channel properties. The completed structure was refined by all-atom molecular dynamics (MD) simulations (220 ns total) in an explicit solvent and POPC membrane system. These MD simulations produced an equilibrated structure with a larger minimal pore diameter, which decreased the height of the permeation barrier formed by the N terminus. GCMC/BD simulations of the MD-equilibrated structure yielded more appropriate single-channel conductance and less anion/cation selectivity. However, the simulations much more closely matched experimentally determined I-V relations when the charge effects of specific co- and posttranslational modifications of Cx26 previously identified by mass spectrometry were incorporated. We conclude that the average equilibrated structure obtained after MD simulations more closely represents the open Cx26 hemichannel structure than does the crystal structure, and that co- and posttranslational modifications of Cx26 hemichannels are likely to play an important physiological role by defining the conductance and ion selectivity of Cx26 channels. Furthermore, the simulations and data suggest that experimentally observed heterogeneity in Cx26 I-V relations can be accounted for by variation in co- and posttranslational modifications.

Molecular dynamics studies of ion permeation in VDAC

The voltage-dependent anion channel (VDAC) in the outer membrane of mitochondria serves an essential role in the transport of metabolites and electrolytes between the cell matrix and mitochondria. To examine its structure, dynamics, and the mechanisms underlying its electrophysiological properties, we performed a total of 1.77 μs molecular dynamics simulations of human VDAC isoform 1 in DOPE/DOPC mixed bilayers in 1 M KCl solution with transmembrane potentials of 0, ±25, ±50, ±75, and ±100 mV. The calculated conductance and ion selectivity are in good agreement with the experimental measurements. In addition, ion density distributions inside the channel reveal possible pathways for different ion species. Based on these observations, a mechanism underlying the anion selectivity is proposed; both ion species are transported across the channel, but the rate for K(+) is smaller than that for Cl(-) because of the attractive interactions between K(+) and residues on the channel wall. This difference leads to the anion selectivity of VDAC.