Conductance hysteresis in the voltage-dependent anion channel

Beta-Barrel Channel Response to High Electric Fields: Functional Gating or Reversible Denaturation?

Ion channels exhibit gating behavior, fluctuating between open and closed states, with the transmembrane voltage serving as one of the essential regulators of this process. Voltage gating is a fundamental functional aspect underlying the regulation of ion-selective, mostly α-helical, channels primarily found in excitable cell membranes. In contrast, there exists another group of larger, and less selective, β-barrel channels of a different origin, which are not directly associated with cell excitability. Remarkably, these channels can also undergo closing, or "gating", induced by sufficiently strong electric fields. Once the field is removed, the channels reopen, preserving a memory of the gating process. In this study, we explored the hypothesis that the voltage-induced closure of the β-barrel channels can be seen as a form of reversible protein denaturation by the high electric fields applied in model membranes experiments-typically exceeding twenty million volts per meter-rather than a manifestation of functional gating. Here, we focused on the bacterial outer membrane channel OmpF reconstituted into planar lipid bilayers and analyzed various characteristics of the closing-opening process that support this idea. Specifically, we considered the nearly symmetric response to voltages of both polarities, the presence of multiple closed states, the stabilization of the open conformation in channel clusters, the long-term gating memory, and the Hofmeister effects in closing kinetics. Furthermore, we contemplate the evolutionary aspect of the phenomenon, proposing that the field-induced denaturation of membrane proteins might have served as a starting point for their development into amazing molecular machines such as voltage-gated channels of nerve and muscle cells.

Modulation of Voltage-Gating and Hysteresis of Lysenin Channels by Cu2+ Ions

The intricate voltage regulation presented by lysenin channels reconstituted in artificial lipid membranes leads to a strong hysteresis in conductance, bistability, and memory. Prior investigations on lysenin channels indicate that the hysteresis is modulated by multivalent cations which are also capable of eliciting single-step conformational changes and transitions to stable closed or sub-conducting states. However, the influence on voltage regulation of Cu2+ ions, capable of completely closing the lysenin channels in a two-step process, was not sufficiently addressed. In this respect, we employed electrophysiology approaches to investigate the response of lysenin channels to variable voltage stimuli in the presence of small concentrations of Cu2+ ions. Our experimental results showed that the hysteretic behavior, recorded in response to variable voltage ramps, is accentuated in the presence of Cu2+ ions. Using simultaneous AC/DC stimulation, we were able to determine that Cu2+ prevents the reopening of channels previously closed by depolarizing potentials and the channels remain in the closed state even in the absence of a transmembrane voltage. In addition, we showed that Cu2+ addition reinstates the voltage gating and hysteretic behavior of lysenin channels reconstituted in neutral lipid membranes in which lysenin channels lose their voltage-regulating properties. In the presence of Cu2+ ions, lysenin not only regained the voltage gating but also behaved like a long-term molecular memory controlled by electrical potentials.

Gating of β-Barrel Protein Pores, Porins, and Channels: An Old Problem with New Facets

β barrels are ubiquitous proteins in the outer membranes of mitochondria, chloroplasts, and Gram-negative bacteria. These transmembrane proteins (TMPs) execute a wide variety of tasks. For example, they can serve as transporters, receptors, membrane-bound enzymes, as well as adhesion, structural, and signaling elements. In addition, multimeric β barrels are common structural scaffolds among many pore-forming toxins. Significant progress has been made in understanding the functional, structural, biochemical, and biophysical features of these robust and versatile proteins. One frequently encountered fundamental trait of all β barrels is their voltage-dependent gating. This process consists of reversible or permanent conformational transitions between a large-conductance, highly permeable open state and a low-conductance, solute-restrictive closed state. Several intrinsic molecular mechanisms and environmental factors modulate this universal property of β barrels. This review article outlines the typical signatures of voltage-dependent gating. Moreover, we discuss recent developments leading to a better qualitative understanding of the closure dynamics of these TMPs.

Rufinamide, a Triazole-Derived Antiepileptic Drug, Stimulates Ca2+-Activated K+ Currents While Inhibiting Voltage-Gated Na+ Currents

