Reorientation dynamics of nanoconfined water: Power-law decay, hydrogen-bond jumps, and test of a two-state model

Universality in the Structure and Dynamics of Water under Lipidic Mesophase Soft Nanoconfinement

Water under soft nanoconfinement features physical and chemical properties fundamentally different from bulk water; yet, the multitude and specificity of confining systems and geometries mask any of its potentially universal traits. Here, we advance in this quest by resorting to lipidic mesophases as an ideal nanoconfinement system, allowing inspecting the behavior of water under systematic changes in the topological and geometrical properties of the confining medium, without altering the chemical nature of the interfaces. By combining Terahertz absorption spectroscopy experiments and molecular dynamics simulations, we unveil the presence of universal laws governing the physics of nanoconfined water, recapitulating the data collected at varying levels of hydration and nanoconfinement topologies. This geometry-independent universality is evidenced by the existence of master curves characterizing both the structure and dynamics of simulated water as a function of the distance from the lipid-water interface. Based on our theoretical findings, we predict a parameter-free law describing the amount of interfacial water against the structural dimension of the system (i.e., the lattice parameter), which captures both the experimental and numerical results within the same curve, without any fitting. Our results offer insight into the fundamental physics of water under soft nanoconfinement and provide a practical tool for accurately estimating the amount of nonbulk water based on structural experimental data.

Effects of Nanoconfinement on Dynamics in Concentrated Aqueous Magnesium Chloride Solutions

Water behavior in various natural and manufactured settings is influenced by confinement in organic or inorganic frameworks and the presence of solutes. Here, the effects on dynamics from both confinement and the addition of solutes are examined. Specifically, water and ion dynamics in concentrated (2.5-4.2 m) aqueous magnesium chloride solutions confined in mesoporous silica (2.8 nm pore diameter) were investigated using polarization selective pump-probe and 2D infrared spectroscopies. Fitting the rotational and spectral diffusion dynamics measured by the vibrational probe, selenocyanate, with a previously developed two-state model revealed distinct behaviors at the interior of the silica pores (core state) and near the wall of the confining framework (shell state). The shell dynamics are noticeably slower than the bulk, or core, dynamics. The concentration-dependent slowing of the dynamics aligns with behavior in the bulk solutions, but the spectrally separated water-associated and Mg2+-associated forms of the selenocyanate probe exhibit different responses to confinement. The disparity in the complete reorientation times is larger upon confinement, but the spectral diffusion dynamics become more similar near the silica surface. The length scales that characterize the transition from surface-influenced to bulk-like behavior for the salt solutions in the pores are discussed and compared to those of pure water and an organic solvent confined in the same pores. These comparisons offer insights into how confinement modulates the properties of different liquids.

Spin-lattice relaxation time in water/graphene-oxide dispersion

We present the results of the calculations of the spin-lattice relaxation time of water in contact with graphene oxide by means of all-atom molecular dynamics simulations. We fully characterized the water-graphene oxide interaction through the calculation of the relaxation properties of bulk water and of the contact angle as a function of graphene oxide oxidation state and comparing them with the available experimental data. We then extended the calculation to investigate how graphene oxide alters the dynamical and relaxation properties of water in different conditions and concentrations. We show that, despite the diamagnetic nature of the graphene oxide, the confining effects of the bilayers strongly affect the longitudinal relaxation properties of interfacial water, which presents a reduced dynamics due to hydrogen bonds with oxygen groups on graphene oxide. This property makes graphene oxide an interesting platform to investigate water dynamics in confined geometries and an alternative contrast-agent for magnetic resonance imaging applications, especially in view of the possibility to functionalize graphene oxide from theranostic perspectives.

Fluids and Electrolytes under Confinement in Single-Digit Nanopores

