Electric Field Induced Wetting of a Hydrophobic Gate in a Model Nanopore Based on the 5-HT3 Receptor Channel

In this study we examined the influence of a transmembrane voltage on the hydrophobic gating of nanopores using molecular dynamics simulations. We observed electric field induced wetting of a hydrophobic gate in a biologically inspired model nanopore based on the 5-HT3 receptor in its closed state, with a field of at least ∼100 mV nm–1 (corresponding to a supra-physiological potential difference of ∼0.85 V across the membrane) required to hydrate the pore. We also found an unequal distribution of charged residues can generate an electric field intrinsic to the nanopore which, depending on its orientation, can alter the effect of the external field, thus making the wetting response asymmetric. This wetting response could be described by a simple model based on water surface tension, the volumetric energy contribution of the electric field, and the influence of charged amino acids lining the pore. Finally, the electric field response was used to determine time constants characterizing the phase transitions of water confined within the nanopore, revealing liquid–vapor oscillations on a time scale of ∼5 ns. This time scale was largely independent of the water model employed and was similar for different sized pores representative of the open and closed states of the pore. Furthermore, our finding that the threshold voltage required for hydrating a hydrophobic gate depends on the orientation of the electric field provides an attractive perspective for the design of rectifying artificial nanopores.

N anopores enable the permeation of water and small molecules through membranes that separate aqueous compartments. The size of these permeant molecules is typically comparable to the diameter of the pores involved (≤1 nm). Consequently, any interaction of the permeant molecules with the lining of the pores will determine the functional properties of the permeation process. 1−4 In biological membranes, ion channels are structurally dynamic protein nanopores that switch between functionally open and closed conformations to control the rapid permeation of ions and water across the membrane. 5,6 However, when open, a typical channel pore has an internal radius of ∼0.5 nm and a length of ∼5 nm that creates a nanoconfined environment where the precise shape and related dynamic physicochemical properties of the pore will also influence the complex behavior of water and of ions, and therefore their permeation across the membrane. For example, we have previously shown how ion permeation can be influenced by both the pore radius and the local hydrophobicity of the pore lining. 7−9 In such cases, the presence of the band of hydrophobic amino acid side chains lining a pore can induce a local liquid-to-vapor phase transition resulting in that section becoming devoid of liquid water. Such pore "dewetting" creates a free energy barrier for ion permeation. Thus, although ion permeation can occur through polar regions only just wider than the radius of the permeating ion, a hydrophobic region of comparable dimensions can prevent ion permeation without complete steric occlusion of the pore. This process has been referred to as hydrophobic gating and has enhanced our understanding of the mechanisms which control ion permeation through biological ion channels and synthetic nanopores. 10−14 The gating or regulation of hydrophobic gates in ion channels is thought to involve either structural changes in pore radii or more subtle changes in hydrophobicity through, e.g., the rotation of helices that contain side chains of different polarities, and examples of such a hydrophobic gating mechanism have now been proposed in a variety of different channels 13,15−28 and synthetic nanopores. 29,30 However, nearly all biological membranes experience a potential difference of between 50 to 200 mV across them and molecular dynamics (MD) simulations of simple model nanopores show that application of a transmembrane electric field can hydrate an otherwise dewetted hydrophobic constriction, which is thus rendered ion permeable. 31,32 This electrowetting effect can be explained by incorporating the electrical energy of the pore volume element which is wetted into a simple thermodynamic model of gating (see eq 1 below 8 ). In addition to this model, which assumes bulk water properties, a number of MD studies have explored the effects of an externally applied electric field on water molecules in nanoconfined hydrophobic environments. 4 For example, simulations of water nanodroplets on hydrophobic surfaces 33 reveal that an electric field modifies the interfacial tension of a nanodroplet by influencing the orientation and hydrogen bonding structure of water molecules located on the droplet surface. Furthermore, MD simulations of water in planar hydrophobic confinement 34 also reveal that interfacial tensions decrease upon application of an electric field, alongside a field-induced change in the average number of hydrogen bonds formed by interfacial water molecules. Electrowetting has also been observed experimentally for several types of artificial biomimetic nanopores. 35 For example, it has been shown that hydrophobic nanopores in polyethylene terephthalate membranes can be wetted and functionally opened through the application of an electric field. 29 For biological ion channels, the role of electrowetting is even less well understood; most studies have been limited to MD simulations, demonstrating that the dewetted hydrophobic gate of the MscS channel can hydrate when a transmembrane voltage of between 0.25 and 1.2 V is applied. 36,37 Comparable electrowetting of a hydrophobic gate in a model protein nanopore has been simulated, 38 and electrowetting has also been seen within the hydrophobic pores formed by carbon nanotubes. 39 However, while it is clear that electrowetting can in principle open a hydrophobic gate, a detailed quantitative relationship between pore hydration probability and the applied electric field has not so far been examined for a biologically realistic model of an ion channel. It therefore remains unclear how such effects of electric fields depend on the structure of the pore and the strength of the applied field or how robust such predictions are to variations in the precise type of water model used in the MD simulations. It is therefore important that these effects are captured by a theoretical model which may subsequently be used in a quantitative predictive fashion.
In this study we have examined the effects of an applied electric field on a hydrophobic gate within a model protein nanopore derived from the transmembrane pore of the 5-HT 3 receptor, a biologically relevant ion channel distributed throughout the central nervous system. 40 Our results demonstrate that an electric field can induce wetting of the hydrophobic gate in this model nanopore, at values that correspond to a (supra-physiological) transmembrane potential difference of ∼1 V. A simple thermodynamic model can be used to quantitatively describe the increase in pore hydration probability with increasing electric field. Interestingly, the observed asymmetry of the hydration probability response curve can also be accounted for by the intrinsic electric field generated by the distribution of charged amino acid residues along the length of the nanopore which provides an attractive perspective for the design of intrinsically rectifying artificial nanopores.