Rufinamide (RFM) is a clinically utilized antiepileptic drug that, as a triazole derivative, has a unique structure. The extent to which this drug affects membrane ionic currents remains incompletely understood. With the aid of patch clamp technology, we investigated the effects of RFM on the amplitude, gating, and hysteresis of ionic currents from pituitary GH₃ lactotrophs. RFM increased the amplitude of Ca²⁺-activated K⁺ currents (I_K(Ca)) in pituitary GH₃ lactotrophs, and the increase was attenuated by the further addition of iberiotoxin or paxilline. The addition of RFM to the cytosolic surface of the detached patch of membrane resulted in the enhanced activity of large-conductance Ca²⁺-activated K⁺ channels (BK_Ca channels), and paxilline reversed this activity. RFM increased the strength of the hysteresis exhibited by the BK_Ca channels and induced by an inverted isosceles-triangular ramp pulse. The peak and late voltage-gated Na⁺ current (I_Na) evoked by rapid step depolarizations were differentially suppressed by RFM. The molecular docking approach suggested that RFM bound to the intracellular domain of K_Ca1.1 channels with amino acid residues, thereby functionally affecting BK_Ca channels' activity. This study is the first to present evidence that, in addition to inhibiting the I_Na, RFM effectively modifies the I_K(Ca), which suggests that it has an impact on neuronal function and excitability.

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.

Characterization of Inhibitory Capability on Hyperpolarization-Activated Cation Current Caused by Lutein (β,ε-Carotene-3,3′-Diol), a Dietary Xanthophyll Carotenoid

Lutein (β,ε-carotene-3,3′-diol), a xanthophyll carotenoid, is found in high concentrations in the macula of the human retina. It has been recognized to exert potential effectiveness in antioxidative and anti-inflammatory properties. However, whether and how its modifications on varying types of plasmalemmal ionic currents occur in electrically excitable cells remain incompletely answered. The current hypothesis is that lutein produces any direct adjustments on ionic currents (e.g., hyperpolarization-activated cation current, I_h [or funny current, I_f]). In the present study, GH₃-cell exposure to lutein resulted in a time-, state- and concentration-dependent reduction in I_h amplitude with an IC₅₀ value of 4.1 μM. There was a hyperpolarizing shift along the voltage axis in the steady-state activation curve of I_h in the presence of this compound, despite being void of changes in the gating charge of the curve. Under continued exposure to lutein (3 μM), further addition of oxaliplatin (10 μM) or ivabradine (3 μM) could be effective at either reversing or further decreasing lutein-induced suppression of hyperpolarization-evoked I_h, respectively. The voltage-dependent anti-clockwise hysteresis of I_h responding to long-lasting inverted isosceles-triangular ramp concentration-dependently became diminished by adding this compound. However, the addition of 10 μM lutein caused a mild but significant suppression in the amplitude of erg-mediated or A-type K⁺ currents. Under current-clamp potential recordings, the sag potential evoked by long-lasting hyperpolarizing current stimulus was reduced under cell exposure to lutein. Altogether, findings from the current observations enabled us to reflect that during cell exposure to lutein used at pharmacologically achievable concentrations, lutein-perturbed inhibition of I_h would be an ionic mechanism underlying its changes in membrane excitability.

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 Effectiveness in Activating M-Type K+ Current Produced by Solifenacin ([(3R)-1-azabicyclo[2.2.2]octan-3-yl] (1S)-1-phenyl-3,4-dihydro-1H-isoquinoline-2-carboxylate): Independent of Its Antimuscarinic Action

Solifenacin (Vesicare^®, SOL), known to be a member of isoquinolines, is a muscarinic antagonist that has anticholinergic effect, and it has been beneficial in treating urinary incontinence and neurogenic detrusor overactivity. However, the information regarding the effects of SOL on membrane ionic currents is largely uncertain, despite its clinically wide use in patients with those disorders. In this study, the whole-cell current recordings revealed that upon membrane depolarization in pituitary GH₃ cells, the exposure to SOL concentration-dependently increased the amplitude of M-type K⁺ current (I_K(M)) with effective EC₅₀ value of 0.34 μM. The activation time constant of I_K(M) was concurrently shortened in the SOL presence, hence yielding the K_D value of 0.55 μM based on minimal reaction scheme. As cells were exposed to SOL, the steady-state activation curve of I_K(M) was shifted along the voltage axis to the left with no change in the gating charge of the current. Upon an isosceles-triangular ramp pulse, the hysteretic area of I_K(M) was increased by adding SOL. As cells were continually exposed to SOL, further application of acetylcholine (1 μM) failed to modify SOL-stimulated I_K(M); however, subsequent addition of thyrotropin releasing hormone (TRH, 1 μM) was able to counteract SOL-induced increase in I_K(M) amplitude. In cell-attached single-channel current recordings, bath addition of SOL led to an increase in the activity of M-type K⁺ (K_M) channels with no change in the single channel conductance; the mean open time of the channel became lengthened. In whole-cell current-clamp recordings, the SOL application reduced the firing of action potentials (APs) in GH₃ cells; however, either subsequent addition of TRH or linopirdine was able to reverse SOL-mediated decrease in AP firing. In hippocampal mHippoE-14 neurons, the I_K(M) was also stimulated by adding SOL. Altogether, findings from this study disclosed for the first time the effectiveness of SOL in interacting with K_M channels and hence in stimulating I_K(M) in electrically excitable cells, and this noticeable action appears to be independent of its antagonistic activity on the canonical binding to muscarinic receptors expressed in GH₃ or mHippoE-14 cells.

