Water, proton, and ion transport: from nanotubes to proteins

Multi-Template-Guided Synthesis of High-Dimensional Molecular Assemblies for Humidity Gradient-Based Power Generators

Systematically orchestrating fundamental building blocks into intricate high-dimensional molecular assemblies at molecular level is imperative for multifunctionality integration. However, this remains a formidable task in crystal engineering due to the dynamic nature of inorganic building blocks. Herein, we develop a multi-template-guided strategy to control building blocks. The coordination modes of ligands and the spatial hindrance of anionic templates are pivotal in dictating the overall structures. Flexible multi-dentate linkers selectively promote the formation of oligomeric assembly ([TeO3(Mo2O2S2)3O2(OH)(C5O2H7)3]4- {TeMo6}) into tetrahedral cages ([(TeO3)4(Mo2O2S2)12(OH)12(C9H9O4P)6]8- {Te4Mo24} and [(AsO4)4(Mo2O2S2)12(OH)12(C9H9O6)4]12- {As4Mo24}), while steric hindrance from anionic templates further assists in assembling cages into an open quadruply twisted Möbius nanobelt ([(C6H5O3P)8(Mo2O2S2)24(OH)24(C8H10O4)12]16- {P8Mo48}). Among these structures, the hydrophilic-hydrophobic hybrid cage {Te4Mo24} emerges as an exemplary molecular model for proton conduction and serves as a prototype for humidity gradient-based power generators (HGPGs). The Te4Mo24-PVDF-based HGPG (PVDF=Poly(vinylidene fluoride)) exhibits notable stability and power generation, yielding an open-circuit voltage of 0.51 V and a current density of 77.8 nA cm-2 at room temperature and 90 % relative humidity (RH). Further insights into the interactions between water molecules and microscale molecules within the generator are achieved through molecular dynamics simulations. This endeavor unveils a universal strategy for synthesizing multifunctional integration molecules.

Investigation of molecular details of a bacterial cationic amino acid transporter (GkApcT) during arginine transportation using molecular dynamics simulation and umbrella sampling techniques

CONTEXT

Cationic amino acid transporters (CATs) facilitate arginine transport across membranes and maintain its levels in various tissues and organs, but their overexpression has been associated with severe cancers. A recent study identified the alternating access mechanism and critical residues involved in arginine transportation in a cationic amino acid transporter from Geobacillus kaustophilus (GkApcT). Here, we used molecular dynamics (MD) simulation methods to investigate the transportation mechanism of arginine (Arg) through GkApcT. The results revealed that arginine strongly interacts with specific binding site residues (Thr43, Asp111, Glu115, Lys191, Phe231, Ile234, and Asp237). Based on the umbrella sampling, the main driving force for arginine transport is the polar interactions of the arginine with channel-lining residues. An in-depth description of the dissociation mechanism and binding energy analysis brings valuable insight into the interactions between arginine and transporter residues, facilitating the design of effective CAT inhibitors in cancer cells.

METHODS

The membrane-protein system was constructed by uploading the prokaryotic CAT (PDB ID: 6F34) to the CHARMM-GUI web server. Molecular dynamics simulations were done using the GROMACS package, version 5.1.4, with the CHARMM36 force field and TIP3P water model. The MM-PBSA approach was performed for determining the arginine binding free energy. Furthermore, the hotspot residues were identified through per-residue decomposition analysis. The characteristics of the channel such as bottleneck radius and channel length were analyzed using the CaverWeb 1.1 web server. The proton wire inside the transporter was investigated based on the classic Grotthuss mechanism. We also investigated the atomistic details of arginine transportation using the path-based free energy umbrella sampling technique (US).

Arginine transportation mechanism through cationic amino acid transporter 1: insights from molecular dynamics studies

