Protection of DNA against low-energy electrons by amino acids: A first-principles molecular dynamics study

Using first-principles molecular dynamics simulations, we have investigated the notion that amino acids can play a protective role when DNA is exposed to excess electrons produced by ionizing radiation. In this study we focus on the interaction of glycine with the DNA nucleobase thymine. We studied thymine-glycine dimers and a condensed phase model consisting of one thymine molecule solvated in amorphous glycine. Our results show that the amino acid acts as a protective agent for the nucleobase in two ways. If the excess electron is initially captured by the thymine, then a proton is transferred in a barrier-less way from a neighboring hydrogen-bonded glycine. This stabilizes the excess electron by reducing the net partial charge on the thymine. In the second mechanism the excess electron is captured by a glycine, which acts as a electron scavenger that prevents electron localization in DNA. Both these mechanisms introduce obstacles to further reactions of the excess electron within a DNA strand, e.g. by raising the free energy barrier associated with strand breaks.

Understanding the Interaction between Low-Energy Electrons and DNA Nucleotides in Aqueous Solution

Reactions that can damage DNA have been simulated using a combination of molecular dynamics and density functional theory. In particular, the damage caused by the attachment of a low energy electron to the nucleobase. Simulations of anionic single nucleotides of DNA in an aqueous environment that was modeled explicitly have been performed. This has allowed us to examine the role played by the water molecules that surround the DNA in radiation damage mechanisms. Our simulations show that hydrogen bonding and protonation of the nucleotide by the water can have a significant effect on the barriers to strand breaking reactions. Furthermore, these effects are not the same for all four of the bases.

A Phase-field Model for Computer Simulation of Lamellar Structure Formation in y-TiAI

In the present paper, a phase-field model is developed to simulate the formation and evolution of lamellar microstructure in γ-TiAl alloys. The mechanism of formation of TiAl lamellae proposed by Denquin and Naka is incorporated into the model. The model describes the formation and evolution of the face-centered cubic (fcc) stacking lamellar zone followed by the subsequent appearance and growth of the γ-phase, involving both the chemical composition change by atom transfer and the ordering of the fcc lattice. The thermodynamics of the model system and the interaction between the displacive and diffusional transformations are described by a non-equilibrium free energy formulated as a function of concentration and structural order parameter fields. The long-range elastic interactions, arising from the lattice misfit between the α, fcc (A1) and the various orientation variants of the γ-phase are taken into account by incorporating of the elastic strain energy into the total free energy. Simulation studies based on the model successfully predicted some essential features of the lamellar structure. It is found that the formation and evolution of the lamellar structure are predominantly controlled by the minimization of the elastic energy of the interfaces between the different fcc stacking groups, low-symmetry product phase γ and the high-symmetry α-phase, as well as between the various orientation variants of the product phase.

Stability of Sr adatom model structures for SrTiO3(001) surface reconstructions

We report results of first-principles calculations on the thermodynamic stability of different Sr adatom structures that have been proposed to explain some of the observed reconstructions of the (001) surface of strontium titanate (Kubo and Nozoye 2003 Surf Sci. 542 177). From surface free energy calculations, a phase diagram is constructed indicating the range of conditions over which each structure is most stable. These results are compared with Kubo and Nozoye's experimental observations. It is concluded that low Sr adatom coverage structures can only be explained if the surface is far from equilibrium. Intermediate coverage structures are stable only if the surface is in or very nearly in equilibrium with the strontium oxide.

Kinetics of phase transformations in a model with metastable fluid-fluid separation: A molecular dynamics study

By molecular dynamics (MD) simulations we study the crystallization process in a model system whose particles interact by a spherical pair potential with a narrow and deep attractive well adjacent to a hard repulsive core. The phase diagram of the model displays a solid-fluid equilibrium, with a metastable fluid-fluid separation. Our computations are restricted to fairly small systems (from 2592 to 10368 particles) and cover long simulation times, with constant energy trajectories extending up to 76x10(6) MD steps. By progressively reducing the system temperature below the solid-fluid line, we first observe the metastable fluid-fluid separation, occurring readily and almost reversibly upon crossing the corresponding line in the phase diagram. The nucleation of the crystal phase takes place when the system is in the two-fluid metastable region. Analysis of the temperature dependence of the nucleation time allows us to estimate directly the nucleation free energy barrier. The results are compared with the predictions of classical nucleation theory. The critical nucleus is identified, and its structure is found to be predominantly fcc. Following nucleation, the solid phase grows steadily across the system, incorporating a large number of localized and extended defects. We discuss the relaxation processes taking place both during and after the crystallization stage. The relevance of our simulation for the kinetics of protein crystallization under normal experimental conditions is discussed. (C) 2002 American Institute of Physics.