Confined fluids and electrolyte solutions in nanopores exhibit rich and surprising physics and chemistry that impact the mass transport and energy efficiency in many important natural systems and industrial applications. Existing theories often fail to predict the exotic effects observed in the narrowest of such pores, called single-digit nanopores (SDNs), which have diameters or conduit widths of less than 10 nm, and have only recently become accessible for experimental measurements. What SDNs reveal has been surprising, including a rapidly increasing number of examples such as extraordinarily fast water transport, distorted fluid-phase boundaries, strong ion-correlation and quantum effects, and dielectric anomalies that are not observed in larger pores. Exploiting these effects presents myriad opportunities in both basic and applied research that stand to impact a host of new technologies at the water-energy nexus, from new membranes for precise separations and water purification to new gas permeable materials for water electrolyzers and energy-storage devices. SDNs also present unique opportunities to achieve ultrasensitive and selective chemical sensing at the single-ion and single-molecule limit. In this review article, we summarize the progress on nanofluidics of SDNs, with a focus on the confinement effects that arise in these extremely narrow nanopores. The recent development of precision model systems, transformative experimental tools, and multiscale theories that have played enabling roles in advancing this frontier are reviewed. We also identify new knowledge gaps in our understanding of nanofluidic transport and provide an outlook for the future challenges and opportunities at this rapidly advancing frontier.

Competing Nuclear Quantum Effects and Hydrogen-Bond Jumps in Hydrated Kaolinite

Recent work has shown that the dynamics of hydrogen bonds in pure clays are affected by nuclear quantum fluctuations, with different effects for the hydrogen bonds holding different layers of the clay together and for those within the same layer. At the clay-water interface there is an even wider range of types of hydrogen bond, suggesting that the quantum effects may be yet more varied. We apply classical and thermostated ring polymer molecular dynamics simulations to show that nuclear quantum effects accelerate hydrogen-bond dynamics to varying degrees. By interpreting the results in terms of the extended jump model of hydrogen-bond switching, we can understand the origins of these effects in terms of changes in the quantum kinetic energy of hydrogen atoms during an exchange. We also show that the extended jump mechanism is applicable not only to the hydrogen bonds involving water, but also those internal to the clay.

Water confinement in small polycyclic aromatic hydrocarbons

The confinement of water molecules is vital in fields from biology to nanotechnology. The conditions allowing confinement in small finite polycyclic aromatic hydrocarbons (PAHs) are unclear, yet are crucial for understanding confinement in larger systems. Here, we report a computational study of water cluster confinement within PAHs dimers. Our results serve as a model for larger carbon allotropes and for understanding molecular interactions in confined systems. We identified size and structural motifs allowing confinement and demonstrated the motifs in various PAHs systems. We show that optimal OH⋯π interactions between water clusters and the PAH dimer permit optimal confinement to occur. However, the lack of such interactions leads to the formation of CH⋯O interactions, resulting in less ideal confinement. Confinement of layered clusters is also possible, provided that the optimal OH⋯π interactions are conserved.

Probing electrolyte–silica interactions through simulations of the infrared spectroscopy of nanoscale pores

The structural and dynamical properties of nanoconfined solutions can differ dramatically from those of the corresponding bulk systems. Understanding the changes induced by confinement is central to controlling the behavior of synthetic nanostructured materials and predicting the characteristics of biological and geochemical systems. A key outstanding issue is how the molecular-level behavior of nanoconfined electrolyte solutions is reflected in different experimental, particularly spectroscopic, measurements. This is addressed here through molecular dynamics simulations of the OH stretching infrared (IR) spectroscopy of NaCl, NaBr, and NaI solutions in isotopically dilute HOD/D2O confined in hydroxylated amorphous silica slit pores of width 1–6 nm and pH [Formula: see text]. In addition, the water reorientation dynamics and spectral diffusion, accessible by pump–probe anisotropy and two-dimensional IR measurements, are investigated. The aim is to elucidate the effect of salt identity, confinement, and salt concentration on the vibrational spectra. It is found that the IR spectra of the electrolyte solutions are only modestly blue-shifted upon confinement in amorphous silica slit pores, with both the size of the shift and linewidth increasing with the halide size, but these effects are suppressed as the salt concentration is increased. This indicates the limitations of linear IR spectroscopy as a probe of confined water. However, the OH reorientational and spectral diffusion dynamics are significantly slowed by confinement even at the lowest concentrations. The retardation of the dynamics eases with increasing salt concentration and pore width, but it exhibits a more complex behavior as a function of halide.

