Water in Nanopores and Biological Channels: A Molecular Simulation Perspective

This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and β-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.


INTRODUCTION
Biological channel proteins form nanoscale pores in cell membranes. 1,2 They are of intrinsic physiological importance and also provide design templates for controllable nanopores in synthetic membranes. Biological channels and nanopores typically have an internal radius of ∼0.5 nm and a length of ∼5 nm. They are filled with water, thus providing permeation pathways across a lipid bilayer membrane for selected ions and/or uncharged low-molecular-weight solutes. Their functional properties (i.e., conductance, selectivity and gating) are therefore dependent on the behavior of water in a nanoconfined environment. To understand the relationship between structure and function of ion channels, and also to aid the design of novel bioinspired nanopores 3,4 which may form components of "smart" membranes for e.g. water treatment and biosensing, 5,6 it is therefore important to understand how such nanoconfined water behaves. Molecular simulations play a key role in understanding the dynamic behavior of water in channels and nanopores. We will review simulation studies alongside key selected experimental investigations of the behavior of water in membrane-spanning nanopores and channels. Our bias will be toward biological or biomimetic pores, especially in lipid bilayer membranes (both synthetic and cellular). We will not cover nanopores in graphene sheets, which have recently been discussed elsewhere (e.g., see He et al. 7 ). Our focus will be largely on the past decade during which advances in molecular simulation have been paralleled by advances in our understanding of the structures of many ion channels and biological pores.
The overall aim of this Review is to explore the relationship between the molecular structure of channels and nanopores, and the dynamic structural behavior of water within them. In particular, we wish to understand the fundamental behavior of water in membrane nanopores and channels, and how this differs from the behavior of bulk liquid water under ambient conditions. This will help us to annotate both existing and newly discovered ion channel and pore structures with respect to their biological function. 8 Furthermore, it will aid and enable the design of novel nanopores such as artificial water channels. 9 Our account is, therefore, complementary to recent accounts of nanopores in a range of advanced materials and biotechnological applications. 10,11 2. NANOPORES AND CHANNELS: STRUCTURES AND SIMULATIONS

Basic Properties
The interactions of water and ions with nanopores and ion channel proteins determine their key functional properties, namely, their permeability and selectivity ( Figure 1). As we shall see, the relationship between nanopore dimensions, the chemical nature of the pore-lining surface and the pore occupancy with water is of particular importance. The dynamic behavior of water within nanopores is such that they can support high, near-diffusion-limited water permeation rates, for example, ∼10 −13 cm 3 /s (i.e., ∼3 water molecules/ns) for a single aquaporin pore. 12 Many biological nanopores exhibit selectivity for water and for other solutes or ions passing through a given pore. Water molecules may also play a key role in "gating" mechanisms of biological nanopores, that is, mechanisms for switching a pore between an open (permeable) and a closed (impermeable) state. This transition may correspond either to changes in the physical dimensions of a pore or to changes in the free energy landscape of water within a pore (e.g., in hydrophobic gatingsee below).

Systems
Nanopores and channels studied via molecular simulations may be classified as follows: (i) simplified models of nanoconfined water which omit most molecular details; (ii) nonbiological nanopores, including carbon nanotubes (CNTs), synthetic (biomimetic) nanopores, synthetic peptide nanopores, and track-etched nanopores in polymer membranes; (iii) porins and related nanopores formed by β-barrel membrane proteins; (iv) aquaporins and related protein Figure 1. Schematic diagram of a nanopore or channel (dark blue) in a membrane (gray; lipid headgroups in dark gray). Water is shown in pale blue, with the orange region inside the nanopore indicating the region of confined water. Flow of ions and water molecules through the pore is indicated by the red arrows.
pores permeable to water and small polar solutes; and (v) ion channels. This classification mirrors the structural complexity of nanopores and also broadly reflects the emphasis and development of simulation studies of behavior within them.

Structures and Experimental Approaches
The main experimental approach for determining highresolution structures of biological nanopores and channels is X-ray crystallography. More recently, cryo-electron microscopy is having an enormous impact on the determination of the structures of channels 13 and related proteins, albeit sometimes at medium resolution (∼3 Å). There are limitations of this resolution when using such structures as the basis for detailed simulation studies, including uncertainties of exact side chain conformations, and lack of experimentally defined initial positions for bound water molecules and ions. It is likely that these limitations will become less of an issue as cryoelectron microscopy of membrane proteins improves and more high-resolution structures (i.e., 2.5 Å and better) become available. For smaller membrane protein channels, NMR structures are also available. 14 In a number of structures of ion channels and pores determined by high (better than 2 Å) resolution crystallography, the low energy positions (temperatures of ∼100 K) of most of the water molecules in pores are revealed. A selection of channel structures (some of which have formed the basis of simulation studies of watersee below) are shown in Figure 2 in order to illustrate the data which structural studies provide, and against which simulation studies may be compared. In a high-resolution (sub-Å) structure of the yeast aquaporin pore Aqy1 (PDB ID 3ZOJ 15 ), there is a well-defined single file of water within the pore (which has a minimum radius R MIN < 0.1 nm) and details of the water−water and water−protein hydrogen bonds (albeit at ∼100 K) can be reliably defined. As discussed below, these interactions play a key role in the high water permeability and selectivity of aquaporins. The 1.5 Å resolution structure of a bacterial porin (see section 6.1 below) Omp32 (PDB ID 2FGR 16 ) corresponds to a wider (R MIN = 0.22 nm) protein nanopore. This structure of an anion-selective porin shows that water molecules in the crystal structure nearly completely fill this wider transbilayer pore. Several structures of potassium channels have been determined at high resolutions and, thus, provide details of water molecules within a highly ion-selective pore. These structures include that of the canonical bacterial potassium channel KcsA at 2.0 Å resolution. 17 The structure of the related bacterial potassium channel MthK at 1.45 Å resolution (PDB ID 3LDC 18 ) is shown in Figure 2C. The positions of both potassium ions and of the water molecules solvating the ions within the pore are thus defined. However, one should recall that all three of the above structures were determined at low (∼100 K) temperatures, as is true for most crystal structures of membrane proteins. Recently, the use of an X-ray free-electron laser (XFEL) source has enabled the structure of the protonselective influenza M2 channel to be determined both at an exceptionally high (1.45 Å) resolution (for a membrane protein) and at room temperature (PDB ID 5JOO 19 ). This structure is notable in that a continuous hydrogen-bonded network of water molecules is seen to span the length of the channel (see below for further discussion of M2) even at room temperature. This helps to support the inference that details of water occupancy of membrane pore and channel proteins seen in low-temperature X-ray structures are relevant to the function of these proteins at room temperature. More generally, from these high-resolution structures of membrane channel proteins we can judge both the strengths and weaknesses of structural studies of water in biological nanopores and channels. High-resolution structures provide details of water−water and water−protein hydrogen bonds within pores. However, they provide at best limited information concerning the dynamics of these interactions and the implications of the dynamic behavior for function. For this, one must turn both to other experimental techniques and, particularly, to molecular simulations.
A number of other experimental techniques have been used to identify bound waters within protein pores, including, for example, neutron diffraction for detecting functionally important wet spots in ion channels. This has been applied to KcsA and influenza M2 channels, yielding one-dimensional water profiles along the length of the pore while embedded in a lipid bilayer. 20 Synchrotron-based radiolytic labeling combined with mass spectroscopy may be used to identify water molecules bound within cavities in pore and channel proteins 21 and have been applied to, for example, a bacterial potassium  Table 1). (A) Aquaporin Aqy1 (PDB ID 3ZOJ) at a resolution of 0.88 Å, (B) the porin Omp32 (PDB ID 2FGR) at resolution 1.5 Å, and (C, D) two ion channels MthK (PDB ID 3LDC) and influenza A M2 (PDB ID 5JOO), resolved at 1.45 and 1.41 Å, respectively. Each protein backbone is in ribbon representation, colored from light to dark gray in the direction of each polypeptide chain. Oxygen atoms of water molecules within the pore and around the mouth region are shown as red spheres. In the case of MthK, potassium ions in the selectivity filter and central cavity of the pore are shown in purple. All other water molecules, ions, and ligands present in the crystal structures are omitted for clarity. Chemical Reviews pubs.acs.org/CR Review channel to explore differences in water exposure of the pore of closed and open states of the channel. 22 Insights into the local conformational dynamics of ion channels and their interactions with water and permeant ions can be obtained by twodimensional infrared spectroscopy, as seen for, for example, gramicidin A and protons 23 and for KcsA and potassium ions. 24,25 When combined with structural studies, these investigations provide a range of experimental data sets against which to compare the outcomes of simulation studies of water within protein nanopores and channels.

Simulations
Levels of molecular simulation range from those employing classical mechanics to quantum mechanics, that is, from modeling the atoms as charges with accompanying short-and long-range interactions to modeling the nuclei as point charges with electronic interactions modeled explicitly through the inclusion of electronic wave functions. These can be used to access different length and time scales (Figure 3; after Cipcigan et al. 26 ). Quantum mechanical simulations, while able to model explicitly a wider range of interactions, including electron transfer, tunnelling and photoexcitations, are currently computationally expensive for biological systems. Instead, MD simulations are favored as these provide a compromise between computational accuracy and efficiency. Hybrid quantum mechanical/molecular mechanical (QM/ MM) techniques have also been used to simulate biological systems. In this case, the region of interest is treated quantum mechanically and embedded in a surrounding classical region. For information on the mathematical formalisms, the associated practical pitfalls, and applications, the reader is directed toward the review of Brunk and Roẗhlisberger. 27 Further reviews on quantum mechanical and hybrid quantumclassical mechanical approaches available for simulating biological molecules include those by Cui 28 and Senn and Thiel. 29 In terms of complete quantum mechanical treatments, a number of large-scale density functional theory (DFT) calculations have been conducted on biological systems, including enzymes and metalloproteins, as discussed in depth by Cole and Hine. 30 However, there have not been many studies of channels and pores. Water is a difficult substance to model because of the competing effects of hydrogen bonding (which tends to open up the molecular structure) and van der Waals interactions (which tend to close the structure). An accurate description of hydrogen bonding involves consideration of the electronic structure of water, including the effects of electrostatics, charge transfer and polarization. 31 Nuclear quantum effects also play an important role and they are often not adequately included in traditional MD force fields. 32 With the use of MD and corresponding water models, there is at best a cancellation of errors which mitigates for the lack of inclusion of quantum effects. 33 A recent paper 34 on a quantum water model (q-SPC/ Fw) concluded that nuclear quantum effects are important, possibly as important as polarization, and should be included in subsequent polarizable water models. When comparing their quantum model to previous nonquantum models, the authors concluded that quantum reparameterization of rigid water models may be sufficient for modeling the quantum effects of water.
The multitude of interactions and methodologies for modeling these have important repercussions for accurately reproducing the water phase diagram both for the liquid phase and the multitude of ice phases (for experimental phase diagrams see, for example, those by Wagner, Feistel, and Soper, 35−37 and for a comparison of simulations to experimental data see, for example, the study by Abascal and Vega 38 ). By extension, this may also have a significant effect on the behavior of water in simulations of biological systems. For example, Anandakrishnan et al. have compared the influence of two explicit solvent water models and one implicit solvent model on protein folding of the CLN025 mini-protein. 39 The free energy landscapes were found to differ substantially, even between the two explicit water models, with the largest difference in contribution to the folding free energy arising from the water−water interactions (in particular from longrange electrostatics). The choice of water model also has a notable effect on the behavior of water in various ion channels (see below).
The majority of recent studies still focus on MD simulations using well-established water models. While these have provided invaluable insights on the structure and function of ion channels and nanopores, their suitability for modeling nanoconfined water has yet to be systematically evaluated. The development of classical water models has been predominantly motivated by the accurate modeling of bulk water. However, it remains unclear how transferrable a well-performing bulk water model is to the modeling of water behavior in novel nanoscale environments. Furthermore, more advanced water models and methodologies, such as polarizable water models and quantum mechanical treatments, potentially offer more detailed and accurate information, in situations where this is required, although these methods remain limited by their significant computational cost. As we progress toward exascale computing, 40,41 large-scale quantum mechanical calculations will become increasingly feasible, but as yet we fall short of being able to treat large protein/membrane systems in this way.

WATER MODELS IN MOLECULAR DYNAMICS
SIMULATIONS Within the scope of MD simulations, the treatment of water can be broadly categorized into implicit solvent, coarse-grained, and explicit atomistic models ( Figure 4). As of 2019, there are over 130 water models (a number of which are summarized in ref 42).
Some water models are developed to better describe a particular physical property, such as the dielectric constant, over a range of temperatures and pressures. Others are more transferable, aiming to provide a compromise between accuracy and applicability over a wide range of properties and conditions. In the case of biological simulations, the latter class of water model is of more use. This is because comparisons with experimental data are not straightforward and so a model designed to be accurately predictive in a wide range of environments and properties is needed. It is important to note that these water models are parametrized with respect to a given MD force field. Therefore, when conducting MD simulations with a given force field, more favorable agreement to experiment might be found with an older water model developed specifically with a given protein force field than with a more recent water model developed with a different force field. Explicit solvent tends to give a more accurate description . Schematic diagram of selected water model geometries. The left column shows the common geometries for rigid fixed charge models; the right is a schematic of the three key polarization schemes applied to a general 3-point geometry. It should be noted that there are various permutations to these, depending on the model. The labels for bond lengths, angles, and charges correspond to those reported in Table 2.  Chemical Reviews pubs.acs.org/CR Review of the system of interest but at significant computational cost compared with coarse-grained or implicit solvent models. An overview of different water models, discussing their uses and limitations, is provided by Onufriev and Izadi. 43

Implicit Solvent
Implicit solvent models provide a continuum approach for treating the presence of water and any solutes as a mean field. The free energy of solvation is typically decomposed into polar and nonpolar components which model the electrostatic and hydrophobic behavior of the solvent, respectively. 44 Formalisms are typically based on the Generalized Born equation 45 or the Poisson−Boltzmann model. 46 There are also extensions to these approaches that include missing elements of the physics, elevating implicit solvent approaches to being closer to explicit solvent methods, sometimes with relatively little increase in computational cost. 43 In the case of calculating the behavior of water in ion channels, the use of implicit solvent models yields at best limited insights into the molecular structure of water inside the channel. There have been some multiscale models of ions in model nanopores focusing on comparison of explicit with implicit water models, for example in the study by Valiskoé t al. 47

