Structure and dynamics of nanoconfined water and aqueous solutions.

This review is devoted to discussing recent progress on the structure, thermodynamic, reactivity, and dynamics of water and aqueous systems confined within different types of nanopores, synthetic and biological. Currently, this is a branch of water science that has attracted enormous attention of researchers from different fields interested to extend the understanding of the anomalous properties of bulk water to the nanoscopic domain. From a fundamental perspective, the interactions of water and solutes with a confining surface dramatically modify the liquid's structure and, consequently, both its thermodynamical and dynamical behaviors, breaking the validity of the classical thermodynamic and phenomenological description of the transport properties of aqueous systems. Additionally, man-made nanopores and porous materials have emerged as promising solutions to challenging problems such as water purification, biosensing, nanofluidic logic and gating, and energy storage and conversion, while aquaporin, ion channels, and nuclear pore complex nanopores regulate many biological functions such as the conduction of water, the generation of action potentials, and the storage of genetic material. In this work, the more recent experimental and molecular simulations advances in this exciting and rapidly evolving field will be reported and critically discussed.

Solvation of Coumarin 480 within nano-confining environments: structure and dynamics.

Equilibrium and dynamical characteristics pertaining to the solvation of the fluorescent probe Coumarin 480 within different confining environments are investigated using molecular dynamics simulations. Three kinds of confining systems are examined: (i) the cetyltrimethylammonium bromide (CTAB)/isooctane/1-hexanol/water; cationic inverse micelle (IM) (ii) a CTAB/water direct micelle (DM), and (iii) a silica-surfactant nanocomposite, comprising a cylindrical silica pore (SP) containing small amounts of water and CTAB species adsorbed at the pore walls. The solvation structures in the three environments differ at a qualitative level: an exchange between bulk- and interface-like solvation states was found in the IM, whereas in the DM, the solvation states of the probe are characterized by its embedding at the interface, trapped among the surfactant heads and tails. Within the SP structure, the coumarin exhibits alternations between internal and interfacial solvation states that occur on a ∼20 ns time scale and operate via 90° rotations of its molecular plane. The solvation responses of the environment following a vertical excitation of the probe are also investigated. Solvation times resulted between 2 and 1000 times longer than those found in bulk water, with a fast-to-slow trend IM→DM→SP, which can be interpreted in terms of the solvation structures that prevail in each case.

Structure and dynamics of nonionic surfactants adsorbed at vacuum/ionic liquid interfaces.

Structural and dynamical properties related to the adsorption of nonionic surfactants at vacuum/ionic liquid interfaces were studied using molecular dynamics simulations. Specifically, the surface activity of pentaethylene glycol monododecyl ether (C12E5) was investigated at the free interface of an imidazolium-based room temperature ionic liquid (RTIL), at different surface densities. At low surface coverages, the incorporation of C12E5 does not produce meaningful changes in the vacuum/RTIL interface: the C12E5 hydrophobic tails remain entangled with those of the RTIL cation groups in the outer shell, whereas the C12E5 hydrophilic heads reside at an inner layer. At high surface coverages, the structure in the substrate-in terms of the features exhibited by the local density profiles-practically vanishes; the interface becomes wider and the surfactant molecules shift toward more external positions. Information about the local structure of the interface at high surface densities can be recovered by performing a tessellation procedure. For the sake of comparison, the surface behavior of two commonly used ionic surfactants, sodium dodecyl sulfate and dodecyl trimethyl ammonium chloride, were also studied. The modifications in the width and structure of the bare vacuum/RTIL interface due to the presence of the ionic surfactants are markedly milder than those observed for the nonionic surfactant. Moreover, the RTIL seemed to behave as a better solvent for the chloride counterions than for sodium ones; which were found to remain bound to the surfactant head groups. An analysis of the dynamics at the surface was also performed. Our results indicate that the presence of increasing amounts of nonionic surfactants leads to a gradual reduction of the mobility of the RTIL species. When ionic surfactants are adsorbed, these retardations are even more severe for the surfactant head groups, where the corresponding diffusion coefficients show reductions of practically 1 order of magnitude.