Ultrafast Dynamics and Liquid Structure in Mesoporous Silica: Propagation of Surface Effects in a Polar Aprotic Solvent

Enhancement of processes ranging from gas sorption to ion conduction in a liquid can be substantial upon nanoconfinement. Here, the dynamics of a polar aprotic solvent, 1-methylimidazole (MeIm), in mesoporous silica (2.8, 5.4, and 8.3 nm pore diameters) were examined using femtosecond infrared vibrational spectroscopy and molecular dynamics simulations of a dilute probe, the selenocyanate (SeCN^-) anion. The long vibrational lifetime and sensitivity of the CN stretch enabled a comprehensive investigation of the relatively slow time scales and subnanometer distance dependences of the confined dynamics. Because MeIm does not readily donate hydrogen bonds, its interactions in the hydrophilic silanol pores differ more from the bulk than those of water confined in the same mesopores, resulting in greater structural order and more dramatic slowing of dynamics. The extent of surface effects was quantified by modified two-state models used to fit three spatially averaged experimental observables: vibrational lifetime, orientational relaxation, and spectral diffusion. The length scales and the models (smoothed step, exponential decay, and simple step) describing the transitions between the distinctive shell behavior at the surface and the bulk-like behavior at the pore interior were compared to those of water. The highly nonuniform distributions of the SeCN^- probe and antiparallel layering of MeIm revealed by the simulations guided the interpretation of the results and development of the analytical models. The results illustrate the importance of electrostatic effects and H-bonding interactions in the behavior of confined liquids.

Advanced Materials for Energy-Water Systems: The Central Role of Water/Solid Interfaces in Adsorption, Reactivity, and Transport

The structure, chemistry, and charge of interfaces between materials and aqueous fluids play a central role in determining properties and performance of numerous water systems. Sensors, membranes, sorbents, and heterogeneous catalysts almost uniformly rely on specific interactions between their surfaces and components dissolved or suspended in the water-and often the water molecules themselves-to detect and mitigate contaminants. Deleterious processes in these systems such as fouling, scaling (inorganic deposits), and corrosion are also governed by interfacial phenomena. Despite the importance of these interfaces, much remains to be learned about their multiscale interactions. Developing a deeper understanding of the molecular- and mesoscale phenomena at water/solid interfaces will be essential to driving innovation to address grand challenges in supplying sufficient fit-for-purpose water in the future. In this Review, we examine the current state of knowledge surrounding adsorption, reactivity, and transport in several key classes of water/solid interfaces, drawing on a synergistic combination of theory, simulation, and experiments, and provide an outlook for prioritizing strategic research directions.

Simulations of the IR and Raman spectra of water confined in amorphous silica slit pores

Water in nano-scale confining environments is a key element in many biological, material, and geological systems. The structure and dynamics of the liquid can be dramatically modified under these conditions. Probing these changes can be challenging, but vibrational spectroscopy has emerged as a powerful tool for investigating their behavior. A critical, evolving component of this approach is a detailed understanding of the connection between spectroscopic features and molecular-level details. In this paper, this issue is addressed by using molecular dynamics simulations to simulate the linear infrared (IR) and Raman spectra for isotopically dilute HOD in D₂O confined in hydroxylated amorphous silica slit pores. The effect of slit-pore width and hydroxyl density on the silica surface on the vibrational spectra is also investigated. The primary effect of confinement is a blueshift in the frequency of OH groups donating a hydrogen bond to the silica surface. This appears as a slight shift in the total (measurable) spectra but is clearly seen in the distance-based IR and Raman spectra. Analysis indicates that these changes upon confinement are associated with the weaker hydrogen-bond accepting properties of silica oxygens compared to water molecules.

Dynamics of water trapped in transition metal oxide-graphene nano-confinement

