Host-guest interactions between xanthones and water: the role of O-H center dot center dot center dot O, C-H center dot center dot center dot O, and pi center dot center dot center dot pi contacts in the channel- and cage-type frameworks

In this article were studied two xanthone derivatives known as 1,5-dihydroxy-8-methoxyxanthone (I) and 1,3,7-trihydroxy-8-methoxyxanthone (II), which show one water molecule into their crystal structures. In xanthone I, there are water wires contributing to build up channel-like cavities along the c axis, whereas in xanthone II the water is surrounded by three xanthone molecules forming a cage-type structure. The geometries of I and II were optimized using the density functional theory method with B3LYP functional, and the results were compared with crystal structure. Both theoretical and experimental investigations reveal a concordance between structural parameters, with the xanthone core presenting an almost flat conformation and substituents adopting the more stable orientations. In the two compounds, the hydroxyl group linked at position 1 is involved in a resonance-assisted hydrogen bond with the carbonyl group. Besides, the supramolecular arrangement of the host/guest systems are stabilized mainly by classical intermolecular hydrogen bonds (O-H center dot center dot center dot O) involving xanthone-to-water and xanthone-to-xanthone. In addition, C-H center dot center dot center dot O weak hydrogen bonds, as well as pi-pi interactions play an important role to stabilize the crystal self-assembly of xanthones I and II. The results reported here underline the role of inclusion of water molecules and their different arrangement into the crystal structure of two xanthone host/guest systems.

Hydrogen bond and the resonance effect on the formamide-water complexes

The interaction of formamide and the two transition states of its amide group rotation with one, two, or three water molecules was studied in vacuum. Great differences between the electronic structure of formamide in its most stable form and the electronic structure of the transition states were noticed. Intermolecular interactions were intense, especially in the cases where the solvent interacted with the amide and the carbonyl groups simultaneously. In the transition states, the interaction between the lone pair of nitrogen and the water molecule becomes important. With the aid of the natural bond orbitals, natural resonance theory, and electron localization function (ELF) analyses an increase in the resonance of planar formamide with the addition of successive water molecules was observed. Such observation suggests that the hydrogen bonds in the formamidewater complexes may have some covalent character. These results are also supported by the quantitative ELF analyses. (C) 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012

A simple analysis of the influence of the solvent-induced electronic polarization on the N-15 magnetic shielding of pyridine in water

Electronic polarization induced by the interaction of a reference molecule with a liquid environment is expected to affect the magnetic shielding constants. Understanding this effect using realistic theoretical models is important for proper use of nuclear magnetic resonance in molecular characterization. In this work, we consider the pyridine molecule in water as a model system to briefly investigate this aspect. Thus, Monte Carlo simulations and quantum mechanics calculations based on the B3LYP/6-311++G (d,p) are used to analyze different aspects of the solvent effects on the N-15 magnetic shielding constant of pyridine in water. This includes in special the geometry relaxation and the electronic polarization of the solute by the solvent. The polarization effect is found to be very important, but, as expected for pyridine, the geometry relaxation contribution is essentially negligible. Using an average electrostatic model of the solvent, the magnetic shielding constant is calculated as -58.7 ppm, in good agreement with the experimental value of -56.3 ppm. The explicit inclusion of hydrogen-bonded water molecules embedded in the electrostatic field of the remaining solvent molecules gives the value of -61.8 ppm.

Atomistic Study of Water Confined in Silica.