RESULTS AND DISCUSSION
Model Nanopore Based on the 5-HT 3 Receptor. To quantify the effects of electrowetting of a hydrophobic gate, we have employed MD simulations of the M2 helix nanopore which corresponds to the pore lining segment of the 5-HT 3 receptor channel. 19 This model nanopore forms an ideal system for studying hydrophobic gating because it embodies many aspects of the wetting/dewetting behavior observed in complex ion channel proteins, 41 while remaining sufficiently simple to allow functional dissection by an extended series of atomistically detailed simulations ( Figure 1). 42 To evaluate our expectation that the model nanopore can capture the hydration state behavior of the larger parent channel protein, we compared simulations (3 × 50 ns) of the complete 5-HT 3 receptor in closed (PDB ID: 4PIR) 43 and open (PDB ID: 6DG8) 19 states with simulations (3 × 150 ns) of the corresponding model nanopores. In the latter simulations, only the pore-lining M2 helices of the channel were included. Both the complete protein channel and the model nanopores were embedded in a phospholipid bilayer and solvated in an ∼0.15 M NaCl electrolyte using the mTIP3P water model.
These simulations show that, in the absence of an electric field, the resulting water density profiles ( Figure 1B) through the transmembrane pore regions are nearly identical for the complete channel and for the nanopore simulation systems. In both cases, the closed state (PDB ID: 4PIR) conformation that had previously exhibited significant dewetting 41 still favors the vapor state even when only the M2 helix nanopore is simulated. In the open conformation (PDB ID: 6DG8), the time-averaged water density in the hydrophobic gate region of the M2 helix nanopore falls slightly below the value observed for the complete channel structure, reflecting the fact that the M2 helix nanopore intermittently dewets. However, the average water density remains close to the bulk value in these simulations, suggesting the influence of the structural model on the hydration equilibrium is limited and not sufficiently strong enough to bias an otherwise hydrated channel toward the vapor state. These results therefore provide confidence these different conformations of the M2 helix nanopore can be used as a suitable model system to quantify the relationship of the pore radius, electric field, and water model in the electrowetting and conductance behavior of a biologically relevant nanopore.
MD Simulations of the Model Nanopore in an Electric Field. We therefore examined the effect of an electric field on the liquid−vapor equilibrium that exists within the hydrophobic barrier of the M2 helix nanopore that is based on the closed state structure of the 5-HT 3 receptor (PDB ID: 4PIR) 43 ( Figure 2A). Without an electric field, this nanopore remained dewetted throughout most of a 150 ns long simulation (as seen previously 41 ) and only entered the wetted state for a small number of very brief periods ( Figure 2B). Correspondingly, the time-averaged water density profile ( Figure 2C) demonstrates that the hydrophobic gate is primarily devoid of bulk ACS Nano www.acsnano.org Article phase water and exists in a vapor state. When an electric field of ±100 mV/nm was imposed (using the methods described by, e.g., refs 44−46) the channel exhibited frequent liquid− vapor oscillations ( Figure 2B). In this situation, the timeaveraged density of water in the hydrophobic gate adopts an intermediate value between the dewetted and the hydrated channels ( Figure 2C). However, it is important to note the relationship between the units of field strength (i.e., mV/nm) used here to describe the external electric field with the values of transmembrane potential commonly measured experimentally (i.e., mV). In our simulation system, a field strength of 100 mV/nm that induces wetting in this pore (see Figure 2) corresponds to a potential difference ΔV of ∼0.85 V across the membrane; this is calculated from ΔV = −EL z , where E is the external field imposed along the z axis of a simulation box of dimension L z , as discussed in ref 47. Therefore, if, e.g., the box dimension L z was doubled the transmembrane voltage ΔV would be doubled. Notably, a value of ΔV = 0.85 V across the membrane is almost an order of magnitude greater than the typical transmembrane voltages experienced physiologically (∼0.1 V). Interestingly, we also observed that a field strength of −100 mV/nm (i.e., one in which the intracellular face of the membrane is at a negative potential relative the extracellular face; Figure 2A) has a greater effect than a field of +100 mV/ nm, indicating that the hydration probability of the pore can depend on the direction of the electric field as well as on its magnitude. At higher magnitude electric fields of ±200 mV/ nm, the channel pore exists almost exclusively in the wetted state irrespective of field direction ( Figure 2B) and the average water density approaches that of bulk water ( Figure 2C).
An external electric field can therefore be used to control the phase behavior of nanoconfined water within a biological nanopore by influencing the relative probability of the liquid and vapor states. The probability distribution of the minimal water density is bimodal ( Figure 2D; SI Figure S1), and hence the underlying free energy landscape has two minima (as discussed previously 4,7,48 ), with one peak near zero and a second peak near the density of bulk water. We note that this figure shows an estimate of the water density in the hydrophobic gate region aggregated over time and the electric field strength. This demonstrates clearly (by aggregating across all field strengths) that the probability distribution of the water density is bimodal, with intermediate density values observed only very rarely. Importantly, this legitimizes our theoretical treatment of the system using a two-state model.
Dependence of the Hydration Probability on the Electric Field. The response of the time averaged pore hydration probability, ⟨ω⟩, to an external electric field, E, was compared for the M2 bundle nanopore in conformations corresponding to closed (PDB ID: 4PIR) 43 and open (PDB ID: 6DG8) 19 states of the 5-HT 3 receptor. The hydration behavior of the nanopore was also compared for four different water models ( Figure 3). For the closed state, the hydration probability vs E response curve is sigmoidal: the hydration probability is near zero in the absence of a field, increasing to a value of one at fields of high magnitude. The threshold for hydration also varies. For example, with the mTIP3P water model 49 (which is widely used) pore wetting is observed at a field strength which is ∼50 mV/nm weaker than with any of the other models (including TIP4P/2005 50 which reproduces well a wide range of properties of water over temperatures from 123 to 573 K and for pressures up to 40 000 bar). However, once pore wetting begins, the transitional regime of intermediate hydration occurs within a further increase of ∼100 mV/nm for all four water models. In the case of the open state structures, the response of hydration probability still resembles a (double) sigmoid, but hydration probability does not start at zero, because the pore is already partially hydrated even in the absence of an electric field.
It is important to point out that in the transitional regime the rate of increase in hydration probability appears similar for both positive and negative electric fields, irrespective of the water model. However, the overall curve is shifted toward the right, which makes it asymmetric with respect to E = 0. This is most clearly seen for the closed state, where the threshold voltage for the onset of pore hydration is ∼25 mV/nm higher at positive fields than for negative fields. Due to the smaller The probability of hydration of a nanopore in the absence of an E-field can be described by where ⟨ω⟩ = the time-averaged hydration probability, ΔΩ 0 = Ω V − Ω L , i.e., the difference between the free energies of the liquid and vapor states in the absence of an E-field, and β = 1/ k B T. 7,11 When ΔΩ 0 < 0, the pore is hydrophobic and favors a vapor state, but when ΔΩ 0 > 0, the liquid (i.e., hydrated) state is favored. To include the effect of the electric field, E, the free energy term 31 can be modified giving where m denotes the strength of the coupling between hydration probability and the magnitude of the field and E INT accounts for the horizontal offset of the response curve in terms of an intrinsic electric field arising from the nanopore structure (see below). The value of m represents the difference in ability to store electrical energy in a water-filled highpermittivity space vs an empty low-permittivity space, and this is related inter alia to the wettable volume of the pore as well as the local dielectric constant of the water model. As can be seen from Figure 3, this model is in good agreement with the simulation data. Model parameters are given Table 1 (with errors estimated from the Bayesian posterior distributions shown in SI Figure  S3). These parameters enable a quantitative analysis of the effects of E-field strength, nanopore conformation, and water model on the hydration probability of the channel. Thus, the zero-field free energy difference is close to zero across all water models for the open state structure with the mTIP3P water model slightly favoring the liquid state by ∼1 k B T, while the   Table 1 and SI Figure S3   Values listed represent the mean and standard deviation of the posterior distributions shown in SI Figure S3. The zero-field free energy difference, ΔΩ 0 , electric field coupling, m, and intrinsic electric field, E INT , are defined in the main text.  Figure S2A). The overall shape of the resultant electrostatic potential profile indicates that there is a substantial potential gradient between the intra-and extracellular openings of the M2 helix bundle, from −160 mV in the intracellular mouth well (IC in SI Figure  S2) to −100 mV in the extracellular well (EC in SI Figure S2, when a lipid bilayer is included in the calculation). This corresponds to an electric field of 20 mV/nm, which is close to the E INT value of ∼15 mV/nm found from the fit ( Table 1). The relative distribution of charged amino acids (SI Figure  S2B) therefore provides a molecular explanation of the resultant electrostatic profile, which could be explored further via channel perturbing mutations.
Ion Conduction through the Nanopore. In addition to determining the influence of the electric field on pore hydration, we also measured the resultant ion conduction through the channel (Figure 4; also SI Figure S4). For the closed (i.e., dehydrated) state the current due to either ion species (i.e., Na + or Cl − ) is approximately zero at low electric fields and begins to increase in magnitude only when the field exceeds the critical value required to hydrate the hydrophobic gate. Following this initial nonlinearity, the current then increases linearly with field strength for large electric fields. This implies that once the pore is hydrated, it acts as an Ohmic resistance. For the open state, the conductance is approximately Ohmic at lower field strengths (>50 mV/nm) reflecting the greater hydration of the pore. Interestingly, if we calculate the slope conductance at the highest field strengths (i.e., G max = dI/dV at E = +200 mV/nm) for both the closed and open states, these values of G max are of comparable magnitude. For example, for the mTIP3P water model the values of G max (for sodium ions) for both the open and closed states are in the range 250 to 300 pS (derived from the data presented in Figure 4; note that there is a relatively high uncertainty associated with the conductance estimate given the error bars on the high field current estimates in this figure). This suggests that once this model nanopore is fully ACS Nano www.acsnano.org Article hydrated and at relatively high E-fields, the overall shape and size of the pore dominates the limiting conductance, as discussed in early theoretical treatments of ion channel conductance. 51 Importantly, these results also indicate that pore dewetting prevents the formation of an ionic current in the closed state. This in turn implies that ion permeation itself does not cause pore hydration, given that the threshold field strength for hydration is dependent on the water model. This is supported by previous simulations of model nanopores, both highly simplified 31 and based on a β-barrel protein architecture. 38 In each case, imposition (via ionic concentration imbalance) of an electrostatic field across a hydrophobic nanopore first resulted in electrowetting of the pore. This was then followed by ion permeation; i.e., water enters the pore first, to be followed at a later stage (in both cases within <1 ns) by ion permeation. Subsequently, while in the electrowetted state, there were periods of several nanoseconds when the pore remained fully wetted, but no ions were present within or passed through the model nanopore. Furthermore, in each case dissipation of the electric field by movement of ions through the pore resulted in subsequent dewetting of the pore. Overall, these observations support the conclusion that the transbilayer E-field itself causes pore hydration, which in turn enables ion permeation, rather than vice versa. The value of G max from our simulations exceeds ∼100 pS ( Figure 4) which is more than 2 orders of magnitude greater than the physiologically observed single channel conductance (∼1 pS). 52 However, this difference likely reflects both the high E-field strength used in these simulations and also the absence of the intracellular and extracellular domains of the protein, which reduce the conductance substantially. Nevertheless, the Na + current is significantly larger than the Cl − current. For example, for the open state with the mTIP3P water model, the G max ratio for Na:Cl is >5; for the closed channel the corresponding ratio is substantially higher. This is in agreement with the cation selectivity of the 5-HT 3 receptor. 53 This is most likely due to the negatively charged residues near the intracellular (D−4′, E−1′) and extracellular (D20′) openings of the M2 helix bundle (SI Figure S2B), which make entry of Cl − ions into the pore highly unfavorable, while remaining attractive for Na + .
Kinetics of Liquid−Vapor Transitions. Having shown that an external E-field influences the hydration probability of a nanopore by changing the free energy difference between its liquid and vapor states, we next explored how the free energy The dependence of the mean survival times of the liquid (τ l , red) and vapor (τ v , black) states on the free energy difference (ΔΩ) between these states is shown. Each data point corresponds to an individual simulation, while the dashed line represents a linear fit to all data points, with confidence bands indicated as a shaded background. The free energy difference was estimated from the mean hydration probability, while the survival times were determined from the time series of hydration probability. Simulations with fewer than two observed transitions are not included in this plot.