Motivated by practical implementation of transition-metal oxide-graphene heterostructures, we use all atom molecular dynamics simulations to study dynamics of water in a nano slit bounded by a transition metal oxide surface, namely, TiO2termination of SrTiO3, and graphene. The resultant asymmetric, strong confinement produces square ice-like crystallites of water pinned at TiO2surface and drives enhanced hydrophobicity of graphene via the proximity effect to the hydrophilic TiO2surface. This importantly brings in dynamic heterogeneity, both in translational and rotational degrees of freedom, due to coupling between the slow relaxing, strongly adsorbed water layer at the hydrophilic oxide surface, and faster relaxation of subsequent water layers. The heterogeneity is signalled in the ruggedness of the effective free energy landscapes. We discuss possible implications of our findings in drug delivery.

Confined Dynamics of Water in Transmembrane Pore of TRPV1 Ion Channel

Solvation effects play a key role in chemical and biological processes. The microscopic properties of water near molecular surfaces are radically different from those in the bulk. Furthermore, the behavior of water in confined volumes of a nanometer scale, including transmembrane pores of ion channels, is especially nontrivial. Knowledge at the molecular level of structural and dynamic parameters of water in such systems is necessary to understand the mechanisms of ion channels functioning. In this work, the results of molecular dynamics (MD) simulations of water in the pore and selectivity filter domains of TRPV1 (Transient Receptor Potential Vanilloid type 1) membrane channel are considered. These domains represent nanoscale volumes with strongly amphiphilic walls, where physical behavior of water radically differs from that of free hydration (e.g., at protein interfaces) or in the bulk. Inside the pore and filter domains, water reveals a very heterogeneous spatial distribution and unusual dynamics: It forms compact areas localized near polar groups of particular residues. Residence time of water molecules in such areas is at least 1.5 to 3 times larger than that observed for similar groups at the protein surface. Presumably, these water “blobs” play an important role in the functional activity of TRPV1. In particular, they take part in hydration of the hydrophobic TRPV1 pore by localizing up to six waters near the so-called “lower gate” of the channel and reducing by this way the free energy barrier for ion and water transport. Although the channel is formed by four identical protein subunits, which are symmetrically packed in the initial experimental 3D structure, in the course of MD simulations, hydration of the same amino acid residues of individual subunits may differ significantly. This greatly affects the microscopic picture of the distribution of water in the channel and, potentially, the mechanism of its functioning. Therefore, reconstruction of the full picture of TRPV1 channel solvation requires thorough atomistic simulations and analysis. It is important that the naturally occurring porous volumes, like ion-conducting protein domains, reveal much more sophisticated and fine-tuned regulation of solvation than, e.g., artificially designed carbon nanotubes.

Comment on "Hierarchically porous, ultra-strong reduced graphene oxide-cellulose nanocrystal sponges for exceptional adsorption of water contaminants" by N. Yousefi, K. K. W. Wong, Z. Hosseinidoust, H. O. Sørensen, S. Bruns, Y. Zheng and N. Tufenkji, Nanoscale, 2018, 10, 7171

The MB-spot test is extensively applied to specific surface area measurement due to its simplicity compared to the classical BET test. However, the MB-spot test is vulnerable to multiple factors, preventing this newly-prevalent technique from accurately characterising the specific surface area. The inaccuracy of the MB-spot test is discussed here at length, which is urgently needed to avoid mistakes and inconsistencies regarding specific surface area characterization for carbon-based materials.

Why is Ice Slippery? Simulations of Shear Viscosity of the Quasi-Liquid Layer on Ice

The temperature and depth dependence of the shear viscosity (η) of the quasi-liquid layer (QLL) of water on ice-I_h crystals was determined using simulations of the TIP4P/Ice model. The crystals display either the basal {0001} or prismatic {101̅0} facets, and we find that the QLL viscosity depends on the presented facet, the distance from the solid/liquid interface, and the undercooling temperature. Structural order parameters provide two distinct estimates of the QLL widths, which are found to range from 6.0 to 7.8 Å, and depend on the facet and undercooling temperature. Above 260 K, the viscosity of the vapor-adjacent water layer is significantly less viscous than the solid-adjacent layer and is also lower than the viscosity of liquid water.

Rotational dynamics and dynamical transition of water inside hydrophobic pores of carbon nanotubes

