Pulsed Electric Fields Can Create Pores in the Voltage Sensors of Voltage-Gated Ion Channels

Biophysica Biophysical Society Article Pulsed Electric Fields Can Create Pores in the Voltage Sensors of Voltage-Gated lon Channels ABSTpACTnpeoeccieasnanocreaengyluednmecneotansentyncaeasahecolmembranepemeabthyva that membrane proteinsare affected as well,particularly ed ion char els(VGICs).The m olecular mechanisms s by take exposing them to fields mimicking electropora that fieldscaninduc pores in domains(VDs)of s and tha CaiedcompiexporesnthSainp anticpate tha such SDs become dysfunctional and unabletorespond ochan es in transmembrane voltage.which is in agree ent tra 1e0 eaoyeeaaoewhakaeetacaEgahaegnscaeae GICs by SIGNIFICANCE Pulsed electric fields are often used for treatment of excitable cells,e.g.,for gene delivery into skeletal muscles.ablation of the heart muscle.or brain tumors.Voltage-gated ion chann elsGiCs)undertiegenerationand the molecular mec ns by whic VGIC estions that have been a by ele hat make the prone to protins.thus future investigations.Finally.we propose a mechanismor ong-ived memrane permeability after pulse treatment,which to date remains poorly understoo INTRODUCTION of ger nhenomenon called electr The integrity of the cell membrane.although essential for the life of any biological cell,presents a barrier that needs tion.Thanks to insights from molecular dynamics (MD) we now und rstand that on to be trans cules 10 in medicine to achie increase in cell membrane permeability.Examples are However.experimental evidence suggests that membrane Proteins.particularly votagon channels GICs) d b drugs and gene therapy techniques for intracel are cla of tran voltage (TMV)with ch Submined November 1.2019.and aceepted for publication May 15.2020. rangements that lead to opening or closure of an ion-selec tive pore.They play crucial roles in the generation and Co espondence:.se propagation c acton potentials in elect rically excitable 16bp2020.05.030 This is an open ess article under the CC BY-NC-ND icese(http/ creanivecommons.org/licenses/by-nc-nd/4.0/) 190 Biophysical Joumal 119.190-205.July 7.2020