In this work, we have used a combined of atomistic simulation methods to explore the effects of confinement of water molecules between silica surfaces. Firstly, the mechanical properties of water severe confined (~3A) between two silica alpha-quartz was determined based on first principles calculations within the density functional theory (DFT). Simulated annealing methods were employed due to the complex potential energry surface, and the difficulties to avoid local minima. Our results suggest that much of the stiffness of the material (46%) remains, even after the insertion of a water monolayer in the silica. Secondly, in order to access typical time scales for confined systems, classical molecular dynamics was used to determine the dynamical properties of water confined in silica cylindrical pores, with diameters varying from 10 to 40A. in this case we have varied the passivation of the silica surface, from 13% to 100% of SiOH, and the other terminations being SiOH2 and SiOH3, the distribution of the different terminations was obtained with a Monte Carlo simulation. The simulations indicates a lowering of the diffusion coefficientes as the diameter decreases, due to the structuration of hydrogen bonds of water molecules; we have also obtained the density profiles of the confined water and the interfacial tension.

Computer simulation study of collective dynamics in the glass former Ca(NO3)(2)center dot 4H(2)O

Time correlation functions of current fluctuations were calculated by molecular dynamics (MD) simulations in order to investigate sound waves of high wavevectors in the glass-forming liquid Ca(NO3)(2)center dot 4H(2)O. Dispersion curves, omega(k), were obtained for longitudinal (LA) and transverse acoustic (TA) modes, and also for longitudinal optic (LO) modes. Spectra of LA modes calculated by MD simulations were modeled by a viscoelastic model within the memory function framework. The viscoelastic model is used to rationalize the change of slope taking place at k similar to 0.3 angstrom(-1) in the omega(k) curve of acoustic modes. For still larger wavevectors, mixing of acoustic and optic modes is observed. Partial time correlation functions of longitudinal mass currents were calculated separately for the ions and the water molecules. The wavevector dependence of excitation energies of the corresponding partial LA modes indicates the coexistence of a relatively stiff subsystem made of cations and anions, and a softer subsystem made of water molecules. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4751548]

Environmental effects in kelvin force microscopy of modified diamond surfaces

We have explored the effects of atmospheric environment on Kelvin force microscopy (KFM) measurements of potential difference between different regions of test polycrystalline diamond surfaces. The diamond films were deposited by microwave plasma-assisted chemical vapor deposition, which naturally produces hydrogen terminations on the surface of the films formed. Selected regions were patterned by electron-beam lithography and chemical terminations of oxygen or fluorine were created by exposure to an oxygen or fluorine plasma source. For KFM imaging, the samples were mounted in a hood with a constant flow of helium gas. Successive images were taken over a 5-h period showing the effect of the environment on KFM imaging. We conclude that the helium flow removes water molecules adsorbed on the surface of the samples, resulting in differences in surface potential between adjacent regions. The degree of water removal is different for surfaces with different terminations. The results highlight the importance of taking into account the atmospheric environment when carrying out KFM analysis. (C) 2012 Wiley Periodicals, Inc.

From pyridoxalrhodanine to a novel 7-azacoumarin complex of dimethylthallium(III) in a one-pot synthesis

The hydrolysis of pyridoxalrhodanine in a basic medium containing the dimethylthallium(III) cation afforded the compound [TlMe2(L)]center dot H2O (1.H2O) [HL = 5-(hydroxymethyl)-8-methyl-3-thiol-7-azacoumarin]. This compound was characterized in solid state by IR spectroscopy and in solution by H-1 and C-13{H-1} NMR spectrometry. X-ray diffraction showed that the crystal consists of associated TlMe2(L) units and hydrogen bonded water molecules. The L- anion is bound to the metal mainly by a bridging S atom [Tl-S = 2.9458(18) angstrom; 2.9616(16) angstrom], although secondary interactions through O atoms (Tl-O: 2.861(5); 2.900(5) angstrom)] are also present. The longer Tl-O interaction and the hydrogen bonds of the water molecules give rise to a tridimensional polymeric structure. (C) 2012 Elsevier B.V. All rights reserved.

Synthesis, crystal structure and magnetic properties of the helical oxalate-bridged copper(II) chain {[(CH3)(4)N](2)[Cu(C2O4)(2)] center dot H2O}(n)