Water in a nanoconfined geometry has attracted great interest from the viewpoint of not only basic science but also nanofluidic applications. Here, the rotational dynamics of water inside single-walled carbon nanotubes (SWCNTs) with mean diameters larger than ca. 1.4 nm were investigated systematically using ²H nuclear magnetic resonance spectroscopy with high-purity SWCNTs and molecular dynamics calculations. The results were compared with those for hydrophilic pores. It was found that faster water dynamics could be achieved by increasing the hydrophobicity of the pore walls and decreasing the pore diameters. These results suggest a strategy that paves the way for emerging high-performance filtration/separation devices. Upon cooling below 220 K, it was found that water undergoes a transition from fast to slow dynamics states. These results strongly suggest that the observed transition is linked to a liquid-liquid crossover or transition proposed in a two-liquid states scenario for bulk water.

Photochemistry and Photophysics in Silica-Based Materials: Ultrafast and Single Molecule Spectroscopy Observation

Silica-based materials (SBMs) are widely used in catalysis, photonics, and drug delivery. Their pores and cavities act as hosts of diverse guests ranging from classical dyes to drugs and quantum dots, allowing changes in the photochemical behavior of the confined guests. The heterogeneity of the guest populations as well as the confinement provided by these hosts affect the behavior of the formed hybrid materials. As a consequence, the observed reaction dynamics becomes significantly different and complex. Studying their photobehavior requires advanced laser-based spectroscopy and microscopy techniques as well as computational methods. Thanks to the development of ultrafast (spectroscopy and imaging) tools, we are witnessing an increasing interest of the scientific community to explore the intimate photobehavior of these composites. Here, we review the recent theoretical and ultrafast experimental studies of their photodynamics and discuss the results in comparison to those in homogeneous media. The discussion of the confined dynamics includes solvation and intra- and intermolecular proton-, electron-, and energy transfer events of the guest within the SBMs. Several examples of applications in photocatalysis, (photo)sensors, photonics, photovoltaics, and drug delivery demonstrate the vast potential of the SBMs in modern science and technology.

Water Dynamics in the Hydration Shells of Biomolecules

The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.

Self-diffusion studies by intra- and inter-molecular spin-lattice relaxometry using field-cycling: Liquids, plastic crystals, porous media, and polymer segments

Field-cycling NMR relaxometry is a well-established technique for probing molecular dynamics in a frequency range from typically a few kHz up to several tens of MHz. For the interpretation of relaxometry data, it is quite often assumed that the spin-lattice relaxation process is of an intra-molecular nature so that rotational fluctuations dominate. However, dipolar interactions as the main type of couplings between protons and other dipolar species without quadrupole moments can imply appreciable inter-molecular contributions. These fluctuate due to translational displacements and to a lesser degree also by rotational reorientations in the short-range limit. The analysis of the inter-molecular proton spin-lattice relaxation rate thus permits one to evaluate self-diffusion variables such as the diffusion coefficient or the mean square displacement on a time scale from nanoseconds to several hundreds of microseconds. Numerous applications to solvents, plastic crystals and polymers will be reviewed. The technique is of particular interest for polymer dynamics since inter-molecular spin-lattice relaxation diffusometry bridges the time scales of quasi-elastic neutron scattering and field-gradient NMR diffusometry. This is just the range where model-specific intra-coil mechanisms are assumed to occur. They are expected to reveal themselves by characteristic power laws for the time-dependence of the mean-square segment displacement. These can be favorably tested on this basis. Results reported in the literature will be compared with theoretical predictions. On the other hand, there is a second way for translational diffusion phenomena to affect the spin-lattice relaxation dispersion. If rotational diffusion of molecules is restricted, translational diffusion properties can be deduced even from molecular reorientation dynamics detected by intra-molecular spin-lattice relaxation. This sort of scenario will be relevant for adsorbates on surfaces or polymer segments under entanglement and chain connectivity constraints. Under such conditions, reorientations will be correlated with translational displacements leading to the so-called RMTD relaxation process (reorientation mediated by translational displacements). Applications to porous glasses, protein solutions, lipid bilayers, and clays will be discussed. Finally, we will address the intriguing fact that the various time limits of the segment mean-square displacement of polymers in some cases perfectly reproduce predictions of the tube/reptation model whereas the reorientation dynamics suggests strongly deviating power laws.