Metabolic and signaling mechanisms in mammalian cells are facilitated by the transportation of L-arginine (Arg) across the plasma membrane through cationic amino acid transporter (CAT) proteins. Due to a lack of argininosuccinate synthase (ASS) activity in various tumor cells such as acute myeloid leukemia, acute lymphocytic leukemia, and chronic lymphocytic leukemia, these tumor entities are arginine-auxotrophic and therefore depend on the uptake of the amino acid arginine. Cationic amino acid transporter-1 (CAT-1) is the leading arginine importer expressed in the aforementioned tumor entities. Hence, in the present study, to investigate the transportation mechanism of arginine in CAT-1, we performed molecular dynamics (MD) simulation methods on the modeled human CAT-1. The MM-PBSA approach was conducted to determine the critical residues interacting with arginine within the corresponding binding site of CAT-1. In addition, we found out that the water molecules have the leading role in forming the transportation channel within CAT-1. The conductive structure of CAT-1 was formed only when the water molecules were continuously distributed across the channel. Steered molecular dynamics (SMD) simulation approach showed various energy barriers against arginine transportation through CAT-1, especially while crossing the bottlenecks of the related channel. These findings at the molecular level might shed light on identifying the crucial amino acids in the binding of arginine to eukaryotic CATs and also provide fundamental insights into the arginine transportation mechanisms through CAT-1. Understanding the transportation mechanism of arginine is essential to developing CAT-1 blockers, which can be potential medications for some types of cancers.Communicated by Ramaswamy H. Sarma.

Reverse Translocation of Nucleotides through a Carbon Nanotube

We report molecular dynamics (MD) simulations of the reverse translocation of single nucleotides through a narrow carbon nanotube (CNT), with a diameter of 1.36 nm, immersed in a 1 M KCl electrolyte solution under an applied electric field along the tube axis. We observe ion selectivity by the narrow CNT, which leads to a high flow of K⁺ ions, in contrast to a negligible and opposing current of Cl^- ions. The K⁺ ions, driven by the electric field, force a negatively charged single nucleotide into the narrow CNT where it is trapped by the incoming K⁺ ions and water molecules, and the nucleotide is driven in the same direction as the K⁺ ions. This illustrates a novel mechanism of nucleotide reverse translocation that is controlled by ion selectivity. An increase in the CNT diameter to 2.71 nm or an increase in nucleotide chain length both lead to translocation in the normal direction of the applied field. The reverse translocation rate of single nucleotides is correlated to the ionic current of K⁺ ions in the narrow tube, unlike translocation in the normal direction in the wider tube.

Molecular simulations of water adsorption and transport in mesopores with varying hydrophilicity arrangements

The adsorption and transport of water in an open cylindrical mesopore with two different inner surface arrangements of hydrophilicities were examined by molecular simulations. The first model has a weak hydrophilic surface at both entrances of the pore and a stronger hydrophilic surface in the mid-section. The second pore has stronger hydrophilic surfaces at the entrances and weaker in the middle region. The simulation results show that the water adsorption isotherms obtained from Grand Canonical Monte Carlo simulations and pore filling curves acquired from Grand Canonical Molecular Dynamics simulations change depending on the arrangement of the strong and weak hydrophilic surfaces. In the first model, water condensation focuses on the mid-section forming a liquid bridge or a film, which creates a concave meniscus accelerating subsequent adsorption within the pore. Two bridges form in the entrance regions, where a cavity naturally occurs in between the films, in the second model. The different filling and emptying mechanisms clearly change the adsorption-desorption characteristics for the two pore types, but the second type generally showed faster transitions overall. Flux and meniscus analysis also reveals a circulating flow at the menisci of the interfaces within the pore. The results are expected to be valuable in understanding the effects of interior surface modification of nanopores in future applications.

Retinal isomerization and water-pore formation in channelrhodopsin-2

Channelrhodopsin-2 (ChR2) is a light-sensitive ion channel widely used in optogenetics. Photoactivation triggers a trans-to-cis isomerization of a covalently bound retinal. Ensuing conformational changes open a cation-selective channel. We explore the structural dynamics in the early photocycle leading to channel opening by classical (MM) and quantum mechanical (QM) molecular simulations. With QM/MM simulations, we generated a protein-adapted force field for the retinal chromophore, which we validated against absorption spectra. In a 4-µs MM simulation of a dark-adapted ChR2 dimer, water entered the vestibules of the closed channel. Retinal all-trans to 13-cis isomerization, simulated with metadynamics, triggered a major restructuring of the charge cluster forming the channel gate. On a microsecond time scale, water penetrated the gate to form a membrane-spanning preopen pore between helices H1, H2, H3, and H7. This influx of water into an ion-impermeable preopen pore is consistent with time-resolved infrared spectroscopy and electrophysiology experiments. In the retinal 13-cis state, D253 emerged as the proton acceptor of the Schiff base. Upon proton transfer from the Schiff base to D253, modeled by QM/MM simulations, we obtained an early-M/P2390-like intermediate. Rapid rotation of the unprotonated Schiff base toward the cytosolic side effectively prevents its reprotonation from the extracellular side. From MM and QM simulations, we gained detailed insight into the mechanism of ChR2 photoactivation and early events in pore formation. By rearranging the network of charges and hydrogen bonds forming the gate, water emerges as a key player in light-driven ChR2 channel opening.

