Single-channel water permeabilities of Escherichia coli aquaporins AqpZ and GlpF

Structural determinants of water permeation through the sodium-galactose transporter vSGLT

Sodium-glucose transporters (SGLTs) facilitate the movement of water across the cell membrane, playing a central role in cellular homeostasis. Here, we present a detailed analysis of the mechanism of water permeation through the inward-facing state of vSGLT based on nearly 10 μs of molecular dynamics simulations. These simulations reveal the transient formation of a continuous water channel through the transporter that permits water to permeate the protein. Trajectories in which spontaneous release of galactose is observed, as well as those in which galactose remains in the binding site, show that the permeation rate, although modulated by substrate occupancy, is not tightly coupled to substrate release. Using a, to our knowledge, novel channel-detection algorithm, we identify the key residues that control water flow through the transporter and show that solvent gating is regulated by side-chain motions in a small number of residues on the extracellular face. A sequence alignment reveals the presence of two insertion sites in mammalian SGLTs that flank these outer-gate residues. We hypothesize that the absence of these sites in vSGLT may account for the high water permeability values for vSGLT determined via simulation compared to the lower experimental estimates for mammalian SGLT1.

Optimizing water permeability through the hourglass shape of aquaporins

The ubiquitous aquaporin channels are able to conduct water across cell membranes, combining the seemingly antagonist functions of a very high selectivity with a remarkable permeability. Whereas molecular details are obvious keys to perform these tasks, the overall efficiency of transport in such nanopores is also strongly limited by viscous dissipation arising at the connection between the nanoconstriction and the nearby bulk reservoirs. In this contribution, we focus on these so-called entrance effects and specifically examine whether the characteristic hourglass shape of aquaporins may arise from a geometrical optimum for such hydrodynamic dissipation. Using a combination of finite-element calculations and analytical modeling, we show that conical entrances with suitable opening angle can indeed provide a large increase of the overall channel permeability. Moreover, the optimal opening angles that maximize the permeability are found to compare well with the angles measured in a large variety of aquaporins. This suggests that the hourglass shape of aquaporins could be the result of a natural selection process toward optimal hydrodynamic transport. Finally, in a biomimetic perspective, these results provide guidelines to design artificial nanopores with optimal performances.

Allosteric mechanism of water-channel gating by Ca2+-calmodulin

Calmodulin (CaM) is a universal regulatory protein that communicates the presence of calcium to its molecular targets and correspondingly modulates their function. This key signaling protein is important for controlling the activity of hundreds of membrane channels and transporters. However, our understanding of the structural mechanisms driving CaM regulation of full-length membrane proteins has remained elusive. In this study, we determined the pseudo-atomic structure of full-length mammalian aquaporin-0 (AQP0, Bos Taurus) in complex with CaM using electron microscopy to understand how this signaling protein modulates water channel function. Molecular dynamics and functional mutation studies reveal how CaM binding inhibits AQP0 water permeability by allosterically closing the cytoplasmic gate of AQP0. Our mechanistic model provides new insight, only possible in the context of the fully assembled channel, into how CaM regulates multimeric channels by facilitating cooperativity between adjacent subunits.

Glycerol modulates water permeation through Escherichia coli aquaglyceroporin GlpF