Effective Perturbations on the Amplitude and Hysteresis of Erg-Mediated Potassium Current Caused by 1-Octylnonyl 8-[(2-hydroxyethyl)[6-oxo-6(undecyloxy)hexyl]amino]-octanoate (SM-102), a Cationic Lipid

SM-102 (1-octylnonyl 8-[(2-hydroxyethyl)[6-oxo-6-(undecyloxy)hexyl]amino]-octanoate) is an amino cationic lipid that has been tailored for the formation of lipid nanoparticles and it is one of the essential ingredients present in the Moderna^TM COVID-19 vaccine. However, to what extent it may modify varying types of plasmalemmal ionic currents remains largely uncertain. In this study, we investigate the effects of SM-102 on ionic currents either in two types of endocrine cells (e.g., rat pituitary tumor (GH₃) cells and mouse Leydig tumor (MA-10) cells) or in microglial (BV2) cells. Hyperpolarization-activated K⁺ currents in these cells bathed in high-K⁺, Ca²⁺-free extracellular solution were examined to assess the effects of SM-102 on the amplitude and hysteresis of the erg-mediated K⁺ current (I_K(erg)). The SM-102 addition was effective at blocking I_K(erg) in a concentration-dependent fashion with a half-maximal concentration (IC₅₀) of 108 μM, a value which is similar to the K_D value (i.e., 134 μM) required for its accentuation of deactivation time constant of the current. The hysteretic strength of I_K(erg) in response to the long-lasting isosceles-triangular ramp pulse was effectively decreased in the presence of SM-102. Cell exposure to TurboFectin^TM 8.0 (0.1%, v/v), a transfection reagent, was able to inhibit hyperpolarization-activated I_K(erg) effectively with an increase in the deactivation time course of the current. Additionally, in GH₃ cells dialyzed with spermine (30 μM), the I_K(erg) amplitude progressively decreased; moreover, a further bath application of SM-102 (100 μM) or TurboFectin (0.1%) diminished the current magnitude further. In MA-10 Leydig cells, the I_K(erg) was also blocked by the presence of SM-102 or TurboFectin. The IC₅₀ value for SM-102-induced inhibition of I_K(erg) in MA-10 cells was 98 μM. In BV2 microglial cells, the amplitude of the inwardly rectifying K⁺ current was inhibited by SM-102. Taken together, the presence of SM-102 concentration-dependently inhibited I_K(erg) in endocrine cells (e.g., GH₃ or MA-10 cells), and such action may contribute to their functional activities, assuming that similar in vivo findings exist.

Exploring lipid-dependent conformations of membrane-bound α-synuclein with the VDAC nanopore

Regulation of VDAC by α-synuclein (αSyn) is a rich and instructive example of protein-protein interactions catalyzed by a lipid membrane surface. αSyn, a peripheral membrane protein involved in Parkinson's disease pathology, is known to bind to membranes in a transient manner. αSyn's negatively charged C-terminal domain is then available to be electromechanically trapped by the VDAC β-barrel, a process that is observed in vitro as the reversible reduction of ion flow through a single voltage-biased VDAC nanopore. Binding of αSyn to the lipid bilayer is a prerequisite of the channel-protein interaction; surprisingly, however, we find that the strength of αSyn binding to the membrane does not correlate in any simple way with its efficiency of blocking VDAC, suggesting that the lipid-dependent conformations of the membrane-bound αSyn control the interaction. Quantitative models of the free energy landscape governing the capture and release processes allow us to discriminate between several αSyn (sub-) conformations on the membrane surface. These results, combined with known structural features of αSyn on anionic lipid membranes, point to a model in which the lipid composition determines the fraction of αSyn molecules for which the charged C terminal domain is constrained to be close, but not tightly bound, to the membrane surface and thus readily captured by the VDAC nanopore. We speculate that changes in the mitochondrial membrane lipid composition may be key regulators of the αSyn-VDAC interaction and consequently of VDAC-facilitated transport of ions and metabolites in and out of mitochondria and, i.e. mitochondrial metabolism.

