Assessment of Hydration Thermodynamics at Protein Interfaces with Grid Cell Theory.

Molecular dynamics simulations have been analyzed with the Grid Cell Theory (GCT) method to spatially resolve the binding enthalpies and entropies of water molecules at the interface of 17 structurally diverse proteins. Correlations between computed energetics and structural descriptors have been sought to facilitate the development of simple models of protein hydration. Little correlation was found between GCT-computed binding enthalpies and continuum electrostatics calculations. A simple count of contacts with functional groups in charged amino acids correlates well with enhanced water stabilization, but the stability of water near hydrophobic and polar residues depends markedly on its coordination environment. The positions of X-ray-resolved water molecules correlate with computed high-density hydration sites, but many unresolved waters are significantly stabilized at the protein surfaces. A defining characteristic of ligand-binding pockets compared to nonbinding pockets was a greater solvent-accessible volume, but average water thermodynamic properties were not distinctive from other interfacial regions. Interfacial water molecules are frequently stabilized by enthalpy and destabilized entropy with respect to bulk, but counter-examples occasionally occur. Overall detailed inspection of the local coordinating environment appears necessary to gauge the thermodynamic stability of water in protein structures.

GPCR structure, function, drug discovery and crystallography: report from Academia-Industry International Conference (UK Royal Society) Chicheley Hall, 1-2 September 2014.

G-protein coupled receptors (GPCRs) are the targets of over half of all prescribed drugs today. The UniProt database has records for about 800 proteins classified as GPCRs, but drugs have only been developed against 50 of these. Thus, there is huge potential in terms of the number of targets for new therapies to be designed. Several breakthroughs in GPCRs biased pharmacology, structural biology, modelling and scoring have resulted in a resurgence of interest in GPCRs as drug targets. Therefore, an international conference, sponsored by the Royal Society, with world-renowned researchers from industry and academia was recently held to discuss recent progress and highlight key areas of future research needed to accelerate GPCR drug discovery. Several key points emerged. Firstly, structures for all three major classes of GPCRs have now been solved and there is increasing coverage across the GPCR phylogenetic tree. This is likely to be substantially enhanced with data from x-ray free electron sources as they move beyond proof of concept. Secondly, the concept of biased signalling or functional selectivity is likely to be prevalent in many GPCRs, and this presents exciting new opportunities for selectivity and the control of side effects, especially when combined with increasing data regarding allosteric modulation. Thirdly, there will almost certainly be some GPCRs that will remain difficult targets because they exhibit complex ligand dependencies and have many metastable states rendering them difficult to resolve by crystallographic methods. Subtle effects within the packing of the transmembrane helices are likely to mask and contribute to this aspect, which may play a role in species dependent behaviour. This is particularly important because it has ramifications for how we interpret pre-clinical data. In summary, collaborative efforts between industry and academia have delivered significant progress in terms of structure and understanding of GPCRs and will be essential for resolving problems associated with the more difficult targets in the future.

Prediction of Small Molecule Hydration Thermodynamics with Grid Cell Theory.

An efficient methodology has been developed to quantify water energetics by analysis of explicit solvent molecular simulations of organic and biomolecular systems. The approach, grid cell theory (GCT), relies on a discretization of the cell theory methodology on a three-dimensional grid to spatially resolve the density, enthalpy, and entropy of water molecules in the vicinity of solute(s) of interest. Entropies of hydration are found to converge more efficiently than enthalpies of hydration. GCT predictions of free energies of hydration on a data set of small molecules are strongly correlated with thermodynamic integration predictions. Agreement with the experiment is comparable for both approaches. A key advantage of GCT is its ability to provide from a single simulation insightful graphical analyses of spatially resolved components of the enthalpies and entropies of hydration.

Evaluation of Host-Guest Binding Thermodynamics of Model Cavities with Grid Cell Theory.