Correlated Rectification Transport in Ultranarrow Charged Nanocones

Using molecular dynamics simulations, we reveal ion rectification in charged nanocones with exit diameters of 1-2 nm. The simulations exhibit an opposite rectification current direction than experiments performed in conical channels with exit diameters larger than 5 nm. This can be understood by the fact that in ultranarrow charged cones screening ions are trapped close to the cone tip at both field directions, which necessitates them to be released from the cone in a correlated multi-ion fashion. Electroosmosis induced by a unidirectional ion flow is also observed.

Structural characteristics of hydrated protons in the conductive channels: effects of confinement and fluorination studied by molecular dynamics simulation

The relationship between the proton conductive channel and the hydrated proton structure is of significant importance for understanding the deformed hydrogen bonding network of the confined protons which matches the nanochannel. In general, the structure of hydrated protons in the nanochannel of the proton exchange membrane is affected by several factors. To investigate the independent effect of each factor, it is necessary to eliminate the interference of other factors. In this paper, a one-dimensional carbon nanotube decorated with fluorine was built to investigate the independent effects of nanoscale confinement and fluorination on the structural properties of hydrated protons in the nanochannel using classical molecular dynamics simulation. In order to characterize the structure of hydrated protons confined in the channel, the hydrogen bonding interaction between water and the hydrated protons has been studied according to suitable hydrogen bond criteria. The hydrogen bond criteria were proposed based on the radial distribution function, angle distribution and pair-potential energy distribution. It was found that fluorination leads to an ordered hydrogen bonding structure of the hydrated protons near the channel surface, and confinement weakens the formation of the bifurcated hydrogen bonds in the radial direction. Besides, fluorination lowers the free energy barrier of hydronium along the nanochannel, but slightly increases the barrier for water. This leads to disintegration of the sequential hydrogen bond network in the fluorinated CNTs with small size. In the fluorinated CNTs with large diameter, the lower degree of confinement produces a spiral-like sequential hydrogen bond network with few bifurcated hydrogen bonds in the central region. This structure might promote unidirectional proton transfer along the channel without random movement. This study provides the cooperative effect of confinement dimension and fluorination on the structure and hydrogen bonding of the slightly acidic water in the nanoscale channel.

Molecular simulation and modeling of complex I

Molecular modeling and molecular dynamics simulations play an important role in the functional characterization of complex I. With its large size and complicated function, linking quinone reduction to proton pumping across a membrane, complex I poses unique modeling challenges. Nonetheless, simulations have already helped in the identification of possible proton transfer pathways. Simulations have also shed light on the coupling between electron and proton transfer, thus pointing the way in the search for the mechanistic principles underlying the proton pump. In addition to reviewing what has already been achieved in complex I modeling, we aim here to identify pressing issues and to provide guidance for future research to harness the power of modeling in the functional characterization of complex I. This article is part of a Special Issue entitled Respiratory complex I, edited by Volker Zickermann and Ulrich Brandt.

Nanofluidic Transport through Isolated Carbon Nanotube Channels: Advances, Controversies, and Challenges