Structure and dynamics of liquid methanol confined within functionalized silica nanopores.

Molecular dynamics simulations have been carried out to investigate the structure and dynamics of liquid methanol confined in 3.3 nm diameter cylindrical silica pores. Three cavities differing in the characteristics of the functional groups at their walls have been examined: (i) smooth hydrophobic pores in which dispersive forces prevail, (ii) hydrophilic cavities with surfaces covered by polar silanol groups, and (iii) a much more rugged pore in which 60% of the previous interfacial hydroxyl groups were replaced by the bulkier trimethylsilyl ones. Confinement promotes a considerable structure at the vicinity of the pore walls which is enhanced in the case of hydroxylated surfaces. Moreover, in the presence of the trimethylsilyl groups, the propagation of this interface-induced spatial ordering extends down to the central region of the pore. Concerning the dynamical modes, we observed an overall slowdown in both the translational and rotational motions. An analysis of these mobilities from a local perspective shows that the largest retardations operate at the vicinity of the interfaces. The gross features of the rotational dynamics were analyzed in terms of contributions arising from bulk and surface states. Compared to the bulk dynamical behavior, the characteristic timescales associated with the rotational motions show the most dramatic increments. A dynamical analysis of hydrogen bond formation and breaking processes is also included.

Encapsulation of small ionic molecules within alpha-cyclodextrins.

Results from molecular dynamics experiments pertaining to the encapsulation of ClO4- within the hydrophobic cavity of an aqueous alpha-cyclodextrin (alpha-CD) are presented. Using a biased sampling procedure, we constructed the Gibbs free energy profile associated with the complexation process. The profile presents a global minimum at the vicinity of the primary hydroxyl groups, where the ion remains tightly coordinated to four water molecules via hydrogen bonds. Our estimate for the global free energy of encapsulation yields DeltaGenc approximately -2.5 kBT. The decomposition of the average forces acting on the trapped ion reveals that the encapsulation is controlled by Coulomb interactions between the ion and OH groups in the CD, with a much smaller contribution from the solvent molecules. Changes in the previous results, arising from the partial methylation of the host CD and modifications in the charge distribution of the guest molecule are also discussed. The global picture that emerges from our results suggests that the stability of the ClO4- encapsulation involves not only the individual ion but also its first solvation shell.

Polar mixtures under nanoconfinement.

We present results from molecular dynamics simulations describing structural and dynamical characteristics of equimolar mixtures of water and acetonitrile, confined between two silica walls separated at interplate distances of d=0.6, 1, and 1.5 nm. Two different environments were investigated: a first one where wall-solvent dispersion forces prevail (hydrophobic confinement) and a second one in which the terminal O atoms at the silica surface are transformed into silanol groups (hydrophilic confinement). For the former case, we found that, at the shortest interplate distance examined, the confined region is devoid of water molecules. At an interplate distance of the order of 1 nm, water moves into the confined region, although, in all cases, there is a clear enhancement of the local concentration of acetonitrile in detriment of that of water. Within hydrophilic environments, we found clear distinctions between a layer of bound water lying in close contact with the silica substrates and a minority of confined water that occupies the inner liquid slab. The bound aqueous layer is fully coordinated to the silanol groups and exhibits minimal hydrogen bonding with the second solvation layer, which exclusively includes acetonitrile molecules. Dynamical characteristics of the solvent mixture are analyzed in terms of diffusive and rotational motions in both environments. Compared to bulk mixtures, we found significant retardations in all dynamical modes, with those ascribed to water molecules bound to the hydrophilic plates being the most dramatic.

Molecular dynamics simulations of AOT-water/formamide reverse micelles: structural and dynamical properties.

