Counterion N.M.R. in heterogeneous aqueous systems

Quadrupolar 23Na+ NMR relaxation as a probe of subpicosecond collective dynamics in aqueous electrolyte solutions

Nuclear magnetic resonance relaxometry represents a powerful tool for extracting dynamic information. Yet, obtaining links to molecular motion is challenging for many ions that relax through the quadrupolar mechanism, which is mediated by electric field gradient fluctuations and lacks a detailed microscopic description. For sodium ions in aqueous electrolytes, we combine ab initio calculations to account for electron cloud effects with classical molecular dynamics to sample long-time fluctuations, and obtain relaxation rates in good agreement with experiments over broad concentration and temperature ranges. We demonstrate that quadrupolar nuclear relaxation is sensitive to subpicosecond dynamics not captured by previous models based on water reorientation or cluster rotation. While ions affect the overall water retardation, experimental trends are mainly explained by dynamics in the first two solvation shells of sodium, which contain mostly water. This work thus paves the way to the quantitative understanding of quadrupolar relaxation in electrolyte and bioelectrolyte systems.

Reverse Micelles As Antioxidant Carriers: An Experimental and Molecular Dynamics Study

Water-in-oil microemulsions with biocompatible components were formulated to be used as carriers of natural antioxidants, such as hydroxytyrosol (HT) and gallic acid (GA). The system was composed of a mixture of natural surfactants, lecithin and monoglycerides, medium chain triglycerides, and aqueous phase. A dual approach was undertaken to study the structure and dynamics of these complicated systems. First, experimental data were collected by using adequate techniques, such as dynamic light scattering (DLS) and electron paramagnetic resonance (EPR) spectroscopy. Following this, a coarse-grained molecular dynamics (CGMD) study based on the experimental composition using the MARTINI force field was conducted. The simulations revealed the spontaneous formation of reverse micelles (RMs) starting from completely random initial conformations, underlying their enhanced thermodynamic stability. The location of the bioactive molecules, as well as the structure of the RM, were in accordance with the experimental findings. Furthermore, GA molecules were found to be located inside the water core, in contrast to the HT ones, which seem to lie at the surfactant interfacial layer. The difference in the antioxidants' molecular location was only revealed in detail from the computational analysis and explains the RM's swelling observed by GA in DLS measurements.

On the microscopic fluctuations driving the NMR relaxation of quadrupolar ions in water

Nuclear Magnetic Resonance (NMR) relaxation is sensitive to the local structure and dynamics around the probed nuclei. The Electric Field Gradient (EFG) is the key microscopic quantity to understand the NMR relaxation of quadrupolar ions, such as (7)Li(+), (23)Na(+), (25)Mg(2+), (35)Cl(-), (39)K(+), or (133)Cs(+). Using molecular dynamics simulations, we investigate the statistical and dynamical properties of the EFG experienced by alkaline, alkaline Earth, and chloride ions at infinite dilution in water. Specifically, we analyze the effect of the ionic charge and size on the distribution of the EFG tensor and on the multi-step decay of its auto-correlation function. The main contribution to the NMR relaxation time arises from the slowest mode, with a characteristic time on the picosecond time scale. The first solvation shell of the ion plays a dominant role in the fluctuations of the EFG, all the more that the ion radius is small and its charge is large. We propose an analysis based on a simplified charge distribution around the ion, which demonstrates that the auto-correlation of the EFG, hence the NMR relaxation time, reflects primarily the collective translational motion of water molecules in the first solvation shell of the cations. Our findings provide a microscopic route to the quantitative interpretation of NMR relaxation measurements and open the way to the design of improved analytical theories for NMR relaxation for small ionic solutes, which should focus on water density fluctuations around the ion.

Signatures of nanoconfinement on the linear and nonlinear vibrational spectroscopy of a model hydrogen-bonded complex dissolved in a polar solvent

The one-dimensional IR (1D-IR) absorption and IR pump-probe spectra of a hydrogen stretch in a model hydrogen-bonded complex dissolved in a polar solvent confined in spherical hydrophobic cavities of different sizes were simulated using ground-state mixed quantum-classical dynamics. Due to a thorough analysis of key properties of the complex and solvent from equilibrium trajectory data, we were able to gain insight into the microscopic details underlying the spectra. Both the 1D-IR and IR pump-probe spectra manifested the effects of confinement on the relative stabilities of the covalent and ionic forms of the complex through pronounced changes in their peak intensities and numbers. However, in contrast to the 1D-IR spectra, the time-resolved pump-probe spectra were found to be uniquely sensitive to the changes in the molecular dynamics as the cavity size is varied. In particular, it was found that the variations in the time evolutions of the peak intensities in the pump-probe spectra reflect the differences in the solvation dynamics associated with the various forms of the complex in different locations within the cavities. The ability to detect these differences underscores the advantage of using pump-probe spectroscopy for studying nanoconfined systems.