ACS Nano
www.acsnano.org Article difference is controlled by the underlying rate constants. The kinetics of liquid−vapor transitions inside a hydrophobic gate can be examined in terms of the mean survival times of the liquid state, τ l , and vapor state, τ v , as they vary with the free energy difference between these two states. In Figure 5 it can be seen that τ l grows exponentially with ΔΩ, while τ v decreases exponentially at approximately the same rate. This trend is seen for both the closed and open states of the nanopore channel and across all four water models. The lines in Figure 5 necessarily cross one another at ΔΩ = 0, as an equal probability for both states implies equal survival times. The crossover point thus defines a characteristic time scale for liquid−vapor transitions, which for the M2 helix nanopore has a value of τ lv = τ l (ΔΩ = 0) = τ v (ΔΩ = 0) ≈ 5 ns. Importantly, this value is nearly identical for both the closed and the open states of the nanopore, despite their different pore radii. Furthermore, ΔΩ = 0 occurs at very different magnitudes of the electric field. This indicates that a time scale of ∼5 ns is likely to characterize liquid−vapor transitions of all pores that share the characteristic dimensions and surface hydrophobicity of the M2 helix nanopore considered here. A systematic survey of a range of nanopore models (see, e.g., refs 7 and 54) could reveal how this time scale may depend on channel properties. Such slow behavior (water has a typical Hbond lifetime on the order of picoseconds) indicates these oscillations arise collectively, reminiscent of, e.g., ultraslow reorientation of water wires with nanotubes 55 or in extended networks at the surface of a protein. 56 There is some slight variation of this oscillation time scale τ lv across water models, but it always falls in the range between 1 and 10 ns, indicating a limited force field sensitivity.
This also suggests that a total simulation time, which is an order of magnitude larger than the characteristic time scale (i.e., 10 τ lv ≈ 50 ns or above), is sufficient to characterize the hydration equilibrium behavior of the M2 helix nanopore systems. It is also worthwhile to compare the duration of the liquid and vapor states with the typical dwell times of the ions inside the channel pore (SI Figure S5). For ions that pass the closed conformation in its partially hydrated state (⟨ω⟩ = 0.25 to 0.75), the average dwell time is comparable to the typical time scale of liquid−vapor transitions (τ lv ≈ 5 ns). This implies that the intermittent periods during which the channel exists in the liquid state last sufficiently long enough to allow the passage of individual ions. For the open state, the typical ion dwell time lies below ∼1 ns. This is less than the average lifetime of the liquid state, which exceeds ∼1 ns. Consequently, even though most water models predict some degree of dehydration of the open state nanopore in the absence of an electric field (Figure 3), this state is still expected to allow the regular passage of ions and so is in agreement with our previous annotation of this structure as a functionally open state. 41