We present results from molecular dynamics simulations performed on reverse micelles immersed in cyclohexane. Three different inner polar phases are considered: water (W), formamide (FM), and an equimolar mixture of the two solvents. In all cases, the surfactant was sodium bis(2-ethylhexyl) sulfosuccinate (usually known as AOT). The initial radii of the micelles were R approximately 15 A, while the corresponding polar solvent-to-surfactant molar ratios were intermediate between w(0)=4.3 for FM and w(0)=7 for W. The resulting overall shapes of the micelles resemble distorted ellipsoids, with average eccentricities of the order of approximately 0.75. Moreover, the pattern of the surfactant layer separating the inner pool from the non-polar phase looks highly irregular, with a roughness characterized by length scales comparable to the micelle radii. Solvent dipole orientation polarization along radial directions exhibit steady growths as one moves from central positions toward head group locations. Local density correlations within the micelles indicate preferential solvation of sodium ionic species by water, in contrast to the behavior found in bulk equimolar mixtures. Still, a sizable fraction of approximately 90% of Na(+) remains associated with the head groups. Compared to bulk results, the translational and rotational modes of the confined solvents exhibit important retardations, most notably those operated in rotational motions where the characteristic time scales may be up to 50 times larger. Modifications of the intramolecular connectivity expressed in terms of the average number of hydrogen bonds and their lifetimes are also discussed.

Influence of Cholesterol on the Dynamics of Hydration in Phospholipid Bilayers.

We investigate the dynamics of interfacial waters in dipalmitoylphosphatidylcholine (DPPC) bilayers upon the addition of cholesterol, by molecular dynamics simulations. Our data reveal that the inclusion of cholesterol modifies the membrane aqueous interfacial dynamics: waters diffuse faster, their rotational decay time is shorter, and the DPPC/water hydrogen bond dynamics relaxes faster than in the pure DPPC membrane. The observed acceleration of the translational water dynamics agrees with recent experimental results, in which, by means of NMR techniques, an increment of the surface water diffusivity is measured upon the addition of cholesterol. A microscopic analysis of the lipid/water hydrogen bond network at the interfacial region suggests that the mechanism underlying the observed water mobility enhancement is given by the rupture of a fraction of interlipid water bridge hydrogen bonds connecting two different DPPC molecules, concomitant to the formation of new lipid/solvent bonds, whose dynamics is faster than that of the former. The consideration of a simple two-state model for the decay of the hydrogen bond correlation function yielded excellent results, obtaining two well-separated characteristic time scales: a slow one (∼250 ps) associated with bonds linking two DPPC molecules, and a fast one (∼15 ps), related to DPPC/solvent bonds.

Molecular dynamics study of polarizability anisotropy relaxation in aromatic liquids and its connection with local structure.

The collective polarizability anisotropy dynamics in a set of three aromatic liquids, benzene (Bz), hexafluorobenzene (HFB), and 1,3,5-trifluorobenzene (TFB), has been studied by molecular dynamics simulation. These liquids have very similar shapes, but different electrostatic interactions due to opposite polarities of C-H and C-F bonds, giving rise to different local intermolecular structures in the liquid phase. We have investigated how these structural arrangements affect polarizability anisotropy dynamics observed in optical Kerr-effect (OKE) spectroscopy. We have modeled the interaction-induced polarizability with the first-order dipole-induced dipole approximation, with the molecular polarizability distributed over the carbon sites. Local contributions to the librational OKE spectrum were computed separately for molecules participating in parallel or perpendicular relative orientations within the first coordination shell. We found that the relative locations of parallel and perpendicular librational bands of the OKE spectra are closely related to the corresponding pair energy distributions of the closest four neighbors of a given molecule, corresponding to a model of a harmonic oscillator in a cage of nearest neighbors. This model predicts higher librational frequencies for more attractive intermolecular interactions, which in all three liquids correspond to parallel local arrangements. On the diffusive orientational time scale, all three liquids exhibit slower relaxation of molecules in parallel arrangements, although the difference in relaxation rates is substantial only in TFB, which has the strongest tendency toward parallel stacking. The analysis of the collective polarizability relaxation was performed using two different approaches, the projection scheme (J. Chem. Phys. 1980, 72, 2801) and the theory developed by Steele (Mol. Phys. 1987, 61, 1031) for the second time derivatives applied to collective time correlations. Both approaches allow the decomposition of the OKE response into contributions from orientational relaxation and other dynamical processes. We find that they lead to different predictions on how the response depends on collective reorientation and processes arising from fluctuations in the interaction-induced polarizability. We discuss the reasons for these differences and the advantages and disadvantages of the two analysis schemes.

