Electronic couplings and on-site energies for hole transfer in DNA: systematic quantum mechanical/molecular dynamic study

The electron hole transfer (HT) properties of DNA are substantially affected by thermal fluctuations of the π stack structure. Depending on the mutual position of neighboring nucleobases, electronic coupling V may change by several orders of magnitude. In the present paper, we report the results of systematic QM/molecular dynamic (MD) calculations of the electronic couplings and on-site energies for the hole transfer. Based on 15 ns MD trajectories for several DNA oligomers, we calculate the average coupling squares 〈 V2 〉 and the energies of basepair triplets X G+ Y and X A+ Y, where X, Y=G, A, T, and C. For each of the 32 systems, 15 000 conformations separated by 1 ps are considered. The three-state generalized Mulliken-Hush method is used to derive electronic couplings for HT between neighboring basepairs. The adiabatic energies and dipole moment matrix elements are computed within the INDO/S method. We compare the rms values of V with the couplings estimated for the idealized B -DNA structure and show that in several important cases the couplings calculated for the idealized B -DNA structure are considerably underestimated. The rms values for intrastrand couplings G-G, A-A, G-A, and A-G are found to be similar, ∼0.07 eV, while the interstrand couplings are quite different. The energies of hole states G+ and A+ in the stack depend on the nature of the neighboring pairs. The X G+ Y are by 0.5 eV more stable than X A+ Y. The thermal fluctuations of the DNA structure facilitate the HT process from guanine to adenine. The tabulated couplings and on-site energies can be used as reference parameters in theoretical and computational studies of HT processes in DNA

Evaluation of the binding ability of a novel dioxatetraazamacrocyclic receptor that contains two phenanthroline units: selective uptake of carboxylate anions

The novel dioxatetraaza macrocycle [26]phen(2)N(4)O(2), which incorporates two phenanthroline units, has been synthesized, and its acid-base behavior has been evaluated by potentiometric and H-1 NMR methods. Six protonation constants were determined, and the protonation sequence was established by NMR. The location of the fifth proton on the phen nitrogen was confirmed by X-ray determinations of the crystal structures of the receptor as bromide and chloride salts. The two compounds have the general molecular formula {(H-5[26]phen(2)N(4)O(2))X-n(H2O)(5-n)}X(n-1)(.)mH(2)O, where X = Cl, n = 3, and m = 6 or X = Br, n = 4, and m = 5.5. In the solid state, the (H-5[26]phen(2)N(4)O(2))(5+) cation adopts a "horseshoe" topology with sufficient room to encapsulate three or four halogen anions through the several N-(HX)-X-... hydrogen-bonding interactions. Two supermolecules {(H-5[26]phen(2)N(4)O(2))X-n(H2O)(5-n)}((5-n)+) form an interpenetrating dimeric species, which was also found by ESI mass spectrum. Binding studies of the protonated macrocycle with aliphatic (ox(2-), mal(2-), suc(2-), cit(3-), cta(3-)) and aromatic (bzc(-), naphc(-), anthc(-), pyrc(-), ph(2-), iph(2-), tph(2-), btc(3-)) anions were determined in water by potentiometric methods. These studies were complemented by H-1 NMR titrations in D2O of the receptor with selected anions. The H-i[26]phen(2)N(4)O(2)(i+) receptor can selectively uptake highly charged or extended aromatic carboxylate anions, such as btc(3-) and pyrc(-), in the pH ranges of 4.0-8.5 and < 4.0, respectively, from aqueous solution that contain the remaining anions as pollutants or contaminants. To obtain further insight into these structural and experimental findings, molecular dynamics (MD) simulations were carried out in water solution.

Computational chemistry for elucidating protein function: energetics and dynamics of myoglobin-ligand systems

State-of-the-art computational methodologies are used to investigate the energetics and dynamics of photodissociated CO and NO in myoglobin (Mb···CO and Mb···NO). This includes the combination of molecular dynamics, ab initio MD, free energy sampling, and effective dynamics methods to compare the results with studies using X-ray crystallography and ultrafast spectroscopy metho ds. It is shown that modern simulation techniques along with careful description of the intermolecular interactions can give quantitative agreement with experiments on complex molecular systems. Based on this agreement predictions for as yet uncharacterized species can be made.

Self-assembly of Fmoc-tetrapeptides based on the RGDS cell adhesion motif