Among aquaglyceroporins that transport both water and glycerol across the cell membrane, Escherichia coli glycerol uptake facilitator (GlpF) is the most thoroughly studied. However, one question remains: Does glycerol modulate water permeation? This study answers this fundamental question by determining the three-dimensional potential of mean force of glycerol along the permeation path through GlpF's conducting pore. There is a deep well near the Asn-Pro-Ala (NPA) motifs (6.5kcal/mol below the bulk level) and a barrier near the selectivity filter (10.1kcal/mol above the well bottom). This profile owes its existence to GlpF's perfect steric arrangement: The glycerol-protein van der Waals interactions are attractive near the NPA but repulsive elsewhere in the conducting pore. In light of the single-file nature of waters and glycerols lining up in GlpF's amphipathic pore, it leads to the following conclusion: Glycerol modulates water permeation in the μM range. At mM concentrations, GlpF is glycerol-saturated and a glycerol residing in the well occludes the conducting pore. Therefore, water permeation is fully correlated to glycerol dissociation that has an Arrhenius activation barrier of 6.5kcal/mol. Validation of this theory is based on the existent in vitro data, some of which have not been given the proper attention they deserved: The Arrhenius activation barriers were found to be 7kcal/mol for water permeation and 9.6kcal/mol for glycerol permeation; The presence of up to 100mM glycerol did not affect the kinetics of water transport with very low permeability, in apparent contradiction with the existent theories that predicted high permeability (0M glycerol).

Glycerol inhibits water permeation through Plasmodium falciparum aquaglyceroporin

Plasmodium falciparum aquaglyceroporin (PfAQP) is a multifunctional membrane protein in the plasma membrane of P. falciparum, the parasite that causes the most severe form of malaria. The current literature has established the science of PfAQP's structure, functions, and hydrogen-bonding interactions but left unanswered the following fundamental question: does glycerol modulate water permeation through aquaglyceroporin that conducts both glycerol and water? This paper provides an affirmative answer to this question of essential importance to the protein's functions. On the basis of the chemical-potential profile of glycerol from the extracellular bulk region, throughout PfAQP's conducting channel, to the cytoplasmic bulk region, this study shows the existence of a bound state of glycerol inside aquaglyceroporin's permeation pore, from which the dissociation constant is approximately 14μM. A glycerol molecule occupying the bound state occludes the conducting pore through which permeating molecules line up in single file by hydrogen-bonding with one another and with the luminal residues of aquaglyceroporin. In this way, glycerol inhibits permeation of water and other permeants through aquaglyceroporin. The biological implications of this theory are discussed and shown to agree with the existent in vitro data. It turns out that the structure of aquaglyceroporin is perfect for the van der Waals interactions between the protein and glycerol to cause the existence of the bound state deep inside the conducting pore and, thus to play an unexpected but significant role in aquaglyceroporin's functions.

Insights into the mechanisms of the selectivity filter of Escherichia coli aquaporin Z

Aquaporin Z (AQPZ) is a tetrameric protein that forms water channels in the cell membrane of Escherichia coli. The histidine residue (residue 174) in the selectivity filter (SF) region plays an important role in the transport of water across the membrane. In this work, we perform equilibrium molecular dynamics (MD) simulations to illustrate the gating mechanism of the SF and the influences of residue 174 in two different protonation states: Hsd174 with the proton at Nδ, and Hse174 with the proton at Nε. We calculate the pore radii in the SF region versus the simulation time. We perform steered MD to compute the free-energy profile, i.e., the potential of mean force (PMF) of a water molecule through the SF region. We conduct a quantum mechanics calculation of the binding energy of one water molecule with the residues in the SF region. The hydrogen bonds formed between the side chain of Hsd174 and the side chain of residue 189 (Arg189) play important roles in the selectivity mechanism of AQPZ. The radii of the pores, the hydrogen-bond analysis, and the free energies show that it is easier for water molecules to permeate through the SF region of AQPZ with residue 174 in the Hse state than in the Hsd state.

Preparation of supported lipid membranes for aquaporin Z incorporation