Molecular dynamics simulations of the dielectric properties of fructose aqueous solutions.

The static dielectric permittivity and dielectric relaxation properties of fructose aqueous solutions of different concentrations ranging from 1.0 to 4.0 mol l(-1) are investigated by means of molecular dynamics simulations. The contributions from intra- and interspecies molecular correlations were computed individually for both the static and frequency-dependent dielectric properties, and the results were compared with the available experimental data. Simulation results in the time- and frequency-domains were analyzed and indicate that the presence of fructose has little effect on the position of the fast, high-frequency (>500 cm(-1)) components of the dielectric response spectrum. The low-frequency (<0.1 cm(-1)) components, however, are markedly influenced by sugar concentration. Our analysis indicates that fructose-fructose and fructose-water interactions strongly affect the rotational-diffusion regime of molecular motions in the solutions. Increasing fructose concentration not only enhances sugar-sugar and sugar-water low frequency contributions to the dielectric loss spectrum but also slows down the reorientational dynamics of water molecules. These results are consistent with previous computer simulations carried out for other disaccharide aqueous solutions.

Effects of molecular association on polarizability relaxation in liquid mixtures of benzene and hexafluorobenzene.

In this work we have studied the relaxation dynamics of the many-body polarizability anisotropy in liquid mixtures of benzene (Bz) and hexafluorobenzene (Hf) at room temperature by femtosecond optical heterodyne-detected Raman-induced Kerr effect spectroscopy (OHD-RIKES) experiments and molecular dynamics (MD) simulations. The computed polarizability response arising from intermolecular interactions was included using the first-order dipole-induced-dipole model with the molecular polarizability distributed over the carbon sites of each molecule. We found good qualitative agreement between experiments and simulations in the features exhibited by the nuclear response function R(t) for pure liquids and mixtures. The long-time diffusive decay of R(t) was observed to vary substantially with composition, slowing down noticeably with dilution of each of the species as compared with that in the corresponding pure liquids. MD simulation shows that the effect on R(t) is due to the formation of strong and localized intermolecular association between Bz and Hf species that hinder the rotational diffusive dynamics. The formation of these Bz-Hf complexes in the liquid mixtures also modifies the rotational diffusive dynamics of the component species in such a way that cannot be explained solely in terms of a viscosity effect. Even though the computed orientational diffusive relaxation times associated with Bz and Hf are larger by a factor of approximately 2 than those from experiments, we found similar trends in experiments and simulations for these characteristic times as a function of composition. Namely, the collective and single-molecule orientational correlation times associated with Bz are observed to grow monotonically with the dilution of Bz, while those corresponding to Hf species exhibit a maximum at the equimolar composition. We attribute the quantitative discrepancy between experiments and simulations to the use of the Williams potential, which seems to overestimate the intermolecular interactions and thus predicts not only a slower translational dynamics but also a slower rotational diffusion dynamics than in real fluids.

Solvation dynamics in acetonitrile: a study incorporating solute electronic response and nuclear relaxation.

The solvent reorganization process after electronic excitation of a polar solute in a polar solvent such as acetonitrile is related mainly to the time evolution of the solute-solvent electrostatic interaction. Modern laser-based techniques have sufficient time resolution to follow this decay in real time, providing information to be confirmed and interpreted by theories and models. We present here a study aimed at the investigation of the different steps involved in the process taking place after a vertical S(0) --> S(1) excitation of a large size chromophore, coumarin 153 (C153), in acetonitrile, from both the solute and the solvent points of view. To do this, we use accurate quantum mechanical calculations for the solute properties within the polarizable continuum model (PCM) and classical molecular dynamics (MD) simulations, both equilibrium and nonequilibrium, for C153 in the presence of the solvent. The geometry of the solute is allowed to change in order to study the role of internal motions in the time-dependent solvation process. The solvent response function has been obtained from the simulation data and compared to experiment, while the comparison between equilibrium and nonequilibrium MD results for the solvation response confirms the validity of the linear response approximation in the C153-acetonitrile system. The MD trajectories have also been used to monitor the structure of the solvation shell and to determine its change in response to the change in the solute partial charges.