A dielectric continuum molecular dynamics method

We introduce a novel method to simulate hydrated macromolecules with a dielectric continuum representation of the surrounding solvent. In our approach, the interaction between the solvent and the molecular degrees of freedom is described by means of a polarization density free energy functional which is minimum at electrostatic equilibrium. After a pseudospectral expansion of the polarization and a discretization of the functional, we construct the equations of motion for the system based on a Car-Parrinello technique. In the limit of the adiabatic evolution of the polarization field variables, our method provides the solution of the dielectric continuum problem "on the fly," while the molecular coordinates are propagated. In this first study, we show how our dielectric continuum molecular dynamics method can be successfully applied to hydrated biomolecules, with low cost compared to free energy simulations with explicit solvent. To our knowledge, this is the first time that stable and conservative molecular dynamic simulations of solutes can be performed for a dielectric continuum model of the solvent. (C) 2001 American Institute of Physics.

Determining the electronic structure and chemical potentials of molecules in solution

A simulation scheme is proposed for determining the excess chemical potential of a substance in solution. First, a Monte Carlo simulation is performed with classical models for solute and solvent molecules. A representative sample of these configurations is then used in a hybrid quantum/classical (QM/MM) calculation, where the solute is treated quantum-mechanically, and the average electronic structure is used to construct an improved classical model. This procedure is iterated to self-consistency in the classical model, which in practice is attained in one or two steps, depending on the quality of the initial guess. The excess free energy of the molecule within the QM/MM approach is determined relative to the classical model using thermodynamic perturbation theory with a cumulant expansion. The procedure provides a method of constructing classical point charge models appropriate for the solution and gives a measure of the importance of solvent fluctuations.

Concertedness and solvent effects in multiple proton transfer reactions: The formic acid dimer in solution

The issue of multiple proton transfer (PT) reactions in solution is addressed by performing molecular dynamics simulations for a formic acid dimer embedded in a water cluster. The reactant species is treated quantum mechanically, within a density functional approach, while the solvent is represented by a classical model. By constraining different distances within the dimer we analyze the PT process in a variety of situations representative of more complex environments. Free energy profiles are presented, and analyzed in terms of typical solvated configurations extracted from the simulations. A decrease in the PT barrier height upon solvation is rationalized in terms of a transition state which is more polarized than the stable states. The dynamics of the double PT process is studied in a low-barrier case and correlated with solvent polarization fluctuations. Cooperative effects in the motion of the two protons are observed in two different situations: when the solvent polarization does not favor the transfer of one of the two protons and when the motion of the two protons is not synchronized. This body of observations is correlated with local structural and dynamical properties of the solvent in the vicinity of the reactant. (C) 2000 American Institute of Physics. [S0021-9606(00)51121-0].

Solvation Structure and Transport of Acidic Protons in Ionic Liquids: A First-principles Simulation Study

Ab initio simulations of a single molecule of HCl in liquid dimethyl imidazolium chloride [dmim][Cl] show that the acidic proton exists as a symmetric, linear ClHCl- species. Details of the solvation structure around this molecule are given. The proton-transfer process was investigated by applying a force along the antisymmetric stretch coordinate until the molecule broke. Changes in the free energy and local solvation structure during this process were investigated. In the reaction mechanism identified, a free chloride approaches the proton from the side. As the original ClHCl- distorts and the incoming chloride forms a new bond to the proton, one of the original chlorine atoms is expelled and a new linear molecule is formed.

Neutral and Charged 1-Butyl-3-methylimidazolium Triflate Clusters: Equilibrium Concentration in the Vapor Phase and Thermal Properties of Nanometric Droplets

