Aqueous imidazole solutions: a structural perspective from simulations with high-rank electrostatic multipole moments.

Imidazole is a small but important molecule occurring as a structure fragment in systems from amino acids, over ionic liquids, to synthetic polymers. Here we focus on the structure and dynamics of imidazole in water at ambient conditions, using both radial and spatial distribution functions. Molecular dynamics simulations were carried out for various imidazole concentrations, using a traditional point-charge potential and a high-rank multipolar potential. The difference in the description of the electrostatics leads to sizable quantitative differences (e.g., the diffusion coefficient) but also qualitative differences in the local structure. In contrast to a point-charge potential, the multipolar potential favors hydrogen-bonded chainlike imidazole dimers over stacked dimers.

Calculation of Raman optical activity spectra of methyl-ß-D-glucose incorporating a full molecular dynamics simulation of hydration effects.

We report calculations of the Raman and Raman optical activity (ROA) spectra of methyl-ß-D-glucose utilizing density functional theory combined with molecular dynamics (MD) simulations to provide an explicit hydration environment. This is the first report of such combination of MD simulations with ROA ab initio calculations. We achieve a significant improvement in accuracy over the more commonly used gas phase and polarizable continuum model (PCM) approaches, resulting in an excellent level of agreement with the experimental spectrum. Modeling the ROA spectra of carbohydrates has until now proven a notoriously difficult challenge due to their sensitivity to the effects of hydration on the molecular vibrations involving each of the chiral centers. The details of the ROA spectrum of methyl-ß-D-glucose are found to be highly sensitive to solvation effects, and these are correctly predicted for the first time including those originating from the highly sensitive low frequency vibrational modes. This work shows that a thorough consideration of the role of water is pivotal for understanding the vibrational structure of carbohydrates and presents a new and powerful tool for characterizing carbohydrate structure and conformational dynamics in solution.

The dynamic behavior of a liquid ethanol-water mixture: a perspective from quantum chemical topology.

Quantum Chemical Topology (QCT) is used to reveal the dynamics of atom-atom interactions in a liquid. A molecular dynamics simulation was carried out on an ethanol-water liquid mixture at its azeotropic concentration (X(ethanol)=0.899), using high-rank multipolar electrostatics. A thousand (ethanol)(9)-water heterodecamers, respecting the water-ethanol ratio of the azeotropic mixture, were extracted from the simulation. Ab initio electron densities were computed at the B3LYP/6-31+G(d) level for these molecular clusters. A video shows the dynamical behavior of a pattern of bond critical points and atomic interaction lines, fluctuating over 1 ns. A bond critical point distribution revealed the fluctuating behavior of water and ethanol molecules in terms of O-H···O, C-H···O and H···H interactions. Interestingly, the water molecule formed one to six C-H···O and one to four O-H···O interactions as a proton acceptor. We found that the more localized a dynamical bond critical point distribution, the higher the average electron density at its bond critical points. The formation of multiple C-H···O interactions affected the shape of the oxygen basin of the water molecule, which is shown in three dimensions. The hydrogen atoms of water strongly preferred to form H···H interactions with ethanol's alkyl hydrogen atoms over its hydroxyl hydrogen.

Simulation of liquid imidazole using a high-rank quantum topological electrostatic potential.

Rigid body molecular dynamics simulations were carried out on pure liquid imidazole at four different temperatures and at 1 atm. Imidazole, which is important both in life science and materials science, is one of the simplest molecules to possess both a lone pair and a π system. These two features are known to benefit from multipolar electrostatics. Here the electrostatic interaction is governed by atomic multipole moments obtained from topologically partitioned ab initio electron densities. The non-electrostatic terms are modeled with Lennard-Jones parameters adjusted to fit the experimental liquid density. All σ values are incrementally increased by one single scaling factor. We report on how the presence of multipolar electrostatics influences the local structure, dynamics and thermodynamics of the liquid compared to electrostatics by atomic point charges. The point charge force field exaggerates the number of π-stacked dimers in the liquid, and underestimates the number of hydrogen-bonded dimers. The effect of the temperature on the local structure of liquid imidazole was analysed using radial and spatial distribution functions.

An investigation of the conductivity of peptide nanotube networks prepared by enzyme-triggered self-assembly.

We demonstrate that nanotubular networks formed by enzyme-triggered self-assembly of Fmoc-L3 (9-fluorenylmethoxycarbonyl-tri-leucine) show significant charge transport. FT-IR, fluorescence spectroscopy and wide angle X-ray scattering (WAXS) data confirm formation of beta-sheets that are locked together viapi-stacking interactions. Molecular dynamics simulations confirmed the pi-pi stacking distance between fluorenyl groups to be 3.6-3.8 A. Impedance spectroscopy demonstrated that the nanotubular xerogel networks possess minimum sheet resistances of 0.1 MOmega/sq in air and 500 MOmega/sq in vacuum (pressure: 1.03 mbar) at room temperature, with the conductivity scaling linearly with the mass of peptide in the network. These materials may provide a platform to interface biological components with electronics.

Properties of liquid water from a systematic refinement of a high-rank multipolar electrostatic potential.

We build on previous work [S. Y. Liem and P. L. A. Popelier, J. Chem. Theory Comput. 4, 353 (2008)], where for the first time, a high-rank multipolar electrostatic potential was used in molecular dynamics simulations of liquid water at a wide range of pressures and temperatures, and using a multipolar Ewald summation. Water is represented as a rigid body, with atomic multipole moments defined by quantum chemical topology partitioning its gas phase electron density. The effect of the level of theory on the local structure of liquid water is systematically addressed. Values for Lennard-Jones (LJ) parameters are optimized, for both oxygen and hydrogen atoms, against bulk properties. The best LJ parameters were then used in a set of simulations at 30 different temperatures (1 atm) and another set at 11 different pressures (at 298 K). Inclusion of the hydrogen LJ parameters significantly increases the self-diffusion coefficient. The behavior of bulk properties was studied and the local water structure analyzed by both radial and spatial distribution functions. Comparisons with familiar point-charge potentials, such as TIP3P, TIP4P, TIP5P, and simple point charge, show the benefits of multipole moments.

Calculation of the relative acidities and oxidation potentials of para-substituted phenols. A model for alpha-tocopherol in solution.

Relative acidities (Delta pK(a)) of phenols and oxidation potentials (Delta E(ox)) of the phenoxide anions have been calculated for nine para-substituted phenols using density functional theory. Solvent effects were incorporated using the conductor-like polarisable continuum method. Using the calculated Delta pK(a) and Delta E(ox) values in a thermodynamic cycle, the DeltaBDE (bond dissociation enthalpy) of the phenols were also determined with all values calculated to within 1.5 kcal mol(-1) of experiment. The Delta pK(a) and Delta E(ox) values were calculated for 6-hydroxy-2,2,5,7,8-pentamethylchroman (HPMC), a model for alpha-tocopherol for which there are no known experimental values. The acidity of this compound is raised by 2.4 pK(a) units and lowered by -0.79 V relative to phenol with a calculated Delta BDE of -14.9 kcal mol(-1). There is a negative correlation (r(2) = 0.86) between the Delta pK(a) and the Delta BDE values. A stronger and positive correlation is found between the Delta E(ox) (r(2) = 0.98) and the Delta BDE values. Using these correlations it is uncovered that hydrogen abstraction of phenols, as measured by the Delta BDE, is driven by electron transfer rather than by proton transfer.