Collective dynamic dipole moment and orientation fluctuations, cooperative hydrogen bond relaxations, and their connections to dielectric relaxation in ionic acetamide deep eutectics: Microscopic insight from simulations

The paper reports a detailed simulation study on collective reorientational relaxation, cooperative hydrogen bond (H-bond) fluctuations, and their connections to dielectric relaxation (DR) in deep eutectic solvents made of acetamide and three uni-univalent electrolytes, lithium nitrate (LiNO3), lithium bromide (LiBr), and lithium perchlorate (LiClO4). Because cooperative H-bond fluctuations and ion migration complicate the straightforward interpretation of measured DR timescales in terms of molecular dipolar rotations for these conducting media which support extensive intra- and inter-species H-bonding, one needs to separate out the individual components from the overall relaxation for examining the microscopic origin of various timescales. The present study does so and finds that reorientation of ion-complexed acetamide molecules generates relaxation timescales that are in sub-nanosecond to nanosecond range. This explains in molecular terms the nanosecond timescales reported by recent giga-Hertz DR measurements. Interestingly, the simulated survival timescale for the acetamide-Li(+) complex has been found to be a few tens of nanosecond, suggesting such a cation-complexed species may be responsible for a similar timescale reported by mega-Hertz DR measurements of acetamide/potassium thiocyanate deep eutectics near room temperature. The issue of collective versus single particle relaxation is discussed, and jump waiting time distributions are determined. Dependence on anion-identity in each of the cases has been examined. In short, the present study demonstrates that assumption of nano-sized domain formation is not required for explaining the DR detected nanosecond and longer timescales in these media.

A new intermolecular potential for simulations of methanol: The OPLS/2016 model

In this work, a new rigid-nonpolarizable model of methanol is proposed. The model has three sites, located at the same positions as those used in the OPLS model previously proposed by Jorgensen [J. Phys. Chem. 90, 1276 (1986)]. However, partial charges and the values of the Lennard-Jones parameters were modified by fitting to an adequately selected set of target properties including solid-fluid experimental data. The new model was denoted as OPLS/2016. The overall performance of this model was evaluated and compared to that obtained with other popular models of methanol using a similar test to that recently proposed for water models. In the test, a certain numerical score is given to each model. It was found that the OPLS/2016 obtained the highest score (7.4 of a maximum of 10) followed by L1 (6.6), L2 (6.4), OPLS (5.8), and H1 (3.5) models. The improvement of OPLS/2016 with respect to L1 and L2 is mainly due to an improvement in the description of fluid-solid equilibria (the melting point is only 14 K higher than the experimental value). In addition, it was found that no methanol model was able to reproduce the static dielectric constant and the isobaric heat capacity, whereas the better global performance was found for models that reproduce the vaporization enthalpy once the so-called polarization term is included. Similar conclusions were suggested previously in the analysis of water models and are confirmed here for methanol.

Hydration in silica based mesoporous materials: a DFT model

In this work, calculable and realistic DFT models of MCM-41 material that follow temperature dependence of silanol density were developed. They can be easily applied in further studies of adsorption or as a support for catalysts.

Exploration of the presence of bulk-like water in AOT reverse micelles and water-in-oil nanodroplets: the role of charged interfaces, confinement size and properties of water

We show that the distance from the interface at which bulk-like properties are recovered strongly depends on the choice of order parameter being probed: translational < tetrahedral ≪ dipolar orientation.

Macro and nano scale modelling of water–water interactions at ambient and low temperature: relaxation and residence times

A method to predict magnitudes in quantitative agreement with experimental data has been devised and applied to model water–water interactions at low temperature.

Exploring the role of hydration and confinement in the aggregation of amyloidogenic peptides Aβ(16-22) and Sup35(7-13) in AOT reverse micelles