Coarse-Grained Water Models
Within coarse-grained simulations, one or several water molecules are grouped together into a single "bead". This increases the computational efficiency of the model but reduces its physical detail. Care must be taken when analyzing time scales, as the diffusivity of each particle is much increased compared with experimental or atomistic MD values. Furthermore, the structure and energetics of the system may not be able to be modeled to a sufficient degree of accuracy using the same coarse-grained method. 48 There are several coarse-grained water models available 49,50 including, for example, the model included within the widely used MARTINI coarse-grained force field, 51 which contains four water molecules per bead, as well as others containing three waters per bead. 52,53 Further adaptive coarse-grained water models have been developed which are able to model polarization effects by combining multiple water molecules into a network of multiple beads. These include the polarizable water model for MARTINI 54 and the big multipole water (BMW) model, 55 which both combine four water molecules into three beads, and the WAT FOUR water model which combines approximately 11 water molecules into a tetrahedral cluster of four beads. 56 The use of polarizable water models comes at the cost of computational efficiency due to the increased number of bead−bead interactions and degrees of freedom. For example, the WAT FOUR model was shown to be ∼50% slower than a water model based on the MARTINI force field. 57 Furthermore, the WAT FOUR water model cannot be used for channels with a small radius, as is the case for many biological systems.
The majority of simulations of channels using coarse-grained water models have focused on the protein or lipids. For example, the MARTINI representation of water was used to investigate the effect of varying the length and diameter of CNT porins on a surrounding POPC lipid bilayer and the tilt of the nanotube. 58 Coarse-grained simulations have also been used to find a permeation pathway along which to determine the free energies of ion conductance through the NanC porin. 59

Atomistic Water Models
These are by far the most widely used group of water models to investigate the effects of solvation in ion channels and nanopores. There are a wide range of types including rigid nonpolarizable models with 3−7 points, upon which charges or Lennard-Jones sites are placed to model the directional and nondirectional intermolecular interactions, respectively. In general, the Lennard-Jones site is placed on the oxygen, and the negative charge is either placed on the oxygen (for 3-point models) or is displaced away from the oxygen (for higher point models). Flexible and polarizable variants have also been developed: these are typically deemed to be more accurate but at increased computational cost. A summary of key nonpolarizable water models, often used in conjunction with ion channels and nanopores, is included in Table 2.
This shows parameters for common water models, along with selected physical properties (at 298 K and 1 atm, unless otherwise stated) of relevance to this review. The σ O and ε O values on the oxygen site are defined according to the Lennard- r OO is the separation between two oxygen sites. q H denotes the charge on the hydrogen site. To maintain charge neutrality, the charge on the oxygen site, q O , is equal to −2q H except in the case of 5-point models (TIP5P above) where q O = −q H . In the case of 4-point models, q O is displaced away from the oxygen site. θ denotes the angle between the hydrogen and oxygen sites. In the case of 5-point models there is an additional angle between the oxygen and displaced charge sites, denoted θ 1 in Figure 4. For TIP5P, θ 1 = 109.47°. l 1 and l 2 denote the oxygen−hydrogen and oxygen−displaced-charge separations, respectively. A graphical summary of these parameters is given in Figure 4. 3.3.1. Rigid Nonpolarizable Water Models. The two most common families of water model are the "simple point charge" models (SPC) 77 and "transferable intermolecular potential" models (TIP), which are based on the geometry of molecules within ice and water vapor, respectively. Within the TIP family, the 3-point model TIP3P and 4-point model TIP4P 78,79 are widely used. SPC was used for the parametrization of the GROMOS force field, 80 whereas TIP3P was used for AMBER 81 and TIP4P for OPLS. 82,83 A modified version of the TIP3P model, mTIP3P (sometimes referred to as TIPS3P), was used for the CHARMM force field and includes Lennard-Jones sites on the hydrogens. 84 Because of the close relationship between water model and force field, it is often the case that the most suitable water model to use for a given force field is the one used to originally parametrize that force field.
There have been various modifications of the water models mentioned above. The SPC model was modified to include a correction to the self-energy interaction of the dipole moments. This "extended simple point charge" model, SPC/ E, produced an improvement in the oxygen−oxygen radial distribution function, the density, and the diffusion constant at 300 K and 1 bar. 67 The TIP4P model was reparameterized to include long-range Lennard-Jones and Coulomb interactions using Ewald summation. The resulting TIP4P-Ew model produced a significant improvement in the bulk density of liquid water and enthalpies of vaporisation across 235.5−400 K at 1 atm. 73 The ability of a water model to accurately reproduce densities close to experimental values is taken as an indication of quality and transferability of the model because it is required for the accurate modeling of hydrophobic effects. 85 Vega and Abascal compared the temperature of maximum density (TMD) for a range of rigid water models. 68 The TMD for TIP4P-Ew was found to be 273 K. 73 By comparison, the TMDs for SPC, SPC/E, TIP4P, and TIP3P were at least 24 K below the experimental value of 277.13 K. 68,86 The 5-point water model TIP5P was specifically parametrized to accurately reproduce liquid densities across 236−373 K and 1−10000 atm, and was found to reproduce a TMD close to the experimental value. 70 (Vega and Abascal similarly obtained a TMD value for TIP5P of 285 K. 68 ) A further modification of TIP5P for use with Ewald summations, TIP5P-E, yields a TMD similar to TIP5P. 68,87 Most recently, the TIP5P/2018 88 and 7-point TIP7P 89 models have been developed which both favorably model the density as a function of temperature.
Several water models have been developed to better reproduce the dielectric constant under ambient conditions. These include the 3-point SPC/DC, 69 SPC/ε, 90 and H 2 O− DC 69 models, as well as the 4-point TIP4P/ε. 91 These were able to provide improvements on some properties but sometimes at the expense of others, providing evidence for the limitations of rigid nonpolarizable models. More recently, the dielectric constants for 19 nonpolarizable water models have been recalculated, showing that these values are particularly sensitive to simulation parameters: 75 this should be taken into account in the reading of Table 2.
In addition to the development of water models which correctly reproduce a particular property (for example the dielectric constant or TMD as discussed above), Abascal and Vega have developed a 4-point rigid water model, which aims to be as predictive as possible across a wide range of properties and conditions. 38 This new model, TIP4P/2005, was parametrized using the TMD, the enthalpy of vaporisation, the density of liquid water at ambient conditions, the densities of two ice polymorphs at particular temperatures and pressures, and the range of temperatures at which another ice polymorph (ice III) is thermodynamically stable at a pressure of 300 MPa.
There have been several studies testing the transferability and predictive quality of different water models. Five rigid nonpolarizable models (TIP3P, SPC/E, TIP4P, TIP4P/2005, and TIP5P) were analyzed using a scoring function which compared 17 simulated properties with experimental values. 64 TIP4P/2005 scored highest with a value of 7.2 out of 10. SPC/ E came second, scoring 5.1 out of 10. Additionally, the surface tension of seven rigid nonpolarizable water models has been evaluated. 92 It was found that SPC/E, TIP4P/Ew, and TIP4P/ 2005 produced values which were in good agreement with experiments. The SPC, TIP3P, and TIP4P models underestimated the surface tension, whereas TIP4P/Ice 93 overestimated it.
Izadi et al. employed a different method for developing new rigid nonpolarizable water models. 71 Rather than deciding on an existing molecular geometry and fixing point charges to this framework, they instead positioned point charges to optimize the dipole, quadrupole, and octupole moments. The rationale for this approach arose from the observation that to optimally mimic the quantum mechanical charge distribution of liquid water, the three classical point charges would counterintuitively have to be placed far away from the oxygen and hydrogen nuclei. The resulting "optimal point charge" (OPC) model has a shorter OH bond length and HOH bond angle compared with the experimental values used for the TIP family of water models (see Table 2). The authors employed the scoring function of Vega et al. 64 to evaluate the quality of OPC across a range of six physical properties. OPC was found to have a higher score (∼8.5) than previous water models. A 3point version of OPC, OPC3, has also been developed using the same optimal point charge method. 66 Using the same scoring system as for OPC, OPC3 obtained a score of 8.1 which was higher than other old and more recent rigid 3-point water models (TIP3P, SPC/E, H 2 O-DC, and TIP3P-FB). It is interesting to note that the OH bond lengths are longer than the experimental gas value; the bond angle was set to the tetrahedral value of 109.47°as a consequence of the fitting of the quadrupole moment.
Wang et al. have developed 3-point and 4-point water models based on the ForceBalance algorithm, which optimizes the model parameters using an adaptive combination of reference data and simulation results. 76 The resulting models TIP3P-FB and TIP4P-FB produce improved densities, dielectric constants, self-diffusion coefficients and shear viscosities over a range of temperatures compared with TIP3P and TIP4P. 76 Initial evaluations of their performance in simulations with amino acids 94 and with lipid bilayers 95 have been performed. The ForceBalance algorithm has been further used to develop water models based on data sets of surface tension. 74 The subsequent 4-point model, TIP4P-ST, performed extremely well over a range of properties. The 3point model, TIP3P-ST, while very accurately reproducing the TMD, leads to an overstructuring of water and a notably low self-diffusion coefficient.
Further reparameterization methods have been utilized. For example, O'Connor and English used the Design of Experiments (DoX) method 96 to systematically investigate how site charge, OH bond length, bond angle, sigma, and epsilon values in SPC affected the self-diffusion coefficient and radial distribution functions (RDFs) at 298 K. 97 It was found that the self-diffusion coefficient was primarily affected by the charge, epsilon, and OH bond length, and that the selfdiffusion coefficient error with respect to experiment could be markedly reduced relative to the original SPC model. As pointed out by Vega and Abascal, 64 any rigid nonpolarizable water model cannot properly describe the dielectric effects of liquid water and tends to underestimate dipole moments. This can be overcome by using polarizable water models.

Polarizable Water
Models. An outline of the theory of polarization and how this can be applied to MD simulations is given by Yu and van Gunsteren. 98 There are three key methods for introducing polarizability: induced dipole, classical Drude oscillator, and fluctuating charge.
1. Induced dipole: Inducible dipoles are placed at specified points in the molecule. The dipole moments are often determined self-consistently within a field of isotropic atomic polarizabilities. Examples of such models include POL3 (developed to better reproduce liquid properties) 99 and the BSV model. 100 A set of water models have been developed for the AMOEBA force field. These use a modified Thole-type interaction 101 in which the charge for one of the dipoles in each interacting pair is smeared over space instead of being treated as acting at a point. The earliest of these models is AMOE-BA03. 102 A more computationally efficient water model, Chemical Reviews pubs.acs.org/CR Review iAMOEBA (inexpensive AMOEBA), was subsequently developed. 103 This uses only one self-consistent cycle and has been parametrized with the ForceBalance algorithm. Although more computationally efficient than the original AMOEBA water model, there are possible limitations to the accuracy of iAMOEBA when modeling systems containing a high degree of heterogeneity and polarity (as discussed in the context of ion solvation 104 ). The AMOEBA03 model has been reparameterized using ForceBalance, yielding the AMOEBA14 model, which gives excellent agreement with experiment, particularly in terms of the TMD and dielectric constant. 105 A "coarse-grained" variant of the AMOEBA family, uAMOEBA, contains only polarizable sites on the oxygen centers, thereby producing a 3-to 5fold computational acceleration with respect to the AMOEBA03 model. 106 Most recently, the aniso-AMOEBA water model has been developed to incorporate the effect of anisotropic polarizability. 107 A further variant on the induced dipole water model, GPCM, has been developed, which models the charges as Gaussian functions. 108 2. Classical Drude oscillator: A mobile point charge is attached to a fixed point charge by means of a harmonic spring. In this way, the changing geometry of the charge ensemble in response to an electric field (E-field) mimics an induced dipole. Models include the SWM4-DP 109 and SWM4-NDP, 110 which incorporate positively and negatively charged Drude particles respectively, and the six-site SWM6 model which has been developed from the SWM4-NDP model. 111 For example, SWM4-DP has been used to explore the properties of the water− graphene interface, 112 while the SWM4-NDP model has been used in development of a Drude polarizable force field for the simulation of lipid bilayer membranes 113 and has been recently employed in simulations of model ion channels (see below, section 3.5). 3. Fluctuating charge: The magnitudes of the charges are allowed to change in response to the E-field. Examples of this class of model include the SPC-FQ and TIP4P-FQ models which are based on the rigid geometries of SPC and TIP4P respectively. 114 An extension of SPC-FQ and TIP4P-FQ models which incorporates an additional coupling between the Lennard-Jones parameters and the fluctuating charges on the oxygen sites has resulted in the SPC-pol and TIP4P-pol models. 115 Some models contain a combination of the above approaches. For example, the POL5 model incorporates both induced dipoles and fluctuating charges, the effects of which are coupled together on a 5-point rigid geometry. 116 Huang et al. have examined the behavior of a classical Drude oscillator model for proteins and water using the Drude-2013 polarizable force field and the SWM4-NDP water model. 104 It was shown that water dipole moments were affected by the surrounding protein environment. The dipole moment was greater than its bulk value of 2.46 D when close to negatively charged amino acid side chains and was smaller when close to positively charged, polar, or hydrophobic side chains. Significantly, the dielectric constant of the hydrophobic interior of the protein was higher with the polarizable Drude model. Polarizable models (TIP4P-FQ and SWM4-NDP) have been compared with standard nonpolarizable models (TIP3P, TIP4P, and SPC/E) of hydrophobically confined water. 117 The two polarizable models show a higher free energy barrier for water entering a dewetted region between two hydrophobic plates than for the nonpolarizable water models (also see discussion of hydrophobic gating below).