Molecular simulation study of water mobility in aerosol-OT reverse micelles

In this work, we present results from molecular dynamics simulations on the single-molecule relaxation of water within reverse micelles (RMs) of different sizes formed by the surfactant aerosol-OT (AOT, sodium bis(2-ethylhexyl)sulfosuccinate) in isooctane. Results are presented for RM water content w(0) = [H(2)O]/[AOT] in the range from 2.0 to 7.5. We show that translational diffusion of water within the RM can, to a good approximation, be decoupled from the translation of the RM through the isooctane solvent. Water translational mobility within the RM is restricted by the water pool dimensions, and thus, the water mean-squared displacements (MSDs) level off in time. Comparison with models of diffusion in confined geometries shows that a version of the Gaussian confinement model with a biexponential decay of correlations provides a good fit to the MSDs, while a model of free diffusion within a sphere agrees less well with simulation results. We find that the local diffusivity is considerably reduced in the interfacial region, especially as w(0) decreases. Molecular orientational relaxation is monitored by examining the behavior of OH and dipole vectors. For both vectors, orientational relaxation slows down close to the interface and as w(0) decreases. For the OH vector, reorientation is strongly affected by the presence of charged species at the RM interface and these effects are especially pronounced for water molecules hydrogen-bonded to surfactant sites that serve as hydrogen-bond acceptors. For the dipole vector, orientational relaxation near the interface slows down more than that for the OH vector due mainly to the influence of ion-dipole interactions with the sodium counterions. We investigate water OH and dipole reorientation mechanisms by studying the w(0) and interfacial shell dependence of orientational time correlations for different Legendre polynomial orders.

Infrared spectra of a model phenol-amine proton transfer complex in nanoconfined CH3Cl

The vibrational spectra of a model phenol-amine proton transfer complex dissolved in CH3Cl solvent confined in a 12 A radius spherical hydrophobic cavity were calculated using mixed quantum-classical molecular dynamics simulations. The reaction free energy of the proton transfer complex was varied in order to explore the contributions to the vibrational absorption band from product and reactant species. The vibrational spectra of the model proton transfer complex resulted in motionally narrowed spectral linewidths with two distinct peaks for products and reactants in cases where the system undergoes chemical exchange. It was found that the n=1 and n=2 vibrational excited states combine to form diabatic states such that the spectra have contributions from both n=0 --> n=1 and n=0 --> n=2 transitions. A strong relationship between the instantaneous vibrational frequency and a collective solvent coordinate was found that assists in understanding the origin of the spectral features.

Dynamics of Water Confined Within Reverse Micelles

We report structural and dynamical properties of water confined within reverse micelles (RMs) ranging in size from R = 10 A to R = 23 A as determined from molecular dynamics simulations. The low-frequency infrared spectra have been calculated using linear response theory and depend linearly on the fraction of bulklike water within the RMs. Furthermore, these spectra show nearly isosbestic behavior in the region near 660 cm(-1). Both of these characteristics are present in previously measured experimental spectra. The single dipole spectra for interfacial trapped, bound, and bulklike water within the RMs have also been calculated and show region-dependent frequency shifts. Specifically, the bound and trapped water spectra have a peak at lower frequencies than that for the inner core water. We therefore assign the low-frequency band in the IR spectra to bound and trapped interfacial water. Finally, region-dependent hydrogen bonding profiles and spatial distribution functions are also presented.

Ion exchange in reverse micelles

The distribution and dynamics of alkali cations inside Na-AOT reverse micelles have been investigated using Monte Carlo and molecular dynamics simulations. Water is modeled using the extended simple point charge (SPC/E) model. Simulations were carried out for alkali salts of Li+, Na+, K+, and Cs+ placed into the aqueous core of the reverse micelle, for situations corresponding to one and three molecules of added salt. In all cases, we observe that the larger K+ and Cs+ ions exchange with the Na+ counterion; however, the smaller Li+ ion prefers to remains solvated within the core of the reverse micelle. Our study reveals that the oil-water interface of the Na-AOT reverse micelle has the greatest selectivity toward Cs+ followed by K+ and Li+. A model based on enthalpic contributions illustrates that the solvation energies of the different cations in water control the ion-exchange process. The hydration number of the first water shell for Li+ situated in the aqueous core of the reverse micelle with radius R = 14.1 A was similar to that observed at infinite dilution in bulk water.