A previously developed cell theory model of liquid water was used to evaluate the excess thermodynamic properties of confined clusters of water molecules. The results are in good agreement with reference thermodynamic integration calculations, suggesting that the model is adequate to probe the thermodynamic properties of water at interfaces or in cavities. Next, the grid cell theory (GCT) method was applied to elucidate the thermodynamic signature of nonpolar association for a range of idealized host-guest systems. Polarity and geometry of the host cavities were systematically varied, and enthalpic and entropic solvent components were spatially resolved for detailed graphical analyses. Perturbations in the thermodynamic properties of water molecules upon guest binding are restricted to the immediate vicinity of the guest in solvent-exposed cavities, whereas longer-ranged perturbations are observed in buried cavities. Depending on the polarity and geometry of the host, water displacement by a nonpolar guest makes a small or large enthalpic or entropic contribution to the free energy of binding. Thus, no assumptions about the thermodynamic signature of the hydrophobic effect can be made in general. Overall the results warrant further applications of GCT to more complex systems such as protein-ligand complexes.

Toward an understanding of agonist binding to human Orexin-1 and Orexin-2 receptors with G-protein-coupled receptor modeling and site-directed mutagenesis.

The class A G-protein-coupled receptors (GPCRs) Orexin-1 (OX1) and Orexin-2 (OX2) are located predominantly in the brain and are linked to a range of different physiological functions, including the control of feeding, energy metabolism, modulation of neuro-endocrine function, and regulation of the sleep-wake cycle. The natural agonists for OX1 and OX2 are two neuropeptides, Orexin-A and Orexin-B, which have activity at both receptors. Site-directed mutagenesis (SDM) has been reported on both the receptors and the peptides and has provided important insight into key features responsible for agonist activity. However, the structural interpretation of how these data are linked together is still lacking. In this work, we produced and used SDM data, homology modeling followed by MD simulation, and ensemble-flexible docking to generate binding poses of the Orexin peptides in the OX receptors to rationalize the SDM data. We also developed a protein pairwise similarity comparing method (ProS) and a GPCR-likeness assessment score (GLAS) to explore the structural data generated within a molecular dynamics simulation and to help distinguish between different GPCR substates. The results demonstrate how these newly developed methods of structural assessment for GPCRs can be used to provide a working model of neuropeptide-Orexin receptor interaction.

Fighting obesity with a sugar-based library: discovery of novel MCH-1R antagonists by a new computational-VAST approach for exploration of GPCR binding sites.

Obesity is an increasingly common disease. While antagonism of the melanin-concentrating hormone-1 receptor (MCH-1R) has been widely reported as a promising therapeutic avenue for obesity treatment, no MCH-1R antagonists have reached the market. Discovery and optimization of new chemical matter targeting MCH-1R is hindered by reduced HTS success rates and a lack of structural information about the MCH-1R binding site. X-ray crystallography and NMR, the major experimental sources of structural information, are very slow processes for membrane proteins and are not currently feasible for every GPCR or GPCR-ligand complex. This situation significantly limits the ability of these methods to impact the drug discovery process for GPCR targets in "real-time", and hence, there is an urgent need for other practical and cost-efficient alternatives. We present here a conceptually pioneering approach that integrates GPCR modeling with design, synthesis, and screening of a diverse library of sugar-based compounds from the VAST technology (versatile assembly on stable templates) to provide structural insights on the MCH-1R binding site. This approach creates a cost-efficient new avenue for structure-based drug discovery (SBDD) against GPCR targets. In our work, a primary VAST hit was used to construct a high-quality MCH-1R model. Following model validation, a structure-based virtual screen yielded a 14% hit rate and 10 novel chemotypes of potent MCH-1R antagonists, including EOAI3367472 (IC50 = 131 nM) and EOAI3367474 (IC50 = 213 nM).

Prediction of cyclin-dependent kinase 2 inhibitor potency using the fragment molecular orbital method.