Self-assembly in aqueous solution has been investigated for two Fmoc [Fmoc ¼ N-(fluorenyl)-9-methoxycarbonyl] tetrapeptides comprising the RGDS cell adhesion motif from fibronectin or the scrambled sequence GRDS. The hydrophobic Fmoc unit confers amphiphilicity on the molecules, and introduces aromatic stacking interactions. Circular dichroism and FTIR spectroscopy show that the self-assembly of both peptides at low concentration is dominated by interactions among Fmoc units, although Fmoc-GRDS shows b-sheet features, at lower concentration than Fmoc-RGDS. Fibre X-ray diffraction indicates b-sheet formation by both peptides at sufficiently high concentration. Strong alignment effects are revealed by linear dichroism experiments for Fmoc-GRDS. Cryo-TEM and smallangle X-ray scattering (SAXS) reveal that both samples form fibrils with a diameter of approximately 10 nm. Both Fmoc-tetrapeptides form self-supporting hydrogels at sufficiently high concentration. Dynamic shear rheometry enabled measurements of the moduli for the Fmoc-GRDS hydrogel, however syneresis was observed for the Fmoc-RGDS hydrogel which was significantly less stable to shear. Molecular dynamics computer simulations were carried out considering parallel and antiparallel b-sheet configurations of systems containing 7 and 21 molecules of Fmoc-RGDS or Fmoc-GRDS, the results being analyzed in terms of both intermolecular structural parameters and energy contributions.

Relaxation of surface tension in the liquid-solid interfaces of Lennard-Jones liquids

We have established the surface tension relaxation time in the liquid-solid interfaces of Lennard-Jones (LJ) liquids by means of direct measurements in molecular dynamics (MD) simulations. The main result is that the relaxation time is found to be almost independent of the molecular structures and viscosity of the liquids (at seventy-fold change) used in our study and lies in such a range that in slow hydrodynamic motion the interfaces are expected to be at equilibrium. The implications of our results for the modelling of dynamic wetting processes and interpretation of dynamic contact angle data are discussed.

Mode-specificity and transition state-specific energy redistribution in the chemisorption of CH4 on Ni{100}

We have investigated methane (CH4) dissociative chemisorption on the Ni{100} surface by first-principles molecular dynamics (MD) simulations. Our results show that this reaction is mode-specific, with the n1 state being the most strongly coupled to efficient energy flow into the reaction coordinate when the molecule reaches the transition state. By performing MD simulations for two different transition state (TS) structures we provide evidence of TS structure-specific energy redistribution in methane chemisorption. Our results are compared with recently reported state-resolved measurement of methane adsorption probability on nickel surfaces, and we find that a strong correlation exists between the highest vibrational efficacy measured on Ni{100} for the n1 state and the calculated highest fractional vibrational energy content in this mode.

Thermodynamic Stabillity of Hydrogen-Bonded Systems in Polar and Nonpolar Environments

The thermodynamic properties of a selected set of benchmark hydrogen-bonded systems (acetic acid dimer and the complexes of acetic acid with acetamide and methanol) was studied with the goal of obtaining detailed information on solvent effects on the hydrogen-bonded interactions using water, chloroform, and n-heptane as representatives for a wide range in the dielectric constant. Solvent effects were investigated using both explicit and implicit solvation models. For the explicit description of the solvent, molecular dynamics and Monte Carlo simulations in the isothermal isobaric (NpT) ensemble combined with the free energy perturbation technique were performed to determine solvation free energies. Within the implicit solvation approach, the polarizable continuum model and the conductor-like screening model were applied. Combination of gas phase results with the results obtained from the different solvation models through an appropriate thermodynamic cycle allows estimation of complexation free energies, enthalpies, and the respective entropic contributions in solution. Owing to the strong solvation effects of water the cyclic acetic acid dimer is not stable in aqueous solution. In less polar solvents the double hydrogen bond structure of the acetic acid dimer remains stable. This finding is in agreement with previous theoretical and experimental results. A similar trend as for the acetic acid dimer is also observed for the acetamide complex. The methanol complex was found to be thermodynamically unstable in gas phase as well as in any of the three solvents. (C) 2010 Wiley Periodicals, Inc. J Comput Chem 31: 2046-2055, 2010

Polarization and Spectral Shift of Benzophenone in Supercritical Water