There has been a recent surge of interest to mimic the performance of natural cellular membranes by incorporating water channel proteins-aquaporins (AQPs) into various ultrathin films for water filtration applications. To make biomimetic membranes one of the most crucial steps is preparing a defect-free platform for AQPs incorporation on a suitable substrate. In this study two methods were used to prepare supported lipid membranes on NF membrane surfaces under a benign pH condition of 7.8. One method was direct vesicle fusion on a hydrophilic membrane NF-270; the other was vesicle fusion facilitated by hydraulic pressure on a modified hydrophilic NF-270 membrane whose surface has been spin-coated with positively charged lipids. Experiments revealed that the supported lipid membrane without AQPs prepared by the spin coating plus vesicle fusion had a much lower defect density than that prepared by vesicle fusion alone. It appears that the surface roughness and charge are the main factors determining the quality of the supported lipid membrane. Aquaporin Z (AqpZ) proteins were successfully incorporated into 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) liposomes and its permeability was measured by the stopped-flow experimental procedure. However, after the proteoliposomes have been fused onto the modified substrate, the AqpZ function in the resultant membrane was not observed and AFM images showed distinct aggregations of unfused proteoliposomes or AqpZ proteins on the substrate surface. It is speculated that the inhibition of AqpZ function may be caused by the low lipid mobility on the NF membrane surface. Further investigations to evaluate and optimize the structure-performance relationship are required.

Intracellular signaling by diffusion: can waves of hydrogen peroxide transmit intracellular information in plant cells?

Amplitude- and frequency-modulated waves of Ca(2+) ions transmit information inside cells. Reactive Oxygen Species (ROS), specifically hydrogen peroxide, have been proposed to have a similar role in plant cells. We consider the feasibility of such an intracellular communication system in view of the physical and biochemical conditions in plant cells. As model system, we use a H(2)O(2) signal originating at the plasma membrane (PM) and spreading through the cytosol. We consider two maximally simple types of signals, isolated pulses and harmonic oscillations. First we consider the basic limits on such signals as regards signal origin, frequency, amplitude, and distance. Then we establish the impact of ROS-removing enzymes on the ability of H(2)O(2) to transmit signals. Finally, we consider to what extent cytoplasmic streaming distorts signals. This modeling allows us to predict the conditions under which diffusion-mediated signaling is possible. We show that purely diffusive transmission of intracellular information by H(2)O(2) over a distance of 1 μm (typical distance between organelles, which may function as relay stations) is possible at frequencies well above 1 Hz, which is the highest frequency observed experimentally. This allows both frequency and amplitude modulation of the signal. Signaling over a distance of 10 μm (typical distance between the PM and the nucleus) may be possible, but requires high signal amplitudes or, equivalently, a very low detection threshold. Furthermore, at this longer distance a high rate of enzymatic degradation is required to make signaling at frequencies above 0.1 Hz possible. In either case, cytoplasmic streaming does not seriously disturb signals. We conclude that although purely diffusion-mediated signaling without relaying stations is theoretically possible, it is unlikely to work in practice, since it requires a much faster enzymatic degradation and a much lower cellular background concentration of H(2)O(2) than observed experimentally.

Molecular dynamics study of the archaeal aquaporin AqpM

Background

Aquaporins are a large family of transmembrane channel proteins that are present throughout all domains of life and are implicated in human disorders. These channels, allow the passive but selective movement of water and other small neutral solutes across cell membranes. Aquaporins have been classified into two sub-families: i) strict aquaporins that only allow the passage of water and ii) the less selective aquaglyceroporins that transport water and other neutral solutes, such as glycerol, CO₂ or urea. Recently, the identification and characterization of a number of archaeal and bacterial aquaporins suggested the existence of a third sub-family; one that is neither a strict aquaporin nor an aquaglyceroporin. The function and phylogeny of this third family is still a matter of debate.

Results

Twenty nanosecond molecular dynamics (MD) simulation of a fully hydrated tetramer of AqpM embedded in a lipid bilayer permitted predictions to be made of key biophysical parameters including: single channel osmotic permeability constant (p_f), single channel diffusive permeability constant (p_d), channel radius, potential water occupancy of the channel and water orientation inside the pore. These properties were compared with those of well characterized representatives of the two main aquaporin sub-families. Results show that changes in the amino acid composition of the aromatic/arginine region affect the size and polarity of the selectivity filter (SF) and could help explain the difference in water permeability between aquaporins. In addition, MD simulation results suggest that AqpM combines characteristics of strict aquaporins, such as the narrow SF and channel radius, with those of aquaglyceroporins, such as a more hydrophobic and less polar SF.

