Dimeric Tubulin Modifies Mechanical Properties of Lipid Bilayer, as Probed Using Gramicidin A Channel.

Using the gramicidin A channel as a molecular probe, we show that tubulin binding to planar lipid membranes changes the channel kinetics-seen as an increase in the lifetime of the channel dimer-and thus points towards modification of the membrane's mechanical properties. The effect is more pronounced in the presence of non-lamellar lipids in the lipid mixture used for membrane formation. To interpret these findings, we propose that tubulin binding redistributes the lateral pressure of lipid packing along the membrane depth, making it closer to the profile expected for lamellar lipids. This redistribution happens because tubulin perturbs the lipid headgroup spacing to reach the membrane's hydrophobic core via its amphiphilic α-helical domain. Specifically, it increases the forces of repulsion between the lipid headgroups and reduces such forces in the hydrophobic region. We suggest that the effect is reciprocal, meaning that alterations in lipid bilayer mechanics caused by membrane remodeling during cell proliferation in disease and development may also modulate tubulin membrane binding, thus exerting regulatory functions. One of those functions includes the regulation of protein-protein interactions at the membrane surface, as exemplified by VDAC complexation with tubulin.

Restricting α-synuclein transport into mitochondria by inhibition of α-synuclein-VDAC complexation as a potential therapeutic target for Parkinson's disease treatment.

Involvement of alpha-synuclein (αSyn) in Parkinson's disease (PD) is complicated and difficult to trace on cellular and molecular levels. Recently, we established that αSyn can regulate mitochondrial function by voltage-activated complexation with the voltage-dependent anion channel (VDAC) on the mitochondrial outer membrane. When complexed with αSyn, the VDAC pore is partially blocked, reducing the transport of ATP/ADP and other metabolites. Further, αSyn can translocate into the mitochondria through VDAC, where it interferes with mitochondrial respiration. Recruitment of αSyn to the VDAC-containing lipid membrane appears to be a crucial prerequisite for both the blockage and translocation processes. Here we report an inhibitory effect of HK2p, a small membrane-binding peptide from the mitochondria-targeting N-terminus of hexokinase 2, on αSyn membrane binding, and hence on αSyn complex formation with VDAC and translocation through it. In electrophysiology experiments, the addition of HK2p at micromolar concentrations to the same side of the membrane as αSyn results in a dramatic reduction of the frequency of blockage events in a concentration-dependent manner, reporting on complexation inhibition. Using two complementary methods of measuring protein-membrane binding, bilayer overtone analysis and fluorescence correlation spectroscopy, we found that HK2p induces detachment of αSyn from lipid membranes. Experiments with HeLa cells using proximity ligation assay confirmed that HK2p impedes αSyn entry into mitochondria. Our results demonstrate that it is possible to regulate αSyn-VDAC complexation by a rationally designed peptide, thus suggesting new avenues in the search for peptide therapeutics to alleviate αSyn mitochondrial toxicity in PD and other synucleinopathies.

Voltage-activated complexation of α-synuclein with three diverse ß-barrel channels: VDAC, MspA, and α-hemolysin.

Voltage-activated complexation is the process by which a transmembrane potential drives complex formation between a membrane-embedded channel and a soluble or membrane-peripheral target protein. Metabolite and calcium flux across the mitochondrial outer membrane was shown to be regulated by voltage-activated complexation of the voltage-dependent anion channel (VDAC) and either dimeric tubulin or α-synuclein (αSyn). However, the roles played by VDAC's characteristic attributes-its anion selectivity and voltage gating behavior-have remained unclear. Here, we compare in vitro measurements of voltage-activated complexation of αSyn with three well-characterized ß-barrel channels-VDAC, MspA, and α-hemolysin-that differ widely in their organism of origin, structure, geometry, charge density distribution, and voltage gating behavior. The voltage dependences of the complexation dynamics for the different channels are observed to differ quantitatively but have similar qualitative features. In each case, energy landscape modeling describes the complexation dynamics in a manner consistent with the known properties of the individual channels, while voltage gating does not appear to play a role. The reaction free energy landscapes thus calculated reveal a non-trivial dependence of the αSyn/channel complex stability on the surface density of αSyn.

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.

Correction to: Molecular mechanism of olesoxime-mediated neuroprotection through targeting α-synuclein interaction with mitochondrial VDAC.

In the published article, an error was noticed and this has been corrected with this erratum publication.