BACKGROUND: The reliable and robust estimation of ligand binding affinity continues to be a challenge in drug design. Many current methods rely on molecular mechanics (MM) calculations which do not fully explain complex molecular interactions. Full quantum mechanical (QM) computation of the electronic state of protein-ligand complexes has recently become possible by the latest advances in the development of linear-scaling QM methods such as the ab initio fragment molecular orbital (FMO) method. This approximate molecular orbital method is sufficiently fast that it can be incorporated into the development cycle during structure-based drug design for the reliable estimation of ligand binding affinity. Additionally, the FMO method can be combined with approximations for entropy and solvation to make it applicable for binding affinity prediction for a broad range of target and chemotypes. RESULTS: We applied this method to examine the binding affinity for a series of published cyclin-dependent kinase 2 (CDK2) inhibitors. We calculated the binding affinity for 28 CDK2 inhibitors using the ab initio FMO method based on a number of X-ray crystal structures. The sum of the pair interaction energies (PIE) was calculated and used to explain the gas-phase enthalpic contribution to binding. The correlation of the ligand potencies to the protein-ligand interaction energies gained from FMO was examined and was seen to give a good correlation which outperformed three MM force field based scoring functions used to appoximate the free energy of binding. Although the FMO calculation allows for the enthalpic component of binding interactions to be understood at the quantum level, as it is an in vacuo single point calculation, the entropic component and solvation terms are neglected. For this reason a more accurate and predictive estimate for binding free energy was desired. Therefore, additional terms used to describe the protein-ligand interactions were then calculated to improve the correlation of the FMO derived values to experimental free energies of binding. These terms were used to account for the polar and non-polar solvation of the molecule estimated by the Poisson-Boltzmann equation and the solvent accessible surface area (SASA), respectively, as well as a correction term for ligand entropy. A quantitative structure-activity relationship (QSAR) model obtained by Partial Least Squares projection to latent structures (PLS) analysis of the ligand potencies and the calculated terms showed a strong correlation (r2 = 0.939, q2 = 0.896) for the 14 molecule test set which had a Pearson rank order correlation of 0.97. A training set of a further 14 molecules was well predicted (r2 = 0.842), and could be used to obtain meaningful estimations of the binding free energy. CONCLUSIONS: Our results show that binding energies calculated with the FMO method correlate well with published data. Analysis of the terms used to derive the FMO energies adds greater understanding to the binding interactions than can be gained by MM methods. Combining this information with additional terms and creating a scaled model to describe the data results in more accurate predictions of ligand potencies than the absolute values obtained by FMO alone.

Compound Design by Fragment-Linking.

The linking together of two fragment compounds that bind to distinct protein sub-sites can lead to a superadditivity of binding affinities, in which the binding free energy of the linked fragments exceeds the simple sum of the binding energies of individual fragments (linking coefficient E<1). However, a review of the literature shows that such events are relatively rare and, in the majority of the cases, linking coefficients are far from optimal being much greater than 1. It is critical to design a linker that does not disturb the original binding poses of each fragment in order to achieve successful linking. However, such an ideal linker is often difficult to design and even more difficult to actually synthesize. We suggest that the chance of achieving successful fragment linking can be significantly improved by choosing a fragment pair that consists of one fragment that binds by strong H-bonds (or non-classical equivalents) and a second fragment that is more tolerant of changes in binding mode (hydrophobic or vdW binders). We also propose that the fragment molecular orbital (FMO) calculations can be used to analyse the nature of the binding interactions of the fragment hits for the selection of fragments for evolution, merging and linking in order to optimize the chance of success.

Using electrophysiology and in silico three-dimensional modeling to reduce human Ether-à-go-go related gene K(+) channel inhibition in a histamine H3 receptor antagonist program.

The histamine H3 receptor (H3R) plays a regulatory role in the presynaptic release of histamine and several other neurotransmitters, and thus, it is an attractive target for central nervous system indications including cognitive disorders, narcolepsy, attention-deficit hyperactivity disorder, and pain. The development of H3R antagonists was complicated by the similarities between the pharmacophores of H3R and human Ether-à-go-go related gene (hERG) channel blockers, a fact that probably prevented promising compounds from being progressed into the clinic. Using a three-dimensional in silico modeling approach complemented with automated and manual patch clamping, we were able to separate these two pharmacophores and to develop highly potent H3R antagonists with reduced risk of hERG liabilities from initial hit series with low selectivity identified in a high-throughput screening campaign.

Modeling neuronal nicotinic and GABA receptors: important interface salt-links and protein dynamics.