Flexible Water Models.
In situations where the effects of intramolecular vibrations need to be taken into account, a set of potentials, typically harmonic functions, can be applied to the OH bond lengths and angles to produce a flexible water molecule. Such situations include the prediction of vibrational spectra. 118,119 A consensus has yet to be reached regarding their effectiveness in nanopores and biological channels. For the flexible SPC/Fw model, the dielectric constant is improved if flexibility is taken into account. 117 Comparison between the surface tension of seven different flexible water models suggests that the surface tension is determined by a combination of bond-stretching and Lennard-Jones interactions. 120,121 Several flexible water models have been developed based on rigid model geometries. These include the SPC/Fw model based on SPC, 122 SPCE-F based on SPCE, 123 a flexible TIP3P, 78 TIP4P/2005f based on TIP4P/2005, 121 and a flexible SPC model (SPC-Fd), which includes an additional harmonic term for H−H interactions. 124 Other flexible models include the F3C model 125 which contains a rapid truncation of LJ and Coulomb interactions.
In the development of SPC/Fw, the influence of equilibrium bond length and angle on key properties was investigated. 122 It was found that the self-diffusion constant and static dielectric constant were primarily affected by the equilibrium OH bond length and bond angle, respectively. The self-diffusion constant was found to decrease with increasing bond length. It was hypothesized that this was due to the hydrogen becoming less shielded and therefore more able to form hydrogen bonds as the OH bond length was increased. Similarly, the static dielectric constant decreased with increasing bond angle. While the latter could be hypothesized to be caused by the decrease in molecular dipole moment with increasing bond angle, an investigation into the correlation between dipole moment and dielectric constant showed that the dominant effect was directly from the changing bond angle, not the dipole moment. 122 This is consistent with the findings of Hochtl et al. 126 In the development of the TIP4P/2005f flexible model the fluctuations in bond angle are represented by a harmonic term, while bond length vibrations are represented by an anharmonic Morse potential. 121 This model was shown to produce improved density−temperature characteristics and a melting point closer to experimental values than does the rigid TIP4P/ 2005 model. The dielectric constant is slightly worse compared with TIP4P/2005, in contrast to the SPC/Fw water model, which produced an improvement relative to the rigid SPC and SPC/E models. 117,122 It should be noted that for flexible models the time step for MD simulations needs to be reduced by a factor of ∼5, which significantly increases the computational cost. 121 3.3.4. Advanced Water Models. While flexibility can partially model the effects of molecular polarizability, it cannot model atomic polarizabilities and their response to the electrostatic environment. A class of models incorporating both flexibility and polarizability have, therefore, been developed. For example, TTM2.1-F, 127 TTM3-F, 128 and TTM4-F 129 are Thole-type 101 flexible and polarizable models. Chemical Reviews pubs.acs.org/CR Review They have been modified from the earlier flexible TTM2-F 130 and rigid TTM2-R 131 models for use in liquid simulations. These TTM models are 4-point models with smeared charges on the hydrogen and M sites and smeared induced dipoles on the hydrogen and oxygen sites. Machine learning techniques have also recently been used to develop a flexible, polarizable, and multipolar electrostatic water model. 132 To account for the quantum mechanical nature of bond flexibility and polarization, as well as for nuclear quantum effects, QM methods have been employed in the development of a further set of flexible water models. These include the q-SPC/Fw model, which was parametrized from SPC/Fw using normal-mode path integral MD and centroid MD simulations. 34 A further flexible quantum model, q-TIP4P/F, has been developed using the rigid TIP4P/2005 as a basis. 33 It was found that quantum effects on the dynamics of liquid water had been previously overestimated. The overall quantum effect is the sum of two mechanisms, which partially cancel each out: (i) the intermolecular zero point fluctuations and tunnelling, which weaken the hydrogen-bond network and result in an increased diffusion coefficient, and (ii) the intramolecular zero point fluctuations, which increases the molecular dipole moment, thereby strengthening intermolecular interactions and decreasing the diffusion coefficient.
More recently, advanced property analysis (for example energy decomposition analysis), automated parametrization methods, such as ForceBalance, and machine learning methods have been used to develop next-generation force fields. Many of these models, while yielding very accurate results, are computationally expensive. For example, MB-pol is regarded as one of the most accurate models, but it is ∼50 times more computationally expensive than AMOEBA, and has not yet been used for simulations beyond water 133−136 and waterhalide (using MB-nrg which is based on MB-pol 137 ) systems. For further information about advanced water models, including parametrization and machine learning methods, the reader is directed to the recent review papers by Cisneros et al. and Demerdash et al. 138,139 Finally, a particularly novel "Quantum Drude Oscillator" (QDO) model for water has been developed. 26,140,141 The electrons are treated as a set of quantum Drude oscillators by including a harmonic oscillator operator into the Hamiltonian operator. In this way, the electrons are effectively treated as "coarse-grained". While this model is able to model long-range effects, such as polarization, dispersion, and mixed interactions to orders higher than the dipole limit and is a potentially promising model for water simulations, it has yet to be developed for use with biological systems.
In conclusion, there are many water models available for MD simulations. The empirical pairwise-additive fixed-charge models, such as those in the SPC and TIP families, are the most commonly used due to being relatively computational inexpensive and transferable. These models have been reparametrized multiple times, yielding significant improvements in the reproducibility of experimental properties. Polarizable and flexible water models are also popular for liquid water simulations, although widespread use of these and other more advanced water models on simulations of nanopores and biological channels has yet to take place.

Quantum Mechanical Descriptions of Water
As discussed in section 2.4, it is possible to take a QM/MM approach to modeling water and biological regions of interest.
For example, Shaw et al. conducted a systematic study of the interaction between a single water molecule, which was treated quantum mechanically, with a surrounding region of classical, that is, MM waters. 142 Three QM theories were investigated: density functional theory (DFT) using the BLYP exchangecorrelation functional, ab initio Hartree−Fock, and secondorder Møller−Plesset perturbation theory. Four different MM water models were considered: TIP3P, TIPS3P, TIP4P, and TIP5P. By calculating the structural and energetic perturbations to the system, TIP4P was found to produce the smallest perturbation and was, therefore, the most compatible water model. By contrast, TIP5P was found to be the least suitable water model with the lone pairs becoming unphysically situated within in the electronic distribution of the quantum mechanical water molecule.
Rather than parametrizing classical water models to include quantum effects and/or using hybrid quantum/classical approaches, an ostensibly more accurate method is to use a full quantum mechanical approach. This does come at a significant computational cost. There are several studies on water using density functional theory (DFT) as reviewed by Gillan et al. 143 Several studies report that ab initio MD (AIMD) simulations of water, in particular with the PBE 144 and BLYP 145,146 exchange correlation functionals, lead to an overstructuring of water as characterized by calculating the pair correlation functions between oxygen atoms. 147−151 There has been speculation that this is due to incomplete convergence of the basis set, although Rempe et al. have demonstrated that this is not so when using the PBE functional. 152 There has also been some discussion about flexible versus rigid water molecules in AIMD simulations. Allesch et al., using the PBE functional, found that fixing the OH bond length and HOH angle produced a structure that was not as overstructured as when using a flexible molecular representation, and that the resulting structural and dynamical properties were much closer to experiment. 153 Leung et al. followed up on this by conducting a study varying the bond length and angles of the water molecule, again using the PBE functional. 154 They found that increasing the bond length by as little as 2−3% had a significant effect on the structural and dynamical properties.
Several AIMD simulations have also been conducted of water in confined environments. 155−157 Mann et al. demonstrated the stabilization of a water wire inside a (6, 6) singlewalled CNT, and that this water wire had an affinity for accommodating an excess proton across the length of the nanotube. 158 Meanwhile, Dellago et al. used both AIMD and an empirical bond valence (EVB) model to investigate how one-dimensional water wires inside CNTs enable proton conduction via a Grotthuss mechanism. 159 It was found that proton mobility was enhanced by 40 times that of bulk water along perfect water wires but that it was reduced in the presence of defects.
In the context of ion channels and other biological nanopores, the channel geometry and electronic structure of gramicidin A has been investigated via large-scale DFT calculations of the entire channel in a vacuum. 160 A largescale DFT study of a complex ion channel in a lipid bilayer with surrounding water has yet to appear. The future of largescale DFT on biological systems clearly points in this direction but current computational resources still delay their feasibility. Chemical Reviews pubs.acs.org/CR Review

Evaluation of Water Models for Simulations of Channel and Nanopore Systems
There have been several investigations into the suitability of different water models for water in confined systems, including channels and nanopores. Prasad et al. compared the effect of water model on water transport through an OH-functionalized graphene nanopore. 161 They found that the water flux between the models evaluated (SPC, SPC/E, SPC/Fw, TIP3P, TIP4P, and TIP4P/2005) varied by up to 84%. They hypothesized that this was due to the differences in partial charges influencing the hydrogen bond network and concluded that overall TIP4P/2005 was the most suitable water model for studying transport in such systems. Guerrero-Avileś et al. conducted DFT calculations and ab initio MD simulations to investigate the transport of water and sodium and chloride ions through "nanoslits" in graphene sheets. 162 Their simulations showed that water molecules passed through the ∼0.8 nm diameter pores at a rate of ∼0.3 ps −1 . This is in favorable agreement with experiments which produced a transport rate of 3 waters per picosecond through pores of ∼1 nm in diameter. Bauer et al. have compared TIP3P, TIP4P, and SPC/ E against the polarizable TIP4P-FQ and SWM4-NDP models for a system of water confined between two hydrophobic plates. 117 It was found that TIP4P-FQ behaved in a markedly different manner to the other water models, with nonbulk-like dynamics being observed at much larger plate separations than for the other models. In the context of hydrophobic gating, it is potentially important to accurately model the effects of induced polarization present at a vapor−liquid interface. The application of polarizable water models (using AMOEBA14) to simulations of ion channels containing a hydrophobic gate has been shown to have small but measurable effects on the wetting/dewetting behavior (ref 163 and also Lynch and Sansom ms. in preparation). Taken together, these studies argue for more widespread evaluation of different water models in nanopore and channel simulations.

SIMPLIFIED MODELS OF NANOPORES
In the early days of simulations of nanopores, the absence of any high-resolution structures dictated the use of simplified, rather abstract models. These enabled a number of initial explorations (see, e.g., the studies by Breed et al. and Sansom et al. 164−166 ), focusing on the diffusional and dielectric properties of water in nanoscale channel-like environments.
Simulations of simplified models have continued; these have proved especially valuable for revealing general principles underlying the variation in water properties in nanoconfined environments compared with bulk water. More general studies of water at interfaces have been treated quite extensively (as reviewed by McCoustra 167 ); simulation studies of water confined between two surfaces (a "nanoslit"; Figure 5A) or within a simple model nanopore ( Figure 5B) will now be reviewed.

Physical Chemistry of Water in Nanopores
There have been a number of studies of the thermodynamics of water in nanopores and, more generally, in hydrophobic nanocavities within proteins and related environments, as reviewed in ref 169. A number of these studies have focused on the thermodynamics of hydration (i.e., of wetting vs. dewetting) of nanopores, which, in the context of gating of nanopores, is discussed more extensively below. Two key studies, of water in a hydrophobic nanopore 168 and in a nonpolar nanocavity, 170 used simulations of nanoconfined waters to estimate free energy landscapes of wetting/dewetting. As shown in Figure 5C, these landscapes display one or two minima: one for the empty pore (n/n BULK = 0) and one in the vicinity of the bulk density (n/n BULK = 1). As the radius of the nanopore is increased, the wetted state becomes globally stable relative to the dewetted state. Similar behavior has been seen in joint experimental and computational studies of the solvation of hydrophobic cavities in the enzyme T4 lysozyme L99A, 171 for which wetting can be favored by increasing pressure, and in the tetrabrachion protein, 172 for which wetting is disfavored at elevated temperatures. Comparable free energy curves are observed for water between two nanoscopic hydrophobic surfaces. For example, Remsing et al. showed that water density fluctuations played a key role in determining the pathways to dewetting for nanoconfined water. 173 Experimentally, the existence of a dehydrated hydrophobic cavity within a protein (of a volume ∼0.3 nm 3 , that is, sufficient to accommodate up to ∼10 waters) has been demonstrated for bovine β-lactoglobulin using water 17 O and 2 H magnetic relaxation dispersion, 174 alongside MD simulations and free energy calculations. More detailed thermodynamic analysis of the hydration of carbon nanotubes 175 has suggested that wetting of hydrophobic CNTs is due mainly to increased entropy of nanoconfined water. This increase in entropy is rotational and translational for small CNTs but mainly translational for larger nanotubes.
In parallel with these studies of the thermodynamics of nanoconfined water, experimental studies based on ultrafast infrared pump−probe spectroscopy of water within reverse micelles 176 have focused on hydrogen bond and orientational Figure 5. Simple idealized models used in simulations of nanoconfined water: (A) water molecules (red) in a nanoslit (the confining surfaces are in blue) and (B) water molecules in a nanopore (shown in cross section). (C) Free energy density for water within a hydrophobic nanopore model (see below) as a function of n/n BULK , that is, the water density in the pore normalized to bulk density. The two minima correspond to the observed two-state (dry vs wet) behavior. The dry (i.e., vapor) state becomes metastable with increasing radius and for R > 0.55 nm the liquid (i.e., wet) state is globally stable. Curves are sketched from data reported in ref 168. Chemical Reviews pubs.acs.org/CR Review dynamics. In small reverse micelles (containing ∼50 water molecules), a single water ensemble is seen, which is dominated by waters interacting with the (hydrophilic) lining of the nanocavity within the reverse micelle. More generally, dynamic interfacial water density fluctuations also occur near convex-shaped nonpolar solutes. 177