CONCLUSIONS
In summary, we have shown that a transbilayer electric field can induce wetting of the hydrophobic gate of a model nanopore derived from the 5-HT 3 receptor. However, a field strength of ∼100 mV nm −1 is required to start to hydrate the pore of the closed conformation of the channel. This corresponds to a potential difference of ∼0.9 V, i.e., substantially greater than the transmembrane voltage an ion channel would normally experience under physiological conditions. We also present a simple model which can quantitatively describe the observed increase in pore hydration probability with increasing electric field, with the asymmetry of the hydration probability response curve accounted for by an offset parameter, E INT , linked to the intrinsic electric field generated by the unequal distribution of charged amino acid residues along the length of the pore.
It is also useful to reflect on possible extensions to the current studies. One possibility would to be extend the comparison of water models to include polarizable models. 42 However, as yet for, e.g., AMOEBA 57,58 external E-fields are not available in the GPU-accelerated OpenMM toolkit, 59,60 because the external-field-induced-dipole interaction is not implemented. Furthermore, it is likely that anisotropic atomic polarizability of the water molecule 61 may be needed for more accurate treatment of electrowetting. Another possible extension to this work would be to explore electrowetting in a complete protein channel. This would be computationally demanding, but would be of interest, e.g., in terms of how a more complex distribution of charge within a protein might influence E INT .
The finding that the threshold voltage required for hydrating a hydrophobic gate depends on the orientation of the electric field also provides an attractive perspective for the design of artificial nanopores. Our results suggest that an excess of only five elementary charges (from amino acid side chains) near one opening of the pore is sufficient to create a substantial intrinsic electric field acting in the transmembrane direction. Therefore, by functionalizing the openings of a hydrophobic pore with static charges, 62 it should be possible to design a nanopore that permits ion transport in only one direction, thus rectifying ionic currents. Importantly, this rectifying property would apply not only to ions but also to water molecules. This not only extends the concept of the electric-field induced smart water gates 63 but may also enable the creation of membrane systems that are semipermeable to water in a switchable fashion.