The preparation, crystal structure and magnetic properties of a new oxalate-containing copper(II) chain of formula {[(CH3)(4)N](2)]Cu(C2O4)(2)] center dot H2O}(n) (1) [(CH3)(4)N+ = tetramethylammonium cation] are reported. The structure of 1 consists of anionic oxalate-bridged copper(II) chains, tetramethylammoniun cations and crystallization water molecules. Each copper(II) ion in 1 is surrounded by three oxalate ligands, one being bidentate and the other two exhibiting bis-bidenate coordination modes. Although all the tris-chelated copper(H) units from a given chain exhibit the same helicity, adjacent chains have opposite helicities and then an achiral structure results. Variable-temperature magnetic susceptibility measurements of 1 show the occurrence of a weak ferromagnetic interaction through the oxalate bridge [J = +1.14(1)cm(-1), the Hamiltonian being defined as H = -J Sigma nm S-i . S-j]. This value is analyzed and discussed in the light of available magnetostructural data for oxalate-bridged copper(H) complexes with the same out-of-plane exchange pathway. (C) 2012 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.

Structure and electronic properties of hydrated mesityl oxide: a sequential quantum mechanics/molecular mechanics approach

The hydration of mesityl oxide (MOx) was investigated through a sequential quantum mechanics/molecular mechanics approach. Emphasis was placed on the analysis of the role played by water in the MOx syn-anti equilibrium and the electronic absorption spectrum. Results for the structure of the MOx-water solution, free energy of solvation and polarization effects are also reported. Our main conclusion was that in gas-phase and in low-polarity solvents, the MOx exists dominantly in syn-form and in aqueous solution in anti-form. This conclusion was supported by Gibbs free energy calculations in gas phase and in-water by quantum mechanical calculations with polarizable continuum model and thermodynamic perturbation theory in Monte Carlo simulations using a polarized MOx model. The consideration of the in-water polarization of the MOx is very important to correctly describe the solute-solvent electrostatic interaction. Our best estimate for the shift of the pi-pi* transition energy of MOx, when it changes from gas-phase to water solvent, shows a red-shift of -2,520 +/- 90 cm(-1), which is only 110 cm(-1) (0.014 eV) below the experimental extrapolation of -2,410 +/- 90 cm(-1). This red-shift of around -2,500 cm(-1) can be divided in two distinct and opposite contributions. One contribution is related to the syn -> anti conformational change leading to a blue-shift of similar to 1,700 cm(-1). Other contribution is the solvent effect on the electronic structure of the MOx leading to a red-shift of around -4,200 cm(-1). Additionally, this red-shift caused by the solvent effect on the electronic structure can by composed by approximately 60 % due to the electrostatic bulk effect, 10 % due to the explicit inclusion of the hydrogen-bonded water molecules and 30 % due to the explicit inclusion of the nearest water molecules.

Salts of the anti-HIV drug lamivudine with phthalic and salicylic acids

Salts of the anti-HIV drug lamivudine, with phthalic acid and salicylic acid as counterions, were investigated in this study. Neither the packing of the (lamivudine)(+)(phthalic acid)(-) ion pairs nor the conformation of the lamivudine moiety itself were similar to those found in other multicomponent molecular salts of the drug, such as hydrogen maleate and saccharinate ones, even though all three salts crystallize in the same P2(1)2(1)2(1) orthorhombic space group with similar unit cell metrics. Lamivudine salicylate assumes a different crystal structure to those of the hydrogen maleate and saccharinate salts, crystallizing in the P2(1) monoclinic space group as a monohydrate whose (lamivudine)(+)(salicylic acid)(-) ion pair is assembled through two hydrogen bonds with cytosine as a dual donor to both oxygens of the carboxylate, such as in the pairing of lamivudine with a phthalic acid counterion. In lamivudine salicylate monohydrate, the drug conformation is related to the hydrogen maleate and saccharinate salts. However, such a conformational similarity is not related to the intermolecular interaction patterns. Lamivudine and water molecules alternate into helical chains in the salicylate salt monohydrate.