Conclusions

MD simulations of AqpM extend previous evidence that this archaeal aquaporin exhibits hybrid features intermediate between the two known aquaporin sub-families, supporting the idea that it may constitute a member of a novel class of aquaporins.

Water wires in atomistic models of the Hv1 proton channel

The voltage-gated proton channel (Hv1) is homologous to the voltage-sensing domain (VSD) of voltage-gated potassium (Kv) channels but lacks a separate pore domain. The Hv1 monomer has dual functions: it gates the proton current and also serves as the proton conduction pathway. To gain insight into the structure and dynamics of the yet unresolved proton permeation pathway, we performed all-atom molecular dynamics simulations of two different Hv1 homology models in a lipid bilayer in excess water. The structure of the Kv1.2-Kv2.1 paddle-chimera VSD was used as template to generate both models, but they differ in the sequence alignment of the S4 segment. In both models, we observe a water wire that extends through the membrane, whereas the corresponding region is dry in simulations of the Kv1.2-Kv2.1 paddle-chimera. We find that the kinetic stability of the water wire is dependent upon the identity and location of the residues lining the permeation pathway, in particular, the S4 arginines. A measurement of water transport kinetics indicates that the water wire is a relatively static feature of the permeation pathway. Taken together, our results suggest that proton conduction in Hv1 may occur via Grotthuss hopping along a robust water wire, with exchange of water molecules between inner and outer ends of the permeation pathway minimized by specific water-protein interactions. This article is part of a Special Issue entitled: Membrane protein structure and function.

Mercury inhibits the L170C mutant of aquaporin Z by making waters clog the water channel

We conduct in silico experiments of the L170C mutant of the Escherichia coli aquaporin Z (AQPZ) with and without mercury bonded to residue Cys 170. We find that bonding mercury to Cys 170 does not induce consequential structural changes to the protein. We further find that mercury does not stick in the middle of the water channel to simply occlude water permeation, but resides on the wall of the water pore. However, we observe that the water permeation coefficient of L170C-Hg(+) (with one mercury ion bonded to Cys 170) is approximately half of that of the mercury-free L170C. We examine the interactions between the mercury ion and the waters in its vicinity and find that five to six waters are strongly attracted by the mercury ion, occluding the space of the water channel. Therefore we conclude that mercury, at low concentration, inhibits AQPZ-L170C mutant by making water molecules clog the water channel.

Exploring transmembrane diffusion pathways with molecular dynamics

Transmembrane exchange of materials is a fundamental process in biology. Molecular dynamics provides a powerful method to investigate in great detail various aspects of the phenomenon, particularly the permeation of small uncharged molecules, which continues to pose a challenge to experimental studies. We will discuss some of the recent simulation studies investigating the role of lipid-mediated and protein-mediated mechanisms in permeation of water and gas molecules across the membrane.

Structural context shapes the aquaporin selectivity filter

Aquaporins are transmembrane channels that facilitate the permeation of water and small, uncharged amphipathic molecules across cellular membranes. One distinct aquaporin subfamily contains pure water channels, whereas a second subfamily contains channels that conduct small alditols such as glycerol, in addition to water. Distinction between these substrates is central to aquaporin function, though the contributions of protein structural motifs required for selectivity are not yet fully characterized. To address this question, we sequentially engineered three signature amino acids of the glycerol-conducting subfamily into the Escherichia coli water channel aquaporin Z (AqpZ). Functional analysis of these mutant channels showed a decrease in water permeability but not the expected increase in glycerol conduction. Using X-ray crystallography, we determined the atomic resolution structures of the mutant channels. The structures revealed a channel surprisingly similar in size to the wild-type AqpZ pore. Comparison with measured rates of transport showed that, as the size of the selectivity filter region of the channel approaches that of water, channel hydrophilicity dominated water conduction energetics. In contrast, the major determinant of selectivity for larger amphipathic molecules such as glycerol was channel cross-section size. Finally, we find that, although the selectivity filter region is indeed central to substrate transport, other structural elements that do not directly interact with the substrates, such as the loop connecting helices M6 and M7, and the C loop between helices C4 and C5, play an essential role in facilitating selectivity.