Ground state energy, structure, and harmonic vibrational modes of 1-butyl-3-methylimidazolium triflate ([bmim][Tf]) clusters have been computed using an all-atom empirical potential model. Neutral and charged species have been considered up to a size (30 [bmim][Tf] pairs) well into the nanometric range. Free energy computations and thermodynamic modeling have been used to predict the equilibrium composition of the vapor phase as a function of temperature and density. The results point to a nonnegligible concentration of very small charged species at pressures (P ~ 0.01 Pa) and temperatures (T 600 K) at the boundary of the stability range of [bmim][Tf]. Thermal properties of nanometric neutral droplets have been investigated in the 0 T 700 K range. A near-continuous transition between a liquidlike phase at high T and a solidlike phase at low T takes place at T ~ 190 K in close correspondence with the bulk glass point Tg ~ 200 K. Solidification is accompanied by a transition in the dielectric properties of the droplet, giving rise to a small permanent dipole embedded into the solid cluster. The simulation results highlight the molecular precursors of several macroscopic properties and phenomena and point to the close competition of Coulomb and dispersion forces as their common origin.

Thermodynamic Studies of Ionic Interactions in Aqueous Solutions of Imidazolium-Based Ionic Liquids [Emim][Br] and [Bmim][Cl]

Experimental measurements of density at different temperatures ranging from 293.15 to 313.15 K, the speed of sound and osmotic coefficients at 298.15 K for aqueous solution of 1-ethyl-3-methylimidazolium bromide ([Emim][Br]), and osmotic coefficients at 298.15 K for aqueous solutions of 1-butyl-3-methylimidazolium chloride ([Bmim][Cl]) in the dilute concentration region are taken. The data are used to obtain compressibilities, expansivity, apparent and limiting molar properties, internal pressure, activity, and activity coefficients for [Emim][Br] in aqueous solutions. Experimental activity coefficient data are compared with that obtained from Debye-Hückel and Pitzer models. The activity data are further used to obtain the hydration number and the osmotic second virial coefficients of ionic liquids. Partial molar entropies of [Bmim][Cl] are also obtained using the free-energy and enthalpy data. The distance of the closest approach of ions is estimated using the activity data for ILs in aqueous solutions and is compared with that of X-ray data analysis in the solid phase. The measured data show that the concentration dependence for aqueous solutions of [Emim][Br] can be accounted for in terms of the hydrophobic hydration of ions and that this IL exhibits Coulombic interactions as well as hydrophobic hydration for both the cations and anions. The small hydration numbers for the studied ILs indicate that the low charge density of cations and their hydrophobic nature is responsible for the formation of the water-structure-enforced ion pairs.

Effect of solution pH, ionic strength, and temperature on adsorption behavior of reactive dyes on activated carbon

The adsorption behavior of C.I. Reactive Blue 2, C.I. Reactive Red 4, and C.I. Reactive Yellow 2 from aqueous solution onto activated carbon was investigated under various experimental conditions. The adsorption capacity of activated carbon for reactive dyes was found to be relatively high. At pH 7.0 and 298 K, the maximum adsorption capacity for C.I. Reactive Blue 2, C.I. Reactive Yellow 2 and C.I. Reactive Red 4 dyes was found to be 0.27, 0.24, and 0.11 mmol/g, respectively. The shape of the adsorption isotherms indicated an L2-type isotherm according to the Giles and Smith classification. The experimental adsorption data showed good correlation with the Langmuir and Ferundlich isotherm models. Further analysis indicated that the formation of a complete monolayer was not achieved, with the fraction of surface coverage found to be 0.45, 0.42, and 0.22 for C.I. Reactive Blue 2, C.I. Reactive Yellow 2 and C.I. Reactive Red 4 dyes, respectively. Experimental data indicated that the adsorption capacity of activated carbon for the dyes was higher in acidic rather than in basic solutions, and further indicated that the removal of dye increased with increase in the ionic strength of solution, this was attributed to aggregation of reactive dyes in solution. Thermodynamic studies indicated that the adsorption of reactive dyes onto activated carbon was an endothermic process. The adsorption enthalpy (?H) for C.I. Reactive Blue 2 and C.I. Reactive Yellow 2 dyes were calculated at 42.2 and 36.2 kJ/mol, respectively. The negative values of free energy (?G) determined for these systems indicated that adsorption of reactive dyes was spontaneous at the temperatures under investigation (298-328 K). © 2007 Elsevier Ltd. All rights reserved.

Scaling of Superdomain Bands in Ferroelectric Dots