Nanoconfined Water
Simplified (i.e., idealized) models of nanotubes and nanopores have been investigated extensively in terms of the behavior of nanoconfined water within these systems. Recent studies have focused largely on two aspects: (i) water flow/transport as a function of the pore dimensions and the chemical properties of the pore-lining surface and (ii) the dielectric properties of water confined within model nanopores.
Simulations of TIP3P water in nanopores crossing a simple membrane model formed by a slab of immobile waters in a rigid ice-like structure 178 focused on the difference between hydrophilic (ice-like water-lining) and hydrophobic (i.e., partial charges on the immobile waters removed) nanopores. The water flux under applied pressure was evaluated as a function of pore radius and geometry (cylindrical vs conical), and the results compared with the predictions of continuum fluid dynamic (CFD) theory via the Hagen−Poiseuille equation. 179 For small (radius 0.6 nm) hydrophobic pores, water flow was greater than predicted by Hagen−Poiseuille, whereas lower water flux was seen for the hydrophilic pores, that is, less than predicted by the Hagen−Poiseuille equation. For the small hydrophilic pores, the structural and dynamical properties of water are strongly influenced by interactions with the pore-lining interfaces. Even in larger (radius 0.9 nm) pores the parabolic distribution of water velocities underlying the Hagen−Poiseuille equation was not observed. Comparable simulations 180 of pressure-driven electrolyte flow through smooth, dipole-lined (i.e., hydrophilic) nanopores 181 of radius 0.6 nm penetrating a model membrane yielded water fluxes which were comparable to those in CNT (see below) simulations. This suggests that reduced water flux in hydrophilic model pores may be related both to hydrogen bonding interactions and to the molecular-scale roughness of such pores. Other studies 182 have employed a highly simplified kinetic generalization of a one-dimensional dipole model of water in a single-file nanopore to study the dielectric response of nanopore water to time-dependent E-fields in the direction of the pore axis (see below). This has helped relate studies of the microscopic properties of single-file water to, for example, data from dielectric spectroscopy experiments. Using a simple nanopore model (the particles of which were able to undergo thermal fluctuations) of radius ∼1 nm and containing a Lennard-Jones fluid 183 in MD simulations under flow conditions, the velocity profile of the fluid particles was approximately parabolic as predicted by CFD theory. These simulations were combined with analytical approaches to provide a framework for exploring the relevance of fluctuation effects, including enhanced diffusion. Simulations of pressuredriven water flow through a silica nanochannel of radius ∼6 nm 184 somewhat surprisingly suggested flow enhancement (relative to the Hagen−Poiseuille flow) similar to that of a CNT (see below). This suggests that the detailed chemical properties of a hydrophilic pore may be important in determining water behavior.
Effects of pore geometry on the pressure-induced flow of SPC/E water flow have been explored 185 using hourglass-shaped (biconical) hydrophobic pores (minimum radius 0.15 nm) in MD simulations. The cone angle α and pore length were varied, revealing that a conical entrance reduces the hydrodynamic entrance hindrance and so leads to higher flow rates for water. Further MD simulations of water transport in these hourglass-shaped model pores 186 explored the effect of the pore length on water energetics and transport, suggesting that water flux increased with pore length, reflecting the decreased importance of entrance effects as the pore length increases. Simulations of a simple hourglass shaped model (designed to mimic an AQP-like pore, see below) 187 indicated that the water flow (albeit in rather short simulations) depended on the hydrophobicity of the pore walls, expressed as a contact angle, with a peak in water flux at θ ∼ 100°, that is, a hydrophobic pore. In addition to pressure-driven flow, MD simulations have been used to explore, for example, the effects of the surface roughness of the pore lining on electro-osmotic flow through nanochannels, 188,189 indicating that surface topography (i.e., local pore shape) plays a significant role in electro-osmotic water flow. Simulations of electro-osmotic flow through idealized (i.e., smooth) nanopores show that both electro-osmotic and Poiseuille (i.e., pressure driven) flow agree with continuum transport descriptions in terms of, for example, comparable velocity profiles, apart from for the first layer of fluid close to the nanopore wall. 190,191 In addition to water flow effects, the behavior of water in nanopores can also differ from the bulk in terms of water dielectric behavior. 166 For example, studies combining MD simulations of SPC/E water with a linear response formalism to obtain a spatially resolved dielectric tensor considered water in generic hydrophobic cylindrical nanoconfinement, including cylindrical nanopores of radius ∼0.5−2.5 nm. 192 Interestingly, the radial component of the dielectric constant showed oscillations near the center of the pore; this would be expected to influence the solvation of ions within such pores.
The studies discussed above have, by and large, used standard fixed point-charge water models (e.g., TIP3P or SPC/ E; see section 3.3.1 above for a discussion of water models) to explore water behavior in nanopores. A number of other studies have either used simplified coarse-grained models or more complex representations of water in nanopores. Coarsegrained simulations 193 using the ELBA water model 194,195 have been used to explore capillary evaporation in simple nanopores (radius ∼2 nm) with either a moderately or a very hydrophilic pore-lining surface. This coarse-grained approach was extended to simple nanopores with different patterns of hydrophobic and hydrophilic walls, demonstrating that pore filling with water depends on the pattern of hydrophilicity/hydrophobicity. 196 As we will see below (section 5.4), this is relevant to possible gating mechanisms for biomimetic pores.
An analytical model within a CFD framework 197 has been used to explore water flow enhancement for nanopores with varying geometries and pore-lining hydrophobicity/hydrophilicity. Internal resistance to flow was reduced by increasing the contact angle (i.e., the hydrophobicity of the nanopore walls). Overall, water flow could be optimized by varying the contact angle and geometrical parameters including those for the conical pore entrance/exit. Earlier studies of idealized model nanopores 198 used simplified finite-element calculations to demonstrate that, compared with a plain cylindrical nanotube, a biconical channel of optimal angle (mimicking AQPssee section 7.1 below) provided increased hydrodynamic permeability. Therefore, the hourglass geometry, Chemical Reviews pubs.acs.org/CR Review associated with small surface friction, optimizes the permeability. Comparison of these results with MD simulations of water flow in a biconical nanopore model based on "fused" CNTs (see section 5.1) 199 yielded good agreement between MD simulations and continuum hydrodynamics predictions of entrance resistance. More recently, a hybrid continuummolecular mechanics model has been used to correct the Hagen−Poiseuille model for water-surface interactions and surface wettability, 200 deriving viscosity functions, which determine factors influencing flow enhancement/inhibition of confined water.
Switching to more complex water models, Zaragoza et al. have explored TIP4P/2005 water in nanoslits and nanopores (graphene and CNTs). 201 As with simpler fixed point-charge models, enhanced diffusion of water is seen in hydrophobic nanotubes. These results and those discussed above agree with a CFD model 202 in which the apparent viscosity of water decreases with an increasing contact angle (i.e., hydrophobicity), and tends to bulk viscosity with an increase of pore dimension. These effects have also been discussed in detail in the context of CNTs and related pores by Wu et al. 203 (see section 5.1) and in a recent review. 204 Moving beyond classical water models, ab initio MD simulations of water confined in model nanoslits between graphene sheets 205 show that confining water in narrow slits (either ∼1 nm to give a water bilayer or ∼0.7 nm to give a water monolayer) changes both the hydrogen bonding and electronic structure of water.
So, in summary, what have we learned from simplified models? A focus of many studies has been on water flow in nanopores. This is faster than Hagen−Poiseuille-predicted flow if the pore-lining is hydrophobic, while it is generally slower if the pore lining is hydrophilic ( Figure 6). Thus, the strength with which water interacts with the pore-lining surface is crucial in determining its local behavior within nanopores. Pore-lining surface roughness and fluctuations as well as the overall pore geometry shape are also important. Therefore, to enhance nanopore design, characterization of the behavior of real "chemical" and "biological" nanopores may be required toward derivation of more general principles. We will now examine what simulations have taught us about the behavior of water in increasingly complex real pores.

Wetting and Dewetting: Hydrophobic Gates in Nanopores
Capillary evaporation or dewetting has been observed in a number of studies of simple model nanopores, and may be related to transient dewetting observed in CNTs (see the study by Hummer et al. 206 and below). More generally, local dewetting has been discussed in terms of hydrophobicity and the folding and ligand-binding by proteins. 207 Here, we focus on simple models of dewetting in nanopores, as this may be related to possible gating mechanisms for both biological and biomimetic nanopores. 208−210 Crucially, simulations have been used to determine that such wetting/dewetting behavior may be modulated by the imposition of an external E-field, such as that which may be formed across a lipid bilayer membrane via an asymmetric distribution of ions.
Early studies of a simple model of a hydrophobic nanopore in a model membrane demonstrated hydrophobic gating, dependent on the radius of a pore of fixed transmembrane length. It was also shown that a closed hydrophobic nanopore could be opened by adding dipoles to its lining. 211 A more extensive study of liquid−vapor oscillations of water in hydrophobic nanopores 168 showed that the free-energy difference between the two states (liquid and vapor, i.e., wet vs dry) depended on the pore radius such that for hydrophobic nanopores with radii above ∼0.55 nm, a liquid-filled pore became the globally stable state (see Figure 7; also see Figure  5C and section 4.1 above). One-dimensional confinement was shown to affect the dynamic behavior of the water molecules, increasing self-diffusion compared with bulk water such that permeabilities for these simple model nanopores were of the same order of magnitude as for biological water pores. Subsequent studies focused on hydrophobic gating from an ionic perspective. 212,213 Both for simple models of hydrophobic nanopores 212 and for (6,6) vs (10,10) CNTs, 214 there was a threshold radius of ∼0.5 nm below which ion permeation was not observed. At the same time comparable studies by Allen et al. showed that water (SPC/E) molecules penetrated a hydrophobic channel of fixed length only beyond a minimum radius and that near the threshold radius water permeation was intermittent and sensitive to, for example, the polarizability of Chemical Reviews pubs.acs.org/CR Review the medium forming the channel. 215 More recent studies have described how hydrophobic gating transitions from closed to open can be explained by a general thermodynamic analysis of nanoconfined water. 216 Building on initial observations of hydrophobic gating in a simplified model pore system, 215 subsequent studies 218,219 demonstrated how an E-field (generated by an ion concentration imbalance resulting in a field of ∼4 V/nm inside the pore) could wet a previously dry pore (radius ∼ 0.3 nm). The ionic charge imbalance across the membrane was shown to induce water permeation of the hydrophobic pore, thus making it permeable to ions. The movement of ions subsequently dissipated the charge imbalance and the pore once again dewetted. This is directly relevant to a number of studies of E-field-induced wetting of simple models of nanopores. Wetting of the narrow hydrophobic pore within a (6,6) CNT has been shown to be favored by an axial E-field of more than ∼0.3 V/nm. The free energy landscape of wetting/ dewetting (see section 4.1 above and Figure 5C) was shown to be tilted in favor of wetting when the E-field was imposed. 220 In a pioneering study, Wallqvist and Berne 221 observed dewetting when two hydrophobic plates were brought closer than 0.95 nm together. Subsequently, in an impressive series of studies of dewetting/wetting of nanopores, Luzar and colleagues have used a simple model consisting of two parallel planar hydrophobic surfaces with water molecules in between. Early studies 222,223 explored capillary evaporation, using a lattice gas model of liquid (water) confined between two planar surfaces. Changing the nature of the surfaces by combining hydrophilic and hydrophobic patches, it was shown that evaporation kinetics depends on this nanoscale structure, allowing estimation of the free energy barrier of vapor tube formation. A similar system was used in atomistic simulations (using both Monte Carlo and MD simulations) 224 to explore capillary evaporation of (SPC) water confined in a smooth hydrocarbon-like slit. At a hydrophobic surface separation of 1.4 nm (i.e., ∼4 waters thick) the barrier to capillary evaporation was ∼19 kT. This corresponds to an evaporation rate of ∼10 5 nm 2 s −1 , indicating that water readily evaporated.
These investigations were extended to studies of "electrowetting" as applied to a simple (nanoslit) model of a nanopore. This has proved to be directly relevant to subsequent experimental studies of (nonbiological) nanopores 225 (see below for a more extended discussion) including hydrophobic nanopores in both silicon nitride 226 and in polyethylene terephthalate. 227 Initially, Daub et al. performed molecular simulations of nanosized water (SPC/E) droplets on a model hydrophobic (graphite) surface in the presence of an imposed E-field of 0.3 V/nm. 228 The contact angle of the water nanodroplet on the hydrophobic surface (normally 96°in the absence of the E-field) was sensitive to the applied E-field direction relative to the liquid/solid interface. A parallel E-field had the largest effect, reducing the contact angle by ∼20°at the leading edge of the nanodrop. Similar studies showed a decrease in contact angle as the E-field was increased up to more than 1 V/nm. 229 Subsequently, simulations were used to explore SPC/E water between parallel hydrophobic plates with a separation, d, of ∼0.9 to ∼4 nm and with applied E-fields up to 4 V/nm. 230 Electrostriction was observed, i.e. an increase in water density, especially when the E-field was parallel to the planes, but no freezing of water was observed. Importantly, for smaller plate separations (d < 1.3 nm) capillary evaporation (which was observed in the absence of an E-field) was suppressed upon application of an E-field. Overall, these studies showed that an E-field induced a transition from strongly hydrophobic to strongly hydrophilic behavior in a system containing water confined in a hydrophobic nanoslit and in equilibrium with zero-field bulk water. Additional simulations 231 showed that water molecules in nanoconfinement could maintain their interactions while polarized by the imposed field. Thus, if a field-exposed nanosized confined system is equilibrated with field-free bulk water, water molecules are attracted to the confined region. Simulations of two molecular hydrophobic plates formed by butyl-functionalized graphene 232 with a separation between the terminal methyl groups of ∼1.1 or ∼1.9 nm, exposed to SPC/E water and ∼1 M NaCl, showed that in an E-field perpendicular to the plane (equivalent to a total voltage across the slab of the order of 0.1 V) the model nanopore wetted within <1 ns. Switching the field off caused a rapid drop in water density (i.e., the reversal of electrostriction) followed by dewetting within a few nanoseconds. Therefore, an applied E-field was shown to reversibly control the wetting and dewetting of a simple hydrophobic nanoslit. It was also shown that imposing an external E-field could switch the state between a narrow channel from salt depletion to salt excess. 233 Electrofreezing of confined water has also been reported, for water confined between two hydrophilic surfaces (separation ∼0.9 nm) with an E-field of 5 V/nm applied parallel to the confinement plane. In related studies, the effect of an E-field on TIP4P/2005 water in a 4 nm channel between two graphene sheets was examined. 234 A field of up to 0.37 V/nm applied perpendicular to the sheets was imposed via placing opposing surface charges on the graphene sheets, and an external force was used to drive water flow. The water velocity profile became more parabolic as the E-field was increased and the apparent viscosity of water increased, reducing the flow rate. Therefore, the walls of the nanochannel became effectively more hydrophilic in the presence of the E-field.
The sensitivity of such effects of imposed E-fields on nanopores to the water model employed has not been widely investigated. The effect of E-fields up to 2 V/nm on "bulk" water (i.e., a box of 512 molecules) was investigated through comparison of a fixed point-charge (SPC/E) and two polarizable (SW4-NDP and BK3) water models. 235 E-field effects were found to be stronger in the polarizable water models, with slower switching of hydrogen bond partners and reduced translational mobility compared to with the nonpolarizable water model. Extending these studies to alternating E-fields 236 showed a nonmonotonic frequency dependence of changes in water hydrogen bonding with peak dynamics at ∼200 GHz, where the period of the E-field oscillation is comparable to the time it takes a water proton to switch between hydrogen bond acceptors. A polarizable model (BK3-AH) of water within either a smooth wall or graphene nanoslit (see above) 237 was explored in response to a perpendicular external field of up to 3 V/nm. Regardless of the E-field, the mean dipoles of interfacial water molecules were lower (by ∼10%) than in the bulk phase with the polarizable model. Simulations of the nanopore-bulk phase equilibrium showed that the polarizable water model was more strongly absorbed to the nanopore in the absence of the applied E-field, while the nonpolarizable water showed greater uptake into the nanopore in response to an E-field.
The importance of such considerations of E-field effects is discussed in a recent review of wettability of nanopores. 10 In particular, wetting/dewetting via changes in externally applied E-fields is an attractive prospect for direct control of gating in novel nanopores. An experimental study 238 has demonstrated the influence of an applied electric potential on water affinity for the interior walls of a CNT channel 14 nm in diameter. Simulation studies of the effect of an external E-field on the wetting of a hydrophobic gate in a simplified model of a membrane protein 239 (discussed in more detail below) demonstrate that these effects of E-field on wetting can also be seen in a protein nanopore environment. However, the studies discussed above on relatively simple systems show a perhaps unexpected complexity of behavior. This in turn suggests that quantitative, predictive design of E-field control of gating for more complex nanopores (whether based on naturally occurring proteins or on biomimetic materials) remains a major challenge and will need careful application to a number of simulation approaches.