Protein motions in the Cys-loop ligand-gated ion receptors that govern the gating mechanism are still not well understood. The details as to how motions in the ligand-binding domain are translated to the transmembrane domain and how subunit rotations are linked to bring about the cooperative movements involved in gating are under investigation. Homology models of the alpha4beta2 nicotinic acetylcholine (nACh) and beta2alpha1gamma2 GABA receptors were constructed based on the torpedo neuromuscular-like nicotinic receptor structure. The template constructed for the full electron microscopy structure must be considered more reliable for structure-function studies due to the preservation of the E45-R209 salt-link. Many other salt-links are seen to transiently form, including switching off of the E45-R209 link, within a network of potential salt-links at the binding domain to the transmembrane domain interface region. Several potentially important intersubunit salt-links form in both the nAChR and GABAR structures during the simulation and appear conserved across many subunit combinations, such as the salt-link between alpha4.E262 and beta2.K255 in nAChR (beta2.E262 and alpha1.K263 in GABAR), at the top of the pore-lining M2 helices, and the intersubunit link of R210 on the M1-linker to E168 on the beta8-sheet of the adjacent subunit in the GABA receptor (E175-K46 being the structurally equivalent link in the nAChR, with reversed polarity). A network of other salt-links may be vital for transmitting the cooperative gating motions between subunits that become biased upon ligand binding. The changes seen in the simulations suggest that this network of salt-links helps to set limits and specific states for the conformational changes involved in gating of the receptor. We hope that these hypotheses will be tested experimentally in the near future.

Cell-free expression for nanolipoprotein particles: building a high-throughput membrane protein solubility platform.

Membrane-associated proteins and protein complexes account for approximately a third or more of the proteins in the cell (1, 2). These complexes mediate essential cellular processes; including signal transduc-tion, transport, recognition, bioenergetics and cell-cell communication. In general, membrane proteins are challenging to study because of their insolubility and tendency to aggregate when removed from their protein lipid bilayer environment. This chapter is focused on describing a novel method for producing and solubilizing membrane proteins that can be easily adapted to high-throughput expression screening. This process is based on cell-free transcription and translation technology coupled with nanolipoprotein par ticles (NLPs), which are lipid bilayers confined within a ring of amphipathic protein of defined diameter. The NLPs act as a platform for inserting, solubilizing and characterizing functional membrane proteins. NLP component proteins (apolipoproteins), as well as membrane proteins can be produced by either traditional cell-based or as discussed here, cell-free expression methodologies.

An ion gating mechanism of gastric H,K-ATPase based on molecular dynamics simulations.

Gastric H,K-ATPase is an electroneutral transmembrane pump that moves protons from the cytoplasm of the parietal cell into the gastric lumen in exchange for potassium ions. The mechanism of transport against the established electrochemical gradients includes intermediate conformations in which the transferred ions are trapped (occluded) within the membrane domain of the pump. The pump cycle involves switching between the E1 and E2P states. Molecular dynamics simulations on homology models of the E2P and E1 states were performed to investigate the mechanism of K(+) movement in this enzyme. We performed separate E2P simulations with one K(+) in the luminal channel, one K(+) ion in the occlusion site, two K(+) ions in the occlusion site, and targeted molecular dynamics from E2P to E1 with two K(+) ions in the occlusion site. The models were inserted into a lipid bilayer system and were stable over the time course of the simulations, and K(+) ions in the channel moved to a consistent location near the center of the membrane domain, thus defining the occlusion site. The backbone carbonyl oxygen from residues 337 through 342 on the nonhelical turn of M4, as well as side-chain oxygen from E343, E795, and E820, participated in the ion occlusion. A single water molecule was stably bound between the two K(+) ions in the occlusion site, providing an additional ligand and partial shielding the positive charges from one another. Targeted molecular dynamics was used to transform the protein from the E2P to the E1 state (two K(+) ions to the cytoplasm). This simulation identified the separation of the water column in the entry channel as the likely gating mechanism on the luminal side. A hydrated exit channel also formed on the cytoplasmic side of the occlusion site during this simulation. Hence, water molecules became available to hydrate the ions. The movement of the M1M2 transmembrane segments, and the displacement of residues Q159, E160, Q110, and T152 during the conformational change, as well as the motions of E343 and L346, acted as the cytoplasmic-side gate.

Gaba receptor insecticide non-competitive antagonists may bind at allosteric modulator sites.

Results from several studies have shown that a series of chemically distinct insecticide compounds (picrotoxin, BIDN, TBPS, fipronil, lindane, EBOB, and alpha-endosulfan) affect GABA A receptor function. In this investigation, docking of this set of insecticides to the GABA receptor identified five potential binding sites. The lowest energy site was found within the base of the transmembrane bundle, interacting with M2 but not in the pore, and includes many of the residues previously experimentally implicated in insecticide binding. Many of the binding modes are suggestive of a non-competitive allosteric mechanism based on interruption of the channel gating mechanism rather than directly blocking the channel. The results also distinguished between isomers of hexachlorohexane (HCH), where gamma-HCH (lindane) binds more favorably than beta-HCH. The results suggest multiple sites for insecticide binding and may suggest further mutagenesis and labeling work to either confirm or rule out these findings.

