Proton transport through water-filled carbon nanotubes

Structural Determinants of Proton Blockage in Aquaporins

Aquaporins are an important class of membrane channels selective for water and linear polyols but impermeable to ions, including protons. Recent computational studies have revealed that the relay of protons through the water-conduction pathway of aquaporin channels is opposed by a substantial free energy barrier peaking at the signature NPA motifs. Here, free-energy simulations and continuum electrostatic calculations are combined to examine the nature and the magnitude of the contribution of specific structural elements to proton blockage in the bacterial glycerol uptake facilitator, GlpF. Potential of mean-force profiles for both hop and turn steps of structural diffusion in the narrow pore are obtained for artificial variants of the GlpF channel in which coulombic interactions between the pore contents and conserved residues Asn68 and Asn203 at the NPA signature motifs, Arg206 at the selectivity filter, and the peptidic backbone of the two half-helices M3 and M7, which are arranged in head-to-head fashion around the NPA motifs, are turned off selectively. A comparison of these results with electrostatic energy profiles for the translocation of a probe cation throughout the water permeation pathway indicates that the free-energy profile for proton movement inside the narrow pore is dominated by static effects arising from the distribution of charged and polar groups of the channel, whereas dielectric effects contribute primarily to opposing the access of H+ to the pore mouths (desolvation penalty). The single most effective way to abolish the free-energy gradients opposing the movement of H+ around the NPA motif is to turn off the dipole moments of helices M3 and M7. Mutation of either of the two NPA Asn residues to Asp compensates for charge-dipole and dipole-dipole effects opposing the hop and turn steps of structural diffusion, respectively, and dramatically reduces the free energy barrier of proton translocation, suggesting that these single mutants could leak protons.

Nucleic acid transport through carbon nanotube membranes

We study the electrophoretic transport of single-stranded RNA molecules through 1.5-nm-wide pores of carbon nanotube membranes by molecular dynamics simulations. From approximately 170 individual RNA translocation events analyzed at full atomic resolution of solvent, membrane, and RNA, we identify key factors in membrane transport of biopolymers. RNA entry into the nanotube pores is controlled by conformational dynamics, and exit by hydrophobic attachment of RNA bases to the pores. Without electric field, RNA remains hydrophobically trapped in the membrane despite large entropic and energetic penalties for confining charged polymers inside nonpolar pores. Differences in RNA conformational flexibility and hydrophobicity result in sequence-dependent rates of translocation, a prerequisite for nanoscale separation devices.

What really prevents proton transport through aquaporin? Charge self-energy versus proton wire proposals

The nature of the control of water/proton selectivity in biological channels is a problem of a fundamental importance. Most studies of this issue have proposed that an interference with the orientational requirements of the so-called proton wire is the source of selectivity. The elucidation of the structures of aquaporins, which have evolved to prevent proton transfer (PT), provided a clear benchmark for exploring the selectivity problem. Previous simulations of this system have not examined, however, the actual issue of PT, but only considered the much simpler task of the transfer of water molecules. Here we take aquaporin as a benchmark and quantify the origin of the water/proton selectivity in this and related systems. This is done by evaluating in a consistent way the free energy profile for transferring a proton along the channel and relating this profile to the relevant PT rate constants. It is found that the water/proton selectivity is controlled by the change in solvation free energy upon moving the charged proton from water to the channel. The reason for the focus on the elegant concept of the proton wire and the related Grotthuss-type mechanism is also considered. It is concluded that these mechanisms are clearly important in cases with flat free energy surfaces (e.g., in bulk water, in gas phase water chains, and in infinitely long channels). However, in cases of biological channels, the actual PT mechanism is much less important than the energetics of transferring the proton charge from water to different regions in the channels.

On the influence of semirigid environments on proton transfer along molecular chains

The dynamics of proton transfer along ammonia chains (chemical composition N(x)H(+)(3x+1), x=2, 4, and 6) in a constraining environment is investigated by ab initio molecular dynamics simulations. A carbon nanotube of defined length and diameter is used as an idealized constraining environment such that the ammonia chain is forced to maintain its quasilinear geometry. It is found that, although the energetics of proton transport shows considerable energetic barriers, proton translocation along the wire is possible at finite temperature for all chain lengths studied. The proton transport involves rotational reorientation of the proton-carrying ammonia molecule. High level ab initio calculations (MP2/aug-cc-pVTZ) yield barriers for internal rotation of 9.1 kcal/mol for NH(4) (+)-NH(3) and 11.7 kcal/mol for OH(3) (+)-OH(2), respectively. The infrared spectrum calculated from the dipole-dipole autocorrelation function shows distinct spectral features in the regions (2000-3000 cm(-1)) where the NHN proton transfer mode is expected to absorb. Assigning moderate opposite total charges between 0.002 and 0.2e to the carbon atoms at the end caps of the nanotube leads to a considerable speedup of the proton transfer.

