Water Nanoconfined in a Hydrophobic Pore: Molecular Dynamics Simulations of Transmembrane Protein 175 and the Influence of Water Models

Water molecules within biological ion channels are in a nano-confined environment and therefore exhibit behaviours which differ from that of bulk water. Here, we investigate the phenomenon of hydrophobic gating, the process by which a nanopore may spontaneously de-wet to form a ‘vapour lock’ if the pore is sufficiently hydrophobic and/or narrow. This occurs without steric occlusion of the pore. Using molecular dynamics simulations with both rigid fixed-charge and polarizable (AMOEBA) force fields, we investigate this wetting/de-wetting behaviour in the TMEM175 ion channel. We examine how a range of rigid fixed-charge and polarizable water models affect wetting/de-wetting in both the wild-type structure and in mutants chosen to cover a range of nanopore radii and pore-lining hydrophobicities. Crucially, we find that the rigid fixed-charge water models lead to similar wetting/de-wetting behaviours, but that the polarizable water model resulted in an increased wettability of the hydrophobic gating region of the pore. This has significant implications for molecular simulations of nano-confined water, as it implies that polarizability may need to be included if we are to gain detailed mechanistic insights into wetting/de-wetting processes. These findings are of importance for the design of functionalised biomimetic nanopores (for e.g. sensing or desalination), as well as for furthering our understanding of the mechanistic processes underlying biological ion channel function.