Article Pulsed Electric Fields Can Create Pores in the Voltage Sensors of Voltage-Gated Ion Channels Lea Rems,1 Marina A. Kasimova,1 Ilaria Testa,1 and Lucie Delemotte1, * 1 Science for Life Laboratory, Department of Applied Physics, KTH Royal Institute of Technology, Solna, Sweden ABSTRACT Pulsed electric fields are increasingly used in medicine to transiently increase the cell membrane permeability via electroporation to deliver therapeutic molecules into the cell. One type of event that contributes to this increase in membrane permeability is the formation of pores in the membrane lipid bilayer. However, electrophysiological measurements suggest that membrane proteins are affected as well, particularly voltage-gated ion channels (VGICs). The molecular mechanisms by which the electric field could affects these molecules remain unidentified. In this study, we used molecular dynamics simulations to unravel the molecular events that take place in different VGICs when exposing them to electric fields mimicking electropora￾tion conditions. We show that electric fields can induce pores in the voltage-sensor domains (VSDs) of different VGICs and that these pores form more easily in some channels than in others. We demonstrate that poration is more likely in VSDs that are more hydrated and are electrostatically more favorable for the entry of ions. We further show that pores in VSDs can expand into so￾called complex pores, which become stabilized by lipid headgroups. Our results suggest that such complex pores are consid￾erably more stable than conventional lipid pores, and their formation can lead to severe unfolding of VSDs from the channel. We anticipate that such VSDs become dysfunctional and unable to respond to changes in transmembrane voltage, which is in agree￾ment with previous electrophysiological measurements showing a decrease in the voltage-dependent transmembrane ionic cur￾rents after pulse treatment. Finally, we discuss the possibility of activation of VGICs by submicrosecond-duration pulses. Overall, our study reveals a new, to our knowledge, mechanism of electroporation through membranes containing VGICs. INTRODUCTION The integrity of the cell membrane, although essential for the life of any biological cell, presents a barrier that needs to be transiently disrupted to deliver therapeutic molecules into the cell. High-intensity pulsed electric fields are increasingly used in medicine to achieve such a transient increase in cell membrane permeability. Examples are cancer treatment for enhanced delivery of chemotherapeutic drugs and gene therapy techniques for intracellular delivery of genetic material (1). The applied electric field induces a phenomenon called electroporation or electropermeabiliza￾tion. Thanks to insights from molecular dynamics (MD) simulations, we now understand that one type of events that takes place in the cell membrane is the formation of pores in the membrane lipid bilayer (2,3). However, experimental evidence suggests that membrane proteins, particularly voltage-gated ion channels (VGICs), could be affected as well (4). VGICs are a class of trans￾membrane proteins that respond to changes in the transmembrane voltage (TMV) with conformational rear￾rangements that lead to opening or closure of an ion-selec￾tive pore. They play crucial roles in the generation and propagation of action potentials in electrically excitable Submitted November 11, 2019, and accepted for publication May 15, 2020. *Correspondence: lucied@kth.se Editor: Vasanthi Jayaraman. SIGNIFICANCE Pulsed electric fields are often used for treatment of excitable cells, e.g., for gene delivery into skeletal muscles, ablation of the heart muscle, or brain tumors. Voltage-gated ion channels (VGICs) underlie generation and propagation of action potentials in these cells, and consequently are essential for their proper function. Our study reveals the molecular mechanisms by which pulsed electric fields directly affect VGICs and addresses questions that have been previously opened by electrophysiologists. We analyze the characteristics of VGICs that make them prone to electroporation, including hydration and electrostatic properties. This analysis is easily transferable to other membrane proteins, thus opening directions for future investigations. Finally, we propose a mechanism for long-lived membrane permeability after pulse treatment, which to date remains poorly understood. 190 Biophysical Journal 119, 190–205, July 7, 2020 https://doi.org/10.1016/j.bpj.2020.05.030
2020 Biophysical Society. This is an open access article under the CC BY-NC-ND license (http:// creativecommons.org/licenses/by-nc-nd/4.0/)

Mechanistic Insights into the Modulation cells,includingneu ons and muscle cells.All VGICs share calcium-independent dow ulation of the Na on archite 一amsncaehothcogPoaae els（eg. caused by alteration of the lipid hilaver Decrease in Nav channel current was also observed by between segments S5 and S6.Segments S1-S4 act as the voltage s r.wher S6.and the d :as ns. m pulses.The esidues and has the not due to a change to respond to changes in TMV.The movement of S4 then or activation of the Nav channels but instead was associ conductance.This ing on the di the e tc n ative effect of the alec c field on Nay cha hundreds of millivolts,far beyond the physiological observed on ky channels.whereas a decrease in conduc resting voltage or voltage generated during action poten tance could also be observed for Cav channels and/or cal ele b cium-depen advanced patch-clamp techniaue rted when using longer millisecond pulses.Chen et al. have demonstrated that high-intensity pulses with submi crosecond duration can decrease the ionic currents medi (10 m rectifier Kv currents after subjecting voltage- clamped frog ated by 20 on all t to hyperpolan g 4 ms Nay and cay channels usin (GH3)and nded on the pulse amplitude,with lower murine neuroblastoma-rat glioma hybrid (NG10)cells. sulting in milder decrease in channel currents Electrophysiological measurements led a decrease Importantly.all 15-2d for pu ampl in me (50 in gested that the decrease in Nay current was not mediated Anart from decreased current through vGlcs ther by Na leakage across the electropermeabilized mem have been other electric-field effects observed as well.Sub- rane or downregulation of the Nav channels by a cal nd pulses have been shown toactivate specifi it remain of th 61s2s3 in thi 51-56 (h)Side vicw lular view of th NavMs HCN1 with S4-S box of th s rib orus toms of lipid h Biophysical Journal 119.190-205.July 7.2020 191

cells, including neurons and muscle cells. All VGICs share a common architecture: each of the four protein domains con￾tains six transmembrane segments (S1–S6) and a pore loop between segments S5 and S6. Segments S1–S4 act as the voltage sensor, whereas S5, S6, and the pore loop serve as the pore-forming module (Fig. 1, a and b; (5)). Segment S4 contains positively charged residues and has the ability to respond to changes in TMV. The movement of S4 then acts on S5 and S6 to open or close the channel pore depend￾ing on the direction of the electric field. Because VGICs are sensitive to changes in TMV and because electroporative pulses induce a TMV of several hundreds of millivolts, far beyond the physiological resting voltage or voltage generated during action poten￾tials, one can speculate that these channels become per￾turbed by pulsed electric fields. Indeed, by using advanced patch-clamp techniques, electrophysiologists have demonstrated that high-intensity pulses with submi￾crosecond duration can decrease the ionic currents medi￾ated by different VGICs during action potentials. Nesin et al. (6) studied the effects of 300 or 600 ns pulses on Nav and Cav channels using murine pituitary (GH3) and murine neuroblastoma-rat glioma hybrid (NG108) cells. Electrophysiological measurements revealed a decrease in Nav and Cav currents for pulse amplitudes above 1.5–2 kV/cm. The results from a follow-up study (7) sug￾gested that the decrease in Nav current was not mediated by Naþ leakage across the electropermeabilized mem￾brane or downregulation of the Nav channels by a cal￾cium-dependent mechanism. The authors thus proposed as the mechanism either an electroconformational change or a calcium-independent downregulation of the Nav channels (e.g., caused by alteration of the lipid bilayer). Decrease in Nav channel current was also observed by Yang et al. (8) when they exposed adrenal chromaffin cells to even shorter 5 ns, 50–100 kV/cm pulses. Their analysis suggested that the decrease in Nav current was not due to a change in either the steady-state inactivation or activation of the Nav channels but instead was associ￾ated with a decrease in maximal Naþ conductance. This decrease could be observed immediately after the pulse (within 0.5 s, earliest time measured), suggesting a direct effect of the electric field on Nav channels. No effect was observed on Kv channels, whereas a decrease in conduc￾tance could also be observed for Cav channels and/or cal￾cium-dependent potassium channels. Although less explored, a decrease in VGIC current has been also re￾ported when using longer millisecond pulses. Chen et al. (9,10) measured up to 40% decrease in Nav and delayed rectifier Kv currents after subjecting voltage-clamped frog skeletal muscle cells to hyperpolarizing 4 ms, 0.5 V pulses. In all the above-discussed studies, the effects de￾pended on the pulse amplitude, with lower amplitude re￾sulting in milder decrease in channel currents. Importantly, all studies also reported that the decrease was not reversible within the observation time (50 min in (10), 10–15 min in (6,8)). Apart from decreased current through VGICs, there have been other electric-field effects observed as well. Sub￾microsecond pulses have been shown to activate specific VGICs (11,12). However, it remains unclear whether the electric field during such short pulses is able to directly + + + + S1 S2 S3 S5 S6 S4 VSD pore NavMs NavPaS HCN1 a b c FIGURE 1 Structure of the voltage-gated ion channels (VGICs) investigated in this study. (a) All VGICs share an overall architecture that includes four subunits, each containing six membrane-span￾ning segments denoted S1–S6. (b) Side view and extracellular view of the NavMs channel, indicating the segments of a single subunit: VSD (S1–S4) colored in red, S5 together with S4–S5 linker colored in green, pore helix colored in blue, and S6 colored in ochre. (c) Simulation box of the investigated sys￾tems: NavMs, NavPaS, and HCN1. Ion channels are shown as ribbons colored as in (a). Water mole￾cules are represented as transparent volume, whereas phosphorus atoms of lipid headgroups, Naþ ions, and Cl ions are shown as gold, yellow, and green spheres, respectively. To see this figure in color, go online. Mechanistic Insights into the Modulation Biophysical Journal 119, 190–205, July 7, 2020 191

Rems eta move the VSDs or whether the channel activation is indi- propensity to be porated with their biophysical rectly triggered by postpulse membrane depolarization characteristics. which is cat sed by no selective ion leakage cross the ele ane.Furne METHODS (157).both of whichare important signals that can initiate Systems preparation The ns with NavMs (Protein Data Bank,PDB:5HVX .depend ing on esting for new applications in wound healing and tissue 150 mM NaCl solution.The CHARMM36 force field (27,2)was engineering. nd it mamembrane.as well as lipid scrambling.which could MD simulations all affect the function of membrane proteins(2).Yet.elec GROMACS 2016(3 using the The trophysiological measurements suggest that pul VGICS Hoo measurements not able to directly det ne whethe 35).D ring the firs 100 and how pulsed electric fields can affect VGICs on the mo 40 lecular le the sing the pe electroporation conditions.Because exp mental studies consistently reported a decrease in Nav cun we per 2 n equ formed simulations of two ed by channels for which x-ray cry Magn cr NavMs)020)and an u Nav from Periplaneta americana(NavPaS)(21).Eukary TMV Na ne composed of single ng Navs,on the other hand,are compr sed of four identical d by apply domains,each being analogous to a single domain of their is li P the act vated cyclic nucleotide-gated (HCN1 channel (24) which is a ctive channe ha TM under hy olarizing TMV.The effects o HCNI interes giving the rapid developm ent of pulsed-elec field-based eardiac ablation for treatment of heart arrhyth that the th spono the nble but did ferences resulting in considerably different hydration and electrostatic profile along their VSD ring si harge i resp nabledu to gain a )The ch nd to harge imng chieved b 192 Biophysical Joumal 119.190-205.July 7.2020

move the VSDs or whether the channel activation is indi￾rectly triggered by postpulse membrane depolarization, which is caused by nonselective ion leakage across the elec￾tropermeabilized membrane. Furthermore, this activation of VGICs has been found to participate in postpulse membrane depolarization (13,14) and intracellular calcium increase (15–17), both of which are important signals that can initiate cell death, proliferation, or differentiation, depending on their spatiotemporal profile (18,19). Understanding the ef￾fects of pulsed electric fields on VGICs is thus also inter￾esting for new applications in wound healing and tissue engineering. Electropermeabilization is associated with a complex set of events, including oxidative lipid damage, disruption of the cytoskeleton network, and its association with the plasma membrane, as well as lipid scrambling, which could all affect the function of membrane proteins (2). Yet, elec￾trophysiological measurements suggest that pulsed electric fields induce some electroconformational change of VGICs during the pulse. However, electrophysiological measurements are not able to directly determine whether and how pulsed electric fields can affect VGICs on the mo￾lecular level. To fill this gap, we used MD simulations and unraveled the molecular events that take place in different VGICs when exposing them to an electric field that mimics electroporation conditions. Because experimental studies consistently reported a decrease in Nav current, we per￾formed simulations of two Nav channels for which the structure had been resolved by x-ray crystallography and cryo-electron microscopy: a bacterial Nav from Magnetococcus marinus (NavMs) (20) and an eukaryotic Nav from Periplaneta americana (NavPaS) (21). Eukary￾otic Nav channels are composed of a single polypep￾tide containing four homologous but nonidentical domains connected by intracellular linkers; bacterial Navs, on the other hand, are comprised of four identical domains, each being analogous to a single domain of their eukaryotic counterparts, and can thus be used as simple models of eukaryotic channels (22,23). In addition, we per￾formed simulations of the human hyperpolarization￾activated cyclic nucleotide-gated (HCN1) channel (24), which is a nonselective voltage-gated cation channel that is responsible for generation of rhythmic activity in heart and brain. Unlike Nav channels, HCN1 activates under hyperpolarizing TMV. The effects on HCN1 are of interest giving the rapid development of pulsed-electric- field-based cardiac ablation for treatment of heart arrhyth￾mias (25). We found that the three tested channels, NavMs, NavPaS, and HCN1, responded differently to electric pulses, which we could relate to the channels’ structural dif￾ferences resulting in considerably different hydration and electrostatic profiles along their VSDs. Comparing their response thus enabled us to gain an atomistic level insight into the biophysical mechanisms governing their interaction with an electric field and to relate their propensity to be porated with their biophysical characteristics. METHODS Systems preparation The computational systems with NavMs (Protein Data Bank, PDB: 5HVX), NavPaS (PDB: 5X0M), and HCN1 (PDB: 5U6O) were built using the CHARMM-GUI webserver (26). Briefly, each protein was embedded into a 1-palmytoyl-2-oleoyl-phosphatidylcholine bilayer and solvated with 150 mM NaCl solution. The CHARMM36 force field (27,28) was used for proteins, lipids, and ions and the TIP3P model (29) for water. The composition of each system is reported in Table S1. MD simulations All simulations were performed in GROMACS 2016 (30,31). Each system was first minimized using the steepest descent algorithm. The equilibra￾tion and production run was then carried out using a leap-frog integrator with a time step of 2.0 fs, Nos
e-Hoover thermostat (t ¼ 0.4 ps, T ¼ 300 K) (32,33), and Parrinello-Rahman barostat (t ¼ 5 ps, P ¼ 1 bar, semi￾isotropic coupling) (34,35). During the first 100 ns, the protein atoms were restrained to their initial positions, after which the simulations were continued for 400–1200 ns without restraints (400 ns for NavMs, 1200 ns for NavPaS, and 700 ns for HCN1). The long-range electrostatic interactions were calculated using the particle mesh Ewald method (36), together with a Fourier grid spacing of 0.15 nm and a cutoff of 1.2 nm. A switching function was used between 0.8 and 1.2 nm to smoothly bring the short-range electrostatic interactions and the van der Waals interac￾tions to 0 at 1.2 nm. The chemical bonds were constrained to their equi￾librium values using the LINCS algorithm (37). Periodic boundary conditions were applied. To simulate the exposure to an electric pulse, we added a force qEz to every atom carrying a charge q (2,3). This method induces a TMV that is approximately equal to the product of the imposed electric field and the simulation box length, TMV ¼ Ez Lz. The Ez was chosen based on the following considerations. Previous simulations on pure lipid bilayers showed that formation of a lipid pore occurs faster with increasing elec￾tric-field magnitude; to observe lipid pores within few tens of nanoseconds, simulations have been typically conducted by applying an electric field re￾sulting in TMV above 3V(38,39). Experimentally, however, the TMV that can be built on the membrane is limited because the membrane starts discharging through pores after it becomes electroporated. According to measurements performed using pulses with duration of 60 ns and longer together with voltage-sensitive dyes (40,41) and microelectrodes (42), the TMV of 1.5 V is the largest that a cell membrane can sustain before dis￾charging. Thus, we chose the value of the electric field such that it resulted in a TMV of about 51.5 V. The chosen value is a compromise between trying to observe an effect of the electric field within reasonable simulation time and at the same time staying within realistic experimental TMV values. These simulations were performed in an NVT ensemble (constant number of atoms, volume, and temperature), which kept the simulation box size constant. This ensured that the voltage applied across the mem￾brane remained constant. We also performed additional control simulations on the NavMs channel by using an electric field in an NPT ensemble but did not observe any considerable differences from NVT simulations that would affect the results and conclusions presented in the work (Supporting Materials and Methods, Section S4.1; Figs. S6 and S10; Table S4). In addition, we performed con￾trol simulations by applying a constant charge imbalance instead of an external electric field (Supporting Materials and Methods, Section S4.2; Fig. S7; Table S5). The charge imbalance was achieved by 1) separating the water bath into two parts by adding a vacuum layer and 2) transferring Rems et al. 192 Biophysical Journal 119, 190–205, July 7, 2020

