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.

Tracing the Chromospheric and Coronal Magnetic Field with AIA, IRIS, IBIS, and ROSA Data

The aim of this study is to explore the suitability of chromospheric images for magnetic modeling of active regions. We use high-resolutionimages (≈0.2"-0.3"), from the Interferometric Bidimensional Spectrometer in the Ca II 8542 Å line, the Rapid Oscillations in the Solar Atmosphere instrument in the Hα 6563Å line, the Interface Region Imaging Spectrograph in the 2796Å line, and compare non-potential magnetic field models obtainedfrom those chromospheric images with those obtained from images of the Atmospheric Imaging Assembly in coronal (171 Å, etc.) and inchromospheric (304 Å) wavelengths. Curvi-linear structures are automatically traced in those images with the OCCULT-2 code, to which we forward-fitted magnetic field lines computed with the Vertical-current Approximation Nonlinear Force Free Field code. We find that the chromospheric images: (1) reveal crisp curvi-linear structures (fibrils, loop segments, spicules) that are extremely well-suited for constraining magnetic modeling; (2) that these curvi-linear structures arefield-aligned with the best-fit solution by a median misalignment angle of μ2 ≈ 4°–7° (3) the free energy computed from coronal data may underestimate that obtained from chromospheric data by a factor of ≈2–4, (4) the height range of chromospheric features is confined to h≲4000 km, while coronal features are detected up to h = 35,000 km; and (5) the plasma-β parameter is β ≈ 10^-5 - 10^-1 for all traced features. We conclude that chromospheric images reveal important magnetic structures that are complementary to coronal images and need to be included in comprehensive magnetic field models, something that is currently not accomodated in standard NLFFF codes.

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].

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 and enthalpies of solution of alcohols CkH2k+1OH (k = 1, 2, 6) in 1-methyl-3-alkylimidazolium bis(trifluoromethylsulfonyl)imide {[C1CnIm][NTf2] n = 2, 4, 6, 8, 10} ionic liquids

We present a study on the effect of the alkyl chain length of the imidazolium ring in 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ionic liquids, [C1CnIm][NTf2] (n = 2 to 10), on the mixing properties of (ionic liquid + alcohol) mixtures (enthalpy and volume). We have measured small excess molar volumes with highly asymmetric curves as a function of mole fraction composition (S-shape) with more negative values in the alcohol-rich regions. The excess molar volumes increase with the increase of the alkyl-chain length of the imidazolium cation of the ionic liquid. The values of the partial molar excess enthalpy and the enthalpy of mixing are positive and, for the case of methanol, do not vary monotonously with the length of the alkyl side-chain of the cation on the ionic liquid – increasing from n = 2 to 6 and then decreasing from n = 8. This non-monotonous variation is explained by a more favourable interaction of methanol with the cation head group of the ionic liquid for alkyl chains longer than eight carbon atoms. It is also observed that the mixing is less favourable for the smaller alcohols, the enthalpy of mixing decreasing to less positive values as the alkyl chain of the alcohol increases. Based on the data from this work and on the knowledge of the vapour pressure of {[C1CnIm][NTf2] + alcohol} binary mixtures at T = 298 K reported in the literature, the excess Gibbs free energy, excess enthalpy and excess entropy could be then calculated and it was observed that these mixtures behave like the ones constituted by a non-associating and a non-polar component, with its solution behaviour being determined by the enthalpy.

Understanding the Performance of Melt-Extruded Poly(ethylene oxide)–Bicalutamide Solid Dispersions: Characterisation of Microstructural Properties Using Thermal, Spectroscopic and Drug Release Methods

In this article, we have prepared hot-melt-extruded solid dispersions of bicalutamide (BL) using poly(ethylene oxide) (PEO) as a matrix platform. Prior to preparation, miscibility of PEO and BL was assessed using differential scanning calorimetry (DSC). The onset of BL melting was signi?cantly depressed in the presence of PEO, and using Flory– Huggins (FH) theory, we identi?ed a negative value of -3.4, con?rming miscibility. Additionally, using FH lattice theory, we estimated the Gibbs free energy of mixing which was shown to be negative, passing through a minimum at a polymer fraction of 0.55. Using these data, solid dispersions at drug-to-polymer ratios of 1:10, 2:10 and 3:10 were prepared via hot-melt extrusion. Using a combination of DSC, powder X-ray diffractometry and scanning electron
microscopy, amorphous dispersions of BL were con?rmed at the lower two drug loadings. At the 3:10 BL to PEO ratio, crystalline BL was detected. The percent crystallinity of PEO was reduced by approximately 10% in all formulations following extrusion. The increased amorphous content within PEO following extrusion accommodated amorphous BL at drug to polymer loadings up to 2:10; however, the increased amorphous domains with PEO following extrusion were not suf?cient to fully accommodate BL at drug-to-polymer ratios of 3:10.