Molecular Dynamics Simulations of Ibuprofen Release from pH-Gated Silica Nanochannels.

The iboprufen delivery process from cylindrical silica pores of 3 nm diameter, with polyamine chains anchored at the pore outlets, was investigated by means of massive molecular dynamics simulations. Effects from pH were introduced by considering polyamine chains with different degrees of protonation. High, low, and intermediate pH environments were investigated. The increment of the acidity of the environment leads to a significant decrease of the pore aperture, yielding an effective diameter, for the lowest pH case, that is 3.5 times smaller than the one associated with the highest pH one. Using a biased sampling procedure, Gibbs free energy profiles for the ibuprofen delivery process were obtained. The joint analysis of the corresponding profiles, time evolution of the ibuprofen position within the channel, orientation of the molecule, and instantaneous effective diameter of the gate suggest a three-step mechanism for ibuprofen delivery. A complementary analysis of the translational mobility of ibuprofen along the axial direction of the channel revealed a subdiffusive dynamics in the low and intermediate pH cases. Deviations from Brownian diffusive dynamics are discussed and compared with direct experimental results.

Equilibrium and Dynamical Characteristics of Imidazole Langmuir Monolayers on Graphite Sheets.

Using molecular dynamics techniques, we examine structural and dynamical characteristics of liquid-like imidazole (Im) monolayers physisorbed onto a planar graphite sheet, at T = 384 K. Our simulations reveal that molecular orientations in the saturated monolayer exhibit a bistable distribution, characterized by an inner parallel arrangement of the molecules in close contact with the substrate and a slanted alignment, in those lying in adjacent, outer locations. Compared to the results found in three-dimensional, bulk phases, the analysis of the spatial correlations between sites participating in hydrogen bonding shows a clear enhancement of the intermolecular interactions, which also leads to stronger dipolar correlations. As a result, the gross structural features of the monolayer can be cast in terms of mesoscopic domains, comprising units articulated via winding hydrogen bonds, that persist along typical time intervals of a few tens of picoseconds. On the dynamical side, a similar comparison of the characteristic decorrelation time for orientational motions shows a 4-fold increment. Contrasting, the reduction of the system dimensionality leads to a larger diffusion constant. Possible substrate-induced anisotropies in the diffusive motions are also investigated.

Purification and some characteristics of a beta-galactoside binding soluble lectin from amphibian ovary.

Soluble extracts of Bufo ovaries agglutinate sialidase-treated rabbit erythrocytes. Unlike other amphibian lectins this agglutination activity does not require the presence of calcium ions. It is specifically inhibited by D-galactose and its derivatives. Thiodi-D-galactoside is the most potent saccharide inhibitor followed by lactose and methyl-beta-D-galactoside, respectively. D-Fucose, D-glucose and D-mannose do not inhibit the activity at concentrations at or above 100 mM. The lectin has been purified 500-fold to apparent homogeneity from the ovaries by salt extraction and affinity chromatography on lactose-aminophenyl-agarose, with a yield of about 0.2%. The molecular mass determined by gel filtration under native conditions was 30 kDa; polyacrylamide gel electrophoresis in SDS gave a molecular mass of 15 kDa, suggesting that the lectin is a dimer. The lectin has an isoelectric point of 40 and contains a high proportion of acidic amino acids.

Animal lectins as self/non-self recognition molecules. Biochemical and genetic approaches to understanding their biological roles and evolution.