NONBIOLOGICAL NANOPORES
Nonbiological nanopores, especially carbon nanotubes (CNTs), are an active area of research on the behavior of nanoconfined water in pores (Figure 8). A recent review, 204 mainly focused on CNTs, compares MD and continuum approaches to water flow in such nanopores. In this section we will first discuss more recent simulation studies of water in CNTs, before examining a range of other nonbiological nanopores.

CNTs
Following on from early suggestions of their wettability 240 and from early MD simulations, 206,241 there have been many computational studies of water in CNTs. 242 Here we will focus on more recent studies, which have gained additional relevance following the experimental demonstration that short CNTs can insert into lipid bilayers to form "CNT porins". These porins can transport water and ions across membranes, 243 thereby providing a biomimetic platform for synthetic nanopores. Further experimental studies have built upon this landmark work: 244 enhanced water permeability has been demonstrated in narrow (of radius 0.4 nm) CNTs which support single-file water, 245 and more recently CNT porins have been incorporated into synthetic membranes formed by amphiphilic block copolymers. The demonstration that the electronic properties of single walled CNTs are sensitive to internal wetting, provides an additional motivation for understanding the behavior of water within CNT nanopores. 246 The earliest simulation studies of water in CNTs 206 showed a pulse-like motion of water through the nanotube, and indicated that small changes in channel hydrophobicity could lead to wet−dry transitions. An extensive review of this and related work 169 concluded that the combination of strong water−water interactions with weak attraction of water to the confining wall resulted in rapid flow of water. Subsequent theoretical investigations 247 focused on single-file water, both in CNTs 248 and in nanopores in general. 249 Simulations of single-file TIP3P water in (6,6) CNTs (of radius ∼0.4 nm) suggested that their filling with water was (just) favorable and was driven by a decrease in energy alongside a loss in water entropy. Recent first-principles MD studies of narrow CNTs have confirmed the one-dimensional single-file hydrogen bonding pattern seen in earlier studies. 245 A number of studies have addressed enhanced water flow in CNTs, as is also seen in studies of idealized hydrophobic   252 showed that water flow was enhanced (relative to the Hagen− Poiseuille equation) and exhibited a nonparabolic radial distribution of velocities. For pore radii from ∼0.4 to 0.6 nm (corresponding to (8,8) to (12,12) CNTs) enhancement factors of ∼10-fold were seen (i.e., somewhat lower than those suggested by experimental studies). 251 A recent comparative study of pressure-induced flow through CNTs with radii between ∼0.4 and ∼0.9 nm predicted by different water models (SPC, SPC/E, TIP3P, TIP4P, and TIP4P/2005) suggested a nearly 3-fold difference between, for example, TIP3P and TIP4P/2005. Pressure-driven flow of water and ions through a (9,9) CNT has been compared for TIP3P and TIP4P/2005 water. 253 TIP3P water flow rates agreed with CFD calculations. In contrast, for TIP4P/2005 water, flow was slowed by the water structure within the CNT. Furthermore, a dependence of the degree of enhancement on the thermostat used in the simulations has been suggested. 254 Thus, while we should not overinterpret quantitative differences, the overall picture remains of significant water flow enhancement within CNTs, as seen by a combined analysis of published MD data and continuum theory to address water flow in CNTs. 255 A decrease in the CNT radius was shown to be accompanied by a decrease in the water−CNT interfacial viscosity and an increase in the water core viscosity. Simulated water flow in CNTs with similar radii but different chirality (obtained by altering the chiral indices (n,m) of the nanotubes) 256 showed that water molecules experienced the lowest resistance to flow in armchair (i.e., n = m) nanotubes. The force exerted by the CNT surface on the water molecules increased as the CNT type was changed from armchair to chiral and to zigzag. Therefore, the details of pore-lining surface structure can influence water flow even in "smooth" CNT nanopores. Given the rapid flow of water through CNTs, it is of interest to explore how this might be regulated or "gated". For example, it has been suggested that water flow might be switched off via lateral deformation of a (6,6) CNT. 257 A deformation of ≥0.25 nm switched off water flux and resulted in local dewetting of the CNT pore. Water flow can also be slowed by mechanically "squashing" a CNT. 258 Other studies have explored the effects of a point charge just outside a (6,6) CNT to increase the flow rate of water 259 and how rotation of a chiral CNT 260 can drive water flux through the pore. Clearly there is considerable potential for biomimetic modification of CNTs to control their ability to enable water flow.
In the context of both water and ionic flow, there have been a range of studies on the effects of an electrostatic field in water in CNTs. The overall picture is complex, with water behavior dependent on the magnitude and direction of the field and on the CNT dimensions. Broadly, lower-strength (≤1 V/nm) fields along the pore axis can enhance water permeation whereas greater (≥2 V/nm) fields can lead to ice-nanotube formation and loss of permeation. A study of a (6,6) CNT with single-file TIP3P water (but no membrane) 220 showed that a field of 0.33 or 1.0 V/nm favored nanotube wetting. Simulations of various CNTs in a phospholipid bilayer (using CHARMM27 and TIP3P waters) with and without an E-field along the nanopore (z) direction of 65 mV/nm 261 revealed that applying an E-field accelerated permeation approximately 2-fold in a (5,5) CNT, but slowed water selfdiffusion in larger CNTs. For narrow CNTs with single-file TIP3P waters, 262 water permeation was substantially reduced for high E-field intensities of 2−3 V/nm perpendicular to the pore axis. Pressure-driven flow of TIP3P water along a (6,6) CNT 263 was shown to be modulated by a 1 V/nm field. By varying the angle of the field relative to the pore, it was shown that a field along z (the nanopore axis) gave a 4-fold increase in water flux, while a field along x gave approximately zero flux. The effect of E-fields on pressure-driven TIP4P/2005 water flux through a wider (15,14) CNT has also been explored. 264 With an axial field of up to 1 V/nm, water flux only increases if the field is applied in the pore mouth region, and not if applied along the whole length of the pore or beyond. This suggests that an applied E-field may reduce the flow of molecules within the CNT but increase flow rates by aligning molecules close to the entry of the CNT. High fields (greater than 1.5 V/nm) were able to change the properties of water from liquid to solid. Fluctuating E-fields may also have an effect: positioning a vibrational charge outside a CNT pore can increase water flux due to changes in the hydrogen bonds. 265 THz E-fields and enhanced water flow have also been explored. 263 Looking beyond single CNTs, Kalra et al. employed MD simulations of a two-compartment system in which, driven by the osmotic imbalance, water flowed from a pure-water compartment into a salt solution through two membranes formed entirely by stacked CNTs. 266 These enabled single-file water flow in an almost frictionless fashion.
A number of studies have suggested that imposed E-fields may induce ice formation within CNTs. For example, in a study 267 of the phase transition of water inside a (4,13) (d = 1.2 nm) CNT, an ice−ice phase transition occurred discontinuously as the strength of the E-field changed. A further study of TIP4P water inside a wider (3,15) singlewalled CNT (d = 1.31 nm) 268 observed populations of three different polygonal ice nanotubes modulated by both temperature and E-field along the pore axis. An exploration of the effect of axial E-fields on the structure of SPC water in CNTs from (7,7) to (10,10) 269 showed that E-fields of 1−2 V/nm along z induced a transition from liquid to ice nanotubes. Similarly, simulations of SPC and of TIP4P waters in a (8,8) CNT 270 showed that an axial E-field of 1 V/nm caused water to form a solid-like structure at all simulation temperatures up to 350 K. This suggested that the E-field induced a phase transition from liquid to ice-nanotube. Applying a high piston pressure (∼500 MPa) can lead to water flow though the CNT even in the presence of an axial electrostatic field. When such a pressure is applied to a water/methanol mixture, along with an electrostatic field along the length of the CNT, only water is able to flow through the CNT due to the ordering of water molecules within the CNT. 271 This results in an effective separation of the mixture.
Much of the motivation for studies of water in CNTs arises from a desire to design novel biomimetic nanopores based on CNT structures, 272,273 especially with a view to potential applications in water desalination. 9,274−276 In particular, such studies have focused on how smaller CNTs support rapid water transport while rejecting ions (as reviewed by Thomas et al. 277 ). Simulation studies have also shown that functionalization of CNTs with charged groups (which is admittedly easier to demonstrate computationally than to achieve experimentally) can modify ion and water interactions with CNTs. 272,275 This is of interest as, for example, recent studies of pressure-Chemical Reviews pubs.acs.org/CR Review induced water flow through (unmodified) CNTs have suggested that increasing ionic strength (from 0 to 2 and to 4 M) decreases flow, most likely due to ions blocking water entry to the nanotube. 278 Taken together, these studies suggest that modifications to CNTs may be used to optimize water flow and ion rejection, with MD simulations playing a key role in designing such modifications.
Other recent studies have focused on more disparate aspects of water in CNTs. Thus, simulations of SPC/Ew water have suggested superpermittivity of nanoconfined water 279 with high (∼700) values of ε Z (i.e., the dielectric constant parallel to the pore axis) near the pore surface arising from fluid structuring and preorientation of water molecules upon confinement. QM/MM simulations (with water treated via QM and the CNT by MM) 280 of water in a (6,6) CNT suggested that the dissociation constant of water K w drops by 3 orders of magnitude with respect to the bulk phase value due to under-coordination of OH − and H 3 O + within the CNT. Of possible relevance to hydrophobic gating (see above), formation of nanobubbles has been studied 281 in wide ( (20,20) of radius ∼ 1.3 nm, and (26,26) of radius ∼ 1.8 nm) CNTs. Stabilization of the nanobubble was seen to be assisted by the adsorption of vapor-like molecules close to the hydrophobic surface and by strong hydrogen bonds at the vapor−liquid interface.
Taken together, this extensive body of research clearly demonstrates that water within CNTs behaves differently from bulk water, especially in terms of the rapid flow of water through CNTs. More recent studies have emphasized that water flow and behavior can be modulated both by local modifications to CNT structure and by external factors such as applied electrostatic fields. It is encouraging that such rich behavior of water can be observed in "simple" nanopores, both for future designs of novel nanopores and for an indication of how water behavior is likely to be "tuned to function" in biological nanopores.

Synthetic Nanopores
Water behavior has been studied in a number of synthetic biomimetic nanopores, including organic and metal−organic structures, as well as synthetic cyclic peptides and pores in polymer membranes. Experimental and computational studies into the effects of water structure and permeability through a range of synthetic nanopore-like structures have been reported, including imidazole I-quartet superstructures, metal−organic nanotubes, zeolites, and metal organic frameworks (MOFs). A useful survey of experimental studies is provided by the papers from a recent Faraday Discussion meeting. 273 Peptide-appended pillar [5]arene synthetic nanopores have been synthesized and their selectivity investigated as a function of pore geometry. 282,283 Most recently, peptide-appended hybrid [4]arene has been synthesized, exhibiting water transport rates of greater than 10 9 waters per second, similar to aquaporins (see below). 284 Experimental studies on water channels formed by stacked imidazole I-quartet superstructures have indicated that water permeability increases with increasing organization in the water structure. 285 In a later study, 286 sum-frequency generation spectroscopy was used to characterize the hydrogen bond network within the channels, revealing the existence of chiral dipolar water wires. This was supported by MD simulations of the structures in lipid bilayers using the CHARMM36 force field for the lipids and the TIP3P water model. Further MD simulations via two different setup procedures suggest that in addition to ordered water wires through the channels, there may be an alternative disordered mechanism of water transport involving perturbations in the surrounding membrane. 287 The diffusion of water into metal organic nanotubes containing UO 2 (UMON) has been measured experimentally. 288 It was found that the rate of diffusion was slower than most porous materials and that the water adopted a structure similar to amorphous ice. Interestingly, the UMON nanotubes demonstrated specificity to water over other common solvents; this is in contrast to many other porous materials. Singlecrystal X-ray diffraction has also revealed structured waters in other novel metal organic nanotubes, for example Zn-(ATIBDC)(bpy)·3H 2 O (MONT-1). 289 More systematic simulation studies of water in metal organic nanotubes have yet to be undertaken.
There is also much interest in water permeation in (nano)porous materials such as zeolites and MOFs, which has been reviewed recently. 290−293 A recent experimental and MD study on zeolites showed that surface defects (which increase the surface hydrophilicity) play a dominant role in decreasing the diffusivity of water through the crystal. 294 This is in agreement with studies of simple nanopores (see above). Simulations of infrared spectra of the MOF MIL-53(Cr) using the TTMF-3 polarizable water model 295 showed that the competing water-framework and water−water interactions considerably reduced water transport rates through the framework compared with in bulk water.
A series of simulation and experimental studies have focused on the behavior of water in biomimetic nanopores based on hydroxylated and functionalized nanoporous silica with pore radii down to <1 nm. 296 A combined quantum (ab initio MD) and MD approach was used to characterize water in 0.6 nm radius pores created in amorphous silica: GCMC was used to determine SPC/E water content and ab initio MD used to evaluate the stability of the solvated model pore. Both hydroxylated (i.e., with SiOH groups lining the pore) and functionalized (with ∼30% of the pore-lining SiOH groups replaced by Si-CH 2 NH 3 + ) pores were simulated, with water behavior and ion exclusion/permeation characterized. Subsequent studies focused on the potential of membranes formed by this class of nanoporous materials for desalination 297 and on hybrid membranes combining nanoporous silica with a porous alumina support, with potential applications for desalination, 297 and for CO 2 separation exploiting nanocontained molecules of the enzyme carbonic anhydrase. 298