Molecular Dynamics
Simulations. MD simulations employing a constant electric field were performed in GROMACS 2018 64 using the protocols described below. Protein structures were obtained from the PDB and missing atoms were added using WHAT IF. 65 A 100 ns long MD simulation using the MARTINI force field 66 of the M2 helix bundle nanopore together with DOPC lipids and water was used to place the nanopore in a lipid bilayer. 67 During this coarse-grained simulation lipid molecules self-assemble into a bilayer structure around the protein. Following the self-assembly simulation, the protein−bilayer system was converted back to an atomistic representation and was solvated in a 0.15 M NaCl solution. The system was then equilibrated through a 10 ns long MD simulation employing the mTIP3P water model 49 and the CHARMM36m force field. 68,69 Constant Electric Field Simulations. Production simulations were carried out using the CHARMM36m protein force field 68 and associated lipid parameters 68,69 in conjunction with the mTIP3P, 49 SPC/E, 70 TIP4P 49 or TIP4P/2005 50 water models. During these simulations, the effect of an external electric field, E, was modeled through an additional force, q i E, acting on each atomic partial charge, q i . The electric field was applied in the z-direction perpendicular to the plane of the membrane. Figure 2A illustrates this simulation setup.
It has been shown 47 that under periodic boundary conditions the application of an external field as described above is equivalent to the application of a transmembrane voltage, such that the voltage drop across the membrane is given by ΔV = −EL z , where L z is the length of ACS Nano www.acsnano.org Article the unit cell in the z-direction. For the system considered here L z ≈ 8.5 nm. The magnitude of the electric field, E, was varied in steps of 25 mV nm −1 between −250 mV nm −1 and 250 mV nm −1 for the closed state structure (PDB ID: 4PIR) and between −150 mV nm −1 and 150 mV nm −1 for the open state structure (PDB ID: 6DG8). In terms of transmembrane potential, this corresponds to a step size of 0.2125 V and maximum transmembrane voltages of ±2.125 V and ±1.275 V, respectively. In each case, the range was chosen to ensure that a complete transition from a (partially) dewetted to a fully hydrated pore could be observed.
Preliminary simulations were carried out to assess the stability of the simulation system under electric fields of this magnitude. In order to assess pore hydration properties under high-field (|E| > 150 mV nm − 1) conditions (under which conditions electroporation may occur 71 ), a harmonic restraining force of 1000 kJ mol −1 nm −2 was placed on the heavy atoms of all lipid molecules (after the initial 10 ns long equilibration period), which prevented electroporation.
Simulation Details. Three independent simulations of duration 150 ns for the closed conformation (PDB ID: 4PIR) and 50 ns for the open conformation (PDB ID: 6DG8) were carried out at each value of the electric field. Time integration was performed using a leapfrog integrator with a step size of 2 fs. Covalent bonds to hydrogen atoms were constrained using the LINCS algorithm, 72 and water model geometry was enforced using the SETTLE method. 73 The conformation of the protein was maintained close to its experimentally determined structure through the application of a harmonic restraining potential of 1000 kJ mol −1 nm −2 to all Cα atoms.
Long-range electrostatic interactions were treated through the smooth PME method 74 employing a real space cutoff of 1.0 nm and a Fourier spacing of 0.12 nm, with charge interpolation onto the grid through fourth order B-splines. Lennard-Jones interactions were switched off between 1.0 and 1.2 nm, and a long-range dispersion correction was applied to energy and pressure. All simulations were carried out in the isothermal isobaric (NPT) ensemble. The temperature was maintained at 310 K through a velocity rescaling thermostat 75 with a coupling constant of 0.1 ps, while a pressure of 1.0 bar was enforced through a semi-isotropic Parinello-Rahman barostat 76 with a coupling constant of 1.0 ps and a compressibility of 4.5 × 10 −5 bar −1 .
Analysis. Unless otherwise specified, MD trajectories were postprocessed using MDAnalysis version 0.19.0 77 while further data analysis and visualization were performed through custom scripts written in the R language (https://www.r-project.org/). Molecular visualizations were created in VMD. 78 As discussed above, the electrowetting of hydrophobic gates can be understood by considering the influence of an external electric field on the free energy difference between the liquid and vapor states of the pore (eq 2 above). Consequently, the hydration probability can be written as which constitutes a nonlinear model with three free parameters, namely, ΔΩ 0 (the free energy difference in the absence of an external electric field), m (which represents how strongly the probability of hydration is coupled to the electric field), and E INT (the horizontal offset accounting for the intrinsic electric field due to the protein; see above for more detailed discussion).
A fit of this equation to the observed hydration probability was obtained within the framework of Bayesian multilevel modeling. 79 Both ΔΩ 0 and m were treated as dependent on both the channel structure and water model, while E INT was considered to only depend on the channel structure. The zero-field free energy difference, ΔΩ 0 , may take both negative and positive values, and a Gaussian prior with mean μ = 0 and standard deviation σ = 10 k B T was used. In contrast to this, the coupling strength, m, is positive definite by definition. A gamma prior with shape α = 2 and rate β = 0.01 was chosen to enforce this constraint, and sampling for this parameter was limited to positive values. While the intrinsic electric field, E INT , can in principle be both positive and negative, allowing this parameter to vary over the entire real range caused divergent transitions in the sampling process. Sampling for this parameter was therefore restricted to positive values (in accordance with the observed data), and a gamma prior with shape α = 2 and rate β = 100 was chosen.
The model was implemented in the probabilistic programming language Stan 80 through the brms 81 library in R. Posterior samples were generated through a Hamiltonian Monte Carlo 82,83 procedure using the No-U-Turn (NUTS) sampler. 84 For each of 10 independent Markov Chains, 10 000 warm-up and 10 000 sampling iterations were carried out. The NUTS sampler used a target average acceptance ratio of 0.9 and a maximum tree depth of 15.

* sı Supporting Information
The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acsnano.0c04387. Figures S1−S5, showing electrowetting of the hydrophobic gate in the 5-HT 3 receptor comparing four water models; electrostatic potential of the M2 helix bundle of the closed state 5-HT 3 receptor; posterior probability of the Bayesian model parameters; cumulative charge flow through the M2 helix bundle of the closed state 5-HT 3 receptor (PDB ID: 4PIR) due to Na + ions; and dwell time of ions inside the M2 helix nanopore as a function of the E-field strength (PDF)