Owing to their simple chemistry and structure, controllable geometry, and a plethora of unusual yet exciting transport properties, carbon nanotubes (CNTs) have emerged as exceptional channels for fundamental nanofluidic studies, as well as building blocks for future fluidic devices that can outperform current technology in many applications. Leveraging the unique fluidic properties of CNTs in advanced systems requires a full understanding of their physical origin. Recent advancements in nanofabrication technology enable nanofluidic devices to be built with a single, nanometer-wide CNT as a fluidic pathway. These novel platforms with isolated CNT nanochannels offer distinct advantages for establishing quantitative structure-transport correlations in comparison with membranes containing many CNT pores. In addition, they are promising components for single-molecule sensors as well as for building nanotube-based circuits wherein fluidics and electronics can be coupled. With such advanced device architecture, molecular and ionic transport can be manipulated with vastly enhanced control for applications in sensing, separation, detection, and therapeutic delivery. Recent achievements in fabricating isolated-CNT nanofluidic platforms are highlighted, along with the most-significant findings each platform enables for water, ion, and molecular transport. The implications of these findings and remaining open questions on the exceptional fluidic properties of CNTs are also discussed.

Non-selective ion channel activity of polymorphic human islet amyloid polypeptide (amylin) double channels

Fundamental understanding of ion channel formation by amyloid peptides, which is strongly linked to cell toxicity, is very critical for (pre)clinical treatment of neurodegenerative diseases. Here, we combine atomistic simulations and experiments to demonstrate a broad range of conformational states of hIAPP double channels in lipid membranes. All individual channels display high selectivity for Cl(-) ions over cations, but the co-existence of polymorphic double channels of different conformations and orientations with different populations determines the non-ionic selectivity nature of the channels, which is different from the typical amyloid-β channels that exhibit Ca(2+) selective ion-permeable characteristics. This work provides a more complete physicochemical mechanism of amyloid-channel-induced toxicity.

Diffusion of water inside carbon nanotubes studied by pulsed field gradient NMR spectroscopy

Diffusion dynamics of guest molecules in nanopores has been studied intensively because diffusion is center on a number of research fields such as separation, drug delivery, chemical reactions, and sensing. In the present work, we report an experimental investigation of the self-diffusion of water inside carbon nanotube (CNT) channels using a pulsed field gradient (PFG) NMR method. The dispersion of CNTs homogeneously in water and cooling to temperatures below the melting point of bulk water allow us to probe the translational motion of confined water molecules. The results demonstrate that the self-diffusion coefficient of water in CNTs is highly dependent on the diffusion time and CNT diameter. In particular, the diffusivity of water in double-walled carbon nanotubes (DWNTs) with an average inner diameter of 2.3 ± 0.3 nm is twice that in multiwalled carbon nanotubes (MWNTs) with an average inner diameter of 6.7 ± 0.8 nm in the temperature range of 263-223 K. In addition, the effective self-diffusion coefficient in DWNTs is 1 order of magnitude higher than that reported for mesoporous silica materials with a similar pore size. The faster diffusivity of water in CNTs could be attributed to the ordered hydrogen bonds formed between water molecules within the confined channels of CNTs and the weak interaction between water and the CNT walls.

The open gate of the K(V)1.2 channel: quantum calculations show the key role of hydration

The open gate of the Kv1.2 voltage-gated potassium channel can just hold a hydrated K(+) ion. Quantum calculations starting from the x-ray coordinates of the channel confirm this, showing little change from the x-ray coordinates for the protein. Water molecules not in the x-ray coordinates, and the ion itself, are placed by the calculation. The water molecules, including their orientation and hydrogen bonding, with and without an ion, are critical for the path of the ion, from the solution to the gate. A sequence of steps is postulated in which the potential experienced by the ion in the pore is influenced by the position of the ion. The gate structure, with and without the ion, has been optimized. The charges on the atoms and bond lengths have been calculated using natural bond orbital calculations, giving K(+) ~0.77 charges, rather than 1.0. The PVPV hinge sequence has been mutated in silico to PVVV (P407V in the 2A79 numbering). The water structure around the ion becomes discontinuous, separated into two sections, above and below the ion. PVPV conservation closely relates to maintaining the water structure. Finally, these results have implications concerning gating.

Hydroxide ion can move faster than an excess proton through one-dimensional water chains in hydrophobic narrow pores