Model Peptide Nanopores
A range of synthetic peptides have been shown to form model nanopores in lipid bilayers. Such peptide nanopores have been studied for a number of years 299 as relatively simple models of more complex ion channel proteins. We will not discuss naturally occurring peptides, for example, gramicidin channels, 300 which have been well reviewed in the past. There are a number of simulations of their stability and conformational dynamics, and of the behavior and conduction properties of the single-file water wire within gramicidin channels. 301 There has also been interest in the behavior of water in channels formed by bundles of amphipathic α-helical peptides, 302−304 especially in the context of proton conduction. 305,306 It should be noted that these simulation studies have been limited by the difficulty of obtaining high-resolution structures. Recent Chemical Reviews pubs.acs.org/CR Review structural and biophysical studies 307−309 may enable such models to be revisited from a computational perspective. There has been considerable interest in transmembrane nanotubes formed by stacks of cyclic peptides. 310−316 Hsieh and Liaw have conducted a recent review of applications, including water and ion transport. 317 Recent experimental and MD studies on cyclic peptides have focused on water and ion permeation to aid future biomimetic nanopore design. 318,319 Ruiz et al. have investigated the effect of different pore linings of cyclic peptide nanotubes (CPNTs) on water transport. 320 Using MD, it was determined that the functional group significantly affected the electrostatic potential and therefore the orientation of water molecules within the channel. The nanotube functionalized with glycine (mCPN-Gly) was found to produce the highest water flux out of the nanotubes studied. A further study employed forward osmosis to investigate water permeation through three types of CPNTs of different degrees of pore-lining hydrophobicity. 321 These simulations showed that modifying the pore-lining with hydrophobic functional groups reduced the pore affinity for water due to the hydrophobic groups preventing entry of water, whereas once inside, water transport through the pore was rapid. This may be compared with simple models where increased hydrophobicity of the pore lining was correlated with higher water permeability.

Track-Etched Nanopores in Polymer Membranes
A number of studies have examined nanopores in polymer membranes and confirmed experimentally the existence of hydrophobic gating, for example, in track-etched pores in polyethylene terephthalate (PET) membranes using ioniccurrent recordings. 227,322 It was found that pores chemically modified to have a more hydrophobic lining could be switched between ion nonconducting and conducting states by application of an E-field (Figure 9). Related studies have been conducted for PET pores in which hydrophobic gating can be achieved by modifying the surface charge density instead of through the application of hydrophobic molecules. 323 Azobenzene-derivative-modified polymer nanochannels have been designed to exhibit hydrophobic gating in the presence of both E-fields and light. 324 Ion selectivity has been investigated in pores etched into PET films. 325 MD experiments using the TIP3P water model and CHARMM36 force field agreed well with experiment. These studies illustrate how design principles derived from computational studies (in this case of hydrophobic gating) may be experimentally tested in biomimetic nanopore systems.

What Have We Learned so Far?
Two major outcomes can be drawn from reviewing the simulation and experimental studies of nonbiological pores, which parallel those from simulations of simplified systems. The first is that there is a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Overall, a wider pore and a more hydrophobic lining favor water flow, provided the pore is wettable. The second is that narrow hydrophobic pores may show dewetting, as suggested by simple models of hydrophobic gating of pores and channels. Some of these concepts have been reviewed recently, drawing together simulation and experimental studies 203 to provide an integrated view of wettability and nanoconfined water flow. A model is presented for confined water flow based on the concept of "effective slip", which is expressed as a linear sum of true slip, depending on contact angle, and apparent slip, caused by a spatial variation of the confined water viscosity as a function of wettability and nanopore dimension. The flow rates for water through CNTs with diameters of 1.3 to 7.0 nm are two to five orders of magnitude greater than those calculated using the no-slip Hagen−Poiseuille equation. This indicates that the behavior of nanoconfined water is dominated by the interaction between water molecules and the lining of the nanopore. The "enhancement factor" (for flow through a nanopore relative to the ideal value from the Hagen−Poiseuille equation) thus depends on both hydrophobicity and the diameter of the nanopore. While this integrative model is formulated in the language of fluid mechanics, the underlying molecular details are embedded within the enhancement factor.

Porins
Porins are a major class of membrane protein nanopores found in the outer membranes of Gram-negative bacteria and mitochondria. They form transmembrane pores based upon a β-barrel architecture ( Figure 2B). They generally have high permeabilities to water, ions, and small solute molecules. Their β-barrel architecture confers structural stability, making them attractive templates for possible biotechnological manipulation and exploitation. The structures of a large number of porins have been determined, mainly by X-ray crystallography. Several have been at high resolution (see Table 1), thus revealing details of bound water molecules within their pores. For example, the high-resolution (1.8 Å) structures of two monomeric porins, Chemical Reviews pubs.acs.org/CR Review CymA (PDB ID 4V3H; see Figure 10) and NanC (PDB ID 2WJR), reveal a substantial number of water molecules bound to the walls of the respective pores. Indeed, the highest resolution porin structure, for the anion-selective porin Omp32 at 1.5 Å (PDB ID 3FGR), shows so many waters that the narrowest region of the pore is (almost) completely occluded by bound waters (albeit in the crystal structure at ∼100 K). This is evidence of favorable interactions between water molecules and the pore-lining residues of porins. There have been a number of simulations of porins, largely focused on their ion and antibiotic permeation properties (see the recent review by Pothula et al. 326 ). Although there has not been a systematic survey of the behavior of water within porins, a number of individual studies of water behavior within these pores provide an emerging overall picture. Early simulations of the canonical porin OmpF 327 showed that water molecules within the pore diffused more slowly than in bulk, and structuring of waters in the narrow charged region of the pore was observed. Subsequent simulations of water diffusion through OmpF 328 (using a modified SPC water model) used a first-passage-time approach and compared water behavior in OmpF with that in a (12,12) CNT. The overall diffusive behavior of water in OmpF was similar to that within a homogeneous CNT channel but a small fraction of the water molecules appeared to be trapped by the protein wall for considerable lengths of time. This is consistent with the structures of porins discussed above. The water diffusion coefficient (D) within the hydrophilic pore was 1.35 nm 2 /ns (compared with D = 4.25 nm 2 /ns for bulk SPC water). This is also lower than in for example the hydrophobic pores of CNTs (D = 2.7 nm 2 /ns). This correlates with structural observations of waters bound to the pore walls in a number of porin X-ray structures (see above and Table 1). Analysis of the effect of the transverse E-field (generated by ionizable residues in the central constriction of the hourglass shaped OmpF pore) on dielectric saturation of water within the pore 329 suggested a space-dependent reduction of the dielectric constant of water within the pore down to values of ∼20. MD simulations (using the TIP3P water model) have been used to compare water behavior in the pores OmpF and OmpC, 330 from which residence times suggest that there are more slow-moving waters in the central constricted region of the respective pores. Simulations have been used to explore a relatively simple cylindrical pore formed by CymA (see Figure 10) from which the N-terminus was removed. 331 Electro-osmotic flow of water associated with ions moving in response to an external transmembrane voltage was observed. MD simulations have also been used to explore water behavior in Wzi, 332 the pore of which is largely occluded by long loops that block the pore.
From these various simulations of porins it is evident that even though the pores are relatively wide (see Table 1), water within the pores behaves differently from bulk water, both in terms of its diffusion rate and dielectric behavior. Bound layers of water are seen both in simulations and structures, and correlate with the reduced diffusion of water within porins. These effects are due to charged side chains lining the pores and hence strong local electrostatic fields experienced within the pore by water molecules.

β-Barrel Nanopores
In addition to studies of porins, there is considerable interest in structurally related nanopores formed by β-barrels. 4,333 These include the heptameric transmembrane pore formed by αhemolysin (αHL), 334 which has been widely investigated in terms of DNA detection 335 and other stochastic sensing 336 technologies based on measurement of ionic currents through single nanopores. There has also been parallel interest in the redesign of nanopores based on simplified β-barrel architectures, 337 both in terms of experimental redesign of bacterial outer membrane proteins for stability 338 and computational redesign for selectivity and other aspects of function 339 as biomimetic protein nanopores.
These studies have been paralleled by computational studies of water within both αHL and in designed β-barrel nanopores. For example, Markosyan et al. used a number of different simulation granularities 340 to investigate the effect of confinement within αHL nanopores on the behavior of water and single-stranded DNA. MD simulations of flexible TIP3P water were used to estimate the one-dimensional profile of the dielectric constant along the length of the αHL nanopore. The water dielectric constant dropped from ∼80 (in bulk water and in the wide extracellular mouth of the nanopore) to ∼55 in the narrower (radius ∼ 0.8 nm) transmembrane pore region. A decrease in water density (∼3-fold) was also shown within the pore in the more constricted regions. This aspect was analyzed in more detail in MD simulations of the behavior of SPC water in pores formed by α-HL molecules embedded in a phospholipid bilayer. 341 These studies revealed fast water transport (66 waters/ns corresponding to D = 2.7 nm 2 /ns) comparable to diffusion of water in a CNT 1.2 nm in diameter (D = 2.7 nm 2 /ns). Analysis of water density along the αHL nanopore axis revealed two drops in density corresponding to two rings of hydrophobic amino acid side chains lining the pore. The free energy profile of water along the pore length revealed barriers of ∼3 to 4 kJ/mol corresponding to these hydrophobic rings. More recently electro-osmotic flow of TIP3P water through an αHL pore has been studied by MD simulations. 342 Perhaps surprisingly, simple continuum models were shown to be able to capture the order of magnitude of the simulated electro-osmotic flow. Simulations of electro-osmotic flow have also been performed for the porin CymA. 331  Table 1). Oxygen atoms of resolved water molecules within the pore are shown as red spheres. In the pore radius profile the broken line indicates the radius of a water molecule. Chemical Reviews pubs.acs.org/CR Review There have been a number of attempts to redesign simplified β-barrel nanopores (recently reviewed by Slusky 343 ). For example, FhuA has been re-engineered 344 to provide a physically robust transmembrane protein for the sampling of biomolecular events within a nanopore. Likewise, a "quiet" OmpG nanopore has been engineered which exhibits negligible spontaneous gating and which can be functionalized using cysteine chemistry. 345 PoreDesigner 346 is a workflow to redesign the β-barrel OmpF as a scaffold for specific pore designs. More purely computational studies of nanopore designs based on porins have also been undertaken. 339,347 The behavior of water has been explored in simple models of β-barrel nanopores, varying the shape and radius of the pores, and the nature of the pore-lining side chains, and focusing on the transport of water through the resultant pores. 347 It was shown that a hydrophobic barrier could be introduced into a βbarrel protein nanopore, and that the height of the energetic barrier to water could be adjusted by modifying the number of consecutive rings of hydrophobic side chains. As with simplified models (see section 4.2 above) a hydrophobic barrier prevented both water and ion permeation even though the pore was sterically unoccluded. It was subsequently shown that a hydrophobic gate in a model β-barrel nanopore could be functionally opened by electro-wetting 239 (see Figure 11). Voltage-induced alignment of water dipoles was produced within the hydrophobic gate region of the nanopore, resulting in water entry followed by the permeation of ions through the opened nanopore. When the ionic imbalance that generated a transbilayer potential was dissipated, water was expelled from the hydrophobic gate (i.e., dewetting took place) and the nanopore functionally reclosed. These studies indicate that the properties observed in simple models of nanopores (see above) can be computationally "transplanted" to biomimetic designs of protein nanopore structures. Similar models of βbarrel nanopores containing hydrophobic barriers may provide template structures for DNA nanosensing. 348