Enzyme Microheterogeneous Hydration and Stabilization in Supercritical Carbon Dioxide

Supercritical carbon dioxide is a promising green-chemistry solvent for many enzyme-catalyzed chemical reactions, yet the striking stability of some enzymes in such unconventional environments is not well understood. Here, we investigate the stabilization of the Candida antarctica Lipase B (CALB) in supercritical carbon dioxide-water biphasic systems using molecular dynamics simulations. The preservation of the enzyme structure and optimal activity depend on the presence of small amounts of water in the supercritical dispersing medium. When the protein is at least partially hydrated, water molecules bind to specific sites on the enzyme surface and prevent carbon dioxide from penetrating its catalytic core. Strikingly, water and supercritical carbon dioxide cover the protein surface quite heterogeneously. In the first solvation layer, the hydrophilic residues at the surface of the protein are able to pin down patches of water, whereas carbon dioxide solvates preferentially hydrophobic surface residues. In the outer solvation shells, water molecules tend to cluster predominantly on top of the larger water patches of the first solvation layer instead of spreading evenly around the remainder of the protein surface. For CALB, this exposes the substrate-binding region of the enzyme to carbon dioxide, possibly facilitating diffusion of nonpolar substrates into the catalytic funnel. Therefore, by means of microheterogeneous solvation, enhanced accessibility of hydrophobic substrates to the active site can be achieved, while preserving the functional structure of the enzyme. Our results provide a molecular picture on the nature of the stability of proteins in nonaqueous media.

Molecular Dynamics Investigations of PRODAN in a DLPC Bilayer

Molecular dynamics computer simulations have been performed to identify preferred positions of the fluorescent probe PRODAN in a fully hydrated DLPC bilayer in the fluid phase. In addition to the intramolecular charge-transfer first vertical excited state, we considered different charge distributions for the electronic ground state of the PRODAN molecule by distinct atomic charge models corresponding to the probe molecule in vacuum as well as polarized in a weak and a strong dielectric solvent (cyclohexane and water). Independent on the charge distribution model of PRODAN, we observed a preferential orientation of this molecule in the bilayer with the dimethylamino group pointing toward the membrane's center and the carbonyl oxygen toward the membrane's interface. However, changing the charge distribution model of PRODAN, independent of its initial position in the equilibrated DLPC membrane, we observed different preferential positions. For the ground state representation without polarization and the in-cyclohexane polarization, the probe maintains its position close to the membrane's center. Considering the in-water polarization model, the probe approaches more of the polar headgroup region of the bilayer, with a strong structural correlation with the choline group, exposing its oxygen atom to water molecules. PRODAN's representation of the first vertical excited state with the in-water polarization also approaches the polar region of the membrane with the oxygen atom exposed to the bilayer's hydration shell. However, this model presents a stronger structural correlation with the phosphate groups than the ground state. Therefore, we conclude that the orientation of the PRODAN molecule inside the DLPC membrane is well-defined, but its position is very sensitive to the effect of the medium polarization included here by different models for the atomic charge distribution of the probe.

Effects of moisture content on structural and dielectric properties of cassava starch films

In this work, it was evaluated the effect of moisture content on the structural and dielectric properties of cassava starch films by means of Fourier Transform infrared spectroscopy (FTIR), impedancimetric, and gravimetric analysis. The film samples were equilibrated in hermetically sealed desiccators, containing different saturated salts in water in order to promote activity of water between 0.11 and 0.85. The position and amplitude of the peaks in the fingerprint region of the FTIR spectra were changed due to the modifications in the interactions between the polymeric chains and water molecules. These effects may be related to the formation of semi-crystalline regions in the film structure. The dielectric properties of the films were also strongly dependent on the moisture content, showing a non-linear and a linear region, which was attributed to the domain of bound and free water in the film, respectively. The gravimetric analyzes showed the typical sigmoidal behavior, attributed to the way the water interacts with the biopolymer. Finally, the flexibility of the films increased with water content increasing.