Plant aquaporin selectivity: where transport assays, computer simulations and physiology meet

Plants contain a large number of aquaporins with different selectivity. These channels generally conduct water, but some additionally conduct NH(3), CO(2) and/or H(2)O(2). The experimental evidence and molecular basis for the transport of a given solute, the validation with molecular dynamics simulations and the physiological impact of the selectivity are reviewed here. The aromatic/arginine (ar/R) constriction is most important for solute selection, but the exact pore requirements for efficient conduction of small solutes remain difficult to predict. Yeast growth assays are valuable for screening substrate selectivity and are explicitly shown for hydrogen peroxide and methylamine, a transport analog of ammonia. Independent assays need to address the relevance of different substrates for each channel in its physiological context. This is emphasized by the fact that several plant NIP channels, which conduct several solutes, are specifically involved in the transport of metalloids, such as silicic acid, arsenite, or boric acid in planta.

Cell-free synthesis of functional aquaporin Z in synthetic liposomes

The challenges involved in producing sufficient quantities of aquaporins for precise biophysical characterization have limited our knowledge of this important class of molecules. This article describes a cell-free protein synthesis method for producing high concentrations of the E. coli water transporter, aquaporin Z (AqpZ), in synthetic liposomes. To our knowledge, this is the first report of in vitro synthesis of a membrane protein directly into synthetic liposomes with verified function, (i.e., transport activity and selectivity). Titration of DOPC lipid vesicles added to the cell-free reaction show that production yields of active AqpZ are dependent on the concentration of DOPC lipid vesicles added to the cell-free reaction, with 224 +/- 24 lipids required per aquaporin monomer. Supplementation of the signal recognition particle receptor (FtsY) to the cell-free reaction increases production of vesicle-associated AqpZ but not active AqpZ. Cell-free reactions using 7 mg/mL lipids that were not supplemented with FtsY produced 507 +/- 11 microg/mL of vesicle-associated AqpZ that exhibited a specific water transport activity of (2.2 +/- 0.3) x 10(-14) cm(3) s(-1) monomer(-1). Proteinase K protection, activation energy determination, and selectivity against glycerol and urea transport also confirmed the production of correctly folded AqpZ. This technique is capable of producing milligram quantities of aquaporin that can be readily assayed for function, facilitating biophysical characterization and high-throughput analysis.

Determinants of water permeability through nanoscopic hydrophilic channels

Naturally occurring pores show a variety of polarities and sizes that are presumably directly linked to their biological function. Many biological channels are selective toward permeants similar or smaller in size than water molecules, and therefore their pores operate in the regime of single-file water pores. Intrinsic factors affecting water permeability through such pores include the channel-membrane match, the structural stability of the channel, the channel geometry and channel-water affinity. We present an extensive molecular dynamics study on the role of the channel geometry and polarity on the water osmotic and diffusive permeability coefficients. We show that the polarity of the naturally occurring peptidic channels is close to optimal for water permeation, and that the water mobility for a wide range of channel polarities is essentially length independent. By systematically varying the geometry and polarity of model hydrophilic pores, based on the fold of gramicidin A, the water density, occupancy, and permeability are studied. Our focus is on the characterization of the transition between different permeation regimes in terms of the structure of water in the pores, the average pore occupancy and the dynamics of the permeating water molecules. We show that a general relationship between osmotic and diffusive water permeability coefficients in the single-file regime accounts for the time averaged pore occupancy, and that the dynamics of the permeating water molecules through narrow non single file channels effectively behaves like independent single-file columns.