Bundles of 90° stripe domains have been observed to form into distinct groups, or bands, in mesoscale BaTiO3 single crystal dots. Vector piezoresponse force microscopy (PFM) shows that each band region, when considered as a single entity, possesses a resolved polarization that lies approximately along the pseudocubic direction; antiparallel alignment of this resultant polarization in adjacent bands means that these regions can be considered as 180° “superdomains.” For dots with sidewall dimensions below ~2 microns, Landau–Kittel like scaling in the width of these superdomains was observed, strongly suggesting that they form in response to lateral depolarizing fields. In larger dot structures, scaling laws break down. We have rationalized these observations by considering changes in the driving force for the adoption of equilibrium superdomain periodicities implied by Landau–Kittel-free energy models; we conclude that the formation of ordered bands of superdomains is a uniquely meso/nanoscale phenomenon. We also note that the superdomain bands found by PFM imaging in air contrast with the quadrant arrangements seen previously by Schilling et al. (Nano Lett., 9, 3359 (2009)) through transmission electron microscopy imaging in vacuum. The importance of the exact nature of the boundary conditions in determining the domain patterns that spontaneously form in nanostructures is therefore clearly implied.

Molecular dynamics and principal components analysis of human telomeric quadruplex multimers.

Guanine-rich DNA repeat sequences located at the terminal ends of chromosomal DNA can fold in a sequence-dependent manner into G-quadruplex structures, notably the terminal 150–200 nucleotides at the 3' end, which occur as a single-stranded DNA overhang. The crystal structures of quadruplexes with two and four human telomeric repeats show an all-parallel-stranded topology that is readily capable of forming extended stacks of such quadruplex structures, with external TTA loops positioned to potentially interact with other macromolecules. This study reports on possible arrangements for these quadruplex dimers and tetramers, which can be formed from 8 or 16 telomeric DNA repeats, and on a methodology for modeling their interactions with small molecules. A series of computational methods including molecular dynamics, free energy calculations, and principal components analysis have been used to characterize the properties of these higher-order G-quadruplex dimers and tetramers with parallel-stranded topology. The results confirm the stability of the central G-tetrads, the individual quadruplexes, and the resulting multimers. Principal components analysis has been carried out to highlight the dominant motions in these G-quadruplex dimer and multimer structures. The TTA loop is the most flexible part of the model and the overall multimer quadruplex becoming more stable with the addition of further G-tetrads. The addition of a ligand to the model confirms the hypothesis that flat planar chromophores stabilize G-quadruplex structures by making them less flexible.

Volumetric properties, viscosity and refractive index of the protic ionic liquid, pyrrolidinium octanoate, in molecular solvents

Densities ([rho]) and viscosities ([eta]) of binary mixtures containing the Protic Ionic Liquid (PIL), pyrrolidinium octanoate with five molecular solvents: water, methanol, ethanol, n-butanol, and acetonitrile are determined at the atmospheric pressure as a function of the temperature and within the whole composition range. The refractive index of all mixtures (nD) is measured at 298.15†K. The excess molar volumes VE and deviation from additivity rules of viscosities [eta]E and refractive index [Delta][phi]n, of pyrrolidinium octanoate solutions were then deduced from the experimental results as well as apparent molar volumes V[phi]i, partial molar volumes and thermal expansion coefficients [alpha]p. The excess molar volumes VE are negative over the entire mole fraction range for mixture with water, acetonitrile, and methanol indicating strong hydrogen-bonding interaction for the entire mole fraction. In the case of longest carbon chain alcohols (such as ethanol and n-butanol)†+†pyrrolidinium octanoate solutions, the VE variation as a function of the composition describes an S shape. The deviation from additivity rules of viscosities is negative over the entire composition range for the acetonitrile, methanol, ethanol, and butanol, and becomes less negative with increasing temperature. Whereas, [eta]E of the {[Pyrr][C7CO2]†+†water} binary mixtures is positive in the whole mole fraction range and decreases with increasing temperature. the excess Gibbs free energies of activation of viscous flow ([Delta]G*E) for these systems were calculated. The deviation from additivity rules of refractive index [Delta][phi]n are positive over the whole composition range and approach a maximum of 0.25 in PIL mole fraction for all systems. The magnitude of deviation for [Delta][phi]n describes the following order: water†>†methanol†>†acetonitrile†>†ethanol. Results have been discussed in terms of molecular interactions and molecular structures in these binary mixtures.