A lower affinity to cytosolic proteins reveals VDAC3 isoform-specific role in mitochondrial biology

Voltage-dependent anion channel (VDAC) is the major pathway for the transport of ions and metabolites across the mitochondrial outer membrane. Among the three known mammalian VDAC isoforms, VDAC3 is the least characterized, but unique functional roles have been proposed in cellular and animal models. Yet, a high-sequence similarity between VDAC1 and VDAC3 is indicative of a similar pore-forming structure. Here, we conclusively show that VDAC3 forms stable, highly conductive voltage-gated channels that, much like VDAC1, are weakly anion selective and facilitate metabolite exchange, but exhibit unique properties when interacting with the cytosolic proteins α-synuclein and tubulin. These two proteins are known to be potent regulators of VDAC1 and induce similar characteristic blockages (on the millisecond time scale) of VDAC3, but with 10- to 100-fold reduced on-rates and altered α-synuclein blocking times, indicative of an isoform-specific function. Through cysteine scanning mutagenesis, we found that VDAC3's cysteine residues regulate its interaction with α-synuclein, demonstrating VDAC3-unique functional properties and further highlighting a general molecular mechanism for VDAC isoform-specific regulation of mitochondrial bioenergetics.

Tunable Electromechanical Nanopore Trap Reveals Populations of Peripheral Membrane Protein Binding Conformations

We demonstrate that a naturally occurring nanopore, the voltage-dependent anion channel (VDAC) of the mitochondrion, can be used to electromechanically trap and interrogate proteins bound to a lipid surface at the single-molecule level. Electromechanically probing α-synuclein (αSyn), an intrinsically disordered neuronal protein intimately associated with Parkinson's pathology, reveals wide variation in the time required for individual proteins to unbind from the same membrane surface. The observed distributions of unbinding times span up to 3 orders of magnitude and depend strongly on the lipid composition of the membrane; surprisingly, lipid membranes to which αSyn binds weakly are most likely to contain subpopulations in which electromechanically driven unbinding is very slow. We conclude that unbinding of αSyn from the membrane surface depends not only on membrane binding affinity but also on the conformation adopted by an individual αSyn molecule on the membrane surface.

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 Q₁₀ = 1.22 ± 0.02, which exceeds that of the bathing electrolyte solution conductivity, Q₁₀ = 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.

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.

Targeting the Multiple Physiologic Roles of VDAC With Steroids and Hydrophobic Drugs

There is accumulating evidence that endogenous steroids and non-polar drugs are involved in the regulation of mitochondrial physiology. Many of these hydrophobic compounds interact with the Voltage Dependent Anion Channel (VDAC). This major metabolite channel in the mitochondrial outer membrane (MOM) regulates the exchange of ions and water-soluble metabolites, such as ATP and ADP, across the MOM, thus governing mitochondrial respiration. Proteomics and biochemical approaches together with molecular dynamics simulations have identified an impressively large number of non-polar compounds, including endogenous, able to bind to VDAC. These findings have sparked speculation that both natural steroids and synthetic hydrophobic drugs regulate mitochondrial physiology by directly affecting VDAC ion channel properties and modulating its metabolite permeability. Here we evaluate recent studies investigating the effect of identified VDAC-binding natural steroids and non-polar drugs on VDAC channel functioning. We argue that while many compounds are found to bind to the VDAC protein, they do not necessarily affect its channel functions in vitro. However, they may modify other aspects of VDAC physiology such as interaction with its cytosolic partner proteins or complex formation with other mitochondrial membrane proteins, thus altering mitochondrial function.

Multiple neurosteroid and cholesterol binding sites in voltage-dependent anion channel-1 determined by photo-affinity labeling

