Search

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/P 2 -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.390 -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., Competing Interests: The authors declare no conflict of interest., (Copyright © 2018 the Author(s). Published by PNAS.)

This special issue presents a series of papers that highlight a new and evolving view of water's molecular role in biological structure and dynamics. Increasingly water is appreciated for playing an active role in biological function rather than just serving as a spectator. The areas represented include molecular interfaces with water, biological self-assembly, conformational changes by macromolecules, and chemical reactivity.

We used molecular dynamics simulations to investigate the thermodynamics of filling of a (6,6) open carbon nanotube (diameter D = 0.806 nm) solvated in TIP3P water over a temperature range from 280 K to 320 K at atmospheric pressure. In simulations of tubes with slightly weakened carbon-water attractive interactions, we observed multiple filling and emptying events. From the water occupancy statistics, we directly obtained the free energy of filling, and from its temperature dependence the entropy of filling. We found a negative entropy of about -1.3 k(B) per molecule for filling the nanotube with a hydrogen-bonded single-file chain of water molecules. The entropy of filling is nearly independent of the strength of the attractive carbon-water interactions over the range studied. In contrast, the energy of transfer depends strongly on the carbon-water attraction strength. These results are in good agreement with entropies of about -0.5 k(B) per water molecule obtained from grand-canonical Monte Carlo calculations of water in quasi-infinite tubes in vacuum under periodic boundary conditions. Overall, for realistic carbon-water interactions we expect that at ambient conditions filling of a (6,6) carbon nanotube open to a water reservoir is driven by a favorable decrease in energy, and opposed by a small loss of water entropy.

Water molecules confined to pores with sub-nanometre diameters form single-file hydrogen-bonded chains. In such nanoscale confinement, water has unusual physical properties that are exploited in biology and hold promise for a wide range of biomimetic and nanotechnological applications. The latter can be realized by carbon and boron nitride nanotubes which confine water in a relatively non-specific way and lend themselves to the study of intrinsic properties of single-file water. As a consequence of strong water-water hydrogen bonds, many characteristics of single-file water are conserved in biological and synthetic pores despite differences in their atomistic structures. Charge transport and orientational order in water chains depend sensitively on and are mainly determined by electrostatic effects. Thus, mimicking functions of biological pores with apolar pores and corresponding external fields gives insight into the structure-function relation of biological pores and allows the development of technical applications beyond the molecular devices found in living systems. In this Perspective, we revisit results for single-file water in apolar pores, and examine the similarities and the differences between these simple systems and water in more complex pores., (This journal is © the Owner Societies 2011)

Proton transfer (pT) reactions in biochemical processes are often mediated by chains of hydrogen-bonded water molecules. We use hybrid density functional calculations to study pT along quasi one-dimensional water arrays that connect an imidazolium-imidazole proton donor-acceptor pair. We characterize the structures of intermediates and transition states, the energetics, and the dynamics of the pT reactions, including vibrational contributions to kinetic isotope effects. In molecular dynamics simulations of pT transition paths, we find that for short water chains with four water molecules, the pT reactions are semi-concerted. The formation of a high-energy hydronium intermediate next to the proton-donating group is avoided by a simultaneous transfer of a proton from the donor to the first water molecule, and from the first water molecule into the water chain. Lowering the dielectric constant of the environment and increasing the water chain length both reduce the barrier for pT. We study the effect of the driving force on the energetics of the pT reaction by changing the proton affinity of the donor and acceptor groups through halogen and methyl substitutions. We find that the barrier of the pT reaction depends linearly on the proton affinity of the donor but is nearly independent of the proton affinity of the acceptor, corresponding to Brønsted slopes of one and zero, respectively., (This journal is © the Owner Societies 2011)