Aquaporins
Aquaporins (AQPs) are α-helical membrane proteins which form water-permeable pores in cell membranes in a wide range of organisms, including animals, plants, and bacteria. They provide a high permeability to water (water permeation through a single aquaporin pore is at rates of ∼10 9 water molecules/s). They are highly selective in that many aquaporins (e.g., AqpZ) are permeable only to water. Others, such as the glycerol-aquaporin GlpF are permeable both to water and to small polar (but uncharged) solutes, such as glycerol. Importantly, despite a high permeability to water, aquaporins are not permeable to protons. Structures of many (∼50) aquaporins have been determined, both by X-ray diffraction and by cryo-electron microscopy. Several of these structures are at high resolution Chemical Reviews pubs.acs.org/CR Review and reveal details of the water molecules within the central aquaporin pore. Currently, the highest resolution structure is for a yeast aquaporin Aqy1, 0.88 Å in resolution (PDB ID 3ZOJ). 15 This structure reveals well-defined positions for water molecules along the pore of the Aqy1 monomer (Figure 2A), albeit at the low temperature (∼100 K) of structure determination by macromolecular X-ray crystallography. The overall structure of all aquaporins is tetrameric, with four identical pore subunits arranged around a central axis. The functional unit is an AQP monomer which is formed by a bundle of six transmembrane α-helices (plus two shorter αhelices each of which spans half of the bilayer) surrounding a central pore. The pore lining is formed by a mixture of peptide backbone atoms and side chains, and contains both polar and hydrophobic regions. Of particular importance is the presence of two asparagine (N) side chains (in a pair of NPA (asparagine-proline-alanine) sequence motifs characteristic of AQPs), which lie close to the middle of the pore and which can hydrogen bond to water molecules. There have been numerous simulation studies of AQPs and water over the past couple of decades since the first structures were determined. We will not discuss all of these but shall rather focus on key earlier studies alongside more recent work in this area. Reviews of some of the older simulation studies are provided by Roux and Schulten 349 and by de Groot and Grubmuller 350 with a particular focus on possible mechanisms of proton exclusion alongside rapid water permeation in the latter review. More recent reviews have focused on, for example, comparison of AQPs with nanopores 351 and in a recent Chemical Reviews article on Aqp0. 352 Early simulations (based on relatively low-resolution AQP structures) focused on the organization of water within the AQP pore in relation to permeation mechanisms 353 and explored pressure-induced water flow. 354 As seen in these and other simulations, the single-file column of water molecules within the central pore is hydrogen-bonded such that there is a bipolar arrangement of the water dipoles, with a switch in their direction at the position of the NPA motifs (see above). The bipolar arrangement of waters appears to be maintained by hydrogen bonds of waters to the protein in the vicinity of the NPA motifs along with the overall protein electrostatic field. At the NPA motifs midway along the pore, water−protein electrostatic interactions facilitate this inversion. Simulations have been used to estimate both the osmotic and diffusion permeability of AQPs, 355 and it has been shown that differences in permeability between different AQP species are not explained simply by (small) differences in pore size. 356 Recent simulations to estimate osmotic permeabilities of aquaporins AQP4, AQP5, and GlpF yielded values in good agreement with experimental estimates. 357 Simulations comparing AqpZ (a highly water-selective AQP) and GlpF (which also allows glycerol to pass through the pore) suggest there is a trade-off between water selectivity and permeability. 358 In analyzing simulations, an osmotic permeability matrix approach has been used to relate permeability to the local structure of the pore lining in different AQP species. 359 Simulations of permeability of H 2 O versus D 2 O and comparison with experiments suggested that, in both cases, water permeability of Aqp1 was reduced by 15−20% in D 2 O, most likely due to the lower self-diffusion coefficient of D 2 O. 360 Simulations have been used to explore the selectivity of AQPs for water. For example, estimation of free energy profiles has been used to explore the selectivity for water over other small polar solutes of AQP1, 361 revealing a negative correlation between permeability and solute hydrophobicity. Of particular importance and interest in terms of the behavior of water within the AQP pore is the mechanism of proton exclusion. Their behavior differs from other single-file water pores (e.g., that within gramicidin A, see above) which conduct protons efficiently 362 via the Grotthuss mechanism. 363 There have been a number of simulation studies of the mechanism of blocking proton permeation by AQPs, resulting in some debate over the relative contribution of different factors (including the water dipolar reversal mentioned above and other electrostatic effects). 364−369 Most of this analysis has focused on MD simulations with "standard" water models (generally TIP3P). Ab initio Car−Parrinello MD simulations 370 suggested that polarization of water molecules varied along the pore axis, and the water wire was broken by the NPA motifs. A more general study of long-range proton transfer in biomolecular pumps by QM/MM 371 revealed that GlpF produced water structure in its interior in agreement with previous MD studies.
Initially it was thought that AQPs were permanently open to water permeation. However, several simulation studies have emphasized the role of the conformational dynamics of certain pore-lining side chains in controlling rates of water permeation. 372,373 Furthermore, higher-resolution experimental structures have revealed alternative local conformations of the NPA motifs 374 in addition to suggesting possible gating at the N-terminal mouth of, for example, the yeast Aqy1. 375 For AqpZ, the possible role of an arginine (R189) side chain conformation in relation to water permeability and gating has been discussed. 376 This may in turn be linked to the protonation state of a histidine side chain. 377,378 Side chain motions and gating have also been explored by simulations of Aqp5 379 and Aqp4. 380 What emerges from these studies is a picture of a pore far from static, with local conformational transitions able to modify rates of water permeation.
AQPs have been the focus of a number of simulation studies of how imposed E-fields may modulate the behavior of water within their pores. For example, relatively large voltages (from a physiological perspective, for example, ± 1.5 V) have been shown to affect water permeability via a switch in the position of an arginine side chain. 381 This indirect effect of water behavior is different from that on water in CNTs where imposing an E-field can align all the water molecules and increase water flow. 382 Thus, effects of fields on both water dipole alignments (as in CNTsalso see above 261 ) and on side chain flipping may have effects on AQPs. 383 Pulsed Efields have been shown to influence the histidine side chain (see above) conformation. 384 Interestingly, THz fields have been shown to influence water permeability in CNTs via the breaking of hydrogen bonds, 263 and E-fields have been shown to have effects on water dipole alignment in cyclic peptide nanotube pores (see above). 385 Brief E-field pulses both perpendicular to 386 and axial to 387 the AQP pore have been shown to have effects on water diffusivity and on pore-lining side chain conformations. 388 It has been suggested that circularly polarized E-fields could yield enhanced water permeation, thus resulting in the "electro-pumping" of water through AQP pores. 389 Overall, these studies indicate how the behavior of water confined within a conformationally dynamic nanopore may be modulated both directly and indirectly by externally applied E-fields. Chemical Reviews pubs.acs.org/CR Review Simulation studies of AQPs have provided considerable insights into the behavior of water confined within a complex and dynamic protein pore. A narrow (radius ∼0.1 nm) pore with a varied lining (containing both hydrophobic and hydrophilic regions) accommodates a single file of water molecules (see Figure 12), which are aligned with a midway dipole switch. This complex pore structure helps to enable proton exclusion. Side chain dynamics may modulate water permeation through the pore, and may in turn be altered by externally applied E-fields. In terms of possible applications, there has been some discussion of the possible uses of AQPcontaining membranes for water purification and desalination. 390,391 In this context, a recent in silico feasibility study 392 has suggested truncation of the AQP structure as a possible route to a "streamlined functional analogue" of the native protein.

Other Protein Pores for Polar Solutes
There are a number of pores which enable the transport of small solutes (e.g., formate, urea, and ammonia) across cell membranes (Table 1). These are generally oligomeric, with the monomeric pore as the functional unit. They include the pentameric formate-nitrite transporter (FNT) family, the trimeric urea transporter (UT) family, hexameric pores for urea and acetate (UreI and SatP), and the Amt/Mep/Rh family of trimeric ammonia transporters. The FNTs have an AQP-like fold. 393 The others have more complex folds, but all have pore-like structures surrounded by α-helices (e.g., for the UT-B fold 394 ). A number of the higher-resolution crystal structures reveal water molecules within the pores (Table 1).
Although simulation studies of these pores have focused on their solute transport mechanisms, a number of these studies have also examined the properties of pore water in passing.
Early simulations of FocA focused on the protonation state of a key histidine residue in the center of the pore. 395 Subsequent, more detailed, studies of FNT pores combined MD, QM, and pK a calculations to explore the energetics of the complete permeation cycle, showing that anion binding into the pore is coupled to protonation of the central histidine. 396 Significantly, water was efficiently translocated only when the central histidine side chain was in a protonated state, indicating that water behavior can be modulated by the charge state of residues within a nanopore.
There have been a number of simulation studies of ureapermeable pores of the UT-B family, mainly focusing on the energetics of urea permeation (e.g., the studies by Levin et al., Azouzi et al. and Ariz-Extreme et al. 394,397,398 ). A combined experimental and computational study of UT-B from red blood cells demonstrated that the per-monomer water permeability of UT-B is similar to that of AQP1. 397 Simulations of water permeability were in good agreement with experiments, and free energy calculations showed a central barrier to water permeation of ∼8 kJ/mol. A study of a bacterial UT showed that the urea selectivity filter in the center of the pore is largely dehydrated due to a hydrophobic constriction. 399 A comparison of UT-B and Aqp1 revealed that despite their different folds, the pore constriction radius in both pores was <0.2 nm, presenting a central energetic barrier to water in UT-B of ∼10 kJ/mol. 398 A detailed simulation study of the hexameric Ure1 pore revealed a low barrier to water (∼4 kJ/mol) and a water transport rate (∼4 × 10 9 s −1 ) comparable to that of AQPs. 400 Simulations of the structurally related SatP acetate channel showed that water hydrates all parts of the pore, including the central constriction site which forms an energetic barrier for acetate. 401 A number of studies have explored water behavior in the ammonia channels of the Amt and Rh proteins. Early MD studies showed that the pore of AmtB can accommodate water molecules. 402 Subsequent studies revealed different hydration patterns in Amt and Rh proteins in relation to their transport mechanisms. 403 Waters were seen to be stable in the AmtB pore but significantly less so in the Rh pore. This was suggested to be due to differences in phenylalanine side chain conformations. Interestingly, free energy studies of water and solutes in Rh50 channels showed a barrier for water of ∼25 kJ/ mol, suggesting that Rh supports NH 3 but not water transport. 404 These simulation studies extend our understanding of water in narrow protein pores beyond just AQPs. In many cases, water permeability in simulations of these "other pores" is comparable to that seen for AQPs. Similarly, the local conformation (for example of phenylalanine in Amt/RH) or charge state (of histidine in UTs) of amino acid side chains within a narrow pore can control water permeability and hydration of the pore. As discussed in a recent simulation study of the SecY pore, 405 water within a narrow protein nanopore can exhibit a slowing down in structural relaxation dynamics, which may be compared to the effects of lowering the temperature of bulk water.

ION CHANNELS
Ion channel proteins function as water-filled pores that control the passage of selected ions across biological membranes, commonly involving conformational transitions in response to intra-and/or intercellular signals. Ion channel pores are  Table 1), with the protein in a semitransparent surface representation, and the broken horizontal lines indicating the approximate location of a bilayer membrane relative to the protein.
Chemical Reviews pubs.acs.org/CR Review typically subnanometer-sized and lined by rings of different amino acid residues whose side chains may be oriented toward or away from the central pore, thus conferring nonuniformity both in pore dimensions and hydrophobicity. Such molecular features of the permeation pathway differentially modulate water behavior in distinct regions of the channel, which in some instances may, in turn, constitute part of its physiological function.

Water in Relation to Potassium Channel Permeation
In a large family of cation-selective channels, the tetrameric pore-loop channels (which include potassium channels as well as voltage-gated sodium and calcium channels), ion selectivity is controlled by a filter region lined by a short polypeptide segment from each of the four subunits on one side of the membrane. There have been extensive simulation studies of ion permeation and selectivity mechanisms of these channels. 24,349,406−409 The exact mechanisms of selectivity and permeation remain unresolved and are areas of active research by a number of groups (see, e.g., ref 410 for an important collection of reviews and ref 411 for a more recent survey of advances and challenges). Here, we will focus more specifically on the behavior of water which remains central to our understanding of permeation mechanisms. MD simulations and electrostatic calculations of TIP3P water inside the canonical potassium channel KcsA showed that the presence of water in the selectivity filter enhanced the depth of potential energy wells occupied by a potassium ion, suggesting a structural role for water in the region. 412 Water molecules (TIP3P and SPC) have also been observed to travel through alternative pathways "behind" the KcsA selectivity filter in the absence of ions, giving rise to a high osmotic permeability of the channel when in a potassium-depleted state. 413 This is consistent with experimental studies of water transport through the KcsA channel. 414,415 Coupling of ion and water flux has also been observed experimentally in other potassium channels. 416,417 More recently, there has been some debate about the interplay of water and ions in conduction through potassium channels. The "canonical" model envisages the narrow (singlefile) selectivity filter being occupied by alternating potassium ions and water molecules. 406 However, recent simulations using computational electrophysiology 418 (i.e., in the presence of an E-field due to an ionic charge imbalance across the membrane) suggest a direct knock-on mechanism in which K + ions move through the selectivity filter in the absence of intervening water molecules. 407 Data from two-dimensional infrared spectroscopy experiments have been interpreted as implicating water in the filter, 24,25 although this has been contested, 419 and solid-state NMR experiments suggest the conduction pathway in the filter may be water-free. 420 Although details remain to be resolved (see ref 411 for a measured review), it is clear that there is a complex role of conformational flexibility within the selectivity filter and elsewhere in the KcsA channel, 25 which determines the relative occupancies of the filter by ions or water in its different conformational states. This is consistent with the more general emerging view of biological nanopores (see above) whereby local changes in conformation and/or electrostatics of narrow protein pores can modulate their water permeability.

Water and Proton Transport in Ion Channels
Water behavior is of particular interest in those ion channels and related membrane proteins which form proton-permeable pores. This has been studied in the past 421,422 in simple model peptide nanopores such as gramicidin A (see section 5.2 above), which contains a single file of water and allows rapid proton flux via a Grotthuss mechanism. This remains a valuable test system for the direct simulation of proton conduction, for example, via multistate EVB simulations 423 or by combining classical MD simulations with a stochastic model for proton hopping. 424 EVB simulations have also been used to explore proton transport through water-filled CNTs, 159 and yield similar results to ab initio MD simulations. Protonconducting water wires are rather more complex in a number of proton-permeable ion channels for which structural data have recently become available, including the voltage-gated proton channel Hv1 and the low-pH-activated proton channel M2 from influenza A.
Early studies of Hv1 function, combined with homology model-based simulations, suggested that a water wire existed within the channel and played a role in proton permeation. 425,426 The determination of the crystal structure 427 and subsequent MD simulations revealed fluctuations in the water density within the pore. 428 These results suggested that formation of a continuous water wire along the entire length of the (dynamic) channel was unlikely to enable a simple Grotthuss-like proton conduction mechanism. Subsequent studies have focused on the interplay of water molecules and charged side chains in the proton conduction mechanism. 429,430 The low-pH-activated M2 proton channel from influenza A virus has been the subject of intensive structural and computational studies of the relationship between channel structure, the nature of water within the channel, and the likely mechanism(s) of proton permeation. Initial simulations of water within M2 were based on plausible models of the tetrameric bundle of M2 helices spanning a membrane 431,432 and included simulations of proton permeation using multistate EVB simulations. 433 Subsequently, structures of the M2 channel were determined using X-ray diffraction, 434 solution NMR of the protein in a detergent micelle, 435 and by solidstate NMR in a lipid bilayer. 436 These and subsequent studies have emphasized the conformational flexibility of this "simple" ion channel. 437 Subsequent simulations based on these structures 438,439 explored the interplay of the protonation state of the key pH-sensing His37 side chain, the channel conformation, and the behavior of water within the transmembrane pore, suggesting that a simple water wire Grotthuss mechanism might not apply. More recent high-resolution structures (see Figure 2D; PDB ID 5JOO resolved at 1.4 Å) obtained using XFELs 19 reveal that at low pH there is a continuous hydrogen-bonded network of water molecules through the length of the channel, consistent with a Grotthuss mechanism for proton transport. Overall these studies reveal the complexities of a dynamic channel structure and its effects on water within a fluctuating charged pore and emphasize the need for further simulations based on recent higher-resolution structures 440 to better understand water behavior in such a complex dynamic pore. An example is provided by recent simulations of a Grotthuss-like mechanism from the Voth group. 441 Proton-permeable channels are also seen in a number of mitochondrial and related proteins, including the proton-Chemical Reviews pubs.acs.org/CR Review translocating transhydrogenase for which structures 442,443 and simulations 442,443 reveal a dry, hydrophobic region in the channel. Wetting is promoted by protonation and the change in conformation of a key histidine residue. In complex I, a redox-driven proton pump, channel hydration has been suggested by simulation to be sensitive to the protonation states of buried lysine residues. 444 In summary, these studies of proton-permeable pores suggest that the influence of pore size, hydrophobicity, and local conformation on water behavior seen in simpler nanopores is exploited by nature/evolution to control water hydrogen bonding and proton permeation through more complex channels and related structures.