Competitive Na(+) and Rb(+) binding in the minor groove of DNA

Sequence-dependent coordination of alkali ions to the nucleotide bases in the minor groove of AT-tract B-DNA has recently been inferred from X-ray crystallography, solution NMR and computer simulations. Here, we present new (23)Na and (87)Rb magnetic relaxation dispersion (MRD) data that demonstrate competitive and long-lived binding of Na(+) and Rb(+) ions in the minor groove of the B-DNA duplex [d(CGCGAATTCGCG)](2). The Na(+)/Rb(+) selectivity of the minor groove is found to be weak, consistent with local structural flexibility. The ion occupancies derived from the MRD data are substantially higher than previously reported, suggesting that groove-bound ions significantly influence the energetics and structural polymorphism of DNA in vivo. For example, in the presence of 0.20 M Na(+) and 0.56 M Rb(+) at 4 degrees C, the ApT site in the minor groove is occupied by a Rb(+) ion, a Na(+) ion, or a water molecule 40, 10, and 50% of the time, respectively. In the absence of Rb(+), the Na(+) occupancy increases to 50%. At 4 degrees C, the mean residence time of groove-bound ions is 0.2 +/- 0.1 micros for Rb(+) and 10 ns to 100 micros for Na(+). A shorter correlation time of 2 ns is attributed to counterions bridging cross-strand phosphate groups.

Sequence-specific binding of counterions to B-DNA

Recent studies by x-ray crystallography, NMR, and molecular simulations have suggested that monovalent counterions can penetrate deeply into the minor groove of B form DNA. Such groove-bound ions potentially could play an important role in AT-tract bending and groove narrowing, thereby modulating DNA function in vivo. To address this issue, we report here (23)Na magnetic relaxation dispersion measurements on oligonucleotides, including difference experiments with the groove-binding drug netropsin. The exquisite sensitivity of this method to ions in long-lived and intimate association with DNA allows us to detect sequence-specific sodium ion binding in the minor groove AT tract of three B-DNA dodecamers. The sodium ion occupancy is only a few percent, however, and therefore is not likely to contribute importantly to the ensemble of B-DNA structures. We also report results of ion competition experiments, indicating that potassium, rubidium, and cesium ions bind to the minor groove with similarly weak affinity as sodium ions, whereas ammonium ion binding is somewhat stronger. The present findings are discussed in the light of previous NMR and diffraction studies of sequence-specific counterion binding to DNA.

Accurate intensities of broad NMR lines from composite pulse experiments

Accurate determination of integral intensities of broad lines is difficult when spin relaxation during the applied pulses cannot be neglected and/or when ringing of the tank circuit interferes with the signal. Here we present an extension of the analytical solution of the generalized Bloch equations (G. A. Morris and P. B. Chilvers, J. Magn. Reson. A 107, 236 (1994)), which is then used to evaluate the signal intensity obtained in a composite pulse experiment designed to cancel ringing effects. Comparing intensities of broad and narrow (81)Br spectral lines tests and proves the accuracy of this approach.

Diffusion of solutes in agarose and alginate gels: 1H and 23Na PFGSE and 23Na TQF NMR studies

Cells immobilized in gels experience potential metabolic restrictions in the form of reduced diffusion rates of metabolites and ions and their possible selective adsorption on the gel matrix. Diffusion and relaxation characteristics of common solutes in agarose and barium alginate gels were investigated at 37 degrees C by using 1H PFGSE and 23Na TQF NMR spectroscopy. Glucose, glycine, alanine, lactate, sodium ions, and HDO were studied. There were no selective interactions between any of the metabolites and the gel materials but the diffusion coefficients were uniformly reduced. The effects of metabolite diffusion and utilization, in gel beads and threads containing cells, were simulated by using a reaction diffusion model incorporating the measured diffusion coefficients. Metabolism is expected to be very significantly limited by diffusion of solutes to and from the cells that are centrally located within gel threads or spheres of radius approximately 2.0 mm, which is a commonly used size.