Dynamics and energetics of permeation through aquaporins. What do we learn from molecular dynamics simulations?

Aquaporins (AQPs) are a family of integral membrane proteins, which facilitate the rapid and yet highly selective flux of water and other small solutes across biological membranes. Molecular dynamics (MD) simulations contributed substantially to the understanding of the molecular mechanisms that underlie this remarkable efficiency and selectivity of aquaporin channels. This chapter reviews the current state of MD simulations of aquaporins and related aquaglyceroporins as well as the insights these simulations have provided. The mechanism of water permeation through AQPs and methods to determine channel permeabilities from simulations are described. Protons are strictly excluded from AQPs by a large electrostatic barrier and not by an interruption of the Grotthuss mechanism inside the pore. Both the protein's electric field and desolvation effects contribute to this barrier. Permeation of apolar gas molecules such as CO(2) through AQPs is accompanied by a large energetic barrier and thus can only be expected in membranes with a low intrinsic gas permeability. Additionally, the insights from simulations into the mechanism of glycerol permeation through the glycerol facilitator GlpF from E. coli are summarized. Finally, MD simulations are discussed that revealed that the aro-matic/arginine constriction region is generally the filter for uncharged solutes, and that AQP selectivity is controlled by a hydrophobic effect and steric restraints.

The roles of pore ring and plug in the SecY protein-conducting channel

The protein-conducting channel, or translocon, is an evolutionarily conserved complex that allows nascent proteins to cross a cellular membrane or integrate into it. The crystal structure of an archaeal translocon, the SecY complex, revealed that two elements contribute to sealing the channel: a small "plug" domain blocking the periplasmic region of the channel, and a pore ring composed of six hydrophobic residues acting as a constriction point at the channel's center. To determine the independent functions of these two elements, we have performed molecular dynamics simulations of the native channel as well as of two recently structurally resolved mutants in which portions of their plugs were deleted. We find that in the mutants, the instability in the plug region leads to a concomitant increase in flexibility of the pore ring. The instability is quantified by the rate of water permeation in each system as well as by the force required for oligopeptide translocation. Through a novel simulation in which the interactions between the plug and water were independently controlled, we find that the role of the plug in stabilizing the pore ring is significantly more important than its role as a purely steric barrier.

Dynamic control of slow water transport by aquaporin 0: implications for hydration and junction stability in the eye lens

Aquaporin 0 (AQP0), the most abundant membrane protein in mammalian lens fiber cells, not only serves as the primary water channel in this tissue but also appears to mediate the formation of thin junctions between fiber cells. AQP0 is remarkably less water permeable than other aquaporins, but the structural basis and biological significance of this low permeability remain uncertain, as does the permeability of the protein in a reported junctional form. To address these issues, we performed molecular dynamics (MD) simulations of water transport through membrane-embedded AQP0 in both its (octameric) junctional and (tetrameric) nonjunctional forms. From our simulations, we measured an osmotic permeability for the nonjunctional form that agrees with experiment and found that the distinct dynamics of the conserved, lumen-protruding side chains of Tyr-23 and Tyr-149 modulate water passage, accounting for the slow permeation. The junctional and nonjunctional forms conducted water equivalently, in contrast to a previous suggestion based on static crystal structures that water conduction is lost on junction formation. Our analysis suggests that the low water permeability of AQP0 may help maintain the mechanical stability of the junction. We hypothesize that the structural features leading to low permeability may have evolved in part to allow AQP0 to form junctions that both conduct water and contribute to the organizational structure of the fiber cell tissue and microcirculation within it, as required to maintain transparency of the lens.

The impact of peptides on lipid membranes