In recent years, the significant contributions from molecular research studies on animal lectins have elucidated structural aspects and provided clues not only to their evolution but also to their multiple biological functions. The experimental evidence has suggested that distinct, and probably unrelated, groups of molecules are included under the term "lectin." Within the invertebrate taxa, major groups of lectins can be identified: One group would include lectins that show significant homology to membrane-integrated or soluble vertebrate C-type lectins. The second would include those beta-galactosyl-specific lectins homologous to the S-type vertebrate lectins. The third group would be constituted by lectins that show homology to vertebrate pentraxins that exhibit lectin-like properties, such as C-reactive protein and serum amyloid P. Finally, there are examples that do not exhibit similarities to any of the aforementioned categories. Moreover, the vast majority of invertebrate lectins described so far cannot yet be placed in one or another group because of the lack of information regarding their primary structure. (See Table 1.) Animal lectins do not express a recombinatorial diversity like that of antibodies, but a limited diversity in recognition capabilities would be accomplished by the occurrence of multiple lectins with distinct specificities, the presence of more than one binding site, specific for different carbohydrates in a single molecule, and by certain "flexibility" of the binding sites that would allow the recognition of a range of structurally related carbohydrates. In order to identify the lectins' "natural" ligands, we have investigated the interactions between those proteins and the putative endogenous or exogenous glycosylated substances or cells that may be relevant to their biological function. Results from these studies, together with information on the biochemical properties of invertebrate and vertebrate lectins, including their structural relationships with other vertebrate recognition molecules, are discussed.

Interaction of ovary lectin with homologous sperm from Bufo arenarum.

A soluble beta-galactoside lectin purified from Bufo arenarum ovary agglutinated homologous neuraminidase-treated spermatozoa. Microscopic observations of sperm clusters showed that spermatozoa agglutinated in a random way, but the head-to-head type of sperm agglutination was the most common (94-98%). The lectin activity was specifically inhibited by D-galactose and its derivatives, thio-digalactoside being the most active saccharide inhibitor.

Purification and partial biochemical characterization of an S-type lectin from blastula embryos of Bufo arenarum.

1. S-type lectin from Bufo arenarum embryos at blastula stage was purified by affinity chromatography. The molecule is a dimer with equal-sized monomers and the apparent subunit molecular weight was found to be 14.5 kDa. 2. Analytical isoelectric focusing of the pure lectin showed an acidic pI of 4.7. 3. Inhibition of the hemagglutination by mono- and oligosaccharides revealed a specificity for sugars bearing a beta-galactoside configuration. 4. Crossreactivity studies between the blastula lectin and the one purified earlier from adult ovary performed by immunodotting, ELISA and immunoblotting showed that these lectins share many epitopes.

A developmentally regulated lectin in Bufo arenarum embryos.

We report the levels of an endogenous beta-galactoside lectin activity from Bufo arenarum whole embryos extracts and specific inhibition by saccharides at different developmental stages. Specific activity measured against trypsinized rabbit red blood cells showed relatively high and fluctuating levels during early stages (up to about 76 h post-fertilization) which fell to significantly lower and more constant values at late stages (77-264 h post-fertilization). Lactose is the most potent inhibitor of this lectin activity, and saccharides having alpha-galactoside configurations are weaker inhibitors. At the last embryonic stage, the agglutinating activity showed a different sugar specificity which suggests either the modification of the preexistent lectin or the synthesis of another type of lectin. The possible physiological roles of these lectins in the blockage of polyspermy or in embryonic cell-cell interactions are discussed.

Galectins: major signaling modulators inside and outside the cell.

Galectins control cell behavior by acting on different signaling pathways. Most of the biological activities ascribed to these molecules rely upon recognition of extracellular glycoconjugates and establishment of multivalente interactions, which trigger adaptive biological responses. However, galectins are also detected within the cell in different compartments, where their regulatory functions still remain poorly understood. A deeper understanding of the entire galectin signalosome and its impact in cell behavior is therefore essential in order to delineate new strategies to specifically manipulate both galectin expression and function. This review summarizes our current knowledge of the signaling pathways activated by galectins, their glycan dependence and the cellular compartment where they become activated and are biologically relevant.