Carbon nanotubes (CNT) are known to provide a hydrophobic, confined environment for water where its structure and dynamics can be very different from those of bulk water. In particular, narrow CNTs of the type (6,6) allow only a single one-dimensional (1D) chain of water molecules inside them, thus providing an idealized scenario to study motion in 1D along water chains. In the present study, we have investigated structural and dynamic behavior of water and also of an excess proton and hydroxide ion in water-filled narrow CNTs by means of ab initio molecular dynamics and combined quantum-classical simulations. The main focus of the present work is on the molecular mechanism and kinetics of hydronium and hydroxide ion migration along 1D water chains of different lengths in confinement. It is found that the hydrogen-bonded structures of water and the excess proton and hydroxide ion in CNTs can be very different from those in bulk, and these altered solvation structures play critical roles in determining the proton-transfer (PT) rates along water chains. For the present 1D chain systems, the hydroxide ion is found to migrate at a slightly faster rate than the excess proton, unlike their relative mobilities in bulk water. This faster migration of the hydroxide ion is found not only in CNTs with periodicity along the tube axis but also in isolated CNTs where the excess proton and the hydroxide ion are allowed to move under the influence of an electric field of an oppositely charged ion. The roles of rotational jumps and hydrogen-bond fluctuations in the PT events are discussed. In addition, the significance of hydrogen-bonding defects on the dynamics of an excess proton and hydroxide ion is also discussed for varying chain lengths.

Confinement by carbon nanotubes drastically alters the boiling and critical behavior of water droplets

Vapor pressure grows rapidly above the boiling temperature, and past the critical point liquid droplets disintegrate. Our atomistic simulations show that this sequence of events is reversed inside carbon nanotubes (CNT). Droplets disintegrate first and at low temperature, while pressure remains low. The droplet disintegration temperature is independent of the CNT diameter. In contrast, depending on CNT diameter, a temperature that is much higher than the bulk boiling temperature is required to raise the internal pressure. The control over pressure by CNT size can be useful for therapeutic drug delivery.

Water transport and purification in nanochannels controlled by asymmetric wettability

Biomimetic asymmetric nanochannels have recently attracted increasing attention from researchers, especially in the aspect of the asymmetric wettability (a hydrophilic-hydrophobic system), which can be utilized to control the wetting behavior of aqueous media and to offer a means for guiding water motion. By using molecular dynamics simulations, a design for a potentially efficient water filter is presented based on (n, n) single-walled carbon nanotubes, where n = 6, 8, 10 and 12, asymmetrically modified with hydrophilic groups (carboxyl, -COOH) at one tip and hydrophobic groups (trifluoromethyl, -CF(3) ) at the other. The reduced water density on the hydrophobic sides of the functionalized nanotubes are observed in both pure water and aqueous electrolyte solution, except for the functionalized (6, 6) tube, due to the change of dipole orientation of the single-file water wire within it. The functionalized (8, 8) tube can significantly maintain the low water density on the hydrophobic side. Both (6, 6) and (8, 8) tubes have relatively high energy barriers at their tips for ion permeation, which can be obtained by calculating the potential of mean force. Such tip functionalization of a nanotube therefore suggests the great possibilities of water transport and filtration, dominated by asymmetric wettability. The functionalized (8, 8) tube could act as a nanofluidic channel for water purification, not only for ion exclusion but also as a stable water column structure.

Control of unidirectional transport of single-file water molecules through carbon nanotubes in an electric field

The transport of water molecules through nanopores is not only crucial to biological activities but also useful for designing novel nanofluidic devices. Despite considerable effort and progress that has been made, a controllable and unidirectional water flow is still difficult to achieve and the underlying mechanism is far from being understood. In this paper, using molecular dynamics simulations, we systematically investigate the effects of an external electric field on the transport of single-file water molecules through a carbon nanotube (CNT). We find that the orientation of water molecules inside the CNT can be well-tuned by the electric field and is strongly coupled to the water flux. This orientation-induced water flux is energetically due to the asymmetrical water-water interaction along the CNT axis. The wavelike water density profiles are disturbed under strong field strengths. The frequency of flipping for the water dipoles will decrease as the field strength is increased, and the flipping events vanish completely for the relatively large field strengths. Most importantly, a critical field strength E(c) related to the water flux is found. The water flux is increased as E is increased for E ≤ E(c), while it is almost unchanged for E > E(c). Thus, the electric field offers a level of governing for unidirectional water flow, which may have some biological applications and provides a route for designing efficient nanopumps.