Knowledge of how intermolecular interactions of amyloid-forming proteins cause protein aggregation and how those interactions are affected by sequence and solution conditions is essential to our understanding of the onset of many degenerative diseases. Of particular interest is the aggregation of the amyloid-β (Aβ) peptide, linked to Alzheimer's disease, and the aggregation of the Sup35 yeast prion peptide, which resembles the mammalian prion protein linked to spongiform encephalopathies. To facilitate the study of these important peptides, experimentalists have identified small peptide congeners of the full-length proteins that exhibit amyloidogenic behavior, including the KLVFFAE sub-sequence, Aβ16-22, and the GNNQQNY subsequence, Sup357-13. In this study, molecular dynamics simulations were used to examine these peptide fragments encapsulated in reverse micelles (RMs) in order to identify the fundamental principles that govern how sequence and solution environment influence peptide aggregation. Aβ16-22 and Sup357-13 are observed to organize into anti-parallel and parallel β-sheet arrangements. Confinement in the sodium bis(2-ethylhexyl) sulfosuccinate (AOT) reverse micelles is shown to stabilize extended peptide conformations and enhance peptide aggregation. Substantial fluctuations in the reverse micelle shape are observed, in agreement with earlier studies. Shape fluctuations are found to facilitate peptide solvation through interactions between the peptide and AOT surfactant, including direct interaction between non-polar peptide residues and the aliphatic surfactant tails. Computed amide I IR spectra are compared with experimental spectra and found to reflect changes in the peptide structures induced by confinement in the RM environment. Furthermore, examination of the rotational anisotropy decay of water in the RM demonstrates that the water dynamics are sensitive to the presence of peptide as well as the peptide sequence. Overall, our results demonstrate that the RM is a complex confining environment where substantial direct interaction between the surfactant and peptides plays an important role in determining the resulting ensemble of peptide conformations. By extension the results suggest that similarly complex sequence-dependent interactions may determine conformational ensembles of amyloid-forming peptides in a cellular environment.

Role of Charge and Solvation in the Structure and Dynamics of Alanine-Rich Peptide AKA2 in AOT Reverse Micelles

The propensity of peptides to form α-helices has been intensely studied using theory, computation, and experiment. Important model peptides for the study of the coil-to-helix transition have been alanine-lysine (AKA) peptides in which the lysine residues are placed on opposite sides of the helix avoiding charge repulsion while enhancing solubility. In this study, the effects of capped versus zwitterionic peptide termini on the secondary structure of alanine-rich peptides in reverse micelles are explored. The reverse micelles are found to undergo substantial shape fluctuations, a property observed in previous studies of AOT reverse micelles in the absence of solvated peptide. The peptides are observed to interact with water, as well as the AOT surfactant, including interactions between the nonpolar residues and the aliphatic surfactant tails. Computation of IR spectra for the amide I band of the peptide allows for direct comparison with experimental spectra. The results demonstrate that capped AKA2 peptides form more stable α helices than zwitterionic AKA2 peptides in reverse micelles. The rotational anisotropy decay of water is found to be distinctly different in the presence or absence of peptide within the reverse micelle, suggesting that the introduction of peptide significantly alters the number of free waters within the reverse micelle nanopool. However, neither the nature of the peptide termini (capped or charged) nor the degree of peptide helicity is found to significantly alter the balance of interactions between the peptides and the environment. Observed changes in the degree of helicity in AKA2 peptides in bulk solution and in reverse micelle environments result from changes in peptide confinement and hydration as well as direct nonpolar and polar interactions with the water-surfactant interface.

Slow dynamics of water confined in Newton black films

Slowdown of translational and reorientational dynamics of water confined in Newton black films (NBFs) is revealed by molecular dynamics simulations. As a film becomes thinner, both translational and reorientational dynamics become slower. The polarization of water molecules in the macroscopic electrostatic field across the NBF and the coordination of Na(+) ions and surfactant anionic groups around water molecules concertedly lead to slowdown of water dynamics. The polarization effect is obvious for water not coordinated by Na(+) ions, which exhibits reorientational dynamics depending on initial dipole orientations. Na(+) ions and surfactant anionic groups retard dynamics of surrounding water by decreasing the hydrogen bond exchange probability and increasing the viscosity of water. The dependences of translational and reorientational dynamics on coordination environments of water are similar. Dynamics of water in positions close to the interfaces of NBFs are mainly retarded by Na(+) ions and surfactant anionic groups, while the macroscopic polarization effect plays the main role in influencing water dynamics in positions far from the interfaces. This study sheds light on the improvement of knowledge about the water dynamics slowdown mechanism in similar environments like reverse micelles and lamellar structures.