Voltage-dependent anion channel-1 (VDAC1) is a mitochondrial porin that is implicated in cellular metabolism and apoptosis, and modulated by numerous small molecules including lipids. VDAC1 binds sterols, including cholesterol and neurosteroids such as allopregnanolone. Biochemical and computational studies suggest that VDAC1 binds multiple cholesterol molecules, but photolabeling studies have identified only a single cholesterol and neurosteroid binding site at E73. To identify all the binding sites of neurosteroids in VDAC1, we apply photo-affinity labeling using two sterol-based photolabeling reagents with complementary photochemistry: 5α-6-AziP which contains an aliphatic diazirine, and KK200 which contains a trifluoromethyl-phenyldiazirine (TPD) group. 5α-6-AziP and KK200 photolabel multiple residues within an E73 pocket confirming the presence of this site and mapping sterol orientation within this pocket. In addition, KK200 photolabels four other sites consistent with the finding that VDAC1 co-purifies with five cholesterol molecules. Both allopregnanolone and cholesterol competitively prevent photolabeling at E73 and three other sites indicating that these are common sterol binding sites shared by both neurosteroids and cholesterol. Binding at the functionally important residue E73 suggests a possible role for sterols in regulating VDAC1 signaling and interaction with partner proteins.

Lipid Headgroup Charge and Acyl Chain Composition Modulate Closure of Bacterial β-Barrel Channels

The outer membrane of Gram-negative bacteria contains β-barrel proteins that form high-conducting ion channels providing a path for hydrophilic molecules, including antibiotics. Traditionally, these proteins have been considered to exist only in an open state so that regulation of outer membrane permeability was accomplished via protein expression. However, electrophysiological recordings show that β-barrel channels respond to transmembrane voltages by characteristically switching from a high-conducting, open state, to a so-called 'closed' state, with reduced permeability and possibly exclusion of large metabolites. Here, we use the bacterial porin OmpF from E. coli as a model system to gain insight on the control of outer membrane permeability by bacterial porins through the modulation of their open state. Using planar bilayer electrophysiology, we perform an extensive study of the role of membrane lipids in the OmpF channel closure by voltage. We pay attention not only to the effects of charges in the hydrophilic lipid heads but also to the contribution of the hydrophobic tails in the lipid-protein interactions. Our results show that gating kinetics is governed by lipid characteristics so that each stage of a sequential closure is different from the previous one, probably because of intra- or intermonomeric rearrangements.

Assessing the role of residue E73 and lipid headgroup charge in VDAC1 voltage gating

The voltage-dependent anion channel (VDAC) is the most abundant protein of the mitochondrial outer membrane (MOM) where it regulates transport of ions and metabolites in and out of the organelle. VDAC function is extensively studied in a lipid bilayer system that allows conductance monitoring of reconstituted channels under applied voltage. The process of switching from a high-conductance state, open to metabolites, to a variety of low-conducting states, which excludes metabolite transport, is termed voltage gating and the mechanism remains poorly understood. Recent studies have implicated the involvement of the membrane-solvated residue E73 in the gating process through β-barrel destabilization. However, there has been no direct experimental evidence of E73 involvement in VDAC1 voltage gating. Here, using electrophysiology measurements, we exclude the involvement of E73 in murine VDAC1 (mVDAC1) voltage gating process. With an established protocol of assessing voltage gating of VDACs reconstituted into planar lipid membranes, we definitively show that mVDAC1 gating properties do not change when E73 is replaced by either a glutamine or an alanine. We further demonstrate that cholesterol has no effect on mVDAC1 gating characteristics, though it was shown that E73 is coordinating residue in the cholesterol binding site. In contrast, we found a pronounced gating effect based on the charge of the phospholipid headgroup, where the positive charge stimulates and negative charge suppresses gating. These findings call for critical evaluation of the existing models of VDAC gating and contribute to our understanding of VDAC's role in control of MOM permeability and regulation of mitochondrial respiration and metabolism.

Structure of voltage-dependent anion channel-tethered bilayer lipid membranes determined using neutron reflectivity