Molecular mechanism of olesoxime-mediated neuroprotection through targeting α-synuclein interaction with mitochondrial VDAC.

An intrinsically disordered neuronal protein α-synuclein (αSyn) is known to cause mitochondrial dysfunction, contributing to loss of dopaminergic neurons in Parkinson's disease. Through yet poorly defined mechanisms, αSyn crosses mitochondrial outer membrane and targets respiratory complexes leading to bioenergetics defects. Here, using neuronally differentiated human cells overexpressing wild-type αSyn, we show that the major metabolite channel of the outer membrane, the voltage-dependent anion channel (VDAC), is a pathway for αSyn translocation into the mitochondria. Importantly, the neuroprotective cholesterol-like synthetic compound olesoxime inhibits this translocation. By applying complementary electrophysiological and biophysical approaches, we provide mechanistic insights into the interplay between αSyn, VDAC, and olesoxime. Our data suggest that olesoxime interacts with VDAC ß-barrel at the lipid-protein interface thus hindering αSyn translocation through the VDAC pore and affecting VDAC voltage gating. We propose that targeting αSyn translocation through VDAC could represent a key mechanism for the development of new neuroprotective strategies.

Regulation of Mitochondrial Respiration by VDAC Is Enhanced by Membrane-Bound Inhibitors with Disordered Polyanionic C-Terminal Domains.

The voltage-dependent anion channel (VDAC) is the primary regulating pathway of water-soluble metabolites and ions across the mitochondrial outer membrane. When reconstituted into lipid membranes, VDAC responds to sufficiently large transmembrane potentials by transitioning to gated states in which ATP/ADP flux is reduced and calcium flux is increased. Two otherwise unrelated cytosolic proteins, tubulin, and α-synuclein (αSyn), dock with VDAC by a novel mechanism in which the transmembrane potential draws their disordered, polyanionic C-terminal domains into and through the VDAC channel, thus physically blocking the pore. For both tubulin and αSyn, the blocked state is observed at much lower transmembrane potentials than VDAC gated states, such that in the presence of these cytosolic docking proteins, VDAC's sensitivity to transmembrane potential is dramatically increased. Remarkably, the features of the VDAC gated states relevant for bioenergetics-reduced metabolite flux and increased calcium flux-are preserved in the blocked state induced by either docking protein. The ability of tubulin and αSyn to modulate mitochondrial potential and ATP production in vivo is now supported by many studies. The common physical origin of the interactions of both tubulin and αSyn with VDAC leads to a general model of a VDAC inhibitor, facilitates predictions of the effect of post-translational modifications of known inhibitors, and points the way toward the development of novel therapeutics targeting VDAC.

VDAC Gating Thermodynamics, but Not Gating Kinetics, Are Virtually Temperature Independent.

The voltage-dependent anion channel (VDAC) is the most abundant protein in the mitochondrial outer membrane and an archetypical ß-barrel channel. Here, we study the effects of temperature on VDAC channels reconstituted in planar lipid membranes at the single- and multichannel levels within the 20°C to 40°C range. The temperature dependence of conductance measured on a single channel in 1 M KCl shows an increase characterized by a 10°C temperature coefficient Q10 = 1.22 ± 0.02, which exceeds that of the bathing electrolyte solution conductivity, Q10 = 1.17 ± 0.01. The rates of voltage-induced channel transition between the open and closed states measured on multichannel membranes also show statistically significant increases, with temperatures that are consistent with activation energy barriers of ∼10 ± 3 kcal/mol. At the same time, the gating thermodynamics, as characterized by the gating charge and voltage of equipartitioning, does not display any measurable temperature dependence. The two parameters stay within 3.2 ± 0.2 elementary charges and 30 ± 2 mV, respectively. Thus, whereas the channel kinetics, specifically its conductance and rates of gating response to voltage steps, demonstrates a clear increase with temperature, the conformational voltage-dependent equilibria are virtually insensitive to temperature. These results, which may be a general feature of ß-barrel channel gating, suggest either an entropy-driven gating mechanism or a role for enthalpy-entropy compensation.

Structural features and lipid binding domain of tubulin on biomimetic mitochondrial membranes.