Inhomogeneity of the amorphous solid water dangling bonds

Amorphous solid water (ASW) is one of the most widely studied molecular systems because of its importance in the physics and chemistry of the interstellar medium and the upper layers of the Earth's atmosphere. Although the global structure of this material, i.e. the bulk and the surface, is well characterised, we are far from having an overall understanding of the changes induced upon chemical or physical perturbation. More specifically, the behaviour of the surface and the immediate sublayers upon mid-infrared irradiation must be understood due to its direct effect on the adsorption capacities of the ASW surface. Small molecules can accrete or form at the surface, adsorbed on the dangling OH groups of surface water molecules. This behaviour allows further reactivity which, in turn, could lead to more complex molecular systems. We have already demonstrated that selective IR irradiations of surface water molecules induce a modification of the surface and the production of a new monomer species which bonds to the surface via its two electronic doublets. However, we did not probe the structure of the dangling bands, namely their homogeneity or inhomogeneity. The structure and orientation of these surface molecules are closely linked to the way the surface can relax its vibrational energy. In this work, we have focussed our attention on the two dH dangling bonds, carrying out a series of selective irradiations which reveal the inhomogeneity of these surface modes. We have also studied the effects of irradiation duration on the surface reorientation, determining that the maximum photoinduced isomerisation yield is ∼15%.

Vibrational spectroscopy of water in hydrated lipid multi-bilayers. III. Water clustering and vibrational energy transfer

Water clustering and connectivity around lipid bilayers strongly influences the properties of membranes and is important for functions such as proton and ion transport. Vibrational anisotropic pump-probe spectroscopy is a powerful tool for understanding such clustering, as the measured anisotropy depends upon the time-scale and degree of intra- and intermolecular vibrational energy transfer. In this article, we use molecular dynamics simulations and theoretical vibrational spectroscopy to help interpret recent experimental measurements of the anisotropy of water in lipid multi-bilayers as a function of both lipid hydration level and isotopic substitution. Our calculations are in satisfactory agreement with the experiments of Piatkowski, Heij, and Bakker, and from our simulations we can directly probe water clustering and connectivity. We find that at low hydration levels, many water molecules are in fact isolated, although up to 70% of hydration water forms small water clusters or chains. At intermediate hydration levels, water forms a wide range of cluster sizes, while at higher hydration levels, the majority of water molecules are part of a large, percolating water cluster. Therefore, the size, number, and nature of water clusters are strongly dependent on lipid hydration level, and the measured anisotropy reflects this through its dependence on intermolecular energy transfer.

Distributions of single-molecule properties as tools for the study of dynamical heterogeneities in nanoconfined water

The explicit trend of the distribution functions of single-molecule rotational relaxation constants and atomic mean-square displacement are used to study the dynamical heterogeneities in nanoconfined water. The trend of the single-molecule properties distributions is related to the dynamic heterogeneities, and to the dynamic crossovers found in water clusters of different shapes and sizes and confined in a variety of zeolites. This was true in all the cases that were considered, in spite of the various shapes and sizes of the clusters. It is confirmed that the high temperature dynamical crossover occurring in the temperature range 200-230 K can be interpreted at a molecular level as the formation of almost translationally rigid clusters, characterized by some rotational freedom, hydrogen bond exchange and translational jumps as cage-to-cage processes. We also suggest a mechanism for the low temperature dynamical crossover (LTDC), falling in the temperature range 150-185 K, through which the adsorbed water clusters are made of nearly rigid sub-clusters, slightly mismatched, and thus permitting a relatively free librational motion at their borders. It appears that the condition required for LTDC to occur is the presence of highly heterogeneous environments for the adsorbed molecules, with some dangling hydrogen bonds or weaker than water-water hydrogen bonds. Under these conditions some dynamics are permitted at very low temperature, although most rotational motion is frozen. Therefore, it is unlikely, though not entirely excluded, that LTDC will be found in supercooled bulk water where no heterogeneous interface is present.