Molecular dynamics of confined glucose solutions

Glucose-Promoted Localization Dynamics of Excess Electrons in Aqueous Glucose Solution Revealed by Ab Initio Molecular Dynamics Simulation

Ab initio molecular dynamics simulations reveal that an excess electron (EE) can be more efficiently localized as a cavity-shaped state in aqueous glucose solution (AGS) than in water. Compared with that (∼1.5 ps) in water, the localization time is shortened by ∼0.7-1.2 ps in three AGSs (0.56, 1.12, and 2.87 M). Although the radii of gyration of the solvated EEs are all close to 2.6 Å in the four solutions, the solvated EE cavities in the AGSs become more compact and can localize ∼80% of an EE, which is considerably larger than that (∼40-60% and occasionally ∼80%) in water. These observations are attributed to a modification of the hydrogen-bonded network by the introduction of glucose molecules into water. The water acts as a promoter and stabilizer, by forming voids around glucose molecules and, in this fashion, favoring the localization of an EE with high efficiency. This study provides important information about EEs in physiological AGSs and suggests a new strategy to efficiently localize an EE in a stable cavity for further exploration of biological function.

Molecular dynamics and neutron scattering study of glucose solutions confined in MCM-41

Glucose aqueous solutions confined in MCM-41 cylindrical pores of diameter 3.2 nm have been studied by molecular dynamics (MD) simulations and quasielastic neutron scattering (QENS). MD simulations reveal a strong preferential interaction of glucose molecules with the silica walls, which induces significant concentration gradients within the pore. The influence of glucose on the structural and dynamical properties of water strongly depends on the region of the pore considered. The distortion of the hydrogen bond network (HBN) and of the tetrahedral organization of interfacial water molecules induced by silica is much stronger than that induced by glucose molecules. The interfacial glucose molecules diffuse about 1 order of magnitude slower than those in the core region. Differences in affinities for silica of the different species in confined hydrogen-bonded mixtures induce significant structural and dynamical heterogeneities not present in bulk solutions.

Confinement effects on monosaccharide transport in nanochannels

Transport theories based on the continuum hypothesis may not be appropriate at the nanoscale in view of surface effects. We employed molecular dynamics simulations to study the effects of confinement and concentration on diffusive transport of glucose in silica nanochannels (10 nm or smaller). We found that glucose modifies the electrical properties of nanochannels and that, below 5 nm in channel height, glucose adsorption and diffusivity are significantly reduced. With increasing concentration, the diffusivity is reduced linearly in the bulk, while it is reduced nonlinearly at the interface. The effective diffusivity reduction is related to the interface thickness, which can be 2-4 nm depending on concentration, and has an unexpected reduction at low concentrations. Results suggest that nanochannels present a one-dimensional cage environment that affects diffusivity in a fashion similar to cage-breaking diffusion. Our simulation results, consistent with the experimental observations presented here, suggest that nanoconfinement is the essential cause of the observed altered fluid diffusive transport, not accounted for by classical theories, because of coupling of confinement and concentration effects.

Translational and rotational dynamics of monosaccharide solutions

Molecular dynamics computer simulations have been carried out on aqueous solutions of glucose at concentrations bracketing those previously measured with quasi-elastic neutron scattering (QENS), in order to investigate the motions and interactions of the sugar and water molecules. In addition, QENS measurements have been carried out on fructose solutions to determine whether the effects previously observed for glucose apply to monosaccharide solutions. The simulations indicate a dynamical analogy between higher solute concentration and lower temperature that could provide a key explanation of the bioprotective phenomena observed in many living organisms. The experimental results on fructose solutions show qualitatively similar behavior to the glucose solutions. The dynamics of the water molecules are essentially the same, while the translational diffusion of the sugar molecules is slightly faster in the fructose solutions.

Characterization of sol-gel-encapsulated proteins using small-angle neutron scattering

Entrapment of biomolecules in silica-derived sol-gels has grown into a vibrant area of research since it was originally demonstrated. However, accessing the consequences of entrapment on biomolecules and the gel structure remains a major challenge in characterizing these biohybrid materials. We present the first demonstration that it is possible with small-angle neutron scattering (SANS) to study the conformation of dilute proteins that are entrapped in transparent and dense sol-gels. Using deuterium-labeled green fluorescent protein (GFP) and SANS with contrast variation, we demonstrate that the scattering signatures of the sol-gel and the protein can be separated. Analysis of the scattering curves of the sol-gels using a mass-fractal model shows that the size of the colloidal silica particles and the fractal dimensions of the gels were similar in the absence and presence of protein, demonstrating that GFP did not influence the reaction pathway for the formation of the gel. The major structural difference in the gels was that the pore size was increased 2-fold in the presence of the protein. At the contrast match point for silica, the scattering signal from GFP inside the gel became distinguishable over a wide q range. Simulated scattering curves representing a monomer, end-to-end dimer, and parallel dimer of the protein were calculated and compared to the experimental data. Our results show that the most likely structure of GFP is that of an end-to-end dimer. This approach can be readily applied and holds great potential for the structural characterization of complex biohybrid and other materials.

Dynamics of trehalose molecules in confined solutions

The dynamics of trehalose molecules in aqueous solutions confined in silica gel have been studied by quasielastic neutron scattering (QENS). Small-angle neutron scattering measurements confirmed the absence of both sugar clustering and matrix deformation of the gels, indicating that the results obtained are representative of homogeneous trehalose solutions confined in a uniform matrix. The pore size in the gel is estimated to be 18 nm, comparable to the distances in cell membranes. For the QENS measurements, the gel was prepared from D2O in order to accentuate the scattering from the trehalose. Values for the translational diffusion constant and effective jump distance were derived from model fits to the scattering function. Comparison with QENS and NMR results in the literature for bulk trehalose shows that confinement on a length scale of 18 nm has no significant effect on the translational diffusion of trehalose molecules.

Carbohydrate Clustering in Aqueous Solutions and the Dynamics of Confined Water

We use molecular dynamics simulations to investigate structure and dynamics of fructose aqueous solutions in the 1-5 M concentration range at ambient conditions. We analyze hydration structures, H-bond statistics, and size distribution of H-bonded carbohydrate clusters as functions of concentration. We find that the local tetrahedral order of water is reasonably well-preserved and that the solute tends to appear as scattered "isolated" molecules at low concentrations and as H-bonded clusters for less diluted solutions. The sugar cluster size distribution exhibits a sharp transition to a percolated cluster between 3.5 and 3.8 M. The percolated cluster forms an intertwined network of H-bonded saccharides that imprisons water. For the dynamics, we find good agreement between simulation and available experimental results for the self-diffusion coefficients. Water librational dynamics is little affected by sugar concentration, whereas reorientational relaxation is described by a concentration-independent bulk-like component attributed to noninterfacial water molecules and a slower component (strongly concentration dependent) that arises from interfacial solvent molecules and, hence, depends on the dynamics of the cluster structure itself. Analysis of H-bonding survival probability functions indicates that the formation of carbohydrate clusters upon increasing concentration enhances the H-bond relaxation time and slows down the entire system dynamics. We find that multiexponential or stretched-exponential fits alone cannot describe the H-bond survival probabilities for the entire postlibrational time span of our data (0.1-100 ps), as opposed to a combined stretched-plus-biexponential function, which provides excellent fits. Our results suggest that water dynamics in concentrated fructose solutions resembles in many ways that of protein hydration water.