Dimeric tubulin, an abundant water-soluble cytosolic protein known primarily for its role in the cytoskeleton, is routinely found to be associated with mitochondrial outer membranes, although the structure and physiological role of mitochondria-bound tubulin are still unknown. There is also no consensus on whether tubulin is a peripheral membrane protein or is integrated into the outer mitochondrial membrane. Here the results of five independent techniques-surface plasmon resonance, electrochemical impedance spectroscopy, bilayer overtone analysis, neutron reflectometry, and molecular dynamics simulations-suggest that α-tubulin's amphipathic helix H10 is responsible for peripheral binding of dimeric tubulin to biomimetic "mitochondrial" membranes in a manner that differentiates between the two primary lipid headgroups found in mitochondrial membranes, phosphatidylethanolamine and phosphatidylcholine. The identification of the tubulin dimer orientation and membrane-binding domain represents an essential step toward our understanding of the complex mechanisms by which tubulin interacts with integral proteins of the mitochondrial outer membrane and is important for the structure-inspired design of tubulin-targeting agents.

Sequence diversity of tubulin isotypes in regulation of the mitochondrial voltage-dependent anion channel.

The microtubule protein tubulin is a heterodimer comprising α/ß subunits, in which each subunit features multiple isotypes in vertebrates. For example, seven α-tubulin and eight ß-tubulin isotypes in the human tubulin gene family vary mostly in the length and primary sequence of the disordered anionic carboxyl-terminal tails (CTTs). The biological reason for such sequence diversity remains a topic of vigorous enquiry. Here, we demonstrate that it may be a key feature of tubulin's role in regulation of the permeability of the mitochondrial outer membrane voltage-dependent anion channel (VDAC). Using recombinant yeast α/ß-tubulin constructs with α-CTTs, ß-CTTs, or both from various human tubulin isotypes, we probed their interactions with VDAC reconstituted into planar lipid bilayers. A comparative study of the blockage kinetics revealed that either α-CTTs or ß-CTTs block the VDAC pore and that the efficiency of blockage by individual CTTs spans 2 orders of magnitude, depending on the CTT isotype. ß-Tubulin constructs, notably ß3, blocked VDAC most effectively. We quantitatively described these experimental results using a physical model that accounted only for the number and distribution of charges in the CTT, and not for the interactions between specific residues on the CTT and VDAC pore. Based on these results, we speculate that the effectiveness of VDAC regulation by tubulin depends on the predominant tubulin isotype in a cell. Consequently, the fluxes of ATP/ADP through the channel could vary significantly, depending on the isotype, thus suggesting an intriguing link between VDAC regulation and the diversity of tubulin isotypes present in vertebrates.

Real-Time Nanopore-Based Recognition of Protein Translocation Success.

A growing number of new technologies are supported by a single- or multi-nanopore architecture for capture, sensing, and delivery of polymeric biomolecules. Nanopore-based single-molecule DNA sequencing is the premier example. This method relies on the uniform linear charge density of DNA, so that each DNA strand is overwhelmingly likely to pass through the nanopore and across the separating membrane. For disordered peptides, folded proteins, or block copolymers with heterogeneous charge densities, by contrast, translocation is not assured, and additional strategies to monitor the progress of the polymer molecule through a nanopore are required. Here, we demonstrate a single-molecule method for direct, model-free, real-time monitoring of the translocation of a disordered, heterogeneously charged polypeptide through a nanopore. The crucial elements are two "selectivity tags"-regions of different but uniform charge density-at the ends of the polypeptide. These affect the selectivity of the nanopore differently and enable discrimination between polypeptide translocation and retraction. Our results demonstrate exquisite sensitivity of polypeptide translocation to applied transmembrane potential and prove the principle that nanopore selectivity reports on biopolymer substructure. We anticipate that the selectivity tag technique will be broadly applicable to nanopore-based protein detection, analysis, and separation technologies, and to the elucidation of protein translocation processes in normal cellular function and in disease.

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.

Evidence of Distinct Channel Conformations and Substrate Binding Affinities for the Mitochondrial Outer Membrane Protein Translocase Pore Tom40.