Theoretic, experimental, clinical bases of the water oscillator hypothesis in near-infrared photobiomodulation

The objective of this review is to propose and document a role for the water oscillator in near-infrared (NIR) photobiomodulation. Greater understanding of the role of the water oscillator may add to a more-coherent description of central effects of NIR light on redox centers and key transmembrane enzymes such as cytochrome c oxidase (CcO). In addition, water provides a complementary pathway for absorption and transportation of NIR energy in photobiomodulation. Because of its unexpected potential, we propose terming it the "water oscillator paradox." Photobiologic mechanisms involved in the treatment of complex diseases are discussed in light of the present state of the art.

pH-tunable ion selectivity in carbon nanotube pores

The selectivity of ion transport in nanochannels is of primary importance for a number of physical, chemical, and biological processes ranging from fluid separation to ion-channel-regulated cellular processes. Fundamental understanding of these phenomena requires model nanochannels with well-defined and controllable structural properties. Carbon nanotubes provide an ideal choice for nanofluidic studies because of their simple chemistry and structure, the atomic scale smoothness and chemical inertness of the graphitic walls, and the tunability of their diameter and length. Here, we investigate the selectivity of single and, for the first time, binary salt mixtures transport through narrow carbon nanotubes that act as the only pores in a silicon nitride membrane. We demonstrate that negatively charged carboxylic groups are responsible for the ion rejection performance of carbon nanotube pores and that ion permeation of small salts can be tuned by varying solution pH. Investigation of the effect of solution composition and ion valences for binary electrolytes with common cation in a pressure-driven flow reveals that the addition of slower diffusing multivalent anions to a solution of faster diffusing monovalent anions favors permeation of the monovalent anion. Larger fractions and valences of the added multivalent anions lower the rejection of the monovalent anion. In some cases, we observe negative rejection at low monovalent ion content.

Simulating proton transport through a simplified model for trans-membrane proteins

Ab initio MD simulations on a polyglycine helix and water-wire expressed under periodic boundary conditions have created a channel that supports proton transfer up to distances of 10.5 A. The effect of varying the density of water molecules in the channel has been investigated. A range of cationic states are identified with widely varying lifetimes. The mechanism of proton transport in this model shares some features with the simulations reported for bulk water, with, e.g., the hydrogen bond distance shortening in the time period leading up to successful proton transfer. However, there are also some important differences such as the observation of a heightened number of proton rattling events. We also observe that the helix plays an important role in directing the behavior of the water wire: the most active proton transport regions of the water-wire are found in areas where the helix is most tightly coiled. Finally, we report on the effects of different DFT functionals to model a water-wire and on the importance of including dispersion corrections to stabilize the alpha-helical structure.

H(aq)+ structures in proton wires inside nanotubes

The hydrated carborane acid H(CHB(11)I(11)).8H(2)O crystallizes in nanometer-diameter tubes of H(aq)(+) enclosed by walls of carborane anions. Three different types of H(aq)(+) clusters are found in these tubes: a symmetrical H(13)O(6)(+) ion with an unusually elongated Zundel-type H(5)O(2)(+) core, two hydrated H(7)O(3)(+) ions, and an unprecedented H(17)O(8)(+) ion having a nearly square core. All of the H(aq)(+) cations show unexpectedly longer O...O separations than in discrete H(aq)(+) ions, indicating greater delocalization of positive charge. The centrosymmetric H(aq)(+) ions are linked via short H bonds, forming a true one-dimensional proton wire.

A one-dimensional dipole lattice model for water in narrow nanopores

We present a recently developed one-dimensional dipole lattice model that accurately captures the key properties of water in narrow nanopores. For this model, we derive three equivalent representations of the Hamiltonian that together yield a transparent physical picture of the energetics of the water chain and permit efficient computer simulations. In the charge representation, the Hamiltonian consists of nearest-neighbor interactions and Coulomb-like interactions of effective charges at the ends of dipole ordered segments. Approximations based on the charge picture shed light on the influence of the Coulomb-like interactions on the structure of nanopore water. We use Monte Carlo simulations to study the system behavior of the full Hamiltonian and its approximations as a function of chemical potential and system size and investigate the bimodal character of the density distribution occurring at small system sizes.

Macroscopically ordered water in nanopores