We review the fundamental strategies used by small peptides to associate with lipid membranes and how the different strategies impact on the structure and dynamics of the lipids. In particular we focus on the binding of amphiphilic peptides by electrostatic and hydrophobic forces, on the anchoring of peptides to the bilayer by acylation and prenylation, and on the incorporation of small peptides that form well-defined channels. The effect of lipid-peptide interactions on the lipids is characterized in terms of lipid acyl-chain order, membrane thickness, membrane elasticity, permeability, lipid-domain and annulus formation, as well as acyl-chain dynamics. The different situations are illustrated by specific cases for which experimental observations can be interpreted and supplemented by theoretical modeling and simulations. A comparison is made with the effect on lipids of trans-membrane proteins. The various cases are discussed in the context of the possible roles played by lipid-peptide interactions for the biological, physiological, and pharmacological function of peptides.

Water transport in aquaporins: osmotic permeability matrix analysis of molecular dynamics simulations

Single-channel osmotic water permeability (p(f)) is a key quantity for investigating the transport capability of the water channel protein, aquaporin. However, the direct connection between the single scalar quantity p(f) and the channel structure remains unclear. In this study, based on molecular dynamics simulations, we propose a p(f)-matrix method, in which p(f) is decomposed into contributions from each local region of the channel. Diagonal elements of the p(f) matrix are equivalent to the local permeability at each region of the channel, and off-diagonal elements represent correlated motions of water molecules in different regions. Averaging both diagonal and off-diagonal elements of the p(f) matrix recovers p(f) for the entire channel; this implies that correlated motions between distantly-separated water molecules, as well as adjacent water molecules, influence the osmotic permeability. The p(f) matrices from molecular dynamics simulations of five aquaporins (AQP0, AQP1, AQP4, AqpZ, and GlpF) indicated that the reduction in the water correlation across the Asn-Pro-Ala region, and the small local permeability around the ar/R region, characterize the transport efficiency of water. These structural determinants in water permeation were confirmed in molecular dynamics simulations of three mutants of AqpZ, which mimic AQP1.

Osmosis in semi-permeable pores: an examination of the basic flow equations based on an experimental and molecular dynamics study

Classically ‘semi-permeable’ pores are generally considered to mediate osmotic flow at a rate dependent upon the hydraulic conductance of the pore and the difference in water potential. The shape or size of the solute molecules is not considered to exert a first-order effect on the flow rate nor is the hydraulic conductance thought to be solute dependent. By the experimental measurement of osmosis in the biological pore AQP (aquaporin) and hard-sphere molecular dynamics simulation of a model pore, we show here that the solute radius can have a profound effect on the osmotic flow rate, causing it to decline steeply with decreasing solute radius.

Using a simple non-equilibrium thermodynamic theory, we propose that an additional ‘osmotic flow coefficient’ is required to describe flows in semi-permeable structures such as AQPs, and that the fall in flow rate with radius represents a conversion from hydraulic to diffusive water flow due to increasing penetration of the pore by the solute. The interaction between the pore geometry and the solute size cannot, therefore, be overlooked, although for every solute the system obeys the criterion for semi-permeability required by basic thermodynamics. The osmotic pore theory therefore reveals a novel and potentially rich structure that remains to be explored in full.

Does CO2 permeate through aquaporin-1?

Aquaporins facilitate water permeation across biological membranes. Additionally, glycerol and other small neutral solutes are permeated by related aquaglyceroporins. The role of aquaporins in gas permeation has been a long-standing and controversially discussed issue. We present an extensive set of atomistic molecular dynamics simulations that address the question of CO(2) permeation through human aquaporin-1. Free energy profiles derived from the simulations display a barrier of approximately 23 kJ/mol in the aromatic/arginine constriction region of the water pore, whereas a barrier of approximately 4 kJ/mol was observed for a palmitoyloleoylphosphatidylethanolamine lipid bilayer membrane. The results indicate that significant aquaporin-1-mediated CO(2) permeation is to be expected only in membranes with a low intrinsic CO(2) permeability.