Monte Carlo simulation and quantum mechanics calculations based on the INDO/CIS and TD-DFT methods were utilized to study the solvatochromic shift of benzophenone when changing the environment from normal water to supercritical (P = 340.2 atm and T = 673 K) condition. Solute polarization increases the dipole moment of benzophenone, compared to gas phase, by 88 and 35% in normal and supercritical conditions, giving the in-solvent dipole value of 5.8 and 4.2 D, respectively. The average number of solute-solvent hydrogen bonds was analyzed, and a large decrease of 2.3 in normal water to only 0.8 in the supercritical environment was found. By using these polarized models of benzophenone in the two different conditions of water, we performed MC simulations to generate statistically uncorrelated configurations of the solute surrounded by the solvent molecules and subsequent quantum mechanics calculations on these configurations. When changing from normal to supercritical water environment, INDO/CIS calculations explicitly considering all valence electrons of the 235 solvent water molecules resulted in a solvatochromic shift of 1425 cm(-1) for the most intense transition of benzophenone, that is, slightly underestimated in comparison with the experimentally inferred result of 1700 cm(-1). TD-B3LYP/6-311+G(2d,p) calculations on the same configurations but with benzophenone electrostatically embedded in the 320 water molecules resulted in a solvatochromic shift of 1715 cm(-1) for this transition, in very good agreement with the experimental result. When using the unpolarized model of the benzophenone, this calculated solvatochromic shift was only 640 cm(-1). Additional calculations were also made by using BHandHLYP/6-311+G(2d,p) to analyze the effect of the asymptotic decay of the exchange functional. This study indicates that, contrary to the general expectation, there is a sizable solute polarization even in the low-density regime of supercritical condition and that the inclusion of this polarization is important for a reliable description of the spectral shifts considered here.

Role of Surfactants in Carbon Nanotubes Density Gradient Separation

Several strategies aimed at sorting single-walled carbon nanotubes (SWNT) by diameter and/or electronic structure have been developed in recent years. A nondestructive sorting method was recently proposed in which nanotube bundles are dispersed in water-surfactant solutions and submitted to ultracentrifugation in a density gradient. By this method, SWNTs of different diameters are distributed according to their densities along the centrifuge tube. A mixture of two anionic amphiphiles, namely sodium dodecylsulfate (SIDS) and sodium cholate (SC), presented the best performance in discriminating nanotubes by diameter. We present molecular dynamics studies of the water-surfactant-SWNT system. The simulations revealed one aspect of the discriminating power of surfactants: they can actually be attracted toward the interior of the nanotube cage. The binding energies of SDS and SC on the outer nanotube surface are very similar and depend weakly on diameter. The binding inside the tubes, on the contrary, is strongly diameter dependent: SDS fits best inside tubes with diameters ranging from 8 to 9 angstrom, while SC is best accommodated in larger tubes, with diameters in the range 10.5-12 angstrom. The dynamics at room temperature showed that, as the amphiphile moves to the hollow cage, water molecules are dragged together, thereby promoting the nanotube filling. The resulting densities of filled SWNT are in agreement with measured densities.

A Multitechnique Study of Structure and Dynamics of Polyfluorene Cast Films and the Influence on Their Photoluminescence