Water confined into the interior channels of narrow carbon nanotubes or transmembrane proteins forms collectively oriented molecular wires held together by tight hydrogen bonds. Here, we explore the thermodynamic stability and dipolar orientation of such 1D water chains from nanoscopic to macroscopic dimensions. We show that a dipole lattice model accurately recovers key properties of 1D confined water when compared to atomically detailed simulations. In a major reduction in computational complexity, we represent the dipole model in terms of effective Coulombic charges, which allows us to study pores of macroscopic lengths in equilibrium with a water bath (or vapor). We find that at ambient conditions, the water chains filling the tube are essentially continuous up to macroscopic dimensions. At reduced water vapor pressure, we observe a 1D Ising-like filling/emptying transition without a true phase transition in the thermodynamic limit. In the filled state, the chains of water molecules in the tube remain dipole-ordered up to macroscopic lengths of approximately 0.1 mm, and the dipolar order is estimated to persist for times up to approximately 0.1 s. The observed dipolar order in continuous water chains is a precondition for the use of nanoconfined 1D water as mediator of fast long-range proton transport, e.g., in fuel cells. For water-filled nanotube bundles and membranes, we expect anti-ferroelectric behavior, resulting in a rich phase diagram similar to that of a 2D Coulomb gas.

Ion exclusion by sub-2-nm carbon nanotube pores

Biological pores regulate the cellular traffic of a large variety of solutes, often with high selectivity and fast flow rates. These pores share several common structural features: the inner surface of the pore is frequently lined with hydrophobic residues, and the selectivity filter regions often contain charged functional groups. Hydrophobic, narrow-diameter carbon nanotubes can provide a simplified model of membrane channels by reproducing these critical features in a simpler and more robust platform. Previous studies demonstrated that carbon nanotube pores can support a water flux comparable to natural aquaporin channels. Here, we investigate ion transport through these pores using a sub-2-nm, aligned carbon nanotube membrane nanofluidic platform. To mimic the charged groups at the selectivity region, we introduce negatively charged groups at the opening of the carbon nanotubes by plasma treatment. Pressure-driven filtration experiments, coupled with capillary electrophoresis analysis of the permeate and feed, are used to quantify ion exclusion in these membranes as a function of solution ionic strength, pH, and ion valence. We show that carbon nanotube membranes exhibit significant ion exclusion that can be as high as 98% under certain conditions. Our results strongly support a Donnan-type rejection mechanism, dominated by electrostatic interactions between fixed membrane charges and mobile ions, whereas steric and hydrodynamic effects appear to be less important.

Computation of binding free energy with molecular dynamics and grand canonical Monte Carlo simulations

The binding of a ligand to a receptor is often associated with the displacement of a number of bound water molecules. When the binding site is exposed to the bulk region, this process may be sampled adequately by standard unbiased molecular dynamics trajectories. However, when the binding site is deeply buried and the exchange of water molecules with the bulk region may be difficult to sample, the convergence and accuracy in free energy perturbation (FEP) calculations can be severely compromised. These problems are further compounded when a reduced system including only the region surrounding the binding site is simulated. To address these issues, we couple molecular dynamics (MD) with grand canonical Monte Carlo (GCMC) simulations to allow the number of water to fluctuate during an alchemical FEP calculation. The atoms in a spherical inner region around the binding pocket are treated explicitly while the influence of the outer region is approximated using the generalized solvent boundary potential (GSBP). At each step during thermodynamic integration, the number of water in the inner region is equilibrated with GCMC and energy data generated with MD is collected. Free energy calculations on camphor binding to a deeply buried pocket in cytochrome P450cam, which causes about seven water molecules to be expelled, are used to test the method. It concluded that solvation free energy calculations with the GCMC/MD method can greatly improve the accuracy of the computed binding free energy compared to simulations with fixed number of water.

Oriented nanostructures for energy conversion and storage