Mechanistic Insights into the Modulation more hydrated and expands enough to conduct ion e charge im ance was kept cons the first h a VSD.we call such an ion-conductive pathway the "VSD on by r setu pore cribed in( stabili 2)AVSD pore forms,begins to expand,and becomes zed by lipic groups cal ned in VSDs,inclu ng thei ionic c ty VSD and then induces mieration of a le wo lipid head groups toward the membrane center.Field-induced expan .Ine type pore that fo Analysis the re ultsobtained for all channels (different type e visualiz are indicated by different markers as shown in the o the p bars in ind ion pa iven n nd the be ith y ius to R ng both ntially More detailed is given in the following sections.For all simulations,the eader can also find in Figs.S1-S4 SD ng sition of t e,the pro helical structure.and the change in the protein's solvent-accessible surface area. 10h 05.19 Electroporation of the NavMs channel d at the We first performed a simulation of NavMs channel under hyperpolarizing TMV of-1.5 V.The VSDs of NavMs ly.we I the average and standard devia already hy rated in the abs ence f an exteral electric of RESULTS ins de ted as VSD4.the S3 helix om S General observations S2.and 4.which further promoted hydration of the do 3 c).Compared to the other helice to the res of the volt larized.whereas the part facing the negative electrode can form and second it has only one salt-bridge intera becomes depolarized Thus,for each of the three investi with S4.unlike.for instance.S2.which has two such sal gated channels(NavMs,NavPaS,and HCN1).we generated bridges we that this s under hyperpo zIng and d The simulations revealed that electric field promotes cre by the first Cl-ion (Fi which was followed by the pas ation of different types of ion-conducting transmembrane sage of the first Na+ion at~130 ns.We will further refer to aqueous pathways. which are schematic such condu thway as the pore During the lined by lipid s the we y hich f d the lipid bilayer next to the rotein and has been previousl dle of the membrane ( observed also in pure lipid bilayers(2.3).We now report that gestions in the literature (4950).we will refer to such lipid-st zed structure as th With expan ion of the complex pore,the Its secondary Biophysical Journal 119.190-205.July 7.2020 193