This article describes the microstructure and dynamics in the solid state of polyfluorene-based polymers, poly(9,)-dioctylfluorenyl-2,7-diyl) (PFO), a semicrystalline polymer, and poly [(9,9-dioctyl- 2,7-divinylene-fluorenylene)-alt-co-{2-methoxy-5-(2-ethyl-hexyloxy)- 1,4-phenylene vinylene}, a copolymer with mesomorphic phase properties. These Structures were determined by wide-angle X-ray scattering (WAXS) measurements, Assuming a packing model for the copolymer structure, where the planes of the phenyl rings are stacked and separated by an average distance of similar to 4.5 angstrom and laterally spaced by about similar to 16 angstrom, we followed the evolution of these distances as a function of temperature using WAXS and associated the changes observed to the polymer relaxation processes identified by dynamical mechanical thermal analysis. Specific molecular motions were studied by solid-state nuclear magnetic resonance. The onset of the side-chain motion at about 213 K (beta-relaxation) produced a small increase in the lateral spacing and in the stacking distance of the phenyl rings in them aggregated Structures, Besides, at about 383 K (alpha-relaxation) there occurs a significant increase in the amplitude of the torsion motion in the backbone, producing a greater increase in the stacking distance of the phenyl rings. Similar results were observed in the semicrystalline phase of PFO, but in this case the presence of the crystalline structure affects considerably the overall dynamics, which tends to be more hindered. Put together, Our data explain many features of the temperature dependence of the photoluminescence of these two polymers.

Effect of SO(2) on the Transport Properties of an Imidazolium Ionic Liquid and Its Lithium Solution

Transport coefficients have been measured as a function of the concentration of sulfur dioxide, SO(2), dissolved in 1-butyl-2,3-dimethylimidazolium bis(trifluoromethylsulfonyl)-imide, [BMMI][Tf(2)N], as well as in its lithium salt solution, Li[Tf(2)N]. The SO(2) reduces viscosity and density and increases conductivity and diffusion coefficients in both the neat [BMMI] [Tf(2)N] and the [BMMI][Tf(2)N]-Li[Tf(2)N] solution. The conductivity enhancement is not assigned to a simple viscosity effect; the weakening of ionic interactions upon SO(2) addition also plays a role. Microscopic details of the SO(2) effect were unraveled using Raman spectroscopy and molecular dynamics (MD) simulations. The Raman spectra suggest that the Li(+)-[Tf(2)N] interaction is barely affected by SO(2), and the SO(2)-[Tf(2)N] interaction is weaker than previously observed in an investigation of an ionic liquid containing the bromide anion. Transport coefficients calculated by MD simulations show the same trend as the experimental data with respect to SO(2) content. The MD simulations provide structural information on SO(2) molecules around [Tf(2)N], in particular the interaction of the sulfur atom of SO(2) with oxygen and fluorine atoms of the anion. The SO(2)-[BMMI] interaction is also important because the [BMMI] cations with above-average mobility have a larger number of nearest-neighbor SO(2) molecules.

Unraveling Dynamical Heterogeneity in the Ionic Liquid 1-Butyl-3-methylimidazolium Chloride

Heterogeneous dynamics within a time range of nanoseconds was investigated by molecular dynamics (MD) simulations of 1 -butyl-3-methylimidazolium chloride ([bmim]Cl). After identifying groups of fast and slow ions, it was shown that the separation between the location of the center of mass and the center of charge of cations, d(CMCC), is a signature of such difference in ionic mobility. The distance d(CMCC) can be used as a signature because it relaxes in the time window of the dynamical heterogeneity. The relationship between the parameter dcmcc and conformations of the side alkyl chain in [bmim] is discussed. Since the relatively slow relaxation of dcmcc is a consequence of coexisting polar and nonpolar domains in the bulk, the MD simulations reveal a subtle interplay between structural and dynamical heterogeneity in ionic liquids.

Host-guest system of 4-nerolidylcatechol in 2-hydroxypropyl-beta-cyclodextrin: preparation, characterization and molecular modeling

The interaction of 4-nerolidylcatechol (4-NRC), a potent antioxidant agent, and 2-hydroxypropyl-beta-cyclodextrin (HP-beta-CD) was investigated by the solubility method using Fourier transform infrared (FTIR) methods in addition to UV-Vis, (1)H-nuclear magnetic resonance (NMR) spectroscopy and molecular modeling. The inclusion complexes were prepared using grinding, kneading and freeze-drying methods. According to phase solubility studies in water a B(S)-type diagram was found, displaying a stoichiometry complexation of 2:1 (drug:host) and stability constant of 6494 +/- A 837 M(-1). Stoichiometry was established by the UV spectrophotometer using Job`s plot method and, also confirmed by molecular modeling. Data from (1)H-NMR, and FTIR, experiments also provided formation evidence of an inclusion complex between 4-NRC and HP-beta-CD. 4-NRC complexation indeed led to higher drug solubility and stability which could probably be useful to improve its biological properties and make it available to oral administration and topical formulations.

Shielding of ionic interactions by sulfur dioxide in an ionic liquid

The effect of adding SO(2) on the structure and dynamics of 1-butyl-3-methylimidazolium bromide (BMIBr) was investigated by low-frequency Raman spectroscopy and molecular dynamics (MD) simulations. The MD simulations indicate that the long-range structure of neat BMIBr is disrupted resulting in a liquid with relatively low viscosity and high conductivity, but strong correlation of ionic motion persists in the BMIBr-SO(2) mixture due to ionic pairing. Raman spectra within the 5 < omega < 200 cm(-1) range at low temperature reveal the short-time dynamics, which is consistent with the vibrational density of states calculated by MD simulations. Several time correlation functions calculated by MD simulations give further insights on the structural relaxation of BMIBr-SO(2).

Transport coefficients, Raman spectroscopy, and computer simulation of lithium salt solutions in an ionic liquid

Lithium salt solutions of Li(CF3SO2)(2)N, LiTFSI, in a room-temperature ionic liquid (RTIL), 1-butyl-2,3-dimethyl-imidazolium cation, BMMI, and the (CF3SO2)(2)N-, bis(trifluoromethanesulfonyl)imide anion, [BMMI][TFSI], were prepared in different concentrations. Thermal properties, density, viscosity, ionic conductivity, and self-diffusion coefficients were determined at different temperatures for pure [BMMI][TFSI] and the lithium solutions. Raman spectroscopy measurements and computer simulations were also carried out in order to understand the microscopic origin of the observed changes in transport coefficients. Slopes of Walden plots for conductivity and fluidity, and the ratio between the actual conductivity and the Nernst-Einstein estimate for conductivity, decrease with increasing LiTFSI content. All of these studies indicated the formation of aggregates of different chemical nature, as it is corroborated by the Raman spectra. In addition, molecular dynamics (MD) simulations showed that the coordination of Li+ by oxygen atoms of TFSI anions changes with Li+ concentration producing a remarkable change of the RTIL structure with a concomitant reduction of diffusion coefficients of all species in the solutions.