Analysis of the gastric H,K ATPase for ion pathways and inhibitor binding sites.

New models of the gastric H,K ATPase in the E1K and E2P states are presented as the first structures of a K+ counter-transport P2-type ATPase exhibiting ion entry and exit paths. Homology modeling was first used to generate a starting conformation from the srCa ATPase E2P form (PDB code 1wpg) that contains bound MgADP. Energy minimization of the model showed a conserved adenosine site but nonconserved polyphosphate contacts compared to the srCa ATPase. Molecular dynamics was then employed to expand the luminal entry sufficiently to allow access of the rigid K+ competitive naphthyridine inhibitor, Byk99, to its binding site within the membrane domain. The new E2P model had increased separation between transmembrane segments M3 through M8, and addition of water in this space showed not only an inhibitor entry path to the luminal vestibule but also a channel leading to the ion binding site. Addition of K+ to the hydrated channel with molecular dynamics modeling of ion movement identified a pathway for K+ from the lumen to the ion binding site to give E2K. A K+ exit path to the cytoplasm operating during the normal catalytic cycle is also proposed on the basis of an E1K homology model derived from the E12Ca2+ form of the srCa ATPase (PDB code 1su4). Autodock analyses of the new E2P model now correctly discriminate between high- and low-affinity K+ competitive inhibitors. Finally, the expanded luminal vestibule of the E2P model explains high-affinity ouabain binding in a mutant of the H,K ATPase [Qiu et al. (2005) J. Biol. Chem. 280, 32349-32355].

Channel opening motion of alpha7 nicotinic acetylcholine receptor as suggested by normal mode analysis.

The gating motion of the human nicotinic acetylcholine receptor (nAChR) alpha7 was investigated with normal mode analysis (NMA) of two homology models. The first model, referred to as model I, was built from both the Lymnaea stagnalis acetylcholine binding protein (AChBP) and the transmembrane (TM) domain of the Torpedo marmorata nAChR. The second model, referred to as model C, was based solely on the recent electron microscopy structure of the T. marmorata nAChR. Despite structural differences, both models exhibit nearly identical patterns of flexibility and correlated motions. In addition, both models show a similar global twisting motion that may represent channel gating. The similar results obtained for the two models indicate that NMA is most sensitive to the contact topology of the structure rather than its finer detail. The major difference between the low-frequency motions sampled for the two models is that a symmetrical pore-breathing motion, favoring channel opening, is present as the second most dominant motion in model I, whilst largely absent from model C. The absence of this mode in model C can be attributed to its less symmetrical architecture. Finally, as a further goal of the present study, an approximate open channel model, consistent with many experimental findings, has been produced.

Membrane protein structure quality in molecular dynamics simulation.

Our goal was to assess the relationship between membrane protein quality, output from protein quality checkers and output from molecular dynamics (MD) simulations. Membrane transport proteins are essential for a wide range of cellular processes. Structural features of integral membrane proteins are still under-explored due to experimental limitations in structure determination. Computational techniques can be used to exploit biochemical and medium resolution structural data, as well as sequence homology to known structures, and enable us to explore the structure-function relationships in several transmembrane proteins. The quality of the models produced is vitally important to obtain reliable predictions. An examination of the relationship between model stability in molecular dynamics (MD) simulations derived from RMSD (root mean squared deviation) and structure quality assessment from various protein quality checkers was undertaken. The results were compared to membrane protein structures, solved at various resolution, by either X-ray or electron diffraction techniques. The checking programs could predict the potential success of MD in making functional conclusions. MD stability was shown to be a good indicator for the quality of structures. The quality was also shown to be dependent on the resolution at which the structures were determined.

Structure and function of two-pore-domain K+ channels: contributions from genetic model organisms.

K+ channels that possess two pore domains in each channel subunit are common in many animal tissues. Such channels are generated from large families of subunits and are implicated in several functions, including temperature sensation, responses to ischaemia, K+ homeostasis and setting the resting potential of the cell. Their activity can be modulated by polyunsaturated fatty acids, pH and oxygen, and some are candidate targets of volatile anaesthetics. However, despite their potential as targets for novel drugs for human health, comparatively little is known about the molecular basis of their diverse physiological and pharmacological properties. Genetic model organisms have considerable potential for improving our understanding of these channels. In this article, we review the contributions of some of these genetic model organisms to recent advances in our knowledge of two-pore-domain K+ channels.