a certain number of positive and negative ions across the membrane, which created the desired initial TMV. The charge imbalance was kept constant despite the transport of ions through transmembrane pores throughout the simulation by replenishing the charge imbalance using a Monte Carlo setup (43). The box size was allowed to expand in the x and y directions by using surface-tension coupling, as described in (44). Simulations in which we characterized the complex pores formed in VSDs, including their growth dynamics, ionic conductance, stability in the absence of an applied electric field, and perturbation of the VSD second￾ary structure (Figs. 8 and 9; Supporting Materials and Methods, Section S8), were carried out in the NPT ensemble. Analysis Trajectories were visualized with VMD (45). Ion passage through pores formed in the system was determined by a custom MATLAB (The Math￾Works, Natick, MA) code based on the positions of ions extracted from the trajectories. The code is available at https://github.com/ delemottelab. The radius of pores was determined by the CHAP tool (46), which uses the same concept as HOLE (47), i.e., using a Monte Carlo simulated annealing procedure to find the best route for a sphere with variable radius to squeeze through the pore cavity. Secondary struc￾ture and solvent-accessible surface area were determined with the GRO￾MACS functions gmx do_dssp and gmx sasa, respectively. Free energy estimates for water molecules along VSDs were determined based on kernel density estimate of the probability distribution of the positions of water molecules extracted from the trajectories at 0 V (see also Support￾ing Materials and Methods, Section S5). To calculate electrostatic profiles along the VSDs, additional 2-ns-long simulations under an electric field were performed while keeping the position of the protein heavy atoms restrained. These short simulations were carried out starting from 10 different configurations extracted from the last 100 ns of the trajectories without applied electric field. For each of these short trajectories, we determined the three-dimensional electrostatic potential using the VMD tool PMEpot (options: ewaldfactor 0.5, grid 1.0 A˚ ) (48). From each of the 10 three-dimensional profiles, we determined a one-dimensional (1D) electrostatic profile along the VSD by averaging the potential along a cylinder with a diameter of 0.5 nm centered at the center of mass of the VSD. Finally, we computed the average and standard deviation of the 10 1D profiles. RESULTS General observations When a cell is exposed to an electric field, the part of its membrane facing the positive electrode becomes hyperpo￾larized, whereas the part facing the negative electrode becomes depolarized. Thus, for each of the three investi￾gated channels (NavMs, NavPaS, and HCN1), we generated up to 600-ns-long trajectories under hyperpolarizing and de￾polarizing TMV of 51.5 V. The simulations revealed that electric field promotes cre￾ation of different types of ion-conducting transmembrane aqueous pathways, which are schematically depicted in Fig. 2 a. One of possible pathways is the well-known hydro￾philic lipid pore, lined by lipid headgroups, which forms in the lipid bilayer next to the protein and has been previously observed also in pure lipid bilayers (2,3). We now report that ion-conductive pathways can also form within the VSDs of VGICs, whereby two scenarios can occur: 1) A VSD be￾comes more hydrated and expands enough to conduct ions but not much more. As soon as the first ion passes through a VSD, we call such an ion-conductive pathway the ‘‘VSD pore.’’ 2) AVSD pore forms, begins to expand, and becomes stabilized by lipid headgroups. We call this pathway a ‘‘complex pore’’ and define it as a pore that initiates in a VSD and then induces migration of at least two lipid head￾groups toward the membrane center. Field-induced expan￾sion of such complex pore can lead to severe unfolding of the VSD. The type of pore that forms preferentially depends on the type of ion channel and the hydration and electro￾static profile of the VSDs. The graphs in Fig. 2 b summarize the results obtained for all channels (different types of pores are indicated by different markers as shown in the figure legend). The horizontal bars in Fig. 2 b indicate the length of each simulation. The markers show the time of the first ion passage through a given pore (orange color for Na, green color for Cl, circles for lipid pores, triangles for VSD pores). Note that unlike lipid pores, which start con￾ducting both Na and Cl ions practically simultaneously, VSD pores usually start conducting one type of ions prefer￾entially. More detailed description of individual simulations is given in the following sections. For all simulations, the reader can also find additional results in Figs. S1–S4, including the number of ions that passed each pore, the pro￾gression of the pore radius, the change in the number of pro￾tein residues in helical structure, and the change in the protein’s solvent-accessible surface area. Electroporation of the NavMs channel We first performed a simulation of NavMs channel under a hyperpolarizing TMV of 1.5 V. The VSDs of NavMs were already hydrated in the absence of an external electric field (Fig. 3 a). After the onset of the electric field, more wa￾ter molecules entered the VSDs (Fig. 3 b). In one of these do￾mains, denoted as VSD4, the S3 helix moved away from S1, S2, and S4, which further promoted hydration of the do￾main’s interior (Fig. 3 c). Compared to the other helices, S3 is weakly connected to the rest of the voltage sensor. First, it is a short helix, which limits the number of interactions it can form, and second, it has only one salt-bridge interaction with S4, unlike, for instance, S2, which has two such salt bridges; we anticipate that this weak connection resulted in the S3 helix being detached first under an excessive electric field. Shortly after S3 moved, at 90 ns, VSD4 was crossed by the first Cl ion (Fig. 3 d), which was followed by the pas￾sage of the first Naþ ion at 130 ns. We will further refer to such conductive pathway as the ‘‘VSD pore.’’ During the simulation, this VSD pore expanded and eventually became stabilized by lipid headgroups that migrated toward the mid￾dle of the membrane (Fig. 3, e–g). Following previous sug￾gestions in the literature (49,50), we will refer to such lipid-stabilized structure as the ‘‘complex pore.’’ With expan￾sion of the complex pore, the S4 helix lost its secondary Mechanistic Insights into the Modulation Biophysical Journal 119, 190–205, July 7, 2020 193

Rems et al 勇 6 Hyperpolarization(-1.5 V) Depolarization(+1.5 V) Legend: NavMs NavMs * un #2 V VSOR run #2 vsD VSD pore NavPaS NavPas n (cplx pore forms)】 n 7米vso2 hyperpolarization HCN1 HCN1 n 01002003004005006000100200300400500600 depolarization Time (ns) Time(ns) the mu he at least t pid he 0.5 nm fron ddle of f a VSD pore be ex oft VSD the VSD pore was formed.The data p时 ed in Ta on the structure,and VSD4 began to unfold,as shown in Fig.3. the first out of two simulations.VSD pores were formed (and Fig.S1;Videos SI and S2).In this simulation,another in VSD2 and VSD4.In the second simulation.a single VSD pore was formed in VSDI at ~120 ns,which allowed Na the passa of a singl C hat eiaropawodhoghvsDe e In all simulations.we also observed pa sage of nat conditions.Again,we observed formation of a VSD pore through the central pore of the channel (Fig.4:see also re om control s upporing Mate n the f g ed to be in ar tate (0).alth formed in VSD4,suggesting that multiple such p es can ulations we oberved that the channel po e became fully by be formed in the same voltage-gated ion channel drated and conductive only upon application of a strons Fig.S9).I NavMs has four identical VSDs beyond physiological values) electric field.Furthermore be formed rand omly in any c then as although in both ca s.the channel pore still became fully hydrated occasionally.similarly as already shown in Fig. expan at 2 ns By In a 194 Biophysical Joumal 119.190-205.July 7.2020