Electric field and temperature effects on water in the narrow nonpolar pores of carbon nanotubes

Water molecules in the narrow cylindrical pore of a (6,6) carbon nanotube form single-file chains with their dipoles collectively oriented either up or down along the tube axis. We study the interaction of such water chains with homogeneous electric fields for finite closed and infinite periodically replicated tubes. By evaluating the grand-canonical partition function term-by-term, we show that homogeneous electric fields favor the filling of previously empty nanotubes with water from the bulk phase. A two-state description of the collective water dipole orientation in the nanotube provides an excellent approximation for the dependence of the water-chain polarization and the filling equilibrium on the electric field. The energy and entropy contributions to the free energy of filling the nanotube were determined from the temperature dependence of the occupancy probabilities. We find that the energy of transfer depends sensitively on the water-tube interaction potential, and that the entropy of one-dimensionally ordered water chains is comparable to that of bulk water. We also discuss implications for proton transfer reactions in biology.

A coarse-grained model for a nanometer-scale molecular pump

A theory for a nanometer-scale pump based on the ratchet concept is presented. A lattice gas model with a set of moves that satisfy hydrodynamic equations is used to describe an asymmetric nanometer channel connecting two reservoirs of fluid. The channel, which is coupled to an external oscillatory (or stochastic) driving force, pumps fluid from one reservoir to the other. The frequency of the external driving force, the fluid density, and the channel dimensions are used to control the fluid flow. We observe a nonmonotonic behavior of the flow with respect to some model parameters and discuss the efficiency of the device.

The Mechanism of Proton Exclusion in the Aquaporin-1 Water Channel

Aquaporins are efficient, yet strictly selective water channels. Remarkably, proton permeation is fully blocked, in contrast to most other water-filled pores which are known to conduct protons well. Blocking of protons by aquaporins is essential to maintain the electrochemical gradient across cellular and subcellular membranes. We studied the mechanism of proton exclusion in aquaporin-1 by multiple non-equilibrium molecular dynamics simulations that also allow proton transfer reactions. From the simulations, an effective free energy profile for the proton motion along the channel was determined with a maximum-likelihood approach. The results indicate that the main barrier is not, as had previously been speculated, caused by the interruption of the hydrogen-bonded water chain, but rather by an electrostatic field centered around the fingerprint Asn-Pro-Ala (NPA) motif. Hydrogen bond interruption only forms a secondary barrier located at the ar/R constriction region. The calculated main barrier height of 25-30 kJ mol(-1) matches the barrier height for the passage of protons across pure lipid bilayers and, therefore, suffices to prevent major leakage of protons through aquaporins. Conventional molecular dynamics simulations additionally showed that negatively charged hydroxide ions are prevented from being trapped within the NPA region by two adjacent electrostatic barriers of opposite polarity.

Water-gated mechanism of proton translocation by cytochrome c oxidase

Cytochrome c oxidase is essential for aerobic life as a membrane-bound energy transducer. O(2) reduction at the haem a(3)-Cu(B) centre consumes electrons transferred via haem a from cytochrome c outside the membrane. Protons are taken up from the inside, both to form water and to be pumped across the membrane (M.K.F. Wikström, Nature 266 (1977) 271; M. Wikström, K. Krab, M. Saraste, Cytochrome Oxidase, A Synthesis, Academic Press, London, 1981 ). The resulting electrochemical proton gradient drives ATP synthesis (P. Mitchell, Chemiosmotic Coupling in Oxidative and Photosynthetic Phosphorylation, Glynn Research, Bodmin, UK, 1966 ). Here we present a molecular mechanism for proton pumping coupled to oxygen reduction that is based on the unique properties of water in hydrophobic cavities. An array of water molecules conducts protons from a conserved glutamic acid, either to the Delta-propionate of haem a(3) (pumping), or to haem a(3)-Cu(B) (water formation). Switching between these pathways is controlled by the redox-state-dependent electric field between haem a and haem a(3)-Cu(B), which determines the water-dipole orientation, and therefore the proton transfer direction. Proton transfer via the propionate provides a gate to O(2) reduction. This pumping mechanism explains the unique arrangement of the metal cofactors in the structure. It is consistent with the large body of biochemical data, and is shown to be plausible by molecular dynamics simulations.