A gating mechanism proposed from a simulation of a human alpha7 nicotinic acetylcholine receptor.

The nicotinic acetylcholine receptor is a well characterized ligand-gated ion channel, yet a proper description of the mechanisms involved in gating, opening, closing, ligand binding, and desensitization does not exist. Until recently, atomic-resolution structural information on the protein was limited, but with the production of the x-ray crystal structure of the Lymnea stagnalis acetylcholine binding protein and the EM image of the transmembrane domain of the torpedo electric ray nicotinic channel, we were provided with a window to examine the mechanism by which this channel operates. A 15-ns all-atom simulation of a homology model of the homomeric human alpha7 form of the receptor was conducted in a solvated palmitoyl-2-oleoyl-sn-glycerol-phosphatidylcholine bilayer and examined in detail. The receptor was unliganded. The structure undergoes a twist-to-close motion that correlates movements of the C loop in the ligand binding domain, via the beta10-strand that connects the two, with the 10 degrees rotation and inward movement of two nonadjacent subunits. The Cys loop appears to act as a stator around which the alpha-helical transmembrane domain can pivot and rotate relative to the rigid beta-sheet binding domain. The M2-M3 loop may have a role in controlling the extent or kinetics of these relative movements. All of this motion, along with essential dynamics analysis, is suggestive of the direction of larger motions involved in gating of the channel.

Homology modelling and molecular dynamics simulations: comparative studies of human aquaporin-1.

The structures of the mammalian water transport protein Aqp1 and of its bacterial homologue GlpF enables us to test whether homology models can be used to explore relationships between structure, dynamics and function in mammalian transport proteins. Molecular dynamics simulations (totalling almost 40 ns) were performed starting from: the X-ray structure of Aqp1; a homology model of Aqp1 based on the GlpF structure; and intermediate resolution structures of Aqp1 derived from electron microscopy. Comparisons of protein RMSDs vs. time suggest that the homology models are of comparable conformational stability to the X-ray structure, whereas the intermediate resolution structures exhibit significant conformation drift. For simulations based on the X-ray structure and on homology models, the flexibility profile vs. residue number correlates well with the crystallographic B-values for each residue. In the simulations based on intermediate resolution structures, mobility of the highly conserved NPA loops is substantially higher than in the simulations based on the X-ray structure or the homology models. Pore radius profiles remained relatively constant in the X-ray and homology model simulations but showed substantial fluctuations (reflecting the higher NPA loop mobility) in the intermediate resolution simulations. The orientation of the dipoles of water molecules within the pore is of key importance in maintaining low proton permeability through Aqp1. This property seems to be quite robust to the starting model used in the simulation. These simulations suggest that homology models based on bacterial homologues may be used to derive functionally relevant information on the structural dynamics of mammalian transport proteins.

Water in ion channels and pores--simulation studies.

The microscopic properties of water in narrow pores are relevant to the function of ion channels and related membrane transport proteins. The emergence of several high-resolution structures allows one to perform molecular dynamics simulation studies of water in such pores. Simulations of bundles of parallel alpha-helical peptides (e.g. alamethicin) have enabled development of methodologies and concepts appropriate to such investigations. In the narrow channels formed by such bundles, water molecules exhibit reduced rotational and translation motion. This reduction in water mobility may be a general property of narrow pores. We have used simplified channel models to explore the role of hydrophobicity/hydrophilicity in the entry of water into pores. Narrow pores with a hydrophobic lining, although physically open, may not admit water molecules, acting as a 'hydrophobic gate' that prevents water and ion permeation. Such a gate can be opened either by widening the pore or making its lining more polar. Simulations have been used to explore the behaviour of water in GlpF, a member of the aquaporin family of water pores, and OmpA, a bacterial outer membrane protein. Preliminary results suggest that a continuous water wire is not formed within the amphipathic GlpF pore. Simulations of OmpA, in which polar residues line the channel, indicate that a small conformational change in one of the channel lining side chains may open the channel. In summary, comparison of the behaviour of water in different narrow transmembrane pores suggests that an amphipathic pore is ideal for water permeation, and that either a highly hydrophobic pore lining or a charged pore-lining region can act as a gate.