Neutron reflectivity (NR) has emerged as a powerful technique to study the structure and behavior of membrane proteins at planar lipid interfaces. Integral membrane proteins (IMPs) remain a significant challenge for NR owing to the difficulty of forming complete bilayers with sufficient protein density for scattering techniques. One strategy to achieve high protein density on a solid substrate is the capture of detergent-stabilized, affinity-tagged IMPs on a nitrilotriacetic acid (NTA)-functionalized self-assembled monolayer (SAM), followed by reconstitution into the lipids of interest. Such protein-tethered bilayer lipid membranes (ptBLMs) have the notable advantage of a uniform IMP orientation on the substrate. Here, NR is used to provide a structural characterization of the ptBLM process from formation of the SAM to capture of the detergent-stabilized IMP and lipid reconstitution. The mitochondrial outer-membrane voltage-dependent anion channel (VDAC), which controls the exchange of bioenergetic metabolites between mitochondria and the cytosol, was used as a model β-barrel IMP. Molecular dynamics simulations were used for comparison with the experimental results and to inform the parameters of the physical models describing the NR data. The detailed structure of the SAM is shown to depend on the density of the NTA chelating groups. The relative content of detergent and protein in surface-immobilized, detergent-stabilized VDAC is measured, while the reconstituted lipid bilayer is shown to be complete to within a few percent, using the known atomic structure of VDAC. Finally, excess lipid above the reconstituted bilayer, which is of consequence for more indirect structural and functional studies, is shown to be present.

Insights into the Voltage Regulation Mechanism of the Pore-Forming Toxin Lysenin

Lysenin, a pore forming toxin (PFT) extracted from Eisenia fetida, inserts voltage-regulated channels into artificial lipid membranes containing sphingomyelin. The voltage-induced gating leads to a strong static hysteresis in conductance, which endows lysenin with molecular memory capabilities. To explain this history-dependent behavior, we hypothesized a gating mechanism that implies the movement of a voltage domain sensor from an aqueous environment into the hydrophobic core of the membrane under the influence of an external electric field. In this work, we employed electrophysiology approaches to investigate the effects of ionic screening elicited by metal cations on the voltage-induced gating and hysteresis in conductance of lysenin channels exposed to oscillatory voltage stimuli. Our experimental data show that screening of the voltage sensor domain strongly affects the voltage regulation only during inactivation (channel closing). In contrast, channel reactivation (reopening) presents a more stable, almost invariant voltage dependency. Additionally, in the presence of anionic Adenosine 5′-triphosphate (ATP), which binds at a different site in the channel’s structure and occludes the conducting pathway, both inactivation and reactivation pathways are significantly affected. Therefore, the movement of the voltage domain sensor into a physically different environment that precludes electrostatically bound ions may be an integral part of the gating mechanism.

Computational membrane biophysics: From ion channel interactions with drugs to cellular function

The rapid development of experimental and computational techniques has changed fundamentally our understanding of cellular-membrane transport. The advent of powerful computers and refined force-fields for proteins, ions, and lipids has expanded the applicability of Molecular Dynamics (MD) simulations. A myriad of cellular responses is modulated through the binding of endogenous and exogenous ligands (e.g. neurotransmitters and drugs, respectively) to ion channels. Deciphering the thermodynamics and kinetics of the ligand binding processes to these membrane proteins is at the heart of modern drug development. The ever-increasing computational power has already provided insightful data on the thermodynamics and kinetics of drug-target interactions, free energies of solvation, and partitioning into lipid bilayers for drugs. This review aims to provide a brief summary about modeling approaches to map out crucial binding pathways with intermediate conformations and free-energy surfaces for drug-ion channel binding mechanisms that are responsible for multiple effects on cellular functions. We will discuss post-processing analysis of simulation-generated data, which are then transformed to kinetic models to better understand the molecular underpinning of the experimental observables under the influence of drugs or mutations in ion channels. This review highlights crucial mathematical frameworks and perspectives on bridging different well-established computational techniques to connect the dynamics and timescales from all-atom MD and free energy simulations of ion channels to the physiology of action potentials in cellular models. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

Nonequilibrium response of a voltage gated sodium ion channel and biophysical characterization of dynamic hysteresis

Here we have studied the dynamic as well as the non-equilibrium thermodynamic response properties of voltage-gated Na-ion channel. Using sinusoidally oscillating external voltage protocol we have both kinetically and energetically studied the non-equilibrium steady state properties of dynamic hysteresis in details. We have introduced a method of estimating the work done associated with the dynamic memory due to a cycle of oscillating voltage. We have quantitatively characterised the loop area of ionic current which gives information about the work done to sustain the dynamic memory only for ion conduction, while the loop area of total entropy production rate gives the estimate of work done for overall gating dynamics. The maximum dynamic memory of Na-channel not only depends on the frequency and amplitude but it also depends sensitively on the mean of the oscillating voltage and here we have shown how the system optimize the dynamic memory itself in the biophysical range of field parameters. The relation between the average ionic current with increasing frequency corresponds to the nature of the average dissipative work done at steady state. It is also important to understand that the utilization of the energy from the external field can not be directly obtained only from the measurement of ionic current but also requires nonequilibrium thermodynamic study.