structure, and VSD4 began to unfold, as shown in Fig. 3, h–i (and Fig. S1; Videos S1 and S2). In this simulation, another VSD pore was formed in VSD1 at 120 ns, which allowed the passage of a single Cl ion. To verify that these results are reproducible, we per￾formed a second 200-ns-long simulation under the same conditions. Again, we observed formation of a VSD pore followed by its expansion into a complex pore, but this time in VSD3, with the first passage of a Cl ion occurring at 110 ns. Subsequently, a second complex pore was formed in VSD4, suggesting that multiple such pores can be formed in the same voltage-gated ion channel (Fig. S9). The NavMs has four identical VSDs; therefore, it is expected that pores can be formed randomly in any of them with equal probability. We further performed simulations under a depolarizing TMV of 1.5 V. In this case, we also observed formation of VSD pores, but they did not expand into complex ones, even though we prolonged the simulations to 400 ns. In the first out of two simulations, VSD pores were formed in VSD2 and VSD4. In the second simulation, a single Naþ and a single Cl ion passed through VSD2, confirming that the formation of complex pores is more difficult under depolarizing compared to hyperpolarizing TMV. In all simulations, we also observed passage of Naþ through the central pore of the channel (Fig. 4; see also re￾sults from control simulations in Supporting Materials and Methods, Section S4.3; Fig. S8). The structure of NavMs was presumed to be in an open state (20), although in sim￾ulations, we observed that the channel pore became fully hy￾drated and conductive only upon application of a strong (beyond physiological values) electric field. Furthermore, the pore stopped conducting as soon as we either lowered the electric field by 50% or completely turned it off, although in both cases, the channel pore still became fully hydrated occasionally, similarly as already shown in Fig. 4 a at 2 ns. By counting the number of Naþ ions that passed through the channel pore within 150 ns at 51.5 V, we a b FIGURE 2 Summary of the simulation results. (a) A schematic representation of the types of ion-conductive transmembrane pathways induced by an elec￾tric field that were observed in the simulations. (b) Graphs show the time of the first Naþ and Cl ion passages through a given pore (see legend) in sim￾ulations for all channels. Black dots and black stars show the time at which a lipid pore and a VSD pore, respectively, became stabilized by lipid headgroups (i.e., the time when at least two lipid headgroups were found within 0.5 nm from the middle of the membrane). If a VSD pore became stabilized by lipid headgroups, such pore was called a ‘‘complex pore.’’ The horizontal bars correspond to the simulation length. The text next to each horizontal bar gives the index of the VSD in which the VSD pore or complex pore was formed. The data presented in this figure are also tabulated in Table S2. The scheme on the bottom right shows how ions redistribute around the cell membrane upon exposure to an electric field; on the side facing the positive electrode (anode), the membrane becomes hyperpolarized, whereas on the opposite side, it becomes depolarized. To see this figure in color, go online. Rems et al. 194 Biophysical Journal 119, 190–205, July 7, 2020

Mechanistic Insights into the Modulation 42n5 100 Time (ns) VSD and this VSD by lipi .(g)As mplex por in gold,and s vely.B 实9 ine TMV (all at fomulation under depopi in the lipid bilav 27C)CTable S2).The estimated conductance rs lowe nel.as illustrated in.In the second simulation.a than the experimental value (33 pS at 2C)(51).Overall. that the channel str we used in ou y in d the pore expand hat the led by the high in the electric field,which promoted the hydration of the hydro perturbation are expected.Ion passage through the centra phobic gate along the channel pore. at the functional stat ctu to electrophysiologi Electroporation of the NavPaS channel To verify that the formation of VSD pores and complex pores is not obs in NavMs we studied the cukaryo Electroporation of the HCN1 channe chanf two under depolarizn TMV.In the first simula prone to poration than those of NavMs and NavPaS (Figs and S )In one ut of two simulatio did not th VSD2 exp in both simulations.we also observed formation of a lipid of the VSD from the channel,similarly as in NavMs com- pore in the lipid bilayer surrounding the protein,suggesting plex pores (Fig. s S3 and S4).In the second sim that formation of pores within the ntical ing TMV.a VSD por e was Ior han latic pore.However,the complex pore did not continue to expand TMV,we observed no pores within the entire 600 ns run. within the VSD.Instead,the pore expanded into the bilayer. The structure of this channel is expected to be in the closed secondary structure (Fig.S2).In the first Biophysical Journal 119.190-205.July 7.2020 195

estimated its conductance to be 17 and 2 pS in the two simulations under hyperpolarizing TMV and 17 and 15 pS in the two simulations under depolarizing TMV (all at 27C) (Table S2). The estimated conductance appears lower than the experimental value (33 pS at 22C) (51). Overall, we speculate that the channel structure that we used in our simulations does not correspond to a fully open state and that the passage of Naþ ions was only enabled by the high electric field, which promoted the hydration of the hydro￾phobic gate along the channel pore. Electroporation of the NavPaS channel To verify that the formation of VSD pores and complex pores is not observable only in NavMs, we studied also the eukaryotic NavPaS channel. The same as for NavMs, we performed two simulations for NavPaS under hyperpola￾rizing and two under depolarizing TMV. In the first simula￾tion under a hyperpolarizing field, VSD pores were formed in VSD1 and VSD2 at 233 ns and 249 ns, respectively. Whereas the pore in VSD1 did not expand considerably and mainly enabled the transport of Cl ions, that in VSD2 expanded into a complex pore that led to unfolding of the VSD from the channel, similarly as in NavMs com￾plex pores (Fig. 5 a; Videos S3 and S4). In the second simu￾lation under hyperpolarizing TMV, a VSD pore was formed in VSD2 at 120 ns, which further expanded into a complex pore. However, the complex pore did not continue to expand within the VSD. Instead, the pore expanded into the bilayer, as shown in Fig. 5 b. This limited the number of residues that lost their helical secondary structure (Fig. S2). In the first simulation under depolarizing TMV, a VSD pore was formed in VSD1 at 80 ns. In addition, at 190 ns a lipid pore was formed in the lipid bilayer surrounding the chan￾nel, as illustrated in Fig. 5 c. In the second simulation, a VSD pore was again formed in VSD1 and also in VSD2, but only in VSD2 did the pore expand into a complex one, similarly as in Fig. 5 a. Importantly, unlike NavMs, NavPaS has four different VSDs; therefore, differences in their perturbation are expected. Ion passage through the central pore was not observed in any of the simulations. It should be noted that the functional state of the NavPaS structure is unclear because it was not possible to electrophysiologi￾cally characterize this channel (21). Electroporation of the HCN1 channel In addition to Nav channels, we performed simulations on HCN1. The VSDs of HCN1 turned out to be much less prone to poration than those of NavMs and NavPaS (Figs. 2 and S3). In one out of two simulations performed under hyperpolarizing TMV, we observed the passage of two Naþ ions through VSD4 at 497 and 525 ns. However, in both simulations, we also observed formation of a lipid pore in the lipid bilayer surrounding the protein, suggesting that formation of pores within the identical VSDs is energet￾ically less favorable than that of lipid pores (Videos S5 and S6). In the single simulation performed at depolarizing TMV, we observed no pores within the entire 600 ns run. The structure of this channel is expected to be in the closed state, and in none of the simulations did we observe passage of ions through the channel pore. ab c d i ef g h FIGURE 3 Formation of a complex pore in VSD4 of NavMs channel. (a) VSD was hydrated already in the absence of electric field. (b and c) Upon electric field application, more water molecules entered the VSD. (d) The first Cl ion passed through the VSD at 90 ns after the onset of the electric field. (e and f) More water and ions entered the VSD, and this VSD pore became stabilized by lipid headgroups forming the so-called complex pore. (g) As the complex pore expanded, the VSD began to unfold from the channel. (h) Unfolded VSD viewed from the extracellular side. In (a)–(h) the VSD is colored in red, water in cyan, lipid phosphorus atoms in gold, and sodium and chloride ions in yellow and green, respectively. Black arrows mark the first Cl ion within VSD and the first lipid headgroup moving into the pore. (i) Disruption of the VSD’s secondary structure. To see this figure in color, go online. Mechanistic Insights into the Modulation Biophysical Journal 119, 190–205, July 7, 2020 195