Nearly all mitochondrial proteins are coded by the nuclear genome and must be transported into mitochondria by the translocase of the outer membrane complex. Tom40 is the central subunit of the translocase complex and forms a pore in the mitochondrial outer membrane. To date, the mechanism it utilizes for protein transport remains unclear. Tom40 is predicted to comprise a membrane-spanning ß-barrel domain with conserved α-helical domains at both the N and C termini. To investigate Tom40 function, including the role of the N- and C-terminal domains, recombinant forms of the Tom40 protein from the yeast Candida glabrata, and truncated constructs lacking the N- and/or C-terminal domains, were functionally characterized in planar lipid membranes. Our results demonstrate that each of these Tom40 constructs exhibits at least four distinct conductive levels and that full-length and truncated Tom40 constructs specifically interact with a presequence peptide in a concentration- and voltage-dependent manner. Therefore, neither the first 51 amino acids of the N terminus nor the last 13 amino acids of the C terminus are required for Tom40 channel formation or for the interaction with a presequence peptide. Unexpectedly, substrate binding affinity was dependent upon the Tom40 state corresponding to a particular conductive level. A model where two Tom40 pores act in concert as a dimeric protein complex best accounts for the observed biochemical and electrophysiological data. These results provide the first evidence for structurally distinct Tom40 conformations playing a role in substrate recognition and therefore in transport function.

α-Synuclein Shows High Affinity Interaction with Voltage-dependent Anion Channel, Suggesting Mechanisms of Mitochondrial Regulation and Toxicity in Parkinson Disease.

Participation of the small, intrinsically disordered protein α-synuclein (α-syn) in Parkinson disease (PD) pathogenesis has been well documented. Although recent research demonstrates the involvement of α-syn in mitochondrial dysfunction in neurodegeneration and suggests direct interaction of α-syn with mitochondria, the molecular mechanism(s) of α-syn toxicity and its effect on neuronal mitochondria remain vague. Here we report that at nanomolar concentrations, α-syn reversibly blocks the voltage-dependent anion channel (VDAC), the major channel of the mitochondrial outer membrane that controls most of the metabolite fluxes in and out of the mitochondria. Detailed analysis of the blockage kinetics of VDAC reconstituted into planar lipid membranes suggests that α-syn is able to translocate through the channel and thus target complexes of the mitochondrial respiratory chain in the inner mitochondrial membrane. Supporting our in vitro experiments, a yeast model of PD shows that α-syn toxicity in yeast depends on VDAC. The functional interactions between VDAC and α-syn, revealed by the present study, point toward the long sought after physiological and pathophysiological roles for monomeric α-syn in PD and in other α-synucleinopathies.

Magic angle spinning nuclear magnetic resonance characterization of voltage-dependent anion channel gating in two-dimensional lipid crystalline bilayers.

The N-terminus of the voltage-dependent anion channel (VDAC) has been proposed to contain the mechanistically important gating helices that modulate channel opening and closing. In this study, we utilize magic angle spinning nuclear magnetic resonance (MAS NMR) to determine the location and structure of the N-terminus for functional channels in lipid bilayers by measuring long-range (13)C-(13)C distances between residues in the N-terminus and other domains of VDAC reconstituted into DMPC lipid bilayers. Our structural studies show that the distance between A14 Cß in the N-terminal helix and S193 Cß is ∼4-6 Å. Furthermore, VDAC phosphorylation by a mitochondrial kinase at residue S193 has been claimed to delay mitochondrial cell death by causing a conformational change that closes the channel, and a VDAC-Ser193Glu mutant has been reported to show properties very similar to those of phosphorylated VDAC in a cellular context. We expressed VDAC-S193E and reconstituted it into DMPC lipid bilayers. Two-dimensional (13)C-(13)C correlation experiments showed chemical shift perturbations for residues located in the N-terminus, indicating possible structural perturbations to that region. However, electrophysiological data recorded on VDAC-S193E showed that channel characteristics were identical to those of wild type samples, indicating that phosphorylation of S193 does not directly affect channel gating. The combination of NMR and electrophysiological results allows us to discuss the validity of proposed gating models.

Acidification asymmetrically affects voltage-dependent anion channel implicating the involvement of salt bridges.

The voltage-dependent anion channel (VDAC) is the major pathway for ATP, ADP, and other respiratory substrates through the mitochondrial outer membrane, constituting a crucial point of mitochondrial metabolism regulation. VDAC is characterized by its ability to "gate" between an open and several "closed" states under applied voltage. In the early stages of tumorigenesis or during ischemia, partial or total absence of oxygen supply to cells results in cytosolic acidification. Motivated by these facts, we investigated the effects of pH variations on VDAC gating properties. We reconstituted VDAC into planar lipid membranes and found that acidification reversibly increases its voltage-dependent gating. Furthermore, both VDAC anion selectivity and single channel conductance increased with acidification, in agreement with the titration of the negatively charged VDAC residues at low pH values. Analysis of the pH dependences of the gating and open channel parameters yielded similar pKa values close to 4.0. We also found that the response of VDAC gating to acidification was highly asymmetric. The presumably cytosolic (cis) side of the channel was the most sensitive to acidification, whereas the mitochondrial intermembrane space (trans) side barely responded to pH changes. Molecular dynamic simulations suggested that stable salt bridges at the cis side, which are susceptible to disruption upon acidification, contribute to this asymmetry. The pronounced sensitivity of the cis side to pH variations found here in vitro might provide helpful insights into the regulatory role of VDAC in the protective effect of cytosolic acidification during ischemia in vivo.