Electron collisions with the HCOOH...(H2O)n (n=1, 2 and 3) complexes in liquid phase: the influence on the π* shape resonance of formic acid

In this conference we report cross sections for elastic collisions of low-energy electrons with the HCOOH…(H2O)n complexes, with n = 1, 2 and 3. The scattering cross sections were computed with the Schwinger multichannel method [K. Takatsuka and V. McKoy, Phys. Rev. A 24 , 2473 (1981); Phys. Rev. A 30 , 1734 (1984)] with pseudopotentials [M. H. F. Bettega, L. G. Ferreira, and M. A. P. Lima, Phys. Rev. A 47, 1111 (1993)] in the static-exchange and static-exchange plus polarization approximations, for energies from 0.5 eV to 6 eV. We considered some diÆerent hydrogen-bonded structures for the complexes that were generated with classical Monte Carlo simulations [K. Coutinho and S. Canuto, J. Chem. Phys. 113, 9132, (2000)]. The aim of this work is to investigate the effect of the surrounding water molecules on the π* shape resonance of the solute. Previous theoretical and experimental studies carried out in the gas phase reported a π* state for HCOOH at around 1.9 eV. For the n = 1 case and for all complexes, the stabilization of the resonance was observed (it appears at lower energy compared to the value obtained in the gas phase), as reported previously for the CH2O…H2O complexes [T. C. Freitas, M. A. P. Lima, S. Canuto, and M. H. F. Bettega, Phys. Rev. A 80, 062710 (2009)]. This result indicates that the presence of the solvent may affect the processes related to the π* state, such as the molecular dissociation by electron impact. For the n = 2 case we have observed both stabilization and destabilization of the π* resonance, that is associated with the hydrogen bond donor or acceptor role of the water molecules in the complexes. For the n = 3 case, preliminary static-exchange results show the stabilization of the π* state. We propose an explanation of the stabilization/destabilization of the π* state in terms of the polarization of the solute due to the surrounding water molecules and the net charge in the solute.

Applying many-body expansion aiming to make ab in initio molecular dynamics more efficient

In this present work we present a methodology that aims to apply the many-body expansion to decrease the computational cost of ab initio molecular dynamics, keeping acceptable accuracy on the results. We implemented this methodology in a program which we called ManBo. In the many-body expansion approach, we partitioned the total energy E of the system in contributions of one body, two bodies, three bodies, etc., until the contribution of the Nth body [1-3]: E = E1 + E2 + E3 + …EN. The E1 term is the sum of the internal energy of the molecules; the term E2 is the energy due to interaction between all pairs of molecules; E3 is the energy due to interaction between all trios of molecules; and so on. In Manbo we chose to truncate the expansion in the contribution of two or three bodies, both for the calculation of the energy and for the calculation of the atomic forces. In order to partially include the many-body interactions neglected when we truncate the expansion, we can include an electrostatic embedding in the electronic structure calculations, instead of considering the monomers, pairs and trios as isolated molecules in space. In simulations we made we chose to simulate water molecules, and use the Gaussian 09 as external program to calculate the atomic forces and energy of the system, as well as reference program for analyzing the accuracy of the results obtained with the ManBo. The results show that the use of the many-body expansion seems to be an interesting approach for reducing the still prohibitive computational cost of ab initio molecular dynamics. The errors introduced on atomic forces in applying such methodology are very small. The inclusion of an embedding electrostatic seems to be a good solution for improving the results with only a small increase in simulation time. As we increase the level of calculation, the simulation time of ManBo tends to largely decrease in relation to a conventional BOMD simulation of Gaussian, due to better scalability of the methodology presented. References [1] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 3, 46 (2007). [2] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 4, 1 (2008). [3] R. Rivelino, P. Chaudhuri and S. Canuto; J. Chem. Phys., 118, 10593 (2003).