Rems et al I pore of t table bridge n color,go online s and omplex pores form more easily in ules VSDs that electrostatically more favorable for the entry of nels,we plotted the height of the free ene rbarrier versus ions or CI)through The results presented above showed he that t bles.c to be perturbed by the electric field ().Consider ate this hypoth simulations ormed,we observed conduc tion In additior s4 VSDs of N n of reported in Tables S4 and S5).In NavPaS.which has four shown in Fig.6a.We then performed simulations of the c conduction only through mutant HCNI under hyperpolarizing TMV.In one out of ons, nd formed a c a VSD pore was for he follow sage of two sodium ions YSD of the mutant HCNIn be more easily observe the th ough VSD4 in one out o the th e simulations,where lipid pore was form our hypothesis thes VSD is and VSD that make it ne to n We thu 6 fomed more poration.The simulations showed that in NavMs.com sized that VSD pores and complex pore bores could only be formed under hyp erpolarizing tmv easily in VSDs more hydrated This hypoth yer the 7 omplex pos th e lipid c52).as well as on th tha t in lipid bilayer,a pore is formed more easily if this bilayer is VSDI was considerably nder depolarizing TMV. We hypoth that the pre-embedded with vater molecules (53) asymr ry in nd con To inve e hyd SD ion of thei ition along the VSDs by lyzing the ability distribution of these molecules in the abse alons one of the VSDs of NavMs (all four VSDs nce of an field (Fig.6 a).The free similar profile:Fig.S19)and VSDI and VSD2 of NavPaS VSD in the e ers in th potential in the VSDs of NavMs and in VSDI and VSD2 of NavPaS ak at th ide tracts Cl but repels Nat The transport of ions is The water probability distribution in the two other VSDs of ded under for NavPaS and the four which we SDs of HCNI has hydrophobic gap ores,bu iplex pore also in con y mg al i energy vPa oth por 196 Biophysical Joumal 119.190-205.July 7.2020

VSD pores and complex pores form more easily in VSDs that are more hydrated and are electrostatically more favorable for the entry of ions The results presented above showed that it is possible to observe conduction of ions through one or more VSDs in all three channels. However, not all VSDs are equally likely to be perturbed by the electric field (Fig. 2). Considering all simulations performed, we observed ionic conduction through all of the identical VSDs of NavMs. In addition, in NavMs we observed the expansion of VSD pores into complex ones in all VSDs (see also the control simulations reported in Tables S4 and S5). In NavPaS, which has four different VSDs, we observed ionic conduction only through VSD1 and VSD2. Moreover, VSD2 formed a complex pore, whereas VSD1 did not. In HCN1, formation of a lipid pore was more favorable compared to a VSD pore, as we could only observe the passage of two sodium ions through VSD4 in one out of the three simulations, whereas a lipid pore was formed in two of these simulations. Overall, these differences suggest that there should be some features of the VSD that make it more prone to porate. We thus hypothe￾sized that VSD pores and complex pores are formed more easily in VSDs that are more hydrated. This hypothesis is based on previous MD simulations of pure lipid bilayers showing the crucial role of water molecules in formation of lipid pores (52), as well as on the observation that in a lipid bilayer, a pore is formed more easily if this bilayer is pre-embedded with water molecules (53). To investigate the hydration profile of individual VSDs, we estimated the free energy of water molecules as a func￾tion of their position along the VSDs by analyzing the prob￾ability distribution of these molecules in the absence of an external electric field (Fig. 6 a). The free energy profiles indeed show the lowest barriers in the VSDs in which we observed formation of VSD pores and complex pores, i.e., in the VSDs of NavMs and in VSD1 and VSD2 of NavPaS. The water probability distribution in the two other VSDs of NavPaS and the four VSDs of HCN1 has hydrophobic gaps, resulting in considerably higher free energy barriers (see also two-dimensional images showing the position of water molecules in Fig. 6 c; Supporting Materials and Methods, Section S5). Furthermore, for all VSDs from all the chan￾nels, we plotted the height of the free energy barrier versus the time of the first ion passage (either Naþ or Cl) through the VSD. The graph in Fig. 6 b shows that there is a positive correlation between these two variables, confirming that hy￾dration of the VSD is an important feature that contributes to VSD’s propensity for poration. To corroborate this hypoth￾esis further, we considered a mutant’s HCN1 VSD that con￾tained three mutations of nonpolar into polar residues (54). The mutations increased the hydration of the VSDs and decreased the free energy barrier for water molecules, as shown in Fig. 6 a. We then performed simulations of the mutant HCN1 under hyperpolarizing TMV. In one out of two simulations, a VSD pore was formed at 250 ns after the onset of the electric field, followed by its transformation into a complex pore; this observation indicates that the VSDs of the mutant HCN1 can be more easily porated compared with the wild-type channel and goes in line with our hypothesis (Fig. S4; Table S2). The electrostatic potential inside the VSD is another candidate feature that contributes to VSDs’ propensity for poration. The simulations showed that in NavMs, complex pores could only be formed under hyperpolarizing TMV. In NavPaS, we observed formation of a complex pore in VSD2 regardless of the TMV polarity, whereas in VSD1 we only observed formation of a VSD pore. Moreover, the pore in VSD1 was considerably more selective for Cl ions under depolarizing TMV. We hypothesized that the observed asymmetry in pore formation and conduction is due to the asymmetric distribution of charges along the VSDs. Fig. 7 compares the electrostatic potential profiles along one of the VSDs of NavMs (all four VSDs show a similar profile; Fig. S19) and VSD1 and VSD2 of NavPaS. Note that the electrostatic potential was determined under an applied electric field before poration. The NavMs VSD has a high positive peak at the extracellular side, which at￾tracts Cl ions but repels Naþ. The transport of ions is thus impeded under depolarizing TMV, for which we observed formation of VSD pores, but not complex pores. In VSD2 of NavPaS, in which a complex pore was formed under both hyperpolarizing and depolarizing TMV, the a b FIGURE 4 Passage of Naþ ions through the cen￾tral pore of the NavMs channel. (a) Before electric- field application, the bottom half of the pore is dehy￾drated. After the application, a stable water bridge is formed in the bottom half of the channel pore, and the number of water molecules hydrating the bottom half gradually increases, allowing the transport of Naþ ions. (b) The cumulative sum of the number of ions that passed through the channel pore in sim￾ulations under hyperpolarizing (blue shades) and depolarizing TMV (red shades). To see this figure in color, go online. Rems et al. 196 Biophysical Journal 119, 190–205, July 7, 2020