Impact of Dendrimer Terminal Group Chemistry on Blockage of the Anthrax Toxin Channel: A Single Molecule Study

Nearly all the cationic molecules tested so far have been shown to reversibly block K⁺ current through the cation-selective PA₆₃ channels of anthrax toxin in a wide nM-mM range of effective concentrations. A significant increase in channel-blocking activity of the cationic compounds was achieved when multiple copies of positively charged ligands were covalently linked to multivalent scaffolds, such as cyclodextrins and dendrimers. Even though multivalent binding can be strong when the individual bonds are relatively weak, for drug discovery purposes we often strive to design multivalent compounds with high individual functional group affinity toward the respective binding site on a multivalent target. Keeping this requirement in mind, here we perform a single-channel/single-molecule study to investigate kinetic parameters of anthrax toxin PA₆₃ channel blockage by second-generation (G2) poly(amido amine) (PAMAM) dendrimers functionalized with different surface ligands, including G2-NH₂, G2-OH, G2-succinamate, and G2-COONa. We found that the previously reported difference in IC₅₀ values of the G2-OH/PA₆₃ and G2-NH₂/PA₆₃ binding was determined by both on- and off-rates of the reversible dendrimer/channel binding reaction. In 1 M KCl, we observed a decrease of about three folds in k o n and a decrease of only about ten times in t r e s with G2-OH compared to G2-NH₂. At the same time for both blockers, k o n and t r e s increased dramatically with transmembrane voltage increase. PAMAM dendrimers functionalized with negatively charged succinamate, but not carboxyl surface groups, still had some residual activity in inhibiting the anthrax toxin channels. At 100 mV, the on-rate of the G2-succinamate binding was comparable with that of G2-OH but showed weaker voltage dependence when compared to G2-OH and G2-NH₂. The residence time of G2-succinamate in the channel exhibited opposite voltage dependence compared to G2-OH and G2-NH₂, increasing with the cis-negative voltage increase. We also describe kinetics of the PA₆₃ ion current modulation by two different types of the "imperfect" PAMAM dendrimers, the mixed-surface G2 75% OH 25% NH₂ dendrimer and G3-NH₂ dendron. At low voltages, both "imperfect" dendrimers show similar rate constants but significantly weaker voltage sensitivity when compared with the intact G2-NH₂ PAMAM dendrimer.

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.

Modeling magnetosensitive ion channels in the viscoelastic environment of living cells

We propose and study a model of hypothetical magnetosensitive ionic channels which are long thought to be a possible candidate to explain the influence of weak magnetic fields on living organisms ranging from magnetotactic bacteria to fishes, birds, rats, bats, and other mammals including humans. The core of the model is provided by a short chain of magnetosomes serving as a sensor, which is coupled by elastic linkers to the gating elements of ion channels forming a small cluster in the cell membrane. The magnetic sensor is fixed by one end on cytoskeleton elements attached to the membrane and is exposed to viscoelastic cytosol. Its free end can reorient stochastically and subdiffusively in viscoelastic cytosol responding to external magnetic field changes and can open the gates of coupled ion channels. The sensor dynamics is generally bistable due to bistability of the gates which can be in two states with probabilities which depend on the sensor orientation. For realistic parameters, it is shown that this model channel can operate in the magnetic field of Earth for a small number (five to seven) of single-domain magnetosomes constituting the sensor rod, each of which has a typical size found in magnetotactic bacteria and other organisms or even just one sufficiently large nanoparticle of a characteristic size also found in nature. It is shown that, due to the viscoelasticity of the medium, the bistable gating dynamics generally exhibits power law and stretched exponential distributions of the residence times of the channels in their open and closed states. This provides a generic physical mechanism for the explanation of the origin of such anomalous kinetics for other ionic channels whose sensors move in a viscoelastic environment provided by either cytosol or biological membrane, in a quite general context, beyond the fascinating hypothesis of magnetosensitive ionic channels we explore.