Conductance hysteresis in the voltage-dependent anion channel.

Hysteresis in the conductance of voltage-sensitive ion channels is observed when the transmembrane voltage is periodically varied with time. Although this phenomenon has been used in studies of gating of the voltage-dependent anion channel, VDAC, from the outer mitochondrial membrane for nearly four decades, full hysteresis curves have never been reported, because the focus was solely on the channel opening branches of the hysteresis loops. We studied the hysteretic response of a multichannel VDAC system to a triangular voltage ramp the frequency of which was varied over three orders of magnitude, from 0.5 mHz to 0.2 Hz. We found that in this wide frequency range the area encircled by the hysteresis curves changes by less than a factor of three, suggesting broad distribution of the characteristic times and strongly non-equilibrium behavior. At the same time, quasi-equilibrium two-state behavior is observed for hysteresis branches corresponding to VDAC opening. This enables calculation of the usual equilibrium gating parameters, gating charge and voltage of equipartitioning, which were found to be almost insensitive to the ramp frequency. To rationalize this peculiarity, we hypothesize that during voltage-induced closure and opening the system explores different regions of the complex free energy landscape, and, in the opening branch, follows quasi-equilibrium paths.

Alpha-synuclein lipid-dependent membrane binding and translocation through the α-hemolysin channel.

Gauging the interactions of a natively unfolded Parkinson disease-related protein, alpha-synuclein (α-syn) with membranes and its pathways between and within cells is important for understanding its pathogenesis. Here, to address these questions, we use a robust ß-barrel channel, α-hemolysin, reconstituted into planar lipid bilayers. Transient, ~95% blockage of the channel current by α-syn was observed when 1), α-syn was added from the membrane side where the shorter (stem) part of the channel is exposed; and 2), the applied potential was lower on the side of α-syn addition. While the on-rate of α-syn binding to the channel strongly increased with the applied field, the off-rate displayed a turnover behavior. Statistical analysis suggests that at voltages >50 mV, a significant fraction of the α-syn molecules bound to the channel undergoes subsequent translocation. The observed on-rate varied by >100 times depending on the bilayer lipid composition. Removal of the last 25 amino acids from the highly negatively charged C-terminal of α-syn resulted in a significant decrease in the binding rates. Taken together, these results demonstrate that ß-barrel channels may serve as sensitive probes of α-syn interactions with membranes as well as model systems for studies of channel-assisted protein transport.

Voltage-dependent anion channels modulate mitochondrial metabolism in cancer cells: regulation by free tubulin and erastin.

Respiratory substrates and adenine nucleotides cross the mitochondrial outer membrane through the voltage-dependent anion channel (VDAC), comprising three isoforms--VDAC1, 2, and 3. We characterized the role of individual isoforms in mitochondrial metabolism by HepG2 human hepatoma cells using siRNA. With VDAC3 to the greatest extent, all VDAC isoforms contributed to the maintenance of mitochondrial membrane potential, but only VDAC3 knockdown decreased ATP, ADP, NAD(P)H, and mitochondrial redox state. Cells expressing predominantly VDAC3 were least sensitive to depolarization induced by increased free tubulin. In planar lipid bilayers, free tubulin inhibited VDAC1 and VDAC2 but not VDAC3. Erastin, a compound that interacts with VDAC, blocked and reversed mitochondrial depolarization after microtubule destabilizers in intact cells and antagonized tubulin-induced VDAC blockage in planar bilayers. In conclusion, free tubulin inhibits VDAC1/2 and limits mitochondrial metabolism in HepG2 cells, contributing to the Warburg phenomenon. Reversal of tubulin-VDAC interaction by erastin antagonizes Warburg metabolism and restores oxidative mitochondrial metabolism.