Mechanistic Insights into the Modulation FIGURE S Formation of ores in the NavPas (a)For of n of th d pore se o. om row sho ws the of t electrostatic profile is similar at either side of the membrane. tion.The NavPaS complex pore was crossed by 104 Na DI of NavPaS the ions under hyperpolarization and by 9 Na sitive pea f CI un ization. comple VSD under denolarizing tmy allowing passa A similar prefe ce for cl-ions has heen ohserved in ions compared with only I Naion(see Table S2).Interest- lipid pores.in which this selectivity could be explained yet comp pores wer only in which nects the ectrosltaicfoy VSDs that influence nore formation such as for instane in the direction of Clion movement (56). the salt-bridge connec salt-bridge The size of complex pores depends on the TM cou asily the 8gnascaneonauatleaspretereniayn Comiebabedicieheapietelectncieid (and th eby the TMV).This is simila to what has been for poro nlex n s in NavMs and NavPaS.The pores were eated at TMVof: at a TMV to prevent the 0 ized ionic transport through such stabilized po eg535 TMV of +15 V the As expected,the conduction of ions through NavMs com crease.For a TMVof1.0 V,the pore radius also pore w to incr but at a ugh thec ndent on TMV polarity (Fig.S20).More specifically helices and linid headgrouns and have an irregular shape within 100 ns of the simulation.the NavMs complex pore The pore radius.determined here as the radius of the largest sphere than can be pushed through the pore,is thus a simpli and 52 Cl ions under ied rep or the pore s12 Biophysical Journal 119.190-205.July 7.2020 197

electrostatic profile is similar at either side of the membrane. In VSD1 of NavPaS, the electrostatic potential also shows a positive peak at the extracellular side. Accordingly, we observed selective conduction of Cl ions through this VSD under depolarizing TMV, allowing passage of 9 Cl ions compared with only 1 Naþ ion (see Table S2). Interest￾ingly, the height of the peaks for VSD1 and VSD2 are similar, yet complex pores were formed only in VSD2. This suggests that there are additional features of the VSDs that influence pore formation such as, for instance, the salt-bridge connections. VSD2 has fewer salt-bridge connections between S1–S4 helices than VSD1 (3 vs. 4; see Table S7), which could be the reason why it unfolds more easily. Complex pores can conduct ions preferentially in one direction To investigate further the asymmetric ion conduction through complex pores, we performed additional simula￾tions of NavMs and NavPaS with these pores under a 3 lower electric field, i.e., resulting in TMV z Ez Lz ¼ 0.5 V (55). We reduced the electric field to prevent the complex pores from further expansion, and we character￾ized ionic transport through such stabilized pores (53,54). As expected, the conduction of ions through NavMs com￾plex pore was considerably larger under hyperpolarizing TMV, especially for Cl ions, whereas the conduction through the NavPaS complex pore (in VSD2) was less dependent on TMV polarity (Fig. S20). More specifically, within 100 ns of the simulation, the NavMs complex pore was crossed by 37 Naþ and 195 Cl ions under hyperpolar￾ization and by 20 Naþ and 52 Cl ions under depolariza￾tion. The NavPaS complex pore was crossed by 104 Naþ and 276 Cl ions under hyperpolarization and by 99 Naþ and 189 Cl ions under depolarization. Both complex pores were less conductive to Naþ compared to Cl ions. A similar preference for Cl ions has been observed in lipid pores, in which this selectivity could be explained by lower bulk mobility of Naþ ions and their binding to lipid headgroups, which also affects the electrostatic envi￾ronment inside the pore and induces an electroosmotic flow in the direction of Cl ion movement (56). The size of complex pores depends on the TMV When maintaining the TMV at 1.5 V, complex pores continue to expand in size. However, the size of the pores can be stabilized by reducing the applied electric field (and thereby the TMV). This is similar to what has been observed before for pores in pure lipid bilayers (44,57). Fig. 8 shows the progression of the radius of example com￾plex pores in NavMs and NavPaS. The pores were created at a TMV of 51.5 V, and then the simulation was continued either at a TMV of 51.5 Vor at a lower electric field, result￾ing in TMV of 51.0 V or 50.5 V. The behavior of both NavMs and NavPaS complex pores was similar. For a TMV of 51.5 V, the pore radius continued to rapidly in￾crease. For a TMVof 51.0 V, the pore radius also continued to increase but at a slower rate. For a TMV of 50.5 V, the pore radius somewhat decreased and became stabilized. Note that complex pores are lined both by transmembrane helices and lipid headgroups and have an irregular shape. The pore radius, determined here as the radius of the largest sphere than can be pushed through the pore, is thus a simpli- fied representation of the pore size. ab c FIGURE 5 Formation of pores in the NavPaS system. (a) Formation of a complex pore in VSD2. Expansion of the pore continues to unfold the VSD. (b) Formation of a complex pore initiates in the VSD and then expands into the bilayer. (c) Formation of a lipid pore close to, but not associated with, the protein. Top row shows the side view of the VSD. Bottom row shows the extracellular view of the protein. Representation of atoms is the same as in Figs. 1 and 3. To see this figure in color, go on￾line. Mechanistic Insights into the Modulation Biophysical Journal 119, 190–205, July 7, 2020 197