Introduction
The transport behaviour of water and other liquids in nanoscale pores is of both fundamental and technological importance. 1 Nanoscale pores across membranes are of especial interest. 2 In biological systems, channels are protein nanopores which enable the flow of water molecules and/or ions across cell membranes. 3 The behaviour of water molecules within nanoscale pores differs in several respects from that of bulk water. [4][5][6][7] In particular, ion channel proteins and nanopores may possess a hydrophobic gate (Fig. 1A). This is a constricted region of the pore which is lined with hydrophobic residues in which de-wetting may occur to form a "vapour lock". [8][9] This in turn presents an energetic barrier to ion permeation and thereby functionally closes the channel. Therefore, in these ion channels, the hydrophobic gate closes the pore to the passage of water and ions without the need for steric occlusion of the pore. 10 The presence of hydrophobic gates in several ion channel species has been inferred from both computational and experimental studies. Hydrophobic gating has also been designed into synthetic nanopores e.g. 11 . Thus, a detailed mechanistic understanding of hydrophobic gating is relevant both to studies of ion channel structure/ function relationships, 12 and also to the design of gating functionality 13 into synthetic nanopores. [14][15] The behaviour of a hydrophobic gate has been shown to depend critically on the radius of the pore and on the hydrophobicity of the pore lining. [16][17][18][19][20][21][22] Furthermore, in addition to these two structural parameters, the application of an electric field [23][24] or pressure [25][26] can cause the pore to wet, thereby opening an otherwise closed channel to the passage of water and ions.
Molecular dynamics (MD) simulations have been widely applied to explore the parameters affecting hydrophobic gating and the wetting/de-wetting behaviour both of ion channels and of simplified models of nanopores (see e.g. 6,12 for two recent reviews). By measuring the number density of water molecules within the gating region, an associated energetic barrier to wetting can be estimated. 27 This has been used to characterise the likely functional state of nearly 200 experimentally determined ion channel structures, 22 leading to a more global analysis of the dependence of hydrophobic gating behaviour (in terms of the energetic barrier to wetting) on channel structure. MD simulations have also been used to explore the effect of pore lining mutations on hydrophobic gating in e.g. the bestrophin-1 chloride channel. 28 When conducting MD simulations of water behaviour in nanopores and channels, a crucial but not yet widely investigated consideration is the choice of water model. For example, recent simulations have examined the effect of varying water models on the behaviour of water and ions within cyclic peptide nanopores 29 and within ligand gated ion channels. 30 Due to the competing effects of hydrogen bonding, van der Waals and other intermolecular interactions, water is a difficult molecule to model accurately. To date, there are over 130 different water models for use in MD simulations and related calculations, each with a range of physical properties which are comparable to experimental values. 6 Typically, these water models are developed to accurately reproduce one or more experimental properties of bulk water. In the case of water in hydrophobic gates, however, water molecules are confined within a region of ~1 nm diameter and ~2 nm in length. This means that there are likely to be only a small number of water molecules in the hydrophobic gating region, with a maximum of ~50 for a cylinder of this volume. Water confined within such a nanoscale channel is unlikely to behave as it does in bulk, e.g. in terms of its dielectric behaviour. 7 Consequently, it is not certain that water models developed to reproduce bulk behaviour will be directly transferable to such nanoconfined environments. Initial studies on the pore domain of the 5-HT 3 (5-hydroxytryptamine) receptor suggest that the choice of water model can have an effect on the wetting/de-wetting behaviour, particularly in the case of structures which are in an intermediate open/closed state. 30 (It should be noted that the 5-HT 3 receptor system corresponds to a modest degree of nanoconfinement, with ~17 water molecules in the key hydrophobic region of the pore). In particular, there are indications that the use of more sophisticated water models such as polarizable water models (the dipole moment of which will be responsive to changes in local environment) may be required to capture more adequately the interactions within hydrophobic gates. Although polarizable models are more computationally expensive, the development of more powerful computing facilities and of improved software enables the use of such models for complex biological systems. [31][32] Thus, recent studies have explored the influence of polarizable models on e.g. simulations of the interactions of ions with lipid bilayer membranes, 33 simple model ion channels e.g. 34 and on electric field strength within the active site of an enzyme. 35 The effect of polarizable water models within synthetic nanoconfined environments under an electric field has also been investigated, with the conclusion that polarizability can affect the structure and dynamics of the systems. [36][37][38] Various rigid fixed-charge water models have also been assessed for ligand binding in host/guest systems 39 and for water behaviour in peptide nanotubes. 29 Wetting/de-wetting of the well-defined hydrophobic gate in the TMEM175 (transmembrane protein 175) ion channel 28,[40][41][42][43] provides an extreme case of non-bulk behaviour of water when nanoconfined in a pore. In this structure the hydrophobic gate is particularly small, with only 1-4 waters able to be accommodated at any one time. We have therefore used this hydrophobic gate to assess how different rigid fixed-charge water models perform compared to the AMOEBA 44-45 polarizable water model in a highly nanoconfined environment.
In the current study we use the prokaryotic CmTMEM175 structure. 40 We investigate the effects of changing both the pore-lining residues and the water model employed in the simulations on wetting/de-wetting behaviour of this hydrophobic gate. This region of the channel is formed by three hydrophobic rings (each with four-fold symmetry): one ring of isoleucine residues (residue Ile23) followed by two of leucine residues (Leu27 and Leu30 respectively). We explore how mutating this hydrophobic gate to either three rings of valines (slightly smaller hydrophobic sidechains), or alanines (hydrophobic but much smaller sidechains), or asparagines (hydrophilic sidechains) affects the water behaviour. We also investigate the transferability of different rigid fixed-charge (i.e. additive) and polarizable water models for simulating water behaviour within the hydrophobic gate, i.e. we evaluate whether the different water models produce consistent wetting/de-wetting behaviours. Our results provide mechanistic insights into the behaviour of water within a Europe PMC Funders Author Manuscripts highly nanoconfined pore. This is of direct relevance to the ongoing interest in nanoconfined water as studied either experimentally (e.g. [46][47] or computationally (e.g. 7,48 ).

Channel System and Simulations
In this study we used the prokaryotic CmTMEM175 (PBD: 5VRE; 3.3 Å resolution) 40 as a model for a pore with a hydrophobic gate which can be wetted/de-wetted. 28 In this channel structure the first transmembrane helix from each of four subunits line a central pore with a constriction (with a minimum radius of < 0.15 nm on the crystal structure) formed by three rings of hydrophobic sidechains: Ile23, Leu27 and Leu30 (see Fig. 2). In two other TMEM175 channels these rings of residues are either conserved (in hTMEM175, PDB: 6WC9) or replaced by three rings of leucines (in MtTMEM175, PDB: 6HD8). However, the structural conservation of this region is less clear-cut, [40][41][42] as might be expected for a hydrophobic gating region.
In order to explore systematically some of the parameters affecting hydrophobic gating (e.g. pore radius and hydrophobicity), the CmTMEM175 structure (PDB code 5VRE) was embedded within a phosphatidylcholine (PC) bilayer and solvated with a ~0.15 M NaCl aqueous solution. In these simulations, backbone restraints were applied to maintain the experimentally observed polypeptide backbone conformation whilst allowing a degree of sidechain flexibility/mobility. In initial simulations of the wild type (i.e. native) protein the gating region was seen to be sufficiently narrow and hydrophobic that the pore spontaneously de-wetted early on in simulations ( Fig. 1B and SI Fig. S1A).
As can be seen (Fig. 2), the hydrophobic gate in the CmTMEM175 structure is formed by a constricted region of three rings of hydrophobic sidechains, forming an 'ILL' motif.
In order to explore the sensitivity of wetting/de-wetting behaviour to the nature of this motif, these residues were in silico mutated to either three rings of valines (VVV), i.e.
hydrophobic sidechains but smaller than isoleucine or leucine, or three rings of alanines (still hydrophobic, but substantially smaller sidechains), or three asparagines (the sidechain of which is the same size as that of leucine, but is polar and able to form H-bonds): see Fig.  2A. For both the wild type and the three mutant channels we conducted simulations with various water models (see below). Simulations were of 100 ns duration with 6 repeats for fixed-charge water models and of up to 50 ns duration with 3 repeats for the polarizable AMOEBA model.

Effect of Mutations on Wetting/De-wetting
In simulations of the wild-type channel, the hydrophobic gate is associated with a constriction of radius down to 0.1-0.2 nm (Fig. 3A). This is comparable to the radius of a water molecule (0.14 nm). The sidechains of the isoleucines and leucines result in a highly hydrophobic lining. Thus, in simulations using the TIP4P/2005 water model and OPLS all-atom protein force field (see Methods below), the gating region de-wets within the first ~10 ns for most of the repeats (SI Fig. S1B). The gating region remains relatively dry for the remainder of each of the simulations. However, in some repeats of the simulations, a water molecule is temporarily trapped within the gate, due to the tight constrictions of ~0.1 nm at each end of the gate. Consequently, the time averaged water density ( Fig. 3B) does not completely fall to 0 nm -3 within the hydrophobic gate of the wild-type channel (this will be discussed in more detail below). The energetic barrier to water of this region ( Fig. 3C) is 13.5 (± 4.8) k B T (averaged over the 6 repeats, ± standard deviation). This conformation of the pore therefore clearly corresponds to a functionally closed state.
The constriction of the pore in the 'VVV' mutant remained hydrophobic but was a little wider (by ~0.1 nm at the lower/extracellular (i.e. s = -1.3 nm) end of the gate (Fig. 3A).
Despite this small increase in radius, the pore remained de-wetted (Fig. 3B), with an energetic barrier comparable to that of the wild-type channel (Fig. 3C).
In contrast, although the AAA mutant preserved the hydrophobic nature of the pore lining, it was wider, with two minima (Fig. 3A) of ~0.2 and ~0.3 nm radius. This was sufficient for the pore to remain fully wetted (Fig. 3B), thus presenting no barrier to ion permeation (Fig.  3C). The more polar/hydrophilic mutant, NNN, had a pore radius profile comparable to that of the VVV mutant ( Fig. 3A). However, the hydrophilic lining provided by the Asn (NNN) sidechains (to which water can form H-bonds) was such that the pore remained fully wetted ( Fig. 3B) with no energetic barrier to water permeation (Fig. 3C).

Effect of Rigid Fixed-charge Water Model on Wetting/De-wetting in the Wild-type Channel
The results discussed so far employed the rigid fixed-charge (i.e. additive) TIP4P/2005 water model, which has been shown to perform well over a range of (bulk) water properties [49][50] and has been used in recent studies of nanoconfined water e.g. 51 . To explore the robustness of these results to variations in rigid fixed-charge water models, we investigated the effect of five different models on the de-wetting behaviour in the hydrophobic gate of the wild-type CmTMEM175 structure (Fig. 4). These models were chosen to cover a broad range of model types (see e.g. 6 For all five rigid fixed-charge models, the pore de-wetted as evidenced by the water density profiles (Fig. 4A) and the associated free energy barriers (Fig. 4B). In the case of TIP4P/2005, some water molecules were trapped by the constrictions at the ends of the hydrophobic gate, as discussed above, resulting in a decrease in the free energy barrier. This difference is however just within the limit of the standard deviations (averaged over repeats) for the five water models (see SI Figs. S2, S3 and S4, and therefore we attribute this to stochastic variation rather than a significant difference due to the water model. We therefore Europe PMC Funders Author Manuscripts conclude that the wetting/de-wetting behaviour is not substantially affected by changes in additive water model, and that these water models are equally valid for investigating the effects of hydrophobic gating. In the absence of direct experimental data, it is difficult to assess whether the water models used are transferrable from bulk to nano-confined (i.e. channel) systems. However, recent experimental studies have attempted to estimate water profiles within the hydrophobic gate of an ion channel 58 (but not yet for TMEM175). In the meantime, we can state that the five additive water models explored appear to be equally transferrable when moving from bulk water to water confined in a nanopore.
Comparing the five rigid fixed-charge water models in the mutants also produced consistent results. For example, in the case of the AAA (

Effect of a Polarizable Water Model on Wetting/De-wetting
We are especially interested in how polarizable water would behave in the TMEM175 hydrophobic gate structure given the stochastic observation of a 'trapped' water (or waters) in this extreme example of nanoconfinement, as seen in our simulations using e.g. the TIP4P/2005 water model (see above and SI Fig. S1B). Experimentally, water molecules in hydrophobic nanoslits exhibit an anomalously low dielectric constant largely reflecting their rotational immobilisation. 46 Simulation studies have been used to compare a rigid fixed-charge (SPC/E) water model and a flexible polarizable (ABEEM-7P) model for water molecules confined in graphene nanocapillaries with a width ranging from 0.59 to 1.1 nm. In the latter case the polarizable model resulted in more ordered water than with SPC/E. 59 Our simulations of the extreme nanoconfinement within the TMEM175 gate (the radius in the middle of the hydrophobic gate is ~0.16 nm for the wild-type (WT) channel, comparable to that of a single water molecule) reveal a difference in the de-wetting behaviour when comparing the polarizable AMOEBA14 water model 45 to the five additive water models. This is most clear cut for the WT channel, for which the pore radius profiles within the hydrophobic gate are identical for the five non-polarizable and for the polarizable (AMOEBA) models (SI Fig. S2A). For this channel, the number density of water in the nanocavity in the centre of the hydrophobic gate (Fig. 6A) was ~0 nm -3 for TIP4P, ~6 nm -3 for TIP4P/2005 and ~15 nm -3 for AMOEBA, corresponding to free energies along s of ~18 (±4), ~13 (±5) and ~7 (±4) kBT respectively, where the numbers in parentheses are ± 1 standard deviation averaged over repeat simulations (Fig. 6B). Thus, water within this extreme nanoconfinement ( Fig. 7A;

Europe PMC Funders Author Manuscripts
For the VVV simulations (SI Fig. S9 and S10) the situation was less clear. With this (modelled) mutant pore there was a degree of dependence of the pore radius profile on the water model employed which complicated such detailed comparison.
Overall, this suggests that the use of a polarizable water model has a significant effect on the de-wetting/wetting behaviour of a 'tight' hydrophobic gate, and that future simulations of hydrophobic gating in general would benefit from a systematic investigation into the effects of including molecular polarizability. It is also worth noting that, whilst the protein/ lipid/ion force field was consistent for the rigid fixed-charge water models (OPLS with united-atom lipids), in the case of our polarizable simulations all atoms were modelled with the AMOEBA force field (see Methods). Therefore, it cannot be ruled out that the induced polarizability of the protein also affected the wetting/de-wetting behaviour. This could be assessed by running simulations using a combination of polarizable and non-polarizable force fields within a simulation. Daub et al. have investigated the combination of a polarizable water model with non-polarizable ions in a silica nanopore, and their results suggest that the polarizability of the water has a greater effect over that of the ions. 36 We can make partial comparisons with other studies. Previous simulation studies of the effect of including polarizability focussed on a less extreme case of nanoconfinement 30 involving the gating region of a M25 nanopore derived from the (partially) open state of the 5HT 3 R channel (PDB 6DG8). For this channel system the time-averaged water density was respectively ~22 nm -3 (for TIP3P water), ~16 nm -3 (TIP4P/2005) and ~35 nm -3 (AMOEBA). This corresponds to an approximate estimate via the Boltzmann equation (where ρ p and ρ NP are the water densities for the polarizable and non-polarizable models respectively): ΔU = − RT ln ρ P /ρ NP of stabilisation of the nanoconfined water within the hydrophobic gate of ~-1.6 kJ/mol by including polarization. A comparable estimate for TMEM175 yields a stabilisation of ~-4.2 kJ/mol by including polarization. An approximate calculation of the relative potential energies of a water dipole in a single molecule sized cavity in dielectrics of 1 (i.e. in vacuo) vs. 4 (for a hydrophobic cavity in a protein) yields about -7 kJ/mol stabilisation of the latter case. 60 Thus, the same trend is seen for the two nanopore systems but the stabilisation by including polarization is more marked in the extreme nanoconfinement case of the closed hydrophobic gate of TMEM175, which can accommodate ~1-4 water molecules compared with the (partially) open hydrophobic gate in the 5HT3R 6DG8 structure which can more readily accommodate multiple water molecules.
Comparisons with experimental data remain difficult but we note that measurements of the refractive index of nanoconfined water (corresponding to its dielectric constant at optical frequencies 47 ) suggest a decrease in molecular polarizability corresponding to a reduction in refractive index from n = 1.33 (bulk) to n = ~1.25 for water confined in ~10 nm gaps of hydrophobic interfaces. It would therefore be of great interest to know how the refractive index of water behaves on nanoconfinement scales of < 1 nm.

Trajectories of De-wetting
We examined the de-wetting behaviour of trapped waters in the hydrophobic nanocavity using the polarizable AMOEBA14 model in a little more detail ( Fig. 7 and SI Fig. S11). We tracked the trajectories projected onto the z axis of all water molecules initially in the hydrophobic gate region. In two of the three simulations, 3-4 waters were initially present but 1-2 of these exited within the first ~40 ns. In the third repeat the hydrophobic gate was fully de-wetted at the start and remained so for the duration (50 ns) of the simulation (SI Fig.  S11).
We also examined the first repeat of these simulations in more detail (Fig. 7). At the start of the simulation three water molecules were present within the central hydrophobic nanocavity. This is divided into two 'sub-cavities', one formed by the I and L sidechain rings and one by the L and L sidechain rings. Just before 40 ns, one water molecule exits from the hydrophobic gate via the wider L side. The remaining two water molecules then continue to switch between the two sub-cavities: sometimes one in each (e.g. Fig. 7D), sometimes both water molecules in one sub-cavity (e.g. Fig. 7E), with the I-L sub-cavity being much more commonly occupied. In the latter case the two waters form a hydrogen bond. Similar behaviour was seen in the second repeat simulation (SI Fig. S11B), during which two water molecules leave early on (< 10 ns) and then two remaining waters switch back and forth dynamically between the same or different nanocavities.
By contrast, in the TIP4P simulations (SI Fig. S1A), in 5 (out of 6) simulations trapped waters were present but were expelled in <15 ns in each case. In the TIP4P/2005 simulations (SI Fig. S1B), in 4 (out of 6) simulations initial trapped water molecules were expelled, whereas in one repeat they mostly remained trapped for the entire simulation. Taken together, these results are consistent with the overall instability of the wetted state of the hydrophobic nanocavity. A metastable state with one or more trapped waters may exist, the exact stability of which is sensitive to the water model employed.

Conclusions
From the results described above, it can be seen that the broad picture of hydrophobic gating (explored for a closed state of a hydrophobic gate in an ion channel protein) is robust to changes in additive (rigid fixed-charge) water models. However, quantitative details alter when electronic polarizability is included in the simulations. This in turn suggests the need to include consideration of polarizability (see e.g. 61 for a recent example) when designing hydrophobic gates or comparable structures into nanopores. 14,62 Relatively small changes in detailed nanoscale structure (and hence in the local dielectric/polarizability properties) could be used to 'fine tune' a hydrophobic gate to e.g. electrowetting.
It would be interesting to consider how an applied electric field or pressure would affect differences in water model on hydrophobic gating in TMEM175. Under physiological conditions, with a membrane potential of ~50-200 mV, this would equate to applied electric fields of ~0.005-0.02 V/nm or pressures of ~2-9 MPa. In a previous study on the 5-HT 3 receptor channel, the effect of supraphysiological electric fields was investigated using different rigid fixed-charge water models; 24 this implied that different water models did Europe PMC Funders Author Manuscripts affect the wetting/de-wetting behaviour. Under high pressures of 110-220 MPa, it has been found that water flow rates through a carbon nanotube significantly depend on the rigid fixed-charge water model employed. [63][64] However, to our knowledge, there have not yet been any studies into the effect of water models on hydrophobic gating in ion channels under physiologically relevant electric fields or pressures. Based on our previous experience, it is possible that under an applied electric field or pressure, the choice of water model could have an effect on hydrophobic gating, particularly if additive and polarizable water models were compared. This would be of importance to the development of designing electrowetting pores, particularly in the case of synthetic pores.
In this study we have used an X-ray structure of a prokaryotic TMEM175 structure (CmTMEM175 at 3.3 Å resolution; PBD 5VRE 40 ) as a model for investigating the effect of mutating the pore lining and water model on wetting/de-wetting behaviour. CryoEM studies of the human TMEM175 channel 41 reveal two possible open and closed structures. Further TMEM175 structures (MtTMEM175) have been solved by X-ray diffraction to 2.4 Å resolution (PDB 6HD8 42 ). These differ from the 5VRE CmTMEM175 structure in some respects. In particular, the pore-lining isoleucine (Ile23) in CmTMEM175 is replaced by a leucine (Leu35) in MtTMEM175. Together these structures indicate conservation of a hydrophobic nanoscale constriction in the transmembrane pore.
The current studies could be extended in the future by applying polarizable (AMOEBA) simulations to characterize water behaviour in a range of channel structures selected to have hydrophobic gates of differing dimensions. 22 This could enable us to explore the effect of the degree of nanoconfinement on the extent of anomalous water behaviour.
Of course, one must also reflect on whether polarizable models are better able to replicate experimental data than additive models. The AMOEBA14 water model is able to accurately reproduce experimental bulk water properties across a range of temperatures. 45 In the context of simulations involving several chemical entities, we note that e.g. local accumulation of anions at water/graphene interfaces as demonstrated experimentally 65 requires a simple model of polarisation (namely charge scaling 66 ) to be reproduced in simulations. Furthermore, inclusion of polarizability (using AMOEBA) seems to be necessary to provide a consistent and accurate explanation of ion permeation through potassium channels. 67 Combining recent computational (this study) and experimental (e.g. 58 ) approaches would therefore allow us to exploit biological nanopores/channels to characterise the effects of confinement on water behaviour on smaller scales than those currently addressed (down to ~1.5 nm) by studies of nanofabricated devices. 46 In the limiting case, this may enable us to use biological systems to examine the behaviour of water under extreme nanoconfinement with just one or two water molecules present. This is of relevance to a number of areas of enquiry, especially given the continued interest in the behaviour of water confined in nanoporous environments. [68][69] Furthermore, in a recent experimental study 70 it has been suggested that nanocavities may be present within other membrane proteins (e.g. GPCRs) and may contribute to the activation mechanism of these receptors. Thus, it is important to understand membrane protein nanocavities both from the perspective of nanopore design, and also for understanding mechanisms of protein nanoswitches within biological membranes.

Protein structure preparation
The prokaryotic TMEM175 structure from Chamaesiphon minutus (CmTMEM175 PBD: 5VRE) was run through the WHATIF server (https://swift.cmbi.umcn.nl/whatif/) to fix any missing sidechains. 71 The three rings of pore-lining residues forming the hydrophobic gate in the wild-type (Ile23, Leu 27 and Leu30 on each of the four subunits) were then mutated in PyMol (http://www.pymol.org) to three rings of alanines, asparagines or valines to form the three mutant structures ('AAA', 'NNN' and 'VVV' structures).

Non-polarizable simulations
TMEM175 structures were embedded in a POPC (1-palmitoyl-2-oleoyl-sn-glycero-3phosphocholine) bilayer using the MemProtMD and CG2AT protocols. [72][73] This uses a multi-scale methodology to solvate the system with ~0.15 M NaCl, self-assemble the lipid bilayer and equilibrate the system. Production runs involving the additive water models were performed using GROMACS 74 version 2016.3 with the OPLS all-atom protein force field with united-atom lipids. 75 We chose to use the OPLS force field because this was developed to accurately reproduce experimental water properties within biological systems. Each simulation was run for 100 ns with an integration time step of 2 fs. For the CHAP analysis (see below), the first 10 ns were discarded as equilibration. 6 repeats were run for each protein structure and water model combination, resulting in a total of 120 simulations. The temperature was maintained at 310 K using the V-rescale thermostat 76 and a time coupling constant of 0.1 ps. The pressure was maintained at 1 bar using the semi-isotropic Parrinello-Rahman barostat 77 and a time coupling constant of 1 ps. The Verlet cut-off scheme was used, along with the Particle Mesh Ewald method for electrostatics. 78 The LINCS algorithm was used to constrain bonds. 79 The positions of the protein backbone atoms were position restrained with a force constant of 1000 kJ mol -1 nm -2 , however protein side chains were allowed to relax. We chose to apply these position restraints because we did not want our overall simulation protein structure to deviate from the original experimental coordinates. It was important to do this because we specifically wanted to assess how changes in the protein pore lining (i.e. via mutants) and water model affected the wetting/de-wetting behaviour within a given overall protein geometry.

Polarizable simulations
These were performed using a similar method to that described in Klesse et al. 30 (https:// github.com/Inniag/openmm-scripts-amoeba). The wild-type and VVV mutant structures were embedded in a DOPC bilayer and solvated with ~0.15 M NaCl solution using the bilayer self-assembly method outlined for the non-polarizable simulations. A 10 ns equilibration simulation was then run using GROMACS with the non-polarizable CHARMM36m force field. Polarizable production runs were performed using OpenMM, 31 and lasted for 50 ns for the wild-type structure and 35-50 ns for the VVV mutant structure. The AMOEBA force field was used for all atoms in the simulation i.e. for the protein, 80 ions 81 and lipids, 82 together with the AMOEBA14 water model. 45 The r-RESPA algorithm was used with an outer timestep of 2 fs and an inner timestep of 0.25 fs. The pressure was maintained at 1 bar using the Monte Carlo barostat and the temperature was maintained at 310 K using the Andersen thermostat. The Cα backbone atoms were position restrained with a force constant of 1000 kJ mol -1 nm -2 for both the equilibration and production runs. There were 3 repeats for each structure, resulting in a total of 6 simulations.

Analysis
Profiles for pore radius, water density and free energy were obtained using ChAP (www.channotation.org). 27 In all cases, the first 10 ns were discarded from the analysis and treated as equilibration. A bandwidth of 0.14 nm was used for calculating the water density through the pore. The free energy profiles were obtained via: where G(s) is the free energy of wetting/de-wetting, n(s) is the water density, s is the coordinate passing down the central axis of the pore (approximately perpendicular to the bilayer), T is the temperature, k B is the Boltzmann constant and C is an arbitrary constant set so that G(s) = 0 outside the pore. Having first discarded the first 10 ns of the simulation as equilibration, these profiles were averaged over all remaining frames of the simulation, and then averaged over all repeats as appropriate.
The positions of individual water molecules were calculated with the aid of MDAnalysis. [83][84][85] It should be noted that the protein is oriented such that z axis is approximately antiparallel to the local pore s defined by CHAP and so for all practical purposes z = -s. Figures were created using VMD 86

Supplementary Material
Refer to Web version on PubMed Central for supplementary material. A The principle of hydrophobic gating: in the wetted state, both water (medium blue) and ions (red) can pass through the ion channel (pale blue) and the hydrophobic gate (orange); whereas in the de-wetted state, the channel, though not physically occluded, is closed to the passage of ions. B An example of simulated de-wetting of a hydrophobic gate. The trajectories of water molecules (dark blue) around and initially located in a hydrophobic gating region (indicated by the orange arrow, centred at z = 0) are shown projected onto the A The simulation set-up with the TMEM175 protein in blue and the wild-type hydrophobic gating residues highlighted in yellow (LEU) and orange (ILE). The LEU and ILE residues in the hydrophobic gate (LLI) are substituted for ALA, ASN and VAL respectively to create three mutants (AAA, NNN and VVV respectively). The position of the lipid headgroups in the bilayer are denoted by the dashed black line, and s indicates the path direction through the centre of the pore. B A close-up looking down the wild-type pore. C A slice through the pore with water (red spheres for oxygen, white spheres for hydrogen) exiting the hydrophobic gating region of the pore surface (blue).  Profiles of pore radius (A), water density (B) and water free energy (C) as a function of distance along the pore centre line, s. The profiles for the wild-type structure (WT, green) and AAA (black), NNN (red) and VVV (blue) mutants are shown, each averaged over the last 90 ns of the simulation and subsequently averaged over 6 independent repeats. The water model TIP4P/2005 was used for all simulations. The hydrophobic gating region is indicated by the orange band, and the approximate positions of the lipid head groups are denoted by the grey dashed lines.   Europe PMC Funders Author Manuscripts

Europe PMC Funders Author Manuscripts
Europe PMC Funders Author Manuscripts Comparison of water density (A) and free energy (B) profiles for the TIP4P/2005 fixedcharge and the AMOEBA polarizable water models for the wild-type TMEM175 structure. The shaded regions denote ±1 standard deviation calculated over all simulation frames included in the analysis (discarding the first 10 ns as equilibration), and averaged over all repeats. The dashed black line denotes the value of bulk water density.  A Schematic of the 'extreme nanocavity' in the TMEM175 hydrophobic gate, showing the cylindrical approximation to the WT cavity volume. B Snapshot of three 'trapped' water molecules (coloured blue, pink and green) within the nanocavity which partially hydrate the hydrophobic gate in the AMOEBA simulations. C Trajectories of the three trapped waters (coloured as in B) for an AMOEBA simulation. The pink water molecule can be seen to leave after ~40 ns. D and E Snapshots of the simulation in C immediately after and ~5 ns after the pink water molecule has departed from the nanocavity, showing the relaxation of