Water and Hydrophobic Gating of Ion Channels
Most ion channels form a continuous pathway for conduction of their permeating species only when in an active (i.e., open) state, with water and ion flow interrupted at one or more positions (i.e., gates) when at rest or inactivated. In the case of a hydrophobic gate (see section 4.2 above), which favors a shift in equilibrium from the liquid to the vapor phase, closure to conduction is achieved energetically without requiring steric occlusion of the pore. Initially characterized using model nanopores, observations of dewetting in hydrophobic regions of channel pores have been increasingly reported in structural studies of ion channels and implicated as an important feature to consider in the mechanistic elucidation of these proteins. 209 The first direct lines of evidence for hydrophobic gating in a biological ion channel were provided by studies on the smallconductance mechanosensitive channel, MscS, one member of a family of bacterial proteins that protect against hypo-osmotic shock. 445 During atomistic MD simulations of TIP3P water in a crystal structure of the protein, 446 liquid−vapor transitions were observed within the channel in a hydrophobic region where the central pore radius was ≥0.25 nm. 447,448 The dewetting region was surrounded by side chains of leucine residues and considered to form a closed hydrophobic gate in the structure, preventing ion conduction in such a conformation. 447 By contrast, a hydrophilic mutation at the constriction, substituting polar serine amino acids for one of the two leucine rings, yielded stable hydration through the channel under the same simulation conditions. 447 An open conformation was subsequently resolved for the protein, showing a widened channel pore with radius >0.6 nm. 449 Further simulation studies examined pore dewetting and hydration using TIP3P and SPC water in models of both MscS and the related large-conductance mechanosensitive channel, MscL, supporting the role of channel pore hydrophobicity (and consequent wettability) in the gating of both proteins. 450 8.3.1. Pentameric Ligand-Gated Ion Channels (pLGICs). Members of the pLGIC family mediate fast synaptic transmission and are the subject of considerable biomedical interest. The transmembrane portion of their pores is each lined by five M2 helices that contain a highly conserved hydrophobic residue in their 9′ position. 451 This region was first demonstrated to constitute a hydrophobic gate by atomistic simulations of SPC water in the M2 helix bundle of a nicotinic acetylcholine receptor structure, 452 in which narrowing of the permeation pathway to ∼0.3 nm in radius by two consecutive pore-lining rings formed by the 9′ leucine and 13′ valine residues was calculated to present an energetic barrier to water and ion passage. 453 Consistent with such simulation results, effects of mutagenesis of the hydrophobic rings on channel electrophysiology were reported by a number of experimental studies. 454−457 Later studies on two prokaryotic homologues of the channel, ELIC and GLIC, Figure 13. Structure and simulation of water in the transmembrane pore domain of the serotonin 5-HT 3 receptor (PDB ID 4PIR). The overall protein architecture is in ribbon representation, with its transmembrane pore surface determined using CHAP and colored by hydrophobicity, 467 and with pore-lining amino acid side chains shown as sticks. Only two out of the five protein subunits are displayed in the zoomed-in view for clarity. Two snapshots from an equilibrium simulation of water (red-and-white spheres) within the membrane-embedded protein (gray surface) are shown. Initially (at t = 0) the pore is hydrated (waters in red/white) but shortly afterward (t = 1 ns) it has dewetted. Chemical Reviews pubs.acs.org/CR Review respectively, having phenylalanine and isoleucine in the 9′ position, showed dewetting of both proteins around this region in their simulated structures, conferring closure to the putative hydrophobic gate. 458,459 For GLIC, liquid-to-vapor transitions were reported during string-method MD simulations of the isolated transmembrane domain from its putative open conformation to a nonconductive state. 460 Comparable studies have also been performed for the entire GLIC channel. 461,462 For the isolated transmembrane domain, the free-energy cost of hydrating the hydrophobic gate (with TIP3P water) was estimated. 463 The hydration state at such a region was suggested to be a principal determinant of the functional (ion conduction) state of the channel based on calculations from the latter study, whereby pore hydration was associated with an increase in free energy of ∼46 kJ mol −1 , compared with the smaller energetic cost (∼5 kJ mol −1 ) of transferring a sodium ion into a hydrated gate. 463 Functioning of pore-lining hydrophobic side chains located midway through the transmembrane pathway as a gate in putative closed conformations of pLGICs, with conduction prevented even in the absence of steric obstruction, has since also been substantiated by MD simulations of the serotonin 5-HT 3 receptor ( Figure 13). 8,464 Using a structure of a glycine receptor, 465 another member of the pLGIC family, the sensitivity of the analysis to various nonpolarizable water models (SPCE, TIP3P, and TIP4P) was also examined, yielding similar water density and free energy profiles. 8 The comparison of water models was further extended in MD simulations of the 5-HT 3 receptor in multiple (both closed and open) states. 163 The open (i.e., wetted) state conformation of the pore exhibited degrees of hydration ranging from partial wetting using the TIP4P/2005 water model to complete wetting for the polarizable AMOEBA14 model. These results indicate that induced polarization can influence the behavior of water and ions within hydrophobic gates of nanoscale pores. The same 5-HT 3 pore has been used to explore electrowetting (see section 4.3 above), revealing that a (supra-physiological) potential difference of ∼0.85 V across the membrane was required to hydrate the hydrophobic gate of the closed (i.e., dewetted) state of the pore. 466 8.3.2. Hydrophobic Gating in Other Channel Families. Dewetting of a hydrophobic region in the inner pore has been implicated to play a role in gating among members of several potassium channel families (as reviewed by Aryal et al. 209 ). This is supported by atomistic MD simulations of SPC and TIP4P water in the voltage-gated potassium channel, 468,469 of the SPC and TIP4P models in the TWIK-1 channel of the two-pore-domain (K2P) family, 470 and of SPC and TIP3P in a homology model of the large-conductance calcium-activated potassium (BK) channel structure. 471 In the case of TWIK-1, significance of the hydrophobic barrier was further validated by electrophysiological measurements, where a hydrophilic substitution led to enhanced channel activity in whole-cell current recordings. 470 Evidence is also emerging for the presence of hydrophobic residues at the gates of TRP (transient receptor potential) channels, a class of nonselective cation channels whose members fulfill diverse functions in sensory physiology. 472,473 The temperature dependence of the hydration equilibrium within TRPV1 has been suggested to contribute to the biological function of the protein as a temperature sensor. 474 The presence of hydrophobic gates and their functional significance continue to be reported for other cation and anion channels as well. Investigated using TIP3P water, pore wetting events in the CorA magnesium channel, concurrent with small increases in pore diameter, have been suggested to mediate control of ion conduction by regulatory ions. 475 A combination of spectroscopy and electrophysiology have demonstrated coupling of ion conduction to hydration of transmembrane helices in a channelrhodopsin. 476,477 In bestrophin-1, a calcium-activated chloride channel, TIP4P water molecules were expelled from the closed state of the hydrophobic gate 478 but hydrated the pore in simulations of in silico mutants of the protein with polar amino acid substitutions in the gate region. 479 8.3.3. Integrative View of Water and Hydrophobic Gating of Ion Channels. An exhaustive survey of water behavior in biological channels has provided an integrative view of hydrophobic gating. MD simulations using the TIP4P water model were run for nearly 200 ion channel structures. 217 The hydration/dewetting behavior of the channels was analyzed as a function of local radius and hydrophobicity of the pores (Figure 14). This revealed a clear relationship between the radius and local hydrophobicity of the pore and the free energy of water within the pore. Narrow (R < 0.4 nm) hydrophobic pores show a clear barrier to water entry. Fitting the simulation data from this analysis allowed a simple (simulation-free) heuristic model to be developed that enables accurate prediction of the presence of a closed hydrophobic gate in an ion channel structure based solely on measurement of pore radius and hydrophobicity profiles. 467 This in turn enables rapid functional annotation of new channel structures as soon as they are determined.

CONCLUSIONS AND CRITICAL REFLECTIONS
The studies reviewed above range from simulations of water in idealized models of nanopores through to complex biological ion channels. For the idealized models, a focus of many studies has been on water flow in nanopores. This has been shown to be potentially faster than Hagen−Poiseuille-predicted flow if Chemical Reviews pubs.acs.org/CR Review the pore-lining is sufficiently hydrophobic. Simulations of nonbiological nanopores reinforce this view, revealing a complex relationship between pore size and geometry, the nature of the pore lining, and rates of water transport. Overall, a wider pore and a more hydrophobic lining favor water flow, provided the pore is wettable. Thus, the strength with which water interacts with the pore-lining surface is crucial in determining its local behavior within nanopores. Pore-lining surface roughness and dynamic fluctuations, as well as the overall pore geometry, are also important. Biological pores and channels exhibit a range of behaviors. β-Barrel bacterial outer membrane proteins confine water within hydrophilic pores with multiple water binding/interaction sites. Aquaporins are narrow and less hydrophilic, and have high water permeation rates. Ion channels confine water and may exhibit hydrophobic gating, which was first predicted in simple models of nanopores, and has been demonstrated experimentally in nonbiological nanopores. The overall picture emerging from these studies is summarized in Figure 15 where the behavior of water within a nanopore is considered as a function of the hydrophobicity and radius of the pore. There are four regions, as exemplified by (i) wide (radius > 1 nm) and hydrophilic nanopores, which will exhibit behavior comparable to bulk water; (ii) wide (radius > 1 nm) and hydrophobic nanopores, which should enable water to flow faster than is predicted by the Hagen− Poiseuille theory, with enhancement factors substantially >1 possible for nanopores; (iii) many biological pores and (open) channels are narrow (radius < 1 nm) and hydrophilic, and thus contain water which interacts with the pore-lining residues; and (iv) some ion channels exhibit hydrophobic gating, with dewetting in their narrow (radius < 1 nm) and hydrophobic constrictions, which in turn closes the channel to the flow of ions.
The biological relevance of these wide-ranging simulation studies is in helping us to understand the emerging diversity of channel structures. Being able to position a new channel structure into its corresponding region in Figure 15, and to characterize how conformational changes in the pore may move a channel from one region to another, will inform our mechanistic understanding of the regulation of biological channels, for example by hydrophobic gating. Such a schematic view may be used as a "template" to aid the design of novel nanopores such as those functioning as artificial water channels. 480 There is also further work to be done regarding the elucidation of the effects of electric fields on water behavior, when designing more complex pores, as a possible mechanism for the switching on/off of novel nanopores.
It is informative to reflect critically on the progress summarized in this review. Two main questions arise, one concerning optimal systems for simulation studies of water in nanopores, and the other concerning simulation methodologies and models.
Nanopore systems for simulations span a spectrum from highly simplified idealized models, offering more general insights into the behavior of nanoconfined water, through to complex biological pores which are gated, that is, switchable between open and closed states. Simplified models provide general insights into e.g. the relationship between pore radius, hydrophobicity, and flow rate of water. However, it remains uncertain to what extent such generalizations can be used in a quantitative and predictive fashion. The latter is an important consideration if we wish to exploit our general understanding in the design of novel nanopores with desired functional properties, and in annotation of the functional properties of novel biological channels and pores (as an increasing number of structures of these are determined).
A more methodological consideration is whether the simulation methodologies and models used to study water in nanopores are of sufficient granularity and/or accuracy for the task. Many studies to date have used conventional MD with well-established water models (e.g., TIP3P or SPC/E). Perhaps surprisingly, continuum models (modified by insights from atomistic MD simulations) have also provided valuable insights into the behavior of water in simple nanopores. It is likely that for the design of novel nanopores and in order to understand complex biological nanopores, accurate atomistic simulations of water properties are needed. It remains uncertain whether fixed point-charge MD models of water are sufficient, or whether polarizable (or even QM) models are needed in some situations (e.g., hydrophobic gating). This is of particular importance in the context of confined water systems as classical (i.e., MM) water models are typically developed and validated for how accurately they reproduce bulk water properties. It, therefore, may be timely to use more advanced water models and to evaluate the robustness to changes in water model when running simulation studies of water behavior in complex nanopores, particularly as preliminary studies indicate that the nanoscale behavior of water is sensitive to water model. Quantum mechanical approaches offer tantalizingly even higher levels of accuracy. This comes at a significant increase in computational cost, and there are still many levels of theory and parametrizations to evaluate. Substantial progress is being made in rendering quantum mechanical methods more computationally feasible for biological systems. There is still some way to go before it will be possible to perform quantum mechanical calculations on Figure 15. Schematic summary of the behavior of water within nanopores as a function of the hydrophobicity and radius of the pores. There are four regions, as exemplified in (i) wide (R > 1 nm) hydrophilic nanopores, exhibiting behavior comparable to bulk water; (ii) wide (R > 1 nm) hydrophobic nanopores, which enable water to flow faster than predicted by the Hagen−Poiseuille theory; (iii) narrow (R < 1 nm) hydrophilic nanopores, corresponding to many biological pores; and (iv) biological channels which exhibit hydrophobic gating, with dewetting in their narrow (R < 1 nm) and hydrophobic constrictions. In this figure, ε is the ratio of the water flow compared to that predicted by the Hagen−Poiseuille equation for bulk water viscosity and no-slip boundary conditions. Chemical Reviews pubs.acs.org/CR Review system sizes and time scales comparable to those accessed via classical methodologies.
Returning to the schematic summary ( Figure 15) in light of these reflections, we are currently able to design a nanopore which promotes wetting over dewetting or vice versa. But can we design a novel nanopore which enables fast, selective flow of water, that is, can we match aquaporins by de novo design? And can we design such a nanopore that it is also switchable in a controlled fashion? Fully answering these two key questions will require more detailed simulations alongside experimental evaluation, as exemplified by, for example, a recent study of gating of hydrophobic nanopores with large anions. 481

AUTHOR INFORMATION
Corresponding Author