Recently, the role of nanostructured materials in addressing the challenges in energy and natural resources has attracted wide attention. In particular, oriented nanostructures demonstrate promising properties for energy harvesting, conversion, and storage. In this Review, we highlight the synthesis and application of oriented nanostructures in a few key areas of energy technologies, namely photovoltaics, batteries, supercapacitors, and thermoelectrics. Although the applications differ from field to field, a common fundamental challenge is to improve the generation and transport of electrons and ions. We highlight the role of high surface area to maximize the surface activity and discuss the importance of optimum dimension and architecture, controlled pore channels, and alignment of the nanocrystalline phase to optimize the transport of electrons and ions. Finally, we discuss the challenges in attaining integrated architectures to achieve the desired performance. Brief background information is provided for the relevant technologies, but the emphasis is focused mainly on the nanoscale effects of mostly inorganic-based materials and devices.

Generalized gradient-augmented harmonic Fourier beads method with multiple atomic and/or center-of-mass positional restraints

We describe a generalization of the gradient-augmented harmonic Fourier beads method for finding minimum free-energy transition path ensembles and similarly minimum potential energy paths to allow positional restraints on the centers of mass of selected atoms. The generalized gradient-augmented harmonic Fourier beads (ggaHFB) method further extends the scope of the HFB methodology to studying molecule transport across various mobile phases such as lipid membranes. Furthermore, the new implementation improves the applicability of the HFB method to studies of ligand binding, protein folding, and enzyme catalysis as well as modeling equilibrium pulling experiments. Like its predecessor, the ggaHFB method provides accurate energy profiles along the specified paths and in certain simple cases avoids the need for path optimization. The utility of the ggaHFB method is demonstrated with an application to the water permeation through a single-wall (5,5) carbon nanotube with a diameter of 6.78 A and length of 16.0 A. We provide a simple rationale as to why water enters the hydrophobic nanotube and why it does so in pulses and in wire assembly.

Dynamics of the glutamic acid 242 side chain in cytochrome c oxidase

In many cytochrome c oxidases glutamic acid 242 is required for proton transfer to the binuclear heme a(3)/Cu(B) site, and for proton pumping. When present, the side chain of Glu-242 is orientated "down" towards the proton-transferring D-pathway in all available crystal structures. A nonpolar cavity "above" Glu-242 is empty in these structures. Yet, proton transfer from Glu-242 to the binuclear site, and for proton-pumping, is well established, and the cavity has been proposed to at least transiently contain water molecules that would mediate proton transfer. Such proton transfer has been proposed to require isomerisation of the Glu-242 side chain into an "up" position pointing towards the cavity. Here, we have explored the molecular dynamics of the protonated Glu-242 side chain. We find that the "up" position is preferred energetically when the cavity contains four water molecules, but the "down" position is favoured with less water. We conclude that the cavity might be deficient in water in the crystal structures, possibly reflecting the "resting" state of the enzyme, and that the "up/down" equilibrium of Glu-242 may be coupled to the presence of active-site water molecules produced by O(2) reduction.

Metastable water clusters in the nonpolar cavities of the thermostable protein tetrabrachion

Water expulsion from the protein core is a key step in protein folding. Nevertheless, unusually large water clusters confined into the nonpolar cavities have been observed in the X-ray crystal structures of tetrabrachion, a bacterial protein that is thermostable up to at least 403 K (130 degrees C). Here, we use molecular dynamics (MD) simulations to investigate the structure and thermodynamics of water filling the largest cavity of the right-handed coiled-coil stalk of tetrabrachion at 365 K (92 degrees C), the temperature of optimal bacterial growth, and at room temperature (298 K). Hydrogen-bonded water clusters of seven to nine water molecules are found to be thermodynamically stable in this cavity at both temperatures, confirming the X-ray studies. Stability, as measured by the transfer free energy of the optimal size cluster, decreases with increasing temperature. Water filling is thus driven by the energy of transfer and opposed by the transfer entropy, both depending only weakly on temperature. Our calculations suggest that cluster formation becomes unfavorable at approximately 384 K (110 degrees C), signaling the onset of drying just slightly above the temperature of optimal growth. "Drying" thus precedes protein denaturation. At room temperature, the second largest cavity in tetrabrachion accommodates a five water molecule cluster, as reported in the X-ray studies. However, the simulations show that at 365 K the cluster is unstable and breaks up. We suggest that the large hydrophobic cavities may act as binding sites for two proteases, possibly explaining the unusual thermostability of the resulting protease-stalk complexes (up to approximately 393 K, 120 degrees C).