Rems et al. Estimated free energy of water along VSDs NavMs NavPaS HCN1 wt HCN1 mut 人 Correlation ition of water along VSDs of NavPaS VSD1 VSD2 VSD () r(nm r(nm) height of the gy barrier and the f the first ion passag either Na*or for this last 200ns of trajctory.To this figure Complex pores are more stable than lipid pores Pores created by electric fields in pure 1-palmytoyl-2-oleoyl- phosphatidylcholine bilayers have been reported to close on ulations did we observe complete recovery of the VSD's sec ondary structure:all simulations suggested that formationof complex pores ex pore we cond in Fig. 3.in the of the pore at .100.and 1000ns after switching off the ele Salt-bridge reorganization in VSDs e but remainec by pass thup The ability of VSDs to respond to cha in TMV is tested the sbility of complex pores ented in Fis These pores also remained opened and stabilized by lipid dues are primary elements to sense an applied electric field ation (Sup- pect that the time for whicha ever, n eectric-ficlda open depends on the size that it reached during the elc move along its direction,causing the disruption of existing comple pore,which ions,application o lectric field me some ex es wer idges wer 198 Biophysical Joumal 119.190-205.July 7.2020

Complex pores are more stable than lipid pores Pores created by electric fields in pure 1-palmytoyl-2-oleoyl￾phosphatidylcholine bilayers have been reported to close on average within 50 ns (38). To investigate the postpulse sta￾bility of complex pores, we conducted a 1-ms-long simulation of the NavMs complex pore, presented in Fig. 3, in the absence of an electric field. Fig. 9 a shows the configuration of the pore at 0, 100, and 1000 ns after switching off the elec￾tric field. The pore reduced in size but remained stabilized by lipid headgroups even after 1000 ns. Ions were able to enter and pass through the pore by diffusion (Fig. 9 b). We further tested the stability of complex pores presented in Fig. 8. These pores also remained opened and stabilized by lipid headgroups at the end of the 400-ns-long simulation (Sup￾porting Materials and Methods, Section S8). However, we expect that the time for which a complex pore remains open depends on the size that it reached during the elec￾tric-field exposure and the extent by which it deformed the VSD. We observed that in a NavMs complex pore, which barely met our criterion for a complex pore and became sta￾bilized by only two lipids, the lipids returned to their default orientation 200 ns after turning off the electric field (Fig. S22). Nevertheless, in none of the abovementioned sim￾ulations did we observe complete recovery of the VSD’s sec￾ondary structure; all simulations suggested that formation of complex pores leads to persistent perturbation of the VSD structure, especially the S4 helix (Figs. S21–S23). Salt-bridge reorganization in VSDs The ability of VSDs to respond to changes in TMV is granted by positively charged residues of the S4 segment. Being embedded into a low dielectric medium, these resi￾dues are primary elements to sense an applied electric field. Inside a VSD, they interact with negative counterparts com￾ing from the remaining S1–S3 segments through salt bridges. Upon electric-field application, the S4 residues move along its direction, causing the disruption of existing salt bridges and formation of new ones (58–60). In our sim￾ulations, application of the electric field modified the salt￾bridge connections within the investigated VSD: some ex￾isting salt bridges were broken, and new salt bridges were a b c FIGURE 6 Hydration of VSDs. (a) Free energy profiles estimated from probabilities of water distribution along the VSD principal axis (z) of NavMs, NavPaS, wild-type HCN1, and mutant HCN1, averaged over 200 ns. Note that the free energy barriers are cut at 6 kcal/mol. (b) Correlation between the height of the free energy barrier and the time of the first ion passage (either Naþ or Cl) through a VSD. Ion passage times are taken from all simulations reported in Table S2. Gray dots represent VSDs that were not porated in any of the simulations. The outlier at (600 ns, 2.3 kcal/mol) corresponds to VSD2 of HCN1 mutant; note that we performed only two simulations for this channel. (c) Positions of water molecules in each VSD of NavPaS projected along the VSD radius (r) and the VSD principal axis (z). The VSDs’ center of mass is located at (0,0). Blue circles show all positions extracted from 200 frames of the last 200 ns of the equilibration trajectory. To see this figure in color, go online. Rems et al. 198 Biophysical Journal 119, 190–205, July 7, 2020

Mechanistic Insights into the Modulation NavMs VSD4 NavPaS VSD1 NavPaS VSD2 Na +Na 2 Na 05 VSD com Distance from VSD com (nm) side ide s.charge (neg.charge) side is on the left.Tose this figure formed.An example is depicted in Fig.10 a and Video S7, Comparison with the charge imbalance method which show how salt bridges were perturbed in VSD2 of the ore appnc ets of event negatively charged residues on the S1-S3 helices(Fig.10 as in simulations with the external electric-field method a).Upon an ele nd depolai ction with asp81 Afte off th .the VSD pore formed in then expanded into a complex pore.Under depolarizing ble for at least I us (last time tested).Two other connections, and ARG109-ASP81.that existed be TMV,the VSD pore formed in VSD3 at 104 ns,whereas and fo of salt bridge also observed in other NavMs VSDs and in other channels. important difference between the charg external electric-field simulations.In charge imbalance as shown in Tables S6-S8.It is interesting to note that the entry of ions nt IO kage. fe to red Biophysical Journal 119.190-205.July 7.2020 199

formed. An example is depicted in Fig. 10 a and Video S7, which show how salt bridges were perturbed in VSD2 of the NavMs channel under hyperpolarizing TMV (note that VSD2 was not porated in this simulation). Before applica￾tion of the electric field, there were four connections formed by positively charged arginine residues on the S4 helix and negatively charged residues on the S1–S3 helices (Fig. 10 a). Upon an electric-field application, all these salt bridges were broken, and ARG106 shifted such that it formed a new connection with ASP81. After turning off the electric field, this new connection, ARG106-ASP81, remained sta￾ble for at least 1 ms (last time tested). Two other connections, ARG103-ASP49 and ARG109-ASP81, that existed before the electric-field application were reformed (Figs. S24– S27). Similar breakage and formation of salt bridges was also observed in other NavMs VSDs and in other channels, as shown in Tables S6–S8. It is interesting to note that the entry of ions into the VSDs was not necessary for the breakage of salt bridges, but it could facilitate this breakage, as shown in Fig. 10 b. Comparison with the charge imbalance method Finally, we performed additional simulations for the NavMs system, in which we mimicked electroporation conditions by establishing a charge imbalance across the membrane. With these simulations we observed similar sets of events as in simulations with the external electric-field method: both under hyperpolarizing and depolarizing TMV, firstly Naþ ions started to pass through the central channel pore, and then a VSD pore formed (Fig. S11). Under hyperpola￾rizing TMV, the VSD pore formed in VSD1 at 29 ns and then expanded into a complex pore. Under depolarizing TMV, the VSD pore formed in VSD3 at 104 ns, whereas this pore was not able to expand into a complex pore within the 200-ns-long simulation. However, we observed an important difference between the charge imbalance and external electric-field simulations. In charge imbalance, the TMV kept dropping though the simulation, even though the charge imbalance was kept constant. This is due to redis￾tribution of ions during simulation (some of the ions move ab c FIGURE 7 Electrostatic profiles along different VSDs under hyperpolarizing (blue lines) and depolarizing TMV (red lines). (a) VSD4 of NavMs. (b) VSD1 of NavPaS. (c) VSD2 of NavPaS. Corresponding profiles in the absence of applied electric field are given in Fig. S18. The line thickness corresponds to standard deviation of 10 1D profiles. The gray area at the back of each graph shows the mass density profile of the VSD. Images below graphs show the position of charged residues on each VSD: positively charged ARG, HIS, and LYS are colored green, negatively charged residues ASP and GLU are colored magenta. Intracellular side is on the left. To see this figure in color, go online. Mechanistic Insights into the Modulation Biophysical Journal 119, 190–205, July 7, 2020 199