Water in the protein interior serves important structural and functional roles and is also increasingly recognized as a relevant factor in drug binding. The nonpolar cavity in the protein interleukin-1β has been reported to be filled by water on the basis of some experiments and simulations and to be empty on the basis of others. Here we study the thermodynamics of filling the central nonpolar cavity and the four polar cavities of interleukin-1β by molecular dynamics simulation. We use different water models (TIP3P and SPC/E) and protein force fields (amber94 and amber03) to calculate the semigrand partition functions term by term that quantify the hydration equilibria. We consistently find that water in the central nonpolar cavity is thermodynamically unstable, independent of force field and water model. The apparent reason is the relatively small size of the cavity, with a volume less than ∼80 Å(3). Our results are consistent with the most recent X-ray crystallographic and simulation studies but disagree with an earlier interpretation of nuclear magnetic resonance (NMR) experiments probing protein-water interactions. We show that, at least semiquantitatively, the measured nuclear Overhauser effects indicating the proximity of water to the methyl groups lining the nonpolar cavity can, in all likelihood, be attributed to interactions with buried and surface water molecules near the cavity. The same methods applied to determine the occupancy of the polar cavities show that they are filled by the same number of water molecules observed in crystallography, thereby validating the theoretical and simulation methods used to study the water occupancy in the nonpolar protein cavity.

We study the effects of nanoscopic roughness on the interfacial free energy of water confined between solid surfaces. SPC/E water is simulated in confinement between two infinite planar surfaces that differ in their physical topology: one is smooth and the other one is physically rough on a sub-nanometre length scale. The two thermodynamic ensembles considered, with constant pressure either normal or parallel to the walls, correspond to different experimental conditions. We find that molecular-scale surface roughness significantly increases the solid-liquid interfacial free energy compared to the smooth surface. For our surfaces with a water-wall interaction energy minimum of -1.2 kcal mol(-1), we observe a transition from a hydrophilic surface to a hydrophobic surface at a roughness amplitude of about 3 angstroms and a wavelength of 11.6 angstroms, with the interfacial free energy changing sign from negative to positive. In agreement with previous studies of water near hydrophobic surfaces, we find an increase in the isothermal compressibility of water with increasing surface roughness. Interestingly, average measures of the water density and hydrogen-bond number do not contain distinct signatures of increased hydrophobicity. In contrast, a local analysis indicates transient dewetting of water in the valleys of the rough surface, together with a significant loss of hydrogen bonds, and a change in the dipole orientation toward the surface. These microscopic changes in the density, hydrogen bonding, and water orientation contribute to the large increase in the interfacial free energy, and the change from a hydrophilic to a hydrophobic character of the surface.

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.

We study the static and dynamic properties of the water-density fluctuations in the interface of large nonpolar solutes. With the help of extensive molecular dynamics simulations of TIP4P water near smooth spherical solutes, we show that for large solutes, the interfacial density profile is broadened by capillary waves. For purely repulsive solutes, the squared width of the interface increases linearly with the logarithm of the solute size, as predicted by capillary-wave theory. The apparent interfacial tension extracted from the slope agrees with that of a free liquid-vapor interface. The characteristic length of local density fluctuations is approximately 0.5 nm, measured along the arc, again consistent with that of a free liquid-vapor interface. Probed locally, the interfacial density fluctuations exhibit large variances that exceed those expected for an ideal gas. Qualitatively consistent with theories of the free liquid-vapor interface, we find that the water interface near large and strongly nonpolar solutes is flickering, broadened by capillary-wave fluctuations. These fluctuations result in transitions between locally wet and dry regions that are slow on a molecular time scale.

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.

Water molecules confined to nonpolar pores and cavities of nanoscopic dimensions exhibit highly unusual properties. Water filling is strongly cooperative, with the possible coexistence of filled and empty states and sensitivity to small perturbations of the pore polarity and solvent conditions. Confined water molecules form tightly hydrogen-bonded wires or clusters. The weak attractions to the confining wall, combined with strong interactions between water molecules, permit exceptionally rapid water flow, exceeding expectations from macroscopic hydrodynamics by several orders of magnitude. The proton mobility along 1D water wires also substantially exceeds that in the bulk. Proteins appear to exploit these unusual properties of confined water in their biological function (e.g., to ensure rapid water flow in aquaporins or to gate proton flow in proton pumps and enzymes). The unusual properties of water in nonpolar confinement are also relevant to the design of novel nanofluidic and molecular separation devices or fuel cells.

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).

Molecular dynamics simulations on quantum energy surfaces are carried out to study the effects of perturbing electric fields on proton transport (PT) in protonated water chains. As an idealized model of a hydrophobic cavity in the interior of a protein the water molecules are confined into a carbon nanotube (CNT). The water chain connects a hydrated hydronium ion (H3O+) at one end of the CNT and an imidazole molecule at the other end. Without perturbing electric fields PT from the hydronium proton donor to the imidazole acceptor occurs on a picosecond time scale. External perturbations to PT are created by electric fields of varying intensities, normal to the CNT axis, generated by a neutral pair of charges on the nanotube wall. For fields above approximately 0.5 VA, the hydronium ion is effectively trapped at the CNT center, and PT blocked. Fields of comparable strength are generated inside proteins by nearby polar/charged amino acids. At lower fields the system displays a rich dynamic behavior, where the excess charge shuttles back and forth along the water chain before reaching the acceptor group on the picosecond time scale. The effects of the perturbing field on the proton movement are analyzed in terms of structural and dynamic properties of the water chain. The implications of these observations on PT in biomolecular systems and its control by external perturbing fields are discussed.

Formation of a water-expelling nonpolar core is the paradigm of protein folding and stability. Although experiment largely confirms this picture, water buried in "hydrophobic" cavities is required for the function of some proteins. Hydration of the protein core has also been suggested as the mechanism of pressure-induced unfolding. We therefore are led to ask whether even the most nonpolar protein core is truly hydrophobic (i.e., water-repelling). To answer this question we probed the hydration of an approximately 160-A(3), highly hydrophobic cavity created by mutation in T4 lysozyme by using high-pressure crystallography and molecular dynamics simulation. We show that application of modest pressure causes approximately four water molecules to enter the cavity while the protein itself remains essentially unchanged. The highly cooperative filling is primarily due to a small change in bulk water activity, which implies that changing solvent conditions or, equivalently, cavity polarity can dramatically affect interior hydration of proteins and thereby influence both protein activity and folding.

Carbon nanotubes can serve as simple nonpolar water channels. Here we report computer simulations exploring the relationship between the mechanical properties of such channels and their interaction with water. We show that on one hand, increasing the flexibility of the carbon nanotubes increases their apparent hydrophobic character, while on the other hand the presence of water inside the channel makes them more resistant to radial collapse. We quantify the effect of increasing flexibility on the hydrophobicity of the nanotube water channel. We also show that flexibility impedes water transport across the nanotube channel by increasing the free-energy barriers to such motion. Conversely, the presence of water inside the nanotube is shown to affect the energetics of radial collapse in a water nanotube, an ostensibly mechanical property. We quantify the magnitude of the effect and show that it arises from the formation of energetically favorable low-dimensional water structures inside the nanotube such as one-dimensional wires and two-dimensional sheets.

Using a coarse-grained molecular dynamics (CMD) approach we study the apparent nonlinear dynamics of water molecules filling or emptying carbon nanotubes as a function of system parameters. Different levels of the pore hydrophobicity give rise to tubes that are empty, water-filled, or fluctuate between these two long-lived metastable states. The corresponding coarse-grained free-energy surfaces and their hysteretic parameter dependence are explored by linking MD to continuum fixed point and bifurcation algorithms. The results are validated through equilibrium MD simulations.

We explore the structure and thermodynamics of water clusters confined in nonpolar cavities. By calculating the grand-canonical partition function term by term, we show that small nonpolar cavities can be filled at equilibrium with highly structured water clusters. The structural and thermodynamic properties of these encapsulated water clusters are similar to those observed experimentally in the gas phase. Water filling is highly sensitive to the size of the cavity and the strength of the interactions with the cavity wall. Water penetration into pores can thus be modulated by small changes in the polarity and structure of the cavity. Implications on water penetration into proteins are discussed.

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.

We use molecular dynamics simulations to study osmotically driven transport of water molecules through hexagonally packed carbon nanotube membranes. Our simulation setup comprises two such semipermeable membranes separating compartments of pure water and salt solution. The osmotic force drives water flow from the pure-water to the salt-solution compartment. Monitoring the flow at molecular resolution reveals several distinct features of nanoscale flows. In particular, thermal fluctuations become significant at the nanoscopic length scales, and as a result, the flow is stochastic in nature. Further, the flow appears frictionless and is limited primarily by the barriers at the entry and exit of the nanotube pore. The observed flow rates are high (5.8 water molecules per nanosecond and nanotube), comparable to those through the transmembrane protein aquaporin-1, and are practically independent of the length of the nanotube, in contrast to predictions of macroscopic hydrodynamics. All of these distinct characteristics of nanoscopic water flow can be modeled quantitatively by a 1D continuous-time random walk. At long times, the pure-water compartment is drained, and the net flow of water is interrupted by the formation of structured solvation layers of water sandwiched between two nanotube membranes. Structural and thermodynamic aspects of confined water monolayers are studied.

Proton-transfer reactions form an integral part of bioenergetics and enzymatic catalysis. The identification of proton-conducting pathways inside a protein is a key to understanding the mechanisms of biomolecular proton transfer. Proton pathways are modeled here as hydrogen bonded networks of proton-conducting groups, including proton-exchanging groups of amino acid side chains and bound water molecules. We focus on the identification of potential proton-conducting pathways inside a protein of known structure. However, consideration of the static structure alone is often not sufficient to detect suitable proton-transfer paths, leading, for example, from the protein surface to the active site buried inside the protein. We include dynamic fluctuations of amino acid side chains and water molecules into our analysis. To illustrate the method, proton transfer into the active site of cytochrome P450cam is studied. The cooperative rotation of amino acids and motion of water molecules are found to connect the protein surface to the molecular oxygen. Our observations emphasize the intrinsic dynamical nature of proton pathways where critical connections in the network may be transiently provided by mobile groups.

Proton transfer along 1D chains of water molecules inside carbon nanotubes is studied by simulations. Ab initio molecular dynamics and an empirical valence bond model yield similar structures and time scales. The proton mobility along 1D water chains exceeds that in bulk water by a factor of 40, but is reduced if orientational defects are present. Excess protons interact with hydrogen-bonding defects through long-range electrostatics, resulting in coupled motion of protons and defects.

Recent molecular dynamics simulations of water transport through the interior channel of a carbon nanotube in contact with an aqueous reservoir showed that conduction occurred in bursts with collective water motion. A continuous-time random-walk model is used to describe concerted transport through channels densely filled with molecules in a single-file arrangement, as also found in zeolites, as well as ion channels and aquaporins in biological membranes. Theoretical predictions for different collective properties of the single-file transport agree with the simulation results.

The kinetics of water filling and emptying the interior channel of carbon nanotubes is studied by molecular dynamics simulations. Filling and emptying occur predominantly by sequential addition of water to or removal from a single-file chain inside the nanotube. Advancing and receding water chains are orientationally ordered. This precludes simultaneous filling from both tube ends, and forces chain rupturing to occur at the tube end where a water molecule donates a hydrogen bond to the bulk fluid. We use transition path concepts and a Bayesian approach to identify a transition state ensemble that we characterize by its commitment probability distribution. At the transition state, the tube is filled with all but one water molecule. Filling thermodynamics and kinetics depend sensitively on the strength of the attractive nanotube-water interactions. This sensitivity increases with the length of the tubes. [ABSTRACT FROM AUTHOR]

Describes an information theory model of hydrophobic interactions. Probability of observing a molecular-sized cavity in the solvent; Basis of two moments available from the density and radial distribution of oxygen atoms in liquid water; Simplicity and flexibility of the approach.

Abstract: In many cytochrome c oxidases glutamic acid 242 is required for proton transfer to the binuclear heme a 3 /CuB 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 O2 